Dear All,
I am having a real problem exporting data from my simulation. This is it:
y2<-function(n,w)
{
z.s<-matrix(0,2,28)
z.o<-matrix(0,2,28)
z.t<-matrix(0,2,28)
z.f<-matrix(0,2,28)
z.s.w<-matrix(0,2,24)
z.o.w<-matrix(0,2,24)
z.t.w<-matrix(0,2,24)
z.f.w<-matrix(0,2,24)
for(j in 1:2)
{
p.n<-c(1:n) #Patient no
ty<-rep(p.n,4) #sorting variable
a.r<-c(0,cumsum(rexp(n-1,rate=70.6903164))) #recruitment times
a<-rep(a.r,4)
x.1<-matrix(0,n/w,w)
for(i in 1:n/w) x.1[i, ]<-sample(1:w,replace=F)
x.2<-t(x.1)
x.3.2<-as.numeric(x.2 > w/2) #block randomisation
x.3<-rep(x.3.2,4)
e2<-as.vector(rep(0.25,n))
e3<-as.vector(rep(0.5,n))
e4<-as.vector(rep(1,n))
e5<-as.vector(rep(2,n))
e.5<-as.vector(matrix(c(e2,e3,e4,e5),ncol=4,byrow=F))
abs.t<-as.vector(e.5+a)
r.e<-rnorm(4*n,mean=0,sd=1) #random error
ran.eff<-rnorm(n,mean=0,sd=1) #esimates of random effects
ran<-rep(ran.eff,4)
b0<-0.1931404
b1<--0.02629946
b2<-0.1995204
b3<--0.03632164
y<-as.vector(b0+b1*a+b2*e.5+b3*x.3+ran+r.e) #generated responce
assign( "df", data.frame(ty=ty,a=a,abs.t=abs.t,x.3=x.3,e.5=e.5,y=y), 0 )
DF<-as.matrix(df)
i<-sort.list(ty); ty<-ty[i]; DF<- DF[i,]
sf<-as.data.frame(DF)
p.n.s<-as.vector(df[,1][abs.t<0.5])
a.s<-as.vector(df[,2][abs.t<0.5])
abs.t.s<-as.vector(df[,3][abs.t<0.5])
x.3.s<-as.vector(df[,4][abs.t<0.5])
e.5.s<-as.vector(df[,5][abs.t<0.5])
y.s<-as.vector(df[,6][abs.t<0.5])
assign( "df.s",
data.frame(p.n.s=p.n.s,a.s=a.s,abs.t.s=abs.t.s,x.3.s=x.3.s,e.5.s=e.5.s,y.s=y.s),
0 )
DF.s<-as.matrix(df.s)
i<-sort.list(p.n.s); p.n.s<-p.n.s[i]; DF.s<- DF.s[i,]
sf.s<-as.data.frame(DF.s)
f.s <-try(lme(fixed=y.s~a.s+e.5.s+x.3.s, est.method =
"ML",cluster=~p.n.s,random=~1,
data=sf.s,na.action = na.omit), first=T)
f.s.w<-try(lme(fixed=y.s~e.5.s+x.3.s, est.method =
"ML",cluster=~p.n.s,random=~1,
data=sf.s,na.action = na.omit), first=T)
f.s.s<-summary(f.s)
f.s.w.s<-summary(f.s.w)
p.n.o<-as.vector(df[,1][abs.t<1])
a.o<-as.vector(df[,2][abs.t<1])
abs.t.o<-as.vector(df[,3][abs.t<1])
x.3.o<-as.vector(df[,4][abs.t<1])
e.5.o<-as.vector(df[,5][abs.t<1])
y.o<-as.vector(df[,6][abs.t<1])
assign( "df.o",
data.frame(p.n.o=p.n.o,a.o=a.o,abs.t.o=abs.t.o,x.3.o=x.3.o,e.5.o=e.5.o,y.o=y.o),
0 )
DF.o<-as.matrix(df.o)
i<-sort.list(p.n.o); p.n.o<-p.n.o[i]; DF.o<- DF.o[i,]
sf.o<-as.data.frame(DF.o)
dim(sf.o)
f.o<-lme(fixed=y.o~a.o+e.5.o+x.3.o, est.method =
"ML",cluster=~p.n.o,random=~1,
data=sf.o,na.action = na.omit)#laird and ware model
f.o.w<-lme(fixed=y.o~e.5.o+x.3.o, est.method =
"ML",cluster=~p.n.o,random=~1,
data=sf.o,na.action = na.omit)#laird and ware model without
f.o.s<-summary(f.o)
f.o.w.s<-summary(f.o.w)
p.n.t<-as.vector(df[,1][abs.t<2])
a.t<-as.vector(df[,2][abs.t<2])
abs.t.t<-as.vector(df[,3][abs.t<2])
x.3.t<-as.vector(df[,4][abs.t<2])
e.5.t<-as.vector(df[,5][abs.t<2])
y.t<-as.vector(df[,6][abs.t<2])
assign( "df.t",
data.frame(p.n.t=p.n.t,a.t=a.t,abs.t.t=abs.t.t,x.3.t=x.3.t,e.5.t=e.5.t,y.t=y.t),
0 )
DF.t<-as.matrix(df.t)
i<-sort.list(p.n.t); p.n.t<-p.n.t[i]; DF.t<- DF.t[i,]
sf.t<-as.data.frame(DF.t)
xd<-dim(sf.s)
f.t<-lme(fixed=y.t~a.t+e.5.t+x.3.t, est.method =
"ML",cluster=~p.n.t,random=~1,
data=sf.t,na.action = na.omit)#laird and ware model
f.t.w<-lme(fixed=y.t~e.5.t+x.3.t, est.method =
"ML",cluster=~p.n.t,random=~1,
data=sf.t,na.action = na.omit)#laird and ware model
f.t.s<-summary(f.t)
f.t.w.s<-summary(f.t.w)
p.n.f<-as.vector(df[,1][abs.t<4])
a.f<-as.vector(df[,2][abs.t<4])
abs.t.f<-as.vector(df[,3][abs.t<4])
x.3.f<-as.vector(df[,4][abs.t<4])
e.5.f<-as.vector(df[,5][abs.t<4])
y.f<-as.vector(df[,6][abs.t<4])
assign( "df.f",
data.frame(p.n.f=p.n.f,a.f=a.f,abs.t.f=abs.t.f,x.3.f=x.3.f,e.5.f=e.5.f,y.f=y.f),
0 )
DF.f<-as.matrix(df.f)
i<-sort.list(p.n.f); p.n.f<-p.n.f[i]; DF.f<- DF.f[i,]
sf.f<-as.data.frame(DF.f)
f.f<-lme(fixed=y.f~a.f+e.5.f+x.3.f, est.method =
"ML",cluster=~p.n.f,random=~1,
data=sf.f,na.action = na.omit)#laird and ware model
f.f.w<-lme(fixed=y.f~e.5.f+x.3.f, est.method =
"ML",cluster=~p.n.f,random=~1,
data=sf.f,na.action = na.omit)#laird and ware model
f.f.s<-summary(f.f)
f.f.w.s<-summary(f.f.w)
f.s.a<-try(anova(f.s))
f.s.w.a<-try(anova(f.s.w))
f.o.a<-anova(f.o)
f.o.w.a<-anova(f.o.w)
f.t.a<-anova(f.t)
f.t.w.a<-anova(f.t.w)
f.f.a<-anova(f.f)
f.f.w.a<-anova(f.f.w)
qe.s<-as.double(c(f.s$coeff$fixed,f.s$var.ran,f.s$loglik,f.s$sigma,
f.s.s$nobs,f.s.s$nclus,f.s.s$AIC,f.s.s$BIC,f.s.s$residuals,f.s.a$Value,f.s.a$Std.Error,f.s.a$"z
ratio"))
qe.o<-as.double(c(f.o$coeff$fixed,f.o$var.ran,f.o$loglik,f.o$sigma,
f.o.s$nobs,f.o.s$nclus,f.o.s$AIC,f.o.s$BIC,f.o.s$residuals,f.o.a$Value,f.o.a$Std.Error,f.o.a$"z
ratio"))
qe.t<-as.double(c(f.t$coeff$fixed,f.t$var.ran,f.t$loglik,f.t$sigma,
f.t.s$nobs,f.t.s$nclus,f.t.s$AIC,f.t.s$BIC,f.t.s$residuals,f.t.a$Value,f.t.a$Std.Error,f.t.a$"z
ratio" ))
qe.f<-as.double(c(f.f$coeff$fixed,f.f$var.ran,f.f$loglik,f.f$sigma,
f.f.s$nobs,f.f.s$nclus,f.f.s$AIC,f.f.s$BIC,f.f.s$residuals,f.f.a$Value,f.f.a$Std.Error,f.f.a$"z
ratio"))
qe.s.w<-as.double(c(f.s.w$coeff$fixed,f.s.w$var.ran,f.s.w$loglik,f.s.w$sigma,
f.s.w.s$nobs,f.s.w.s$nclus,f.s.w.s$AIC,f.s.w.s$BIC,f.s.w.s$residuals,f.s.w.a$Value,f.s.w.a$Std.Error,f.s.w.a$"z
ratio"))
qe.o.w<-as.double(c(f.o.w$coeff$fixed,f.o.w$var.ran,f.o.w$loglik,f.o.w$sigma,
f.o.w.s$nobs,f.o.w.s$nclus,f.o.w.s$AIC,f.o.w.s$BIC,f.o.w.s$residuals,f.o.w.a$Value,f.o.w.a$Std.Error,f.o.w.a$"z
ratio"))
qe.t.w<-as.double(c(f.t.w$coeff$fixed,f.t.w$var.ran,f.t.w$loglik,f.t.w$sigma,
f.t.w.s$nobs,f.t.w.s$nclus,f.t.w.s$AIC,f.t.w.s$BIC,f.t.w.s$residuals,f.t.w.a$Value,f.t.w.a$Std.Error,f.t.w.a$"z
ratio"))
qe.f.w<-as.double(c(f.f.w$coeff$fixed,f.f.w$var.ran,f.f.w$loglik,f.f.w$sigma,
f.f.w.s$nobs,f.f.w.s$nclus,f.f.w.s$AIC,f.f.w.s$BIC,f.f.w.s$residuals,f.f.w.a$Value,f.f.w.a$Std.Error,f.f.w.a$"z
ratio"))
z.s[j,]<-c(qe.s)
z.o[j,]<-c(qe.o)
z.t[j,]<-c(qe.t)
z.f[j,]<-c(qe.f)
z.s.w[j,]<-c(qe.s.w)
z.o.w[j,]<-c(qe.o.w)
z.t.w[j,]<-c(qe.t.w)
z.f.w[j,]<-c(qe.f.w)
}
assign( "s.m", data.frame(z.s=z.s),frame=0,immediate=T )
assign( "o.y", data.frame(z.o=z.o),frame=0,immediate=T )
assign( "t.y", data.frame(z.t=z.t),frame=0,immediate=T )
assign( "f.y", data.frame(z.f=z.f),frame=0,immediate=T )
assign( "s.m.w",
data.frame(z.s.w=z.s.w),frame=0,immediate=T )
assign( "o.y.w",
data.frame(z.o.w=z.o.w),frame=0,immediate=T )
assign( "t.y.w",
data.frame(z.t.w=z.t.w),frame=0,immediate=T )
assign( "z.f.w",
data.frame(z.f.w=z.f.w),frame=0,immediate=T )
export.data(DataSet="z.f.w",ColumnNames=T,FileName="C:\\simulation\\four
years without.sdd",
FileType="SPLUS")
}
When I run it I get this error:
ERROR:Cannot export DataSet = "z.f.w". The data must be committed to the
database first such as calling assign(...,immediate = T).
As far as I am aware frame=0 is the session frame and so why is the data not
commited to the database? Do you know? I would really appreciate any
guidence you can give me.
Thanx,
Paul.
_________________________________________________________________________
Get Your Private, Free E-mail from MSN Hotmail at http://www.hotmail.com.
|