Displaying 1 result from an estimated 1 matches for "rowx1".
Did you mean:
row1
2005 May 19
1
R 2.1.0 RH Linux Built from Source Segmentation Fault
...s=sum(yrows)/nalt
> if (nrow(X)!=nobs*nalt){print("data dimensions wrong")}
>
> betastore=NULL
> delta=diag(nvar)
> z=read.table("z300.txt",header=TRUE)
> z=as.matrix(z)
> k=ncol(z)
> A1=.1*diag(k)
> nu1=3+nalt
> V=diag(nvar)*nu1
>
> rowx1=rep(0,nphy+1)
> rowx2=rep(0,nphy)
> rowy1=rep(0,nphy+1)
> rowy2=rep(0,nphy)
> rowx1[1]=1
> rowy1[1]=1
> for (i in 1:nphy){
+ # for each physician i draw the relevant X and Y obs using yrows
+ #
+ rowx2[i]=rowx1[i]-1+yrows[i,1]
+ rowy2[i]=rowy1[i]-1+yrows[i,1]/nalt
+ rowx1[i+...