[R] sm library

[BXC (Bendix Carstensen)]

If you instead use

points(X,V2/(V2+V3),col=2)

things make more sense. Cf. the help page for glm.

Bendix
----------------------
Bendix Carstensen
Senior Statistician
Steno Diabetes Centre
Niels Steensens Vej 2
DK-2820 Gentofte
Denmark
tel: +45 44 43 87 38
mob: +45 28 25 87 38
fax: +45 44 43 73 13
bxc at novo.dk
www.biostat.ku.dk/~bxc
----------------------



-----Original Message-----
From: Oscar [mailto:orp at intras.es]
Sent: 20. december 2000 09:28
To: r-help at lists.R-project.org
Subject: Re: Re: Re: [R] sm library


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