search for: blomberg_at_anu

Dear Lister, I understand merge() can be used to join 2 data frames based on 1 variable. But how about merge based on more than 2 variables? Thank you so much! -- WenSui Liu (http://spaces.msn.com/statcompute/blog) Senior Decision Support Analyst Health Policy and Clinical Effectiveness Cincinnati Children Hospital Medical Center [[alternative HTML version deleted]]

...e pub to celebrate.) -- Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat. Centre for Resource and Environmental Studies The Australian National University Canberra ACT 0200 Australia T: +61 2 6125 7800 email: Simon.Blomberg_at_anu.edu.au F: +61 2 6125 0757 CRICOS Provider # 00120C The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. - John Tukey.

...uot;graphics" "grDevices" "utils" "datasets" [7] "base" -- Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat. Centre for Resource and Environmental Studies The Australian National University Canberra ACT 0200 Australia T: +61 2 6125 7800 email: Simon.Blomberg_at_anu.edu.au F: +61 2 6125 0757 CRICOS Provider # 00120C

I had a question about scan in R. For better code readability, I would like to have lines in the block of data to be scanned that are commented - not just lines that have a comment at the end. For example #age, weight, height 33,128,65 34,56,155 instead of having to do something like 33,128,65 #age, weight, height 34,56,155 Is this at all possible?

I need to specify a model similar to this lme.formula(fixed = sqrt(lbPerAc) ~ y + season + y:season, data = cy, random = ~y | observer/set, correlation = corARMA(q = 6)) except that observer and set are actually crossed instead of nested. observer and set are factors y and lbPerAc are numeric If you know how to do it or have suggestions for reading I will be grateful. eal ps I have

Hello Group, What is the meaning of the error. is there any place to look for this. I guess 'atomic' seems to be OOP related concept. thank you srini

Hello, I am trying to create a stripchart for my data, where y axis labels are characters (ie,names of cities). I would like to change the orientation of the y - axis labels, ie perpendicular to y axis. Below is the code i am using: par(srt=90) with(ozone.ne.trim, stripchart(Median~City,main = "stripchart(ozone)")) The par option doesn't seem to working. Kindly help. Thanks

Does someone have code, or point to a source to get it, to model repeated measures zero-inflated poisson data. The data come from a replicated field trial comparing two treatments - a control and a test treatment. Thanks in advance ------------------------------------ Subhash Chandra, DSc Senior Biometrician Primary Industries Research Victoria Department of Primary Industries 1 Ferguson Road

Hello! I just noticed new link on R wiki on R galleries and wanted to share this info with YOU! - R graphical manuals (this is awesome page as there are all help pages of all packages on CRAN and probably even more and all graphics examples are displayed! - more than 8000 images!) http://bg9.imslab.co.jp/Rhelp/ This is a very nice addition to already existing R graph and movies galleries - R

I am using R in a Windows environment. I have a basic question regarding lm(). I have a dataframe ?data1? with ncol=w. I know that my dependent variable is in column1. Is there a way to write the regression formula so that I can use columns 2 thru w as my independent variables? e.g. something like: ? lm(data1[,1] ~ data1[,2:w] ) ? Thanks

Hi all, Does anybody know whether one can nest IF statement in multiple FOR loops ? According to results of my code, only the first iteration of each loop is performed... Thanks, Regards, Natalie [[alternative HTML version deleted]]

Hello, I've just located the illuminating explanation by Douglas Bates on degrees of freedom in mixed models. The take-home message appears to be: don't trust the p-values from lme. Questions: Should I give up hypothesis testing for fixed effects terms in mixed models? Has my time spent reading Pinheiro & Bates been in vain? Is there a publication on this issue? Thanks, Dan Bebber

Can someone please recommend a novice/beginner's reading list for non-programmers learning R? --------------------------------- 8:00? 8:25? 8:40? Find a flick in no time [[alternative HTML version deleted]]

Hello. I am looking for some help using anova's in RGui. My experiment ie data, has a fixed burning treatment (Factor A) 2 levels, unburnt/burnt. Nested within each level of Factor A are 2 random sites (Factor B). All sites are crossed with a fixed temperature treatment (Factor C) 2 levels, 0 degreesC/2 degreesC, with many replicates of these temperature treatments randomly located at

Hi netters, This is probably a silly question,but I can't find the answer after searching the R-help archives online. ok, I have a matrix. I know there is a "10" somewhere in it. Now I want to know the index of the element "10" in this matrix. That is, if X[i,j]=10, I want to know i and j. Is there a R function to do this? Just like the "find" function in

Dear List, Is there a way to extract variance components from lmeObjects or summary.lme objects without using intervals()? For my purposes I don't need the confidence intervals which I'm obtaining using parametric bootstrap. Thanks, Mike [[alternative HTML version deleted]]

Dear All: I would like to center the data in 'x' by 'group'. The following code scale the data and I have not been able to figure out how to change it so I get the centered data. x <- c(1, 2, 3, 4, 5, 6, 7, 8) group <- c(1,1,1,2,2,2,2,2) unsplit(lapply(split(x,group),scale),group) I would appreciate your help. Ashraf

Hi, How do I get R to output the 95% cutoff from a distribution of a continous variable? summary() only displays a few statistics.... Thanks!

I would appreciate help translating the following lme model to an lmer function. lme(lognrms ~ Group*Rotation*muscle*side*support*arms, random=~1|Subject/Stratum2/rep, data=Data) Many thanks Ross Darnell r.darnell at uq.edu.au

Hi there, Wonder if someone who is R-savvy can help me with the following task (see below) that I occasionally do work and gets quite tedious if I do it manually. Thanks in advance! jenny. ----------------- I have a column of data that looks like this: NA18501 NA18502 NA18504 NA18505 NA18507 NA18508 NA18516 NA18517 NA18522 NA18523 And I want to duplicate the values and sort of