search for: wegelin

...not yet found an answer to this question in Paul Murrell's book. Does Deepayan Sarkar's _Lattice_ go into that kind of detail? (I use lattice, not traditional graphics, because the plot will eventually be conditioned on a third categorical variable.) Thanks for any suggestions. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwegelin@vcu.edu URL: http://www.people.vcu.edu/~jwegelin [[alternative HTML version deleted]]

...ot;col" argument has no effect. I'm on a Windows 2000 machine with > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 5.0 year 2002 month 04 day 29 language R Thanks for any advice Jake Jacob A. Wegelin, Ph.D. Adjunct Assistant Professor Epidemiology and Preventive Medicine, Univ. California One Shields Ave, TB-168 Davis CA 95616-8638 USA TEL 530.752.2793 FAX 530.752.3239 http://wegelin.ucdavis.edu/ jawegelin at ucdavis.edu

...av") Error in read.spss("indata/z2EXvideo.sav") : Error reading system-file header. In addition: Warning message: indata/z2EXvideo.sav: File layout code has unexpected value 50331648. Value should be 2, in big-endian or little-endian format. Thanks for any information Jake Wegelin > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 8.0 year 2003 month 10 day 08 language R

...ot( yield ~ nitro , data=Oats , scales=list( x=list( at=unique(Oats$nitro) ) ) , subset=Variety=="Victory" ) # But the following returns an error: xyplot( yield ~ nitro , data=Oats , scales=list( x=list( at=unique(nitro) ) ) ) Thanks for any insight Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwegelin at vcu.edu URL: http://www.people.vcu.edu/~jwegelin

...t the human eye does not interpret area accurately, so that making the area of a circle proportional to sample size may be a flawed approach. Thus I'd also be interested in any comments on how one might best graphically display x, y, and sample size in a single plot.) Thanks for any tips Jake Wegelin

...ch_1.2.8 RColorBrewer_1.0-5 [13] reshape2_1.2.2 scales_0.2.3 stringr_0.6.2 > Incidentally (and for the sake of keyword searches): Although a google search initially led me to posts about expression() and plotmath, those eventually had nothing to do with the solution. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 830 E. Main St., Seventh Floor P. O. Box 980032 Richmond VA 23298-0032 U.S.A. CTSA grant: UL1TR000058 URL: http://www.people.vcu.edu/~jwegelin

...wnloaded 18Kb > library("lme4") > groupedData Error: Object "groupedData" not found > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 9.1 year 2004 month 06 day 21 language R Jake Wegelin

...thematically savvy? For instance, someone who doesn't remember any math beyond algebra and doesn't think in terms of P(A|B)? Or have we done a disservice to our students by choosing to make them learn R, rather than making ourselves learn SAS, Stata or SPSS? Thank you for any ideas Jake Wegelin

...[1] "y ~ x" > deparse(formula(myglm)[2]) [1] "y()" > deparse(formula(myglm)[3]) [1] "x()" Ideally the elegant method would, in this example, return the character strings "x", "y", and "DAT". Thanks for any insights. Jake Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwegelin@vcu.edu URL: http://www.people.vcu.edu/~jwegelin [[alternative HTML version deleted]]

...h=ntimepts) , time=rep(1:(2*nb+1), length=ntimepts) ) DAT<-rbind(girls,boys) DAT$y<-DAT$y + rnorm(nrow(DAT))/5 DAT$sex<-factor(substring( as.character(DAT[,"id"]), 1,1)) row.names(DAT)<-paste( DAT[,"id"], DAT[,"time"], sep=".") Jacob A. Wegelin Assistant Professor Division of Biostatistics, School of Medicine University of California, Davis One Shields Ave, TB-168 Davis CA 95616-8638 USA http://wegelin.ucdavis.edu/ jawegelin at ucdavis.edu

...1 (2009-12-14) x86_64-apple-darwin9.8.0 locale: [1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tools_2.10.1 > Thanks for any insight Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A.

...s( ~ fac1 + fac2, data=stuff) Error in crosstabs(~fac1 + fac2, data = stuff) : Object "NULL" not found The same command does not fail in S-PLUS. Thanks for any info. The R behavior occurred both in Windows and in unix. In both cases the R version was dated 2002/04/29. Jacob A. Wegelin, Ph.D. University of Washington Department of Statistics B-313 Padelford Box 354322 Seattle WA 98195 http://www.stat.washington.edu/wegelin/ wegelin at stat.washington.edu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.t...

...st figure out whether the variable is continuous or categorical, then perform an operation depending on its type--maybe lm() for continuous but wilcox.test() for dichotomous. But the use of paste(), eval(), and parse() seems awkward. Is there a more elegant way to approach this? Thanks Jacob A. Wegelin Department of Biostatistics Virginia Commonwealth University Richmond VA 23298-0032 U.S.A.

...aphrase Douglas Bates, the above approach is like the diner scene in the movie "Five Easy Pieces". You get an order of toast by first ordering a chicken sandwich and then telling the waitress to hold (that is, to subtract) the meat, lettuce, and mayonnaise. Thanks for any insights Jacob Wegelin

...the correlation matrix. One could print it to a file and then use perl to get the numbers back, but that seems a bit circuitous. What am I overlooking? Of course, this correlation matrix is not the sample correlation matrix of the random effects, as is confirmed below. Thanks for any ideas. Jake Wegelin > junk Linear mixed-effects model fit by REML Data: blockmat Log-restricted-likelihood: -426.3231 Fixed: RESP ~ -1 + intS + tS + intC + tC + intA + tA intS tS intC tC intA tA -0.049391372 0.034116310 -0.003453588 0.200868693 0.01939...

...-8/C/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] lattice_0.17-17 loaded via a namespace (and not attached): [1] Matrix_0.999375-21 boot_1.2-34 grid_2.8.1 lme4_0.999375-28 > Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwegelin@vcu.edu URL: http://www.people.vcu.edu/~jwegelin [[alternative HTML version deleted]]

...its = "secs", class = "difftime") > methods("max") no methods were found > methods("difftime") no methods were found Warning message: function 'difftime' appears not to be generic in: methods("difftime") Thanks for any insights Jake Wegelin

...attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] lme4_0.999375-28 Matrix_0.999375-21 lattice_0.17-17 foreign_0.8-29 loaded via a namespace (and not attached): [1] boot_1.2-34 grid_2.8.1 Thanks for any insights. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwegelin at vcu.edu URL: http://www.people.vcu.edu/~jwegelin

...t;II",] junk$Block<-factor(as.character(junk$Block)) xyplot(yield~nitro, group=Block, data=junk, auto.key=T) What is the elegant or "proper R thinking" way to do this? That is, I want to get a key that only mentions the levels of Block that are used in the plot. Thanks Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwegelin at vcu.edu URL: http://www.people.vcu.edu/~jwegelin

...is. What is the right approach? It would be nice to be able to specify both the aspect ratio and the amount of whitespace outside the plot, and then let the software compute the necessary dimensions of the graphic device that are necessary to accomodate these specs. Thanks for any ideas Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwegelin@vcu.edu URL: http://www.people.vcu.edu/~jwegelin [[alternative HTML version deleted]]