[R] sm library

Oscar orp at intras.es
Wed Dec 20 09:27:30 CET 2000


If you could please try this code:

Y<-cbind(V2,V3)
X<-log(V1)
model1<-glm(Y~X,family=binomial(link=cloglog))
plot(X,model1$fitted.values,ylim=c(0,1))
points(X,V2/V3,col=2)

With data

1 6 23
2 9 24
4 20 25
8 24 24
16 21 22
1 1 23
2 14 25
4 17 24
8 21 23
16 24 24
1 11 24
2 14 24
4 20 26
8 19 23
16 22 23
1 4 23
2 9 21
4 15 24
8 23 23
16 24 24
1 7 23
2 16 26
4 23 25
8 17 18
16 25 25
1 5 22
2 8 24
4 23 25
8 23 23
16 26 26

It doesn't produce the right adjust.

Thanks again for your kindness

Oscar

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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