search for: zhongm

Dear All: I need draw some figure through trellis.device and save them as pdf files. How can I close trellis.device (something like dev.off() in nonlattice figure)? Many thanks Zhongming Yang --------------------------------- [[alternative HTML version deleted]]

...: How can I catch R errors and depend on this error call other R functions in Perl? The foolowing is the error message when I use perl call R. Error in nls(modout ~ exp(-b1 * modin)/(b2 + b3 * modin), trace = F, start = c(b1 = 0.005, : singular gradient Segmentation fault Thanks, Zhongming

Hi, I use lm or loess to make smoothing. After smoothing I need "Residual Standard Error" in my script. Could you please tell me how can I get this information? Thanks,

Hi, I am trying to use piecewise linear regression to approximate a nonlinear function. Actually, I don't know how many linear functions I need, therefore, I want build an array of regression models to automate the approximation job. Could you please give me any clue? Attached is ongoing code: rawData = scan("c:/zyang/mass/data/A01/1.PRN", what=list(numeric(),numeric())); len =

Dear All: Could you please tell me how I can draw figure formatted like figure 4.18 of MASS (4th) with the attached data set? Thanks Zhongming Yang --------------------------------- -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: sample.txt Url: https://stat.ethz.ch/pipermail/r-help/attachments/20050316/abfdb85e/sample.txt

Hi, I want to use R handle hundreds data file. So I want use some string class to handle the names of data file. My idea procedure like: for i in 1:20 { filename = "***" + i +".PRN" scan(filename) .... } Do anyone know how to do this? Thanks,

An embedded message was scrubbed... From: "Zhongming Yang" <Zhongming.Yang at cchmc.org> Subject: how to make exprSet Date: Mon, 14 Jul 2003 12:16:33 -0400 Size: 683 Url: https://stat.ethz.ch/pipermail/r-help/attachments/20030714/6e4061d3/attachment.mht

Here is the test code. foo<-function(x) exp(-x) for (i in 1:5) { foo<-function(x) foo(x)*x foo(2) } The error is "evalution nested too deeply". I tried Recall() but it didn't work either. Thanks a lot for your input.

Dear All: I am working on writing some R functions to make statistical reports automatically. Dr. Harrell's Hmisc has all the wonderful stuff. But sometimes I need change some formats, so I want to read through it and make some modifications to fit my project. Ideally, I want proceed as following: 1. change some source of Hmisc 2. compile and install the modified Hmisc 3. debug my

...decided by the prefix of the idno). I want draw a repeat measurement comparison figure. My goal is to list 5 patients from same group on one horizontal line. But xyplot sounds pick them randomly (or I was confused?). Could you please help me modify the following code to accomplish this? Thanks Zhongming Yang library(nlme) library(lattice) md <- read.table('../data/sample.txt', header=T) md <- groupedData(md ~ month | idno, data=md) trellis.device(theme=col.whitebg()) xyplot(md ~ month | idno, data=md, main="title", xlab="x", ylab="y", layout=c(...

...ol(surface='interpolate', statistics='approximate', trace.hat='approximate')); modpre = predict(mod.model); mod.std = sqrt(var(modout-modpre, na.rm=TRUE)); (#mod.std = sqrt(var(modout[1:length(modpre)] - modpre[1:length(modpre)], na.rm=TRUE));) print(mod.std); Thanks Zhongming