[R], Netlogit
hkimscil
hkimscil at scils.rutgers.edu
Wed May 8 23:01:49 CEST 2002
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list