search for: revorder

Greetings R-helpers, I am attempting to fit an lme() while specifying a correlation structure, but I'm getting into trouble long before I get to that point. I am receiving the error: Error in y[revOrder] - Fitted : non-conformable arrays It doesn't seem to matter how simple or complex the model I specify is, it always gives this same error message. This makes me suspect something is wrong with my data set, but I'm not sure where to begin looking. I can run through the Pinheiro & Bate...

...05440 1 1 5 18 2.722942793 2 1 5 19 4.395861278 3 1 5 20 4.849296475 4 1 5 21 3.049616421 1 1 6 22 2.874405962 2 1 6 23 4.359511097 3 1 6 24 6.165419699 4 1 6 This happened: > testLME <- lme(Y~t+D,data=simpledata,random=~1|ID) > predict(testLME, simpledata[1,]) Error in val[revOrder, level + 1] : incorrect number of dimensions This has occurred with other datasets as well. Is this a bug in the code, or am I doing something wrong? (Also, is there a way to parse a formula of a type given to "random"? For example, given ~1+t|ID, I'd like to be able to extract al...

Hi all, I tried to do a variance components using nlme, but I got the following warning mesage ##################################################### not meaningful for factors in: Ops.factor(y[revOrder], Fitted) ###################################################### Can someone point out what is the meaning of this warning message? I tried to look at Ops.factor, but I don't understand it since I'm relatively new to R. Thanks. Yen Lin [[alternative HTML version deleted]]

I'm attempting to run an LME model with random intercepts for location and have received an error message. The model and error message are given below: > mixed.distance <- lme (fixed= log(min5inpnc) ~ i93centerline, random=~1|location, data = nosmokecomp) Error in y[revOrder] - Fitted : non-conformable arrays I believe that my dataset may have dimensions that are not working with the model, but I'm not sure how to diagnose the exact problem. Can someone provide assistance? Thank you, Christina [[alternative HTML version deleted]]

I want to take the first row of each unique ID value from a data frame. For instance > ddTable <- data.frame(Id=c(1,1,2,2),name=c("Paul","Joe","Bob","Larry")) I want a dataset that is Id Name 1 Paul 2 Bob > unique(ddTable) Will give me all 4 rows, and > unique(ddTable$Id) Will give me c(1,2), but not accompanied by the name column.