similar to: multiple graphs

Hi. I'm a student at SFU in Canada. The basic thing I want to do is calculate means of different strata. I have 2 vectors. One has the values I want to take the means from, the other is the four strata I am interested in. So I essentially want to break up the information vector into the four strata and calculate four means, one for each stratum. How can I do this in a reasonable way? Thanks

Is there a way to fit mixed proportional odds models in R? -- Using Opera's revolutionary e-mail client: http://www.opera.com/m2/

Dear All, I have small samples of data (between 6 and 15) for numerious time series points. I am assuming the data for each time point is normally distributed. The problem is that the data arrvies sporadically and I would like to detect the number of outliers after I have six data points for any time period. Essentially, I would like to detect the number of outliers when I have 6 data points then

Hello all, Is AIC calculated incorrectly in lmer? It appears as though it uses AIC = -2*logLik - 2*#parms, instead of -2*LogLik + 2*#parms? Below is output from one of many models I have tried: Generalized linear mixed model fit using PQL Formula: cswa ~ pcov.ess1k + (1 | year) Data: ptct50.5 Family: poisson(log link) AIC BIC logLik deviance 224.8466 219.19 -114.4233 228.8466

Dear list, I would like to fit a negative binomial regression model as described in "Byers AL, Allore H, Gill TM, Peduzzi PN., Application of negative binomial modeling for discrete outcomes: a case study in aging research. J Clin Epidemiol. 2003 Jun;56(6):559-64" to my data in which the response is count data. There are also 10 predictors that are count data, and I have also 3

Hi, I've an issue adding an interaction to a GAMM: My model was of form: gamm1 <- gamm(TOTSR ~ fROT + s(PH) + s(LOI) + s(ASP) + s(SQRT_ELEV) + CANCOV + s(SQRT_TOTCWD) + s(WELLF) + s(WELLN) + s(OLDWDLD) + s(DISTWOOD) + s(Annprec) + s(OLDWDLD:DISTWOOD) + (1|fSITE), family = poisson, data = BIOFOR2) with interaction of s(OLDWDLD:DISTWOOD). However I got this error message below but

dear all, This is a probably a silly question. If I type > grep("x",c("a.x" ,"b.x","a.xx"),value=TRUE) [1] "a.x" "b.x" "a.xx" Instead, I would like to obtain only "a.x" "b.x" How is it possible to get this result with grep()? many thanks for your attention, best, vito --

Dear all, It appears that MASS::polr() and Design::lrm() return the same point estimates but different st.errs when fitting proportional odds models, grade<-c(4,4,2,4,3,2,3,1,3,3,2,2,3,3,2,4,2,4,5,2,1,4,1,2,5,3,4,2,2,1) score<-c(525,533,545,582,581,576,572,609,559,543,576,525,574,582,574,471,595, 557,557,584,599,517,649,584,463,591,488,563,553,549) library(MASS) library(Design)

With version 8 of acrobat reader, it is now possible to have 3D in PDf documents. Does it exist already an R package who manage to produce 3D plots which can be saved as interactive 3D graphs in a PDF file? Best Regards Bruno Cavestro ------------------------------------------------------ Leggi GRATIS le tue mail con il telefonino i-mode? di Wind http://i-mode.wind.it/

Hi all, Could anybody recommend some easy-to-understand and example based notes/tutorials on how to use cubic splines to do filtering on multi-dimension data? I am confused by the 1-dimensional case, and more confused by multi-dimensional case. I found all the books suddenly become very abstract when it comes to this subject. They don't provide examples in R or Splus at all. Specifically,

I have estimated a multiple nonparametric regression using the loess command in R. I have also estimated an additive version of the model using the gam function. Is there a way of using the output of these two models to test the restrictions imposed by the additive model?

Dear all, there are R packages able to simulate or estimate bilinear model for time series? I know it is an open problem, but do exist something for very simplified bilinear models? Alternatively, what kinfd of non linear time series models are performed in R? If R is not able, could someone suggest me for some commercial softwares to deal with bilinear models? i'm afraid of a negative

Dear all, I have tried to consider a correlation structure in lme (package lme4), but without success. I have used something like: > risul<-lme(y~x+ z , data=mydata, random=~ x | g, correlation = corAR1()) but the result is the same as: > risul<-lme(y~x+ z , data=mydata, random=~ x | g). Can anybody help me? Antonella ************************************************** Prof.

Dear all, Is there *any* reason explaining what I describe below? I have the following line myfun(x) If I type them directly in R (or copy/past), it works.. However if I type in R 2.6.1 > source("code.R") ##code.R includes the above line Error in inherits(x, "data.frame") : object "d" not found namely myfun() does not work correctly. In particular the

dear all, Could anyone explain me the behaviour of median() within by()? (I am running R.2.7.0) thanks, vito > H<-cbind(rep(0:1,l=20),matrix(rnorm(20*2),20,2)) > by(H[,-1],H[,1],mean) INDICES: 0 V1 V2 -0.2101069 0.2954377 --------------------------------------------------------------------------------------------------------------------- INDICES: 1 V1

dear all, When I use all.vars(), I am interest in extracting only the variable names.. Here a simple example all.vars(as.formula(y~poly(x,k)+z)) returns [1] "y" "x" "k" "z" and I would like to obtain "y" "x" "z" Where is the trick? many thanks vito -- ==================================== Vito M.R. Muggeo Dip.to Sc

dear all, I would like to get the lme call without fitting the relevant model. library(nlme) data(Orthodont) fm1 <- lme(distance ~ age, random=list(Subject=~age),data = Orthodont) To get fm1$call without fitting the model I use call(): my.cc<-call("lme.formula", fixed= distance ~ age, random = list(Subject = ~age)) However the two calls are not the same (apart from the data

Hi, When using lapply (or sapply) to loop over a list, can I somehow access the index of the list from inside the function? A trivial example: df1 <- split( x=rnorm(n=100, sd=seq(from=1, to=10, each=10)), f=letters[seq(from=1, to=10, each=10)] ) str(df1) #List of 10 # $ a: num [1:10] -0.801 0.418 1.451 -0.554 -0.578 ... # $ b: num [1:10] -2.464 0.279 4.099 -2.483 1.921 ... # $ c: num

Hi everybody, I have the following problem. I have a vector containing character elements, such as: list = c("aa","bb","cc","dd","ee") I want to create an index which identifies the elements that are different from, e.g. "aa" and "bb". When I do the following: jj = list!="aa" & list!="bb" > jj

Dear all, I am dealing with the following (apparently simple problem): For some reasons I am interested in passing variables from a dataframe to a specific environment, and in fitting a standard glm: dati<-data.frame(y=rnorm(10),x1=runif(10),x2=runif(10)) KK<-new.env() for(i in 1:ncol(dati)) assign(names(dati[i]),dati[[i]],envir=KK) #Now the following two lines work correctly: