search for: torder

...odont, 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,pos=c(1,1)) # random coefficient for age carma(data,ccov=~sex,torder=1,pre=rep(1,2),pos=rbind(c(1,1),c(2,2)))...

Dear all, I am working with taxonomic data, represented as a list of classes, orders, families, genera and finally species. > class(mydata) [1] "data.frame" > mode(mydata) [1] "list" > names(mydata) [1] "tclass" "torder" "tfamily" "tgenus" "tspecies" > length(mydata$tclass) [1] 161590 The first 10 rows look like the following: > mydata[1:10,] tclass torder tfamily tgenus 1 Chlorophyta Chlorophyceae Dunaliellaceae Collodictyon 2...

...f 11 variables: $ site : Factor w/ 24 levels "B(W)","BC-1",..: 1 1 2 2 2 1 1 1 ... $ sampdate : Date, format: "2000-07-18" "2000-07-18" ... $ tclass : Factor w/ 8 levels "Annelida","Arachnida",..: 1 5 5 5 5 ... $ torder : Factor w/ 18 levels "Achtinedida",..: 13 5 5 10 13 5 7 ... $ tfamily : Factor w/ 81 levels "","Ameletidae",..: 79 46 46 14 42 ... $ tgenus : Factor w/ 206 levels "","Acentrella",..: 1 10 10 1 140 ... $ tspecies : Factor...

...0 rows some of which would be empty, but instead I get a dataframe with 995 rows. In this case missing 5 rows were omitted. Any suggestions? Thanks, Olga ---------------- result<-list() for(i in 1:1000)try({ hclass<-mydata2$tclass[mydata2$tspecies%in%com[[i]][[5]][[1]]] horder<-mydata2$torder[mydata2$tspecies%in%com[[i]][[5]][[1]]] hfamily<-mydata2$tfamily[mydata2$tspecies%in%com[[i]][[5]][[1]]] hgenus<-mydata2$tgenus[mydata2$tspecies%in%com[[i]][[5]][[1]]] hspecies<-mydata2$tspecies[mydata2$tspecies%in%com[[i]][[5]][[1]]] hier<-cbind(hspecies,hgenus,hfamily,horder,hclass) s...

An executable version 0.6 of my libraries is now available at www.luc.ac.be/~jlindsey/rcode.html This works with MS R0.64.2 and appears possibly to work with R65.0. There is a serious problem with the Fortran compiler as some of the examples for elliptic and carma crash it. These same examples do not crash R63.0 with the library executables of Jan 99. I am releasing this anyway because of the