lfranzen / scwat-st Goto Github PK
View Code? Open in Web Editor NEWVisium (10x Genomics) on human abdominal s.c white adipose tissue.
License: MIT License
Visium (10x Genomics) on human abdominal s.c white adipose tissue.
License: MIT License
Hello @lfranzen, I will bother you again
Block {r clustering_annotations_export}
which starts in line 471 gives me
Error in file(file, ifelse(append, "a", "w")) : cannot open the connection
6.
file(file, ifelse(append, "a", "w"))
5.
write.table(c_anno, file = file.path(DIR_RES, "tables", fname), row.names = F, col.names = TRUE, sep = ",", dec = ".", qmethod = "double")
4.
eval(expr, p)
3.
eval(expr, p)
2.
eval.parent(Call)
1.
write.csv(c_anno, file = file.path(DIR_RES, "tables", fname), row.names = F)
It seems to be specifically in the second line of the block, which is
write.csv(c_anno, file = file.path(DIR_RES, "tables", fname), row.names = F)
It seems to be a problem with the path. I see DIR_RES is defined in the very first block {r setup, include=FALSE}
, as DIR_RES <- file.path(DIR_ROOT, "results" , PROJECT_ID)
. No folder called results
is available and none was created when I ran the script. I decided to create both the results
and figures
folder manually and ran the code from the start, but I still got the error.
Hello
I am trying to reproduce the code and got the following error on block 179-185 in the visium_baseline.main_analysis.Rmd file.
Here is the traceback
Error in open.connection(con, "rb") : cannot open the connection
12.
open.connection(con, "rb")
11.
open(con, "rb")
10.
parse_con(txt, bigint_as_char)
9.
parseJSON(txt, bigint_as_char)
8.
parse_and_simplify(json, simplifyVector = simplifyVector, ...)
7.
parse_json(file(path), simplifyVector = simplifyVector, ...)
6.
read_json(f)
5.
FUN(X[[i]], ...)
4.
lapply(X = X, FUN = FUN, ...)
3.
sapply(infotable[, "json"], function(f) { read_json(f)$tissue_hires_scalef })
2.
scaleVisium %||% sapply(infotable[, "json"], function(f) { read_json(f)$tissue_hires_scalef })
1.
InputFromTable(infotable = infoTable, min.gene.count = 100, min.gene.spots = 5, min.spot.count = 200, platform = "Visium")
I ran into another problem from the block in line 322
#' Read in original se object in case minor random changes in the SCTransform/ICA/UMAP has occurred
se <- readRDS(file = file.path(DIR_DATA, "se-object.visium_baseline_ori.rds"))
se@tools$Staffli@imgs <- infoTable$imgs
The error is
Error in gzfile(file, "rb") : cannot open the connection
2.
gzfile(file, "rb")
1.
readRDS(file = file.path(DIR_DATA, "se-object.visium_baseline_ori.rds"))
I see no rds called se-object.visium_baseline_ori.rds
, in the data folder, nor in the downloadable data. I also see no saveRDS()
command upstream, so I assume it is not created by the script. There were no errors or warnings previous to this. There is a .zip file called se-object.visium_baseline.rds.zip
though, which I am unsure if it is the same. However, when I extract it, it creates a folder instead of a rds file. Is this working as intended?
A declarative, efficient, and flexible JavaScript library for building user interfaces.
๐ Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. ๐๐๐
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google โค๏ธ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.