[R] fitting a model with the nlme package

[Frank Johannes]

Dear all,
I am analyzing some data that requires a mixed model. I have been
reading Pinheiro and Bates' book,
but cannot find the notation to fit the following model:

Suppose I have the dataset below. Here I am fitting variable p as a
fixed effect, variable h
as a random effect and variable t as nested within h.  I would like to
include the variable j as well
as an independent (non-nested) random effect. I can't find the notation
in the book to tell R
to do this. Does anybody know? 

library(nlme)

h<-c(1,1,1,2,2,2,3,3,3)
j<-c(2,3,8,3,4,3,9,5,4)
t<-c(1,2,3,1,2,3,1,2,3)
f<-c(1,2,1,2,1,2,1,2,1)
p<-c(45,32,35,12,23,12,2,9,12)
set<-data.frame(p,h,j,t)

out<-lme(p~ -1 + f,data=set, random=~1|h/t)

Thanks a alot,
frank.

-- 

                          unladen european swallow