search for: rmna

...s <- c(ctimes,list(as.matrix(Soybean[4][Soybean[1]==i]))) plot <- c(plot,Soybean[1][Soybean[1]==i][1]) variety <- c(variety,Soybean[2][Soybean[1]==i][1]) year <- c(year,Soybean[3][Soybean[1]==i][1])} # unfortunately, category names of factor variables were lost here by c() data <- rmna(restovec(soybean), ccov=tcctomat(data.frame(plot,variety=variety-1, year=as.factor(year)),dataframe=F)) # for discrete time AR(1) #data <- rmna(restovec(soybean), # ccov=tcctomat(data.frame(plot,variety=variety-1, # year=as.factor(year)),dataframe=F),tvc=tvctomat(ctimes,name="ctimes"...

...distance ~ age, data = Orthodont)) # random is ~ age summary(fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1)) # set up data object distance <- matrix(Orthodont[1],ncol=4,byrow=T) age <- matrix(Orthodont[2],ncol=4,byrow=T) sex <- Orthodont[4][seq(4,108,by=4)] data <- rmna(restovec(distance,times=age),ccov=tcctomat(sex)) rm(Orthodont,age,sex,distance) #model fm1 carma(data,torder=1,pre=rep(1,3),pos=rbind(c(1,1),c(2,2),c(1,2))) # uncorrelated random coefficients carma(data,torder=1,pre=rep(1,2),pos=rbind(c(1,1),c(2,2))) #model fm2 carma(data,ccov=~sex,torder=1,pre=1,...