231 results • Page 1 of 4
I try to use ConsensusClusterPlus to subgroup my samples. This package will generate a random seed number and put the seed number in the output log file. However, I found I couldn’t repeat the result with the seed number even...other question is about the reproducibility. Using my own data, the output from manually assigned a seed number (seed = as.numeric(Sys.time())) is very different from the …
reading the vignette for DelayedMatrixStats, I'm confused about the mention of two types of RleArray seeds: SolidRleArraySeed & ChunkedRleArraySeed What's the difference between these two seeds and a simple DataFrame...of Rle's as seed? Are there perfomance/memory differences between them? Are there any pro/cons associated with each of them
updated 7.2 years ago • maltethodberg
group. Cannot calculate variance. This problem may be resolved if you try again with a different seed. Is there anyone who can help me to understand what I do wrong? Many thanks for your attention. Best, Stefano </div
updated 17.3 years ago • Stefano Moretti
If &gt; you have an example DCF file I could look at, The BSgenomeForge vignette contains the full seed file that was used to forge the BSgenome.Rnorvegicus.UCSC.rn4 package. If you need more examples of seed files, the vignette...also gives you a list of more than 30 seed files that you can look at. All these seed files are text files in the DCF format, that is, the information in them is o…
Hi, Is it recommended to set a seed before running AUCell build rankings if you want to generate the same output each time you run it? From what I understand
updated 3.2 years ago • Lucy
getting acquainted with the EBImage package and generating an R script to analyze x-ray images of seed. We have about 150 x-ray images that are 5888 x 4620 pixels. Within the image, there are two "reference disks" we 3D printed and...between 0.20 and 0.25. Here is an example image I generated using this package. You can see the 2 reference disks, 10 seeds, and lengths and widths lines drawn on…
updated 2.1 years ago • Brian
Hi, I'm using the 'segment' function on my data. In the manual it's mentioned that a seed should be specified in order to make the results reproducible. However, even when I use a seed ```set.seed(73684)``` the results...5, undo.splits="sdundo", undo.SD=20) ``` Is there any reason why, even when specifying a seed, the outcome of segment can be different when the same input is used? ```…
updated 21 months ago • BID
devel/bioc/vignettes/BSgenome/inst/doc/BSgenomeForge.pdf. However, I have trouble to prepare the seed file. Just wondering whether you have any suggestions on which editor I can use to generate DCF format - Debian Control
updated 5.3 years ago • jiazhou0116
div class="preformatted"> I'm preparing a custom reference genome for use with the MEDIPS package. I see that one field of the seed file, which is apparently not optional, is the
updated 12.7 years ago • Guest User
subset of probes to do the between array normalization. In order to achive reproducible results, the seed needs to be set using set.seed." So I am wondering: 1. Why is it unnecessary to set the seed for intra-array normalization and
updated 10.2 years ago • maden.sean
Chr13.fa Chr14.fa Chr15.fa Chr16.fa Chr17.fa Chr18.fa Chr19.fa.fa: file(s) not found &gt; &gt; My seed file: &gt; &gt; Package: BSgenome.Ptrichocarpa.Phytozome.v3 &gt; Title: Populus trichocarpa full genome (Phytozome version...guest]; bioconductor at r-project.org; MEDIPS Maintainer &gt;&gt; Subject: Re: [BioC] BSgenomeForge seed file - seqnames field &gt;&gt;…
from 2005. Here is my directory with the .fasta files for each chromosomes and the seed file: ```bash ls /media/nmary/DONNEES/Abeille/Temp/ BSgenome.Amellifera.NCBI.Amel-HAv3.1.seed chrLG15.fa.gz chrLG5.fa.gz...chrLG9.fa.gz chrLG14.fa.gz chrLG4.fa.gz ``` Here is my seed file ```bash cat /media/nmary/DONNEES/Abeille/Temp/BSgenome.Amellifera.NCBI.Amel…
updated 5.1 years ago • n.mary
Kelly &gt; Sent: Thursday, December 12, 2013 6:51 AM &gt; To: Hervé Pagès &gt; Subject: Success: RE: Seed file error: RE: BioStrings for current R: RE: [BioC] BSgenomeForge seed file - seqnames field &gt; &gt; Hi Herve, &gt; I have now succeeded...10, 2013 2:38 PM &gt; To: Vining, Kelly &gt; Cc: bioconductor at r-project.org &gt; Subject: Re: Seed file error: R…
ConsensusClusterPlus', which automatically sets the seed to 'as.numeric(Sys.time())' if it is not specified with the 'seed' argument; hence it ignores any seeds set externally with 'set.seed...found the FlowSOM package to be very useful in our CyTOF data analysis pipelines, and having this seed argument in the release version would make things easier for getting reproducible results. In addition…
updated 8.8 years ago • Lukas Weber
I have performed an RNA-seq experiment on leaves (five time points) and maturing seeds (three time points - from early maturation to dry seed) on a non-model plant species, and used the Salmon -&gt; tximport -&gt; DESeq2...data looks pretty good. However, I have come across an interesting issue when it comes to the dry seed data (last seed time point). Although most genes are not DE (base…
updated 7.4 years ago • Poikilo
comparison. I have two treatments (WT and mutant), each treatment in three tissues (root, leaf, seed), each tissue either has replicates or not. I want to compare the expression between WT and mutant in general (considering...treatment=factor(c(rep('WT',8),rep('MUTANT',8))) tissue=factor( c('root','root','root','leaf','seed','seed','seed','seed', 'root','root','root','leaf','seed','seed','seed',…
updated 14.2 years ago • Xiaohui Wu
gave me a warning&nbsp;of " WARNING: Remove set.seed usage in R code" because I used set.seed(seed) in my R function to reproduce&nbsp;results, where "seed" is an argument of&nbsp;my R function and&nbsp;can be specified by the...user. To bypass Bioccheck, I have to remove the set.seed(seed) in my R function, but I want to let the user have the option to select a seed&nbsp;in m…
updated 7.5 years ago • xyluo1991
nmf(V, 3, my.algorithm, scale.factor=10) But I get the following error: Error in metagenes(fit(seed)) &lt;- abs(pca$rotation[, 1:factorization.rank]) : could not find function "metagenes&lt;-" What is wrong here? This is the algorithm...my.algorithm &lt;- function(x, seed, scale.factor=1){ # do something with starting point # for example: # 1. compute principal components pc…
updated 14.8 years ago • Pavel Goldstein
generateSeeds' requires an matrix of integer. Coerced &gt; QUBIC.1 &lt;- quBicluster(seeds=data.seeds, eset=data.exprs, report.no=30) *** caught segfault *** address 0xffffffff1efabdd0, cause 'memory not mapped' Traceback...1: .Call("qubicluster", seeds = seeds, exprs = exp.leveled, sigma_val = sigma, symbols = symbols, report_no = report.no, tolerance_val = tolerance, filter_…
updated 11.7 years ago • Guest User
sce &lt;- cluster(sce, features = "type", xdim = 10, ydim = 10, maxK = maxk, seed = 1) ``` Each time I run the above function in a new R session I get slightly different results (e.g. percentages, clusters) using
updated 2.1 years ago • elefth.pavlos
vignettes/BSgenome/inst/doc/BSgenomeForge.pdf> I have: 1. Genome file (gzipped) of V.radiata. 2. Seed file called BSgenome.Vradiata.UCSC.vr1-seed created referring to the manual (above link) section __2.2 Prepare the...BSgenome data package seed file __at page\#4. Now, I ran following commands: <pre> &gt; library(BSgenome) &gt; forgeBSgenomeDataPkg("path/to/my/seed")</…
updated 9.2 years ago • Vijay Lakhujani
initial package (say, using gff files), or am I limited such that I would have to re-forge a new reference? &gt; &gt; A BSgenome data package is for storing the sequences of a given &gt; genome/assembly. No annotations should...in a BSgenome package would duplicate a lot of DNA &gt; sequences that are already in the reference sequences (those &gt; small sequences…
you generate? Best Rodrigo -- output of sessionInfo(): &gt; forgeBSgenomeDataPkg("path/to/my/seed") Error: no se pudo encontrar la funci'on "forgeBSgenomeDataPkg" &gt; library(BSgenome) Loading required package: GenomicRanges...2: package 'GenomicRanges' was built under R version 2.12.1 &gt; forgeBSgenomeDataPkg("path/to/my/seed") Error en .readSeedFile(x, verbose = verbose) : …
<div class="preformatted">Hi everybody, In forging my custom Bsgenome data package I encountered a problem with gap masks. For my genome of interest, NCBI has 2 different gap masks for each assembled chromosome: the chr?.comp.agp file (chromosome from component AGP) and the chr?.agp file (chromosome from scaffold AGP). There is only one agp file for the unlocalized and one for the unplaced…
updated 12.5 years ago • Ugo Borello
I already have the fa files of hg38 and HIV-1 from NCBI. I keep getting errors when running the seed file. I don't know what to use for the genome. In addition, can you help me find out what's wrong with this seed file? Thanks in
updated 2.2 years ago • 13273976287
for one chromosome) and put it into one folder: <pre> marcin@marcin:/media/marcin/ngs/maizeGenome/seed$ ls chr10.fa chr2.fa chr4.fa chr6.fa chr8.fa maize chr1.fa chr3.fa chr5.fa chr7.fa chr9.fa</pre> I also create seed file...chr4", "chr5", "chr6", "chr7", "chr8", "chr9", "chr10") seqs_srcdir: /media/marcin/ngs/maizeGenome/seed </pre> Unfortunately, after run…
updated 10.3 years ago • Marcin Grzybowski
<div class="preformatted">Pretty much copied the code from http://linusben.net/wiki/index.php/Oligo.bioconductor. Something ala, """ &gt; library(pdInfoBuilder) &gt; filename_ndf &lt;- list.files('.',pattern='.ndf',full.names=T) &gt; filename_xys &lt;- list.files('.',pattern=".xys",full.names=T)[1] &gt; seed &lt;- new("NgsExpressionPDInfoPkgSeed", ndfFile=filen…
updated 15.9 years ago • Chris Fenton
Hey, I have a question regarding the reproducibility of limma. I am using it for differential expression analysis for proteomics experiments. I also use the `decideTests` function with the option 'global' in order to correct for multiple testing over &gt;1 contrasts. I am having some problems regarding the reproducibility of the DE calculations: the results I am seeing in R Studio during …
updated 2.9 years ago • TSauer
I kept getting the Error: could not find function "forgeSeqFiles". I know I need to make the seed files so I gave the command line: forgeSeqlengthsFile(seqnames, prefix="pi1&gt;", suffix=".fa", seqs_srcdir="/Users/Me/Documents...Microarray", seqs_destdir="/Users/Me/Documents/Microarray/Seeds", verbose=TRUE) Error: could not find function "forgeSeqFiles". I would appreciate it if you could…
forging a BSgenome I'm getting this error : ``` &gt; "Error in .readSeedFile(x, verbose = verbose) : seed file &gt; 'SARSCOV2-seed' must have exactly 1 record" ``` This is the code I'm running: ``` &gt; library(BSgenome) myfile &lt;- read.dcf("seed_file.txt
updated 5.6 years ago • Didi
AnnotationDbi, verbatim, and hit an error in 'code chunk 22', where calling: &gt; makeAnnDbPkg(seed, dbfile, dest_dir = tempdir()) Error in matrix(unlist(value, recursive = FALSE, use.names = FALSE), nrow = nr, : negative extents to...number 22: build new package (eval = FALSE) ################################################### &gt; seed &lt;- new("AnnDbPkgSeed", Packa…
Im a beginner in R and I want to forge Casuarina equistifolia genome in BSgenome package as it is not available as inbuilt package. After spending so much time in solving errors but still I couldn't solve one error. Kindly help me to solve the below error. #Seed File I have used ``` Package: BSgenome.Cequisetifolia.NCBI.ASM379533v1 Title: Full genome sequences for Casuarina equisetifolia...tim…
updated 3.4 years ago • Manoj Kumar R
https://support.bioconductor.org/p/72378/ I followed the same steps and created z and seed. I am not able to install it. This is the error that R is throwing me. ``` &gt; install.packages("pd.mogene21stmmense/", type="source...tmp.zip") unzip("tmp.zip") dir() z &lt;- cdf2table("mogene21st_Mm_ENSE.cdf") seed &lt;- new("GenericPDInfoPkgSeed", table=z, author = "me", email = "pr…
updated 6.7 years ago • pbachali
<div class="preformatted"> Hello, I'm trying to build a data package for Sus scrofa with BSgenome (R version 2.13.2 and BSgenome version 1.20.0). At the bottom of this email I've copied my seed file. I've downloaded the sequence files from UCSC and checked the md5sums. I've also downloaded the gap.txt and masks files...BSgenome (R version 2.13.2 and BSgenome version 1.20.0). At the bottom…
pattern = ".ndf", full.names = TRUE) xys &lt;- list.files(pattern = ".xys", full.names = TRUE) seed &lt;- new("NgsExpressionPDInfoPkgSeed", ndfFile = ndf, xysFile = xys, author = "MOZAR Michael", email = "michael.mozar@univ-perp.fr...biocViews = "AnnotationData", genomebuild = "Retroryza", organism = "Rice", species = "Oryza Sativa", url = "none") here my sessionInfo : library(pdInfoBui…
updated 14.7 years ago • Michael MOZAR
I noticed that DiffBind returns different grey lists every time I ran dba.blacklist even if I set seed as follows (set.seed(42)). If there is a hint to solve this problem, please let me know. Thank you. ``` &gt;library("DiffBind") &gt;set.seed
updated 2.3 years ago • r.murakami
NROW) s = generateSeeds(m) eset &lt;- quantileDiscretize(ExpressionSet(m)) x = quBicluster(seeds=s, eset=eset, report.no=30) In the next release I shall consider making generateSeeds and quantileDiscretize internal...cause 'memory not mapped' &gt;&gt;&gt; &gt;&gt;&gt; Traceback: &gt;&gt;&gt; 1: .Call("qubicluster", seeds = seeds, exprs = exp.…
updated 11.5 years ago • Zhang, Jitao David
why are these two so different: `` plotPCA(splatSimulateGroups(groupCells=c(1000,1000), nGenes=21, seed=4)) `` `` plotPCA(splatSimulateGroups(groupCells=c(1000,1000), nGenes=21, seed=5
updated 8.6 years ago • luke.zappia
<div class="preformatted">SZSN Expands To Become 3rd Largest Agricultural Seed Provider! Shandong Zhouyuan Seed and Nursery Co., Ltd (SZSN) $0.24 SZSN is expanding. Recent acquisitions put it in the top 3 Seed providers in China. It is also expanding its distribution chains to 500 new regional agencies. Big news expected! Get on...div class="preformatted">SZSN Expands To Become 3rd Lar…
updated 18.5 years ago • Jack Whaley
Plants, Ayurveda, Bio-fuel, Environment, Herbs, Medicinal Plants, Organic Farming, Pharmaceutical, Seed-technology, Spices. We are selling the above publications not only to make money but also for the development of this...the botanical species, the volatile oils in these plants are stored and preserved in roots, leaves, seeds, bark, fruits, flowers and stem which may be emitting characteristic …
updated 16.4 years ago • Sanjeev Prajapati
Construction fo the package: library(pdInfoBuilder) setwd("/my_bioc_packages/") seed &lt;- new("AffyHTAPDInfoPkgSeed", version = "3.8.0", license = "Artistic-2.0", pgfFile = ".../HTA-2_0.r1.pgf", clfFile = ".../HTA-2_0.r1.clf", probeFile...hg19", organism = "Homo sapiens", species = "…
I am using the rbsurv packaeg to select surival-assosiated probes. My data was pre-processed microarray data with ~7,000 probes and 192 tumor samples. Data was split in half for modeling and validation(train and test data).&nbsp; I also included in my model a risk predictor, tumor stage, to help model fitting. But I am confused as to how to determine a suitable value for the&nbsp;`` max.n…
updated 8.2 years ago • JackieMe
use pdInfoBuilder to make an annotation package. However, I get an error after the makePdInfoPackage(seed, destDir = "./") command. Error in aggregate.data.frame(as.data.frame(x), ...) : no rows to aggregate Has anybody encountered this...TRUE)) [1] "/cchome/qichen2/Nimblegen_methylation/21863/081229_HG18_Promoter_MeDI P_HX1. pos" &gt; seed &lt;- new("NgsTilingPDInfoPkgSeed", ndfFile…
c("X", "Y")) 3: parseNgsPair(object@ndfFile, object@xysFile, verbose = !quiet) 2: makePdInfoPackage(seed, destDir = ".") 1: makePdInfoPackage(seed, destDir = ".") and this is my script : &gt; library(pdInfoBuilder) &gt; baseDir &lt;- "/home/michael...pattern = ".xys", full.names = TRUE)[1]) [1] "/home/michael/Documents/test_2007/nobless1.xys" &gt; seed &lt;- new("NgsExpressi…
<div class="preformatted">Dear list, For R version 2.13.0, I had no problems running the following code (CMA package was used to load golub data, but detach afterwards trying to isolate the problem): library(CMA) data(golub) detach(package:CMA) library(e1071) golubY &lt;- golub[,1] golubX &lt;- as.matrix(golub[,-1]) ratio &lt;- 2/3 golubY2 &lt;- as.factor(golubY) levels(go…
updated 13.8 years ago • Javier Pérez Florido
<div class="preformatted">I offer below a short example of using rGADEM with MotifDb which you may find useful. The example motivates and then demonstrates "seeded" search, in which a candidate TF is provided to rGADEM. (The rGADEM vignette illustrates the simpler task of de novo motif search. - Paul On Feb 27, 2013, at 9:03 AM, Zhu, Lihua (Julie) wrote: &gt; Jose, &gt; &g…
updated 12.8 years ago • Paul Shannon
over 7000 unplaced contigs and scaffolds. This means that the `` seqnames: `` line in the BSgenome seed file must be very very long (136762 characters), and I seem to be running into some line length limitation in R that causes...that actually parses. I'm not sure where to go from here. Can anyone offer some help? [Here is the seed file I'm using](https://gist.github.com/DarwinAwardWinner/329e40…
updated 11.0 years ago • Ryan C. Thompson
the makeAnnDbPkg failed due to an error (see below) Any suggestions? &gt;myPath &lt;- getwd() &gt; seed &lt;- new("AnnDbPkgSeed", Package="ATHPfam", Version="1.0", AnnObjPrefix="ATHPfam") &gt; makeAnnDbPkg(seed, dbfile="Pfam.sqlite", dest_dir
updated 17.6 years ago • Samuel Wuest
number of clusters and the cell-to-cluster assignments are different. This is despite setting the seed at the very beginning (`bp` is what I pass to all the functions that have a `BPPARAM` argument): ```{r} set.seed(16) bp &lt;- BiocParallel...component of the UMAP and Louvain clustering algorithms, though I would have thought setting the seed was enough. Interestingly, I can't quite reprodu…
Make info package MUST have only one POS and ## ## NDF file in the directory ## #################################################### seed &lt;- new("NgsTilingPDInfoPkgSeed", ndfFile = ndf.file, xysFile = xys.files[1], posFile = pos.file, author = "x", email = "x at gmail.com...Yeast", manufact…
getting a warning which suggests they aren't. Is this the case, or have I done something wrong in my seed file? ```r &gt; library(BSgenome) &gt; forgeBSgenomeDataPkg("R://Sophie/Pea_RNASeq/genomes/BSGenomes_seed_PisumSativum.txt...seqs_srcdir = seqs_srcdir, destdir = destdir, : field 'release_name' is deprecated ``` The seed file reads: ``` Package: BSgenome.Psativum.URGI.…
updated 5.0 years ago • sophie.harrington
<div class="preformatted">Hi Peter, Your input file looks OK at 1st glance. Are you trying to make an annotation package? Have you tried following the instructions in the vignette titled "SQLForge: An easy way to create a new annotation package with a standard database schema."? If not, you can read it here: http://www.bioconductor.org/packages/2.11/bioc/html/AnnotationForge.ht ml …
updated 12.8 years ago • Marc Carlson
Hi, I am trying to use ConsensusClusterPlus to do cluster classification. After clustering, how can I get the class ID of patients? For example, results = ConsensusClusterPlus(d,maxK=6,reps=50,pItem=0.8,pFeature=1, title=title,clusterAlg="hc",distance="pearson",seed=1262118388.71279,plot='pdf') Then, I wanna divide patients by 2 groups. How can I get the c…
updated 2.7 years ago • hsin hsiung
motif but in other runs it disappears. Now i am trying to find the motif with some motif as a seed in my DNAstringset object of sequences. My seeded motif is given here is present in motif.txt file &nbsp; &nbsp;&nbsp; A&nbsp...nbsp; &nbsp;0.3224&nbsp;&nbsp; &nbsp;0.2151 Please tell me possible command to get similar motif to seeded motif in DNAstringset. &n…
updated 10.0 years ago • vinod.acear
r forgeBSgenomeDataPkg("cws_analysis/BSGenome_NC_035897/BSgenome.Amexicanus.SIMR.astMexNC2-seed") Error in .make_Seqinfo_from_genome(genome) : "astMexNC" is not a registered NCBI assembly or UCSC genome (use registered_NCBI_assemblies...or UCSC assemblies/genomes currently registered in the GenomeInfoDb package) ``` Here is the seed file: ``` Package: BSgenome.Amexicanus.SIMR.…
updated 4.7 years ago • Chris Seidel
clusters 2015-06-12 14:09:48 calculatePvalues: calculating F-statistics for permutation 1 and seed 20140924 2015-06-12 14:09:51 findRegions: segmenting F-stats information 2015-06-12 14:09:51 findRegions: identifying...regions 2015-06-12 14:09:51 calculatePvalues: calculating F-statistics for permutation 2 and seed 20140925 2015-06-12 14:09:53 findRegions: segmenting F-stats inform…
updated 10.5 years ago • drramki.chop
Array","MANUFACTURERURL"="http://www.affymetrix.com "),baseMapType="eg",outputDir="primeview"); &gt; seed&lt;-new("AnnDbPkgSeed",Package = "primeview.db",Version="1.0.0",PkgTemplate="HUMANCHIP.DB",AnnObjPrefix ="primeview"); &gt; makeAnnDbPkg...seed,file.path("primeview","primeview.sqlite"),dest_dir=" primeview"); Error in match(x, table, nomatch = 0L) : object 'doc_template_names
updated 12.7 years ago • wang peter
Even I tried to change the arguments of the flowpeaks but it didn't make any change. I found the seeds that must be formed within the 4 clusters, they are found at the midst of the 4 clusters which makes them dividing horizontally
updated 8.5 years ago • shadhiksk
Hello, I am trying to forge a genome for a non-model organism. I have generated the following seed file (in dcf format) and using the available Athaliana seed file as an example. Package: BSgenome.Bnapus.NCBI.Bra_napus_v2.0...Hello, I am trying to forge a genome for a non-model organism. I have generated the following seed file (in dcf format) and using the available Athaliana seed…
updated 4.7 years ago • sieminsk
dot edu ############################################################ &gt; B &lt;- 1000 &gt; seed &lt;- 99 &gt; TTBoot &lt;- MTP(X=esetGcrmaExprsFiltered, Y=TT, test = "t.twosamp.unequalvar", alternative = "two.sided", typeone="fdr...method="ss.maxT", fdr.method="conservative", B=B, seed=seed) running bootstrap... iteration = 100 200 300 400 500 600 700 800 900 1000 …
updated 17.9 years ago • Timur Shtatland
231 results • Page 1 of 4
Traffic: 1201 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6