search for: acefix

Hi Ace, You can just read the file first to find out: max_fields<-function(file,sep=" ") { rlines<-readLines(file) return(max(unlist(lapply(sapply(rlines,strsplit,sep),length)))) } nmax<-max_fields(test.txt,"\t") Jim On Wed, Aug 30, 2017 at 2:22 AM, Fix Ace <acefix at rocketmail.com> wrote: > Thank you very much! Looks like I have to know the length of each record > ahead of time. > > Ace > > > On Monday, August 28, 2017 12:56 AM, Jim Lemon <drjimlemon at gmail.com> wrote: > > > Hi Ace, > With tabs as separators: >...

...s as separators: testdf<-read.table("test.txt",header=FALSE,fill=TRUE,sep="\t", col.names=paste("V",1:19,sep=""),stringsAsFactors=FALSE) Also note that I got the number of columns wrong the first time. Jim On Mon, Aug 28, 2017 at 12:56 PM, Fix Ace <acefix at rocketmail.com> wrote: > Hi, Jim, > > Thank you very much for pointing out the format issue. Here is the original > text: > > === > I have a text file (test.txt) with different number of columns: > > 0610007P14Rik%%% Tcf19 Gtf2i > 0610010O12Rik%%% Ivns1abp Etv6 &...

...s as separators: testdf<-read.table("test.txt",header=FALSE,fill=TRUE,sep="\t", col.names=paste("V",1:19,sep=""),stringsAsFactors=FALSE) Also note that I got the number of columns wrong the first time. Jim On Mon, Aug 28, 2017 at 12:56 PM, Fix Ace <acefix at rocketmail.com> wrote: > Hi, Jim, > > Thank you very much for pointing out the format issue. Here is the original > text: > > === > I have a text file (test.txt) with different number of columns: > > 0610007P14Rik%%% Tcf19 Gtf2i > 0610010O12Rik%%% Ivns1abp Etv6 &...

Hi, Jim, Thank you very much for pointing out the format issue. Here is the original text: ===I have a text file (test.txt) with different number of columns: 0610007P14Rik%%% Tcf19 Gtf2i 0610010O12Rik%%% Ivns1abp Etv6 1100001G20Rik%%% Nmi 1500015O10Rik%%% Foxi1 Ascl3 Sirt3 1700003E16Rik%%% Ascl2 Ifnar2 1700028J19Rik%%% Musk Nfe2l3 1810011O10Rik%%% Ppp1r13b Bpnt1 Cdkn2c Foxc1 Sox10 Smarca2

...to the r-sig-mixed-models list instead. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Nov 15, 2017 at 5:09 AM, Fix Ace <acefix at rocketmail.com> wrote: > Hi, Bert, > > Sorry about that! David seemed to be able to read the post since he > replied. Here I just email you the sample code and error message: > > > example.3=data.frame(levels=as.numeric(XXX[,c(4)]), > replicate=rep(c("0",&qu...

...to the r-sig-mixed-models list instead. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Nov 15, 2017 at 5:09 AM, Fix Ace <acefix at rocketmail.com> wrote: Hi, Bert, Sorry about that! David seemed to be able to read the post since he replied. Here I just email you the sample code and error message: > example.3=data.frame(levels= as.numeric(XXX[,c(4)]), replicate=rep(c("0","1","2"," 3&...

Hi Ace, As your example seems to have spaces as separators, testdf<-read.table("test.txt",header=FALSE,fill=TRUE, col.names=paste("V",1:14,sep=""),stringsAsFactors=FALSE) By specifying the number of columns with "col.names" and using "fill=TRUE" you can get a data frame with zero length strings where values are missing in the input file. Jim

Hi, David, Thank you very much for getting back to me! Sorry about the messy code example. I am re-posting here (including the error message): > example.3=data.frame(levels=as.numeric(XXX[,c(4)]),replicate=rep(c("0","1","2","3","4","5"),3),conditions=c(rep("11",6),rep("12",6),rep("13",6)))> example.3? ?

Dear R Community, Is there a way to adjust the line width of the dendrogram in heat map.2 function? I tried "lwd", but it didn't work... Also, is there a way to adjust the general height/width/position of the dendrogram using heatmap.2 function? I feel the portion of the dendrogram is huge compared with the overall graph. I would like to squeeze them a bit. Thank you very much for

Dear R community, I would like to plot a partial hclust output, so I?am looking for a subtree function that would return an tree structure I can plot. I ran the test code of subtree following the instruction on?http://finzi.psych.upenn.edu/library/extracat/html/subtree.html However, an error message popped up: ====> library(extracat) Attaching package: ?extracat? The following object is masked

Dear R Community, I recently read an article and found a plot as attached. It has scatterplot, barplot, and error bar. Could anyone help me to figure out what package I can use in R to generate such plot? Thank you very much for any inputs! Kind regards, Ace -------------- next part -------------- A non-text attachment was scrubbed... Name: Screen Shot 2018-02-04 at 6.38.14 AM.png Type: image/png

Dear R Community, I tried to generate heatmap for a matrix of 1500 columns by 106 rows using the following R script: > pheatmap(tf.vs.DE.1.removeAllZeroCol, fontsize=3,border_color=NA) and got the graph (as attached Fig 1) Since the column labels appear very crowded, I tried to increase the cellwidth to stretch the graph horizontally. The idea was to show the graph section by section, but with

Dear R Community, My data have 3 conditions and each condition has 6 replicates. I am trying to fit my data for a linear mixed model using the lmer function from lme4 package to find the random effects of the replicates; however, I got the error message. Here are the example codes:

Dear R community, I have a data matrix (531X314), and would like to apply the prcomp. However, I got this error Lapack message. I am using R3.2.2. I googled a bit and found that it might be related to converge issue. ?Just wonder if there is a way to get around it? Thank you very much! Ace On Thursday, December 29, 2016 11:44 AM, Ista Zahn <istazahn at gmail.com> wrote: Use

Dear R community, I have a text file (test.txt) with different number of columns: 0610007P14Rik%%% Tcf19 Gtf2i 0610010O12Rik%%% Ivns1abp Etv6 1100001G20Rik%%% Nmi 1500015O10Rik%%% Foxi1 Ascl3 Sirt3 1700003E16Rik%%% Ascl2 Ifnar2 1700028J19Rik%%% Musk Nfe2l3 1810011O10Rik%%% Ppp1r13b Bpnt1 Cdkn2c Foxc1 Sox10 Smarca2 1810019D21Rik%%% Asb8 1810037I17Rik%%% Zfp612 1810055G02Rik%%% Nkx2-3 Maged1 Runx1