search for: yao

...91 3906 0 NULL 275 4219 0 NULL 273 4236 1 49410 281 4163 1 49408 274 4226 1 49406 295 3869 2 49403 49404 287 4113 0 NULL 295 3871 1 49401 292 3895 4 49396 49397 49398 49399 291 3900 3 49392 library (plyr) and library (reshape2) and other good packages are OK for me. Thanks a lot! Yao He ????????????????????????? Master candidate in 2rd year Department of Animal genetics & breeding Room 436,College of Animial Science&Technology, China Agriculture University,Beijing,100193 E-mail: yao.h.1988 at gmail.com ??????????????????????????

...-read.table("03-04.txt",header=T) I don't want to type hundreds of read.table(),so how I read it in one time? I think the core problem is that I can't create different objects' name in the use of loop or sapply() ,but there may be a better way to do what I want. Thanks a lot Yao He Yao He -- ????????????????????????? Master candidate in 2rd year Department of Animal genetics & breeding Room 436,College of Animial Science&Technology, China Agriculture University,Beijing,100193 E-mail: yao.h.1988 at gmail.com ??????????????????????????

...en, the source for RCS should be easily compiled and installed, I used to do it all the time. Except there is no GCC installed. Well, I can get gcc from sunfreeware ... enter a loop. Any suggestions? Thanks. -- /*********************************************************************\ ** ** Joe Yao jsdy@tux.org - Joseph S. D. Yao ** \*********************************************************************/

Could anyone please explain to me why the following writes nothing into "all.Rout" file? If the "for" loop is removed, t.test output can be written into "all.out". Thanks in advance. Minghua Yao ...... zz <- file("all.Rout", open="wt") sink(zz) for(i in 1:n) { Cy3<-X[,2*i-1]; Cy5<-X[,2*i]; t.test(Cy3, Cy5) } sink() close(zz) ......

...ifferent variables in a loop , so I want to add all results to an object then dump() them to an text. But I don't know how to append T-test result to the object? I have already plot the barplot and I want to know an elegant way to report raw result. Can anybody give me some pieces of advice? Yao He ????????????????????????? Master candidate in 2rd year Department of Animal genetics & breeding Room 436,College of Animial Science&Technology, China Agriculture University,Beijing,100193 E-mail: yao.h.1988 at gmail.com ??????????????????????????

...1 2 4633 TWF2H23 44.42 40.89 12.96 0.3169479 I try this code > aggregate(df[,4:7],df[,1],mean) But I couldn't set the agrument na.rm=T in the mean() function,so the results are all NAs Please tell me how to handle NA values in the use of aggregate() Thanks a lot Yao He ????????????????????????? Master candidate in 2rd year Department of Animal genetics & breeding Room 436,College of Animial Science&Technology, China Agriculture University,Beijing,100193 E-mail: yao.h.1988 at gmail.com ??????????????????????????

...e about 180 individuals and 12 variables and 60000 Fixed effect estimates As asreml-R is not free ,is there any packages for my study? I heard nlme or lme4 but I'm not sure whether they could incorporate covariates and what about their computational efficiency? Thanks for you recommendation Yao He ????????????????????????? Master candidate in 2rd year Department of Animal genetics & breeding Room 436,College of Animial Science&Technology, China Agriculture University,Beijing,100193 E-mail: yao.h.1988 at gmail.com ??????????????????????????

...t;,"r") geno_t<-list() repeat{ line<-readLines(c,n=1) if (length(line)==0)break #end of file line<-unlist(strsplit(line,"\t")) geno_t<-cbind(geno_t,line) } write.table(geno_t,"xxx.txt") It works but it is too slow ,how to optimize it??? Thank you Yao He ????????????????????????? Master candidate in 2rd year Department of Animal genetics & breeding Room 436,College of Animial Science&Technology, China Agriculture University,Beijing,100193 E-mail: yao.h.1988 at gmail.com ??????????????????????????

Hello, I have R1.9.0 under Windows XP. My program plots several plots using x11() par(cex = 0.75) ...... x11() par(cex = 0.75) ...... x11() par(cex = 0.75) ...... x11() par(cex = 0.75) ...... Sometimes, one of them generates a small frame only with title area "R graphics: Device X (ACTIVE)". The message in the console window is Error in plot.new() : Figure margins too large

...t;FID", and fill the bar with different colour. I try these codes: p<-ggplot(data,aes(x=FID)); p+geom_bar(aes(x=factor(FID),y=..count..,fill=STATUS)) But how could I exclude "nosperm" or other levels just in the use of ggplot2 without generating another dataframe Thanks a lot Yao He ————————————————————————— Master candidate in 2rd year Department of Animal genetics & breeding Room 436,College of Animial Science&Technology, China Agriculture University,Beijing,100193 E-mail: yao.h.1988@gmail.com <mingxie@vt.edu> —————————————————————————— [[alternative HTML...

Hi, all, I want to use lowess(x, y) where x and y are vectors of length of 4000+. In fact, x and y are log of some vectors. So, some of the elements are NaN. lowess() can not take away those elements then do the fitting. It will give the error message and do nothing. 1. Can anybody tell me how to get rid of those NaN's and use lowess()? 2. How to get the LOWESS fitting values for any

...led catalog" for a specific node (on the master node)? And apply this on the client node? What are the limitations for doing this? Basically we have special handling for the webserver/CA part of puppet. Thank you very much! -Yushu +-------------------------------------------------+ | Yushu Yao | Ph:1-510-486-4690 | | Lawrence Berkeley National Lab | Mailstop 50B-6222 | 1 Cyclotron Road | Berkeley CA 94720-8147 - USA +-------------------------------------------------+ -- You received this message because you are subscribed to the Google Groups "Puppet Users" group. To post to...

Hi, Is there any function that can arrange a character in alphabetic order? Thanks for answer -MY [[alternative HTML version deleted]]

...off virtfs, but also enables ZFS file-based vdevs on virtfs to be used without killing QEMU. Also, special thanks to both Avi Kivity and Alexander Graf for their interpretation of QEMU backtraces. Without their guidence, tracking down this bug would have taken much longer. Signed-off-by: Richard Yao <ryao at gentoo.org> --- net/9p/trans_virtio.c | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/net/9p/trans_virtio.c b/net/9p/trans_virtio.c index 9c5a1aa..5d1d04b 100644 --- a/net/9p/trans_virtio.c +++ b/net/9p/trans_virtio.c @@ -340,7 +340,10 @@ static int p9_get_ma...

...off virtfs, but also enables ZFS file-based vdevs on virtfs to be used without killing QEMU. Also, special thanks to both Avi Kivity and Alexander Graf for their interpretation of QEMU backtraces. Without their guidence, tracking down this bug would have taken much longer. Signed-off-by: Richard Yao <ryao at gentoo.org> --- net/9p/trans_virtio.c | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/net/9p/trans_virtio.c b/net/9p/trans_virtio.c index 9c5a1aa..5d1d04b 100644 --- a/net/9p/trans_virtio.c +++ b/net/9p/trans_virtio.c @@ -340,7 +340,10 @@ static int p9_get_ma...

After upgrading R to 2.9.2, I can't use the anova() fuction. It says "could not find function "Varcov" ". What's wrong with my computer? Help needed, thanks! Yao Zhu Department of Urology Fudan University Shanghai Cancer Center No. 270 Dongan Road, Shanghai, China [[alternative HTML version deleted]]

Hi Richard, On 11/25/2013 04:50 PM, Richard Yao wrote: > I figured out the problem. There is zerocopy IO is being done via DMA to > a buffer allocated with valloc(). Right now, I am running a hack-fix > locally so I can get some other stuff done first. I will propose a > proper fix to the list in a few days. I've also encountere...

Hi Richard, On 11/25/2013 04:50 PM, Richard Yao wrote: > I figured out the problem. There is zerocopy IO is being done via DMA to > a buffer allocated with valloc(). Right now, I am running a hack-fix > locally so I can get some other stuff done first. I will propose a > proper fix to the list in a few days. I've also encountere...

Dear all, I have the following output generated by linear regression. Since there is only one regression intercept and one slope for one set of data, what is the meaning of std. error for intercept and that of slope? Thanks in advance. Sincerely, Minghua > data(thuesen) > attach(thuesen) > lm(short.velocity~blood.glucose) Call: lm(formula = short.velocity ~ blood.glucose)

...frame(x,y) library(rms) fito<-lrm(y~x) preo<-predict(fito) pre<-matrix(NA,nrow=100,ncol=200) for (i in 1:200) { sam<-sample(1:nrow(dat)) sam<-split(sam,1:10) for (j in 1:10) { fit<-lrm(y~x,data=dat[-sam[[j]],]) pre[sam[[j]],i]<-predict(fit,data=dat[sam[[j]],]) } } *Yao Zhu* *Department of Urology Fudan University Shanghai Cancer Center Shanghai, China* [[alternative HTML version deleted]]