[R], Netlogit

[hkimscil]

Hello -

I have a problem of getting output from a logit test (network). The command, "summary" would  not finish its report; because of an obscure error...

I have a stack of network graphs, which are arrayed as follows
== description starts here.
#### graphs structure
> rrel<-array(dim=c(4,69,69))
> rrel[1,,]<-aspRel
> rrel[2,,]<-aspArea
> rrel[3,,]<-aspCrowding
> rrel[4,,]<-aspSharing
### Then do logit test
rrel.logit.qap<-netlogit(rrel[1,,],rrel[c(8,3,4),,],nullhyp="qap",reps=10)
### check the result
summary(rrel.logit.qap)

The summary output is terminated with an error sign as follows

coefficients stuffs ........
Goodness of fitness .........
chi-square ............
         Actual
Predicted      0     1
        0   4564   128
Error in cat("\n        Total Fraction Correct:", (n$ctable[1, 1] + n$ctable[2,  : 
        subscript out of bounds
==== description ends here.

I think it has a problem in rendering the second prediction line? I'd appreciate any help!

-hyo
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._