Federico Bonofiglio
2011-Feb-08 21:14 UTC
[R] intervals {nlme} lower CI greater than upper CI !!!????
Hi folks... check this out..> GLU<-lme(gluc~rt*cd4+sex+age+rf+nadir+pharmac+factor(hcv)+factor(hbs)++ haartd+hivdur+factor(arv), + random= ~rt|id, na.action=na.omit)> intervals(GLU)$fixedlower est. upper (Intercept) 67.3467070345 7.362307e+01 7.989944e+01 rt *0.0148050160* 6.249304e-02 1.101811e-01 cd4 -0.0032133709 -6.566501e-04 1.900071e-03 sexM 0.9601467541 3.565606e+00 6.171066e+00 age *0.2436472425* 3.420025e-01 4.403578e-01 rfOTHER -5.0933264232 -1.198439e+00 2.696449e+00 rfETEROSESSUALE -2.1085096013 1.974725e+00 6.057960e+00 rfOMOSESSUALE -5.8156940466 -1.891870e+00 2.031953e+00 nadir -0.0063302516 8.604450e-04 8.051142e-03 pharmacNRTI -4.2586487617 -1.252531e+00 1.753586e+00 pharmacPI -3.3132784651 -8.829819e-01 1.547315e+00 factor(hcv)1 -0.4338530808 2.374438e+00 5.182729e+00 factor(hbs)1 -0.7670414150 3.274923e+00 7.316888e+00 haartd -0.0837790540 -4.655330e-02 -9.327549e-03 hivdur -0.1351044027 8.417971e-02 3.034638e-01 factor(arv)1 -6.7807591712 -2.075808e+00 2.629143e+00 rt:cd4 -0.0001150701 -4.667572e-05 2.171870e-05 attr(,"label") [1] "Fixed effects:" ...check for instance on row 2 and 5 (bold numbers)....isn't something really weird there ? lower interval greater than upper !!???? Federico -- *Little u can do against ignorance,....it will always disarm u: is the 2nd principle of thermodinamics made manifest, ...entropy in expansion.**....But setting order is the real quest 4 truth, ......and the mission of a (temporally) wise dude. * [[alternative HTML version deleted]]
David Winsemius
2011-Feb-08 21:50 UTC
[R] intervals {nlme} lower CI greater than upper CI !!!????
On Feb 8, 2011, at 4:14 PM, Federico Bonofiglio wrote:> Hi folks... > > check this out.. > >> GLU<-lme(gluc~rt*cd4+sex+age+rf+nadir+pharmac+factor(hcv) >> +factor(hbs)+ > + haartd+hivdur+factor(arv), > + random= ~rt|id, na.action=na.omit) > > >> intervals(GLU)$fixed > lower est. upper > (Intercept) 67.3467070345 7.362307e+01 7.989944e+01 > rt *0.0148050160* 6.249304e-02 1.101811e-01 > cd4 -0.0032133709 -6.566501e-04 1.900071e-03 > sexM 0.9601467541 3.565606e+00 6.171066e+00 > age *0.2436472425* 3.420025e-01 4.403578e-01 > rfOTHER -5.0933264232 -1.198439e+00 2.696449e+00 > rfETEROSESSUALE -2.1085096013 1.974725e+00 6.057960e+00 > rfOMOSESSUALE -5.8156940466 -1.891870e+00 2.031953e+00 > nadir -0.0063302516 8.604450e-04 8.051142e-03 > pharmacNRTI -4.2586487617 -1.252531e+00 1.753586e+00 > pharmacPI -3.3132784651 -8.829819e-01 1.547315e+00 > factor(hcv)1 -0.4338530808 2.374438e+00 5.182729e+00 > factor(hbs)1 -0.7670414150 3.274923e+00 7.316888e+00 > haartd -0.0837790540 -4.655330e-02 -9.327549e-03 > hivdur -0.1351044027 8.417971e-02 3.034638e-01 > factor(arv)1 -6.7807591712 -2.075808e+00 2.629143e+00 > rt:cd4 -0.0001150701 -4.667572e-05 2.171870e-05 > attr(,"label") > [1] "Fixed effects:" > > ...check for instance on row 2 and 5 (bold numbers)....isn't > something > really weird there ? > lower interval greater than upper !!????Perhaps you need to review your notes on scientific notation: > 1.101811e-01 > 0.0148050160 [1] TRUE > 0.2436472425 < 4.403578e-01 [1] TRUE -- David> > > Federico > > -- > *Little u can do against ignorance,....it will always disarm u: > is the 2nd principle of thermodinamics made manifest, ...entropy in > expansion.**....But setting order is the real quest 4 > truth, ......and the > mission of a (temporally) wise dude. > * > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.David Winsemius, MD West Hartford, CT
Ben Bolker
2011-Feb-08 22:00 UTC
[R] intervals {nlme} lower CI greater than upper CI !!!????
Federico Bonofiglio <bonoricus <at> gmail.com> writes:>You're confused by the formatting differences betwen columns. The first column uses 'standard' notation, the second and third use 'scientific' notation -- see comments inline.> > GLU<-lme(gluc~rt*cd4+sex+age+rf+nadir+pharmac+factor(hcv)+factor(hbs)+ > + haartd+hivdur+factor(arv), > + random= ~rt|id, na.action=na.omit) > > > intervals(GLU)$fixed > lower est. upper > (Intercept) 67.3467070345 7.362307e+01 7.989944e+01 > rt *0.0148050160* 6.249304e-02 1.101811e-01this is 0.0148... vs. 0.101811 ...> cd4 -0.0032133709 -6.566501e-04 1.900071e-03 > sexM 0.9601467541 3.565606e+00 6.171066e+00 > age *0.2436472425* 3.420025e-01 4.403578e-01this is 0.2436... vs. 0.4403578> rfOTHER -5.0933264232 -1.198439e+00 2.696449e+00[snip]> ...check for instance on row 2 and 5 (bold numbers)....isn't something > really weird there ? > lower interval greater than upper !!???? >
Eik Vettorazzi
2011-Feb-08 22:00 UTC
[R] intervals {nlme} lower CI greater than upper CI !!!????
Hi, there is nothing to complain?! 0.0148 < 6.24e-02=0.0624 < 1.1e-01 = 0.11 same for the 5th line. so I didn't get the point. Am 08.02.2011 22:14, schrieb Federico Bonofiglio:> Hi folks... > > check this out.. > >> GLU<-lme(gluc~rt*cd4+sex+age+rf+nadir+pharmac+factor(hcv)+factor(hbs)+ > + haartd+hivdur+factor(arv), > + random= ~rt|id, na.action=na.omit) > > >> intervals(GLU)$fixed > lower est. upper > (Intercept) 67.3467070345 7.362307e+01 7.989944e+01 > rt *0.0148050160* 6.249304e-02 1.101811e-01 > cd4 -0.0032133709 -6.566501e-04 1.900071e-03 > sexM 0.9601467541 3.565606e+00 6.171066e+00 > age *0.2436472425* 3.420025e-01 4.403578e-01 > rfOTHER -5.0933264232 -1.198439e+00 2.696449e+00 > rfETEROSESSUALE -2.1085096013 1.974725e+00 6.057960e+00 > rfOMOSESSUALE -5.8156940466 -1.891870e+00 2.031953e+00 > nadir -0.0063302516 8.604450e-04 8.051142e-03 > pharmacNRTI -4.2586487617 -1.252531e+00 1.753586e+00 > pharmacPI -3.3132784651 -8.829819e-01 1.547315e+00 > factor(hcv)1 -0.4338530808 2.374438e+00 5.182729e+00 > factor(hbs)1 -0.7670414150 3.274923e+00 7.316888e+00 > haartd -0.0837790540 -4.655330e-02 -9.327549e-03 > hivdur -0.1351044027 8.417971e-02 3.034638e-01 > factor(arv)1 -6.7807591712 -2.075808e+00 2.629143e+00 > rt:cd4 -0.0001150701 -4.667572e-05 2.171870e-05 > attr(,"label") > [1] "Fixed effects:" > > ...check for instance on row 2 and 5 (bold numbers)....isn't something > really weird there ? > lower interval greater than upper !!???? > > > Federico >-- Eik Vettorazzi Institut f?r Medizinische Biometrie und Epidemiologie Universit?tsklinikum Hamburg-Eppendorf Martinistr. 52 20246 Hamburg T ++49/40/7410-58243 F ++49/40/7410-57790