search for: zad

...1<-a/(2-a) > q2<-(2*a^2)/(1+a) > se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na)) > se > } > > cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){ > se1<-seROC(AUC1,na1,nn1) > se2<-seROC(AUC2,na2,nn2) > > sed<-sqrt(se1^2+se2^2-2*r*se1*se2) > zad<-(AUC1-AUC2)/sed > p<-dnorm(zad) > a<-list(zad,p) > a Could somebody kindly suggest me how to calculate the value of "r" or some ways to calculate the statistical significance measure for the differences of auc for a part of the curve like AUC0.1. Thank You -- Dukka K...

...30:37 EST seROC<-function(AUC,na,nn){ a<-AUC q1<-a/(2-a) q2<-(2*a^2)/(1+a) se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na)) se } cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){ se1<-seROC(AUC1,na1,nn1) se2<-seROC(AUC2,na2,nn2) sed<-sqrt(se1^2+se2^2-2*r*se1*se2) zad<-(AUC1-AUC2)/sed p<-dnorm(zad) a<-list(zad,p) a } The author of this script says: "The first function (seROC) calculate the standard error of ROC curve, the second function (cROC) compare ROC curves." What do you think of this script? Is there any function to do it better in R...

...30:37 EST seROC<-function(AUC,na,nn){ a<-AUC q1<-a/(2-a) q2<-(2*a^2)/(1+a) se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na)) se } cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){ se1<-seROC(AUC1,na1,nn1) se2<-seROC(AUC2,na2,nn2) sed<-sqrt(se1^2+se2^2-2*r*se1*se2) zad<-(AUC1-AUC2)/sed p<-dnorm(zad) a<-list(zad,p) a } The author of this script says: "The first function (seROC) calculate the standard error of ROC curve, the second function (cROC) compare ROC curves." What do you think of this script? Is there any function to do it better in R...

...1<-a/(2-a) > q2<-(2*a^2)/(1+a) > se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na)) > se > } > > cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){ > se1<-seROC(AUC1,na1,nn1) > se2<-seROC(AUC2,na2,nn2) > > sed<-sqrt(se1^2+se2^2-2*r*se1*se2) > zad<-(AUC1-AUC2)/sed > p<-dnorm(zad) > a<-list(zad,p) > a > } > -- Waverley @ Palo Alto

Dear R-helper, I would like to compare the AUC of two logistic regression models (same population). Is it possible with R ? Thank you Roman Rouzier [[alternative HTML version deleted]]

...command to fit(plot) a regression line but I haven't found any command on the html.help. My question is if I am using a regression for instance reg<-lm(y~x) and I would like to make a regression line, to see how well my data is fit. Can you please help me with this?? Yours Sincerely Soheila Zad -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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....

...consonance densities and distributions. And it can utilize bootstrapping (the BCa and t method) to approximate the curves, distributions, and densities. And to top it off, it can export tables of relevant statistics in various formats. We welcome any feedback and bug reports ( https://github.com/zadchow/concurve/issues) Best, Zad R. Chow Andrew D. Vigotsky [[alternative HTML version deleted]] _______________________________________________ R-packages mailing list R-packages at r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

...consonance densities and distributions. And it can utilize bootstrapping (the BCa and t method) to approximate the curves, distributions, and densities. And to top it off, it can export tables of relevant statistics in various formats. We welcome any feedback and bug reports ( https://github.com/zadchow/concurve/issues) Best, Zad R. Chow Andrew D. Vigotsky [[alternative HTML version deleted]] _______________________________________________ R-packages mailing list R-packages at r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

I am a novice user of R and am stumbling on how to order a dataset produced during my session. I have a 1863 row X 14 column dataset that I want to put out to a file. I want the output sorted by the first column and then by the second column both in ascending order. The first column is character and the second is numeric (I hope). I used an "as.numeric" function to assign that

...gp6sLLE8 > dywsGoSAsrqjoJa4UaIouaO5FrOIyMi1lYb519jRqMX72EGi5sLBi0rR0oEY1P4U9Of1ucXd > 4WI+IvYv7PVkVFAF7SKGHwIUJMGhvRuJAPX+a+nQlEYASLrscbpwsBYHKSlMHls7KRlHIQhW > Vldfo1tYLTaHTEwHUWklSwLJRBRQQEI6W103XwFbUlZRRVVLluciSyleIlRQV72UMsAAoZp2 > 5HDhdJsCIat4ek7wQGal6YhiQPcCbtjHYybZYLbVgyp8ACjlg4NeZaDR14AC4AU6YqcKoFnL > eQUAyNkj5ubfNbykMuzZqOSDdBTGkK8cpG6Dq0f21be6m6uBVmJhiPGXpn6w/usAv1W4lZlG > RYQAzf5ZDMhIT2QYTSJ+Sn9AeCpDLWEwbv4N2Itk63YIM9qjoKQjoxz/pRuqplYvjEwETeSP > iR+MkNT1zlr49AhKATWDR376cDtzKEzH2hj+M6rS13oz4A4R+GCuAfTkIbWyaNRnJVTYL6dY > WGgzzaIJ++oxMgDhjezLyUdv1Jn7KrUkOBvDd1W2sDE...

...`JI2c zk1&Qm^R_;J;{#`zImgwdflv|%NFnJc!BL~RaZvK~Ud_#uhU2AWQ9BwT;w(@JkO3h> zbl+0qWwbzH1%dc?Lih+`0B<6=TFbR9`Beu-3Y36EdJ<n at x*c57V=<kPBncNToV&}* z+!u+Bb%Dk^QVW}#hEIOI;a5L$ol3%co+$X0KYJEyVLbC70t$tZiW`%fKe{sJslAfJ z<C4lNlul8?p|l`T3aJHNfNWWDUWC|A6QF~@D*-{^ga_Ioued|daxP}2VKo32apq#A zAd!MXdx}n%I@{oA|0Y|*4bGjr at akR3=TQw#xB6b`w%^^UTW|zdA%$aiV)?aCe}kiC zpSM3+ at GF1zBDHrMR06F^3X^baB>Ci5_xQrIOU84JNdhWKD3kykK`2B3Ap&0dn5}x% zS_dr&LtqgIABe^(2_hz1L4a?+yp$3bLX3ej#1^4<NTra<v#Bhyzd>dSv`X1Ie~!(; z5b${aDp_N=8^jRY3u5sr at j33cMDJqz$_f6D5B~{YSw6z+FKqMIKDdd_72}blClz&o zf...

...lt;z8k+kuq7g6B@4me$ckjyA2`230xDD?!{!p4t={)I44;E~` zc++?+vEXAQFz%TC-kMte;w&Lg(euf_wHeyT^{<E>LPpvtD)qRyk+H60th8nP;^GQ8 z%^?&4Y&iZcH?=gK?Q(dhzou_~mTRF&5w<!gxDU{|33Y=#`O=2k93x$Sxf#-KZP;F& zk9#WD8mW!QIbCl3-g4yVtcBKu^oh>%yJ6mt37S9DY0yU@V`3gSGNojC!U7Fh1x_{I zeOs_sDyhcsC@F6y=Zad!WWtW6T$R?<Ici~3_0$Rk4$g^<{keuN3;Yh2jeJ&eFlFtp z6vsY9s9=!<&plSFyV4fzD3_0V29o6E2kDtIv(fDOdCLW_b{)7{-wbTaOL(jpJ+%Ff zX5!@~o#c1Ay4OoNo?FgYlX!kB>KgMeEcOlIAOFLgnmg5D2swsmxQHNnI$bcOy#bAI zNFeh+1@cCckToEo{WBLl@YD*dp`iUu7sFpb8^ZMjTEzUXBM8E9AA;HNNVeNLCY=f_ z=iPbkgaZ+poaNKQ{Q!ddC)oa$_PrXaB...

...`%eUvgP5J%mCIabj3&;<jxwuEhJ9R%cyjt}<ZM`9Uh#tgzpCovC?Ns^7 z!|N?Ge%v5g^(!Xj%o#vhIS*|8t$x{8`*zIr+90{r{z%EAYnN7lT-B&P;+TsL*WjWO zv6`Hxp1N~L at ld&AYtsqo>A66;|0XzDJs|qs{B^T(hZOG>ms|3RDgc;U;E>$gbpLw7 zs;X~h<z{W%>toKp(7U<66)3+ at 5q+g!_H=ChX6nwYjYET0F1WDiY6O^o-?8;$c7JrQ zad5ATl~1aJ59AJxYMO6LYZ%p-B#!M1a=Yl#(+5{RUp$~gB#>Rt0FGHH at 8uncxbr*k zRmB3o?$&})(`EyST(;=`iu!p(ZmgeR12MMCbP&1GAtX9{Z}IvANT0)zzVke|p*dvT zY7jXBfM(seVa<1eFMe&m|L9KcW5KZUME6HylIHwJcM}YUU%gCKc;`;21MFCrM2r<q ztJ|yay}zPUe33cxt<F1u?}i<~<z=6BgC!G=O?|6c;T`A at I0UcOf0G#n#JBdiuO9cZ...