Logistisk regression
Logistisk regression
Logistisk regression
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
u n i v e r s i t y o f c o p e n h a g e n<br />
Prædiktion<br />
d e p a r t m e n t o f b i o s t a t i s t i c s<br />
Vi ønsker at sammenligne sandsynligheden for CHD for en<br />
ikke-rygende mand på 50 år med sandsynligheden for en ikke-rygende<br />
kvinde på 50 år som funktion af SBP.<br />
Prædiktion i R<br />
> Mdata head(Mdata,3)<br />
mand ryger age sbp<br />
1 1 0 50 90<br />
2 1 0 50 91<br />
3 1 0 50 92<br />
> pM head(pM)<br />
1 2 3 4 5 6<br />
0.06037988 0.06228133 0.06421687 0.06618649 0.06819014 0.07022779<br />
><br />
> Kdata pK<br />
> plot(90:300,pM,xlab="SBP",ylab="Sandsynlighed for CHD",ylim=c(0,1),col="blue",type="l")<br />
> lines(90:300,pK,col="red")<br />
> legend(100,1,c("Mand, ikke-ryger, 50 aar","Kvinde, ikke-ryger, 50 aar"),<br />
> col=c("blue","red"),lty=1)<br />
><br />
41 / 42