I have been having some trouble with a certain function and getting the message
''subscript out of bounds''. The problem is that in the below
function, the first entry is not defined b/c it is singular and the tstat cannot
be calculated (it is defined as ''NA''). R though will not put
an ''NA'' in tstat[1,1] as I would like it to. Is there any way
I can make this happen and therefore allow R to calculate the next entry in
tstat?
Yi<-c(23,54,67,34,56,78)
x1<-c(1,1,1,1,1,1)
x2<-c(1,0,1,1,0,0)
x<-rbind(x1,x2)
kn<-2 #number of rows in x
tstat<-matrix(0,kn,1)
for (k in 1:kn){
lm1<-lm(Yi~x[k,])
cof<-t(as.matrix(lm1$coefficients))
b5<-summary.lm(lm1)
se<-t(as.matrix(b5$coefficient[,2]))
se1<-se[,2]
cof1<-t(as.matrix(cbind(cof[,2],se1)))
cof2a<-t(cof1)
b1<-as.matrix(cof2a[,1])
se11<-cof2a[,2]
se11<-as.matrix(se11)
tstat[k,]<-b1/se11
}
---------------------------------
[[alternative HTML version deleted]]