| To: | s-news@lists.biostat.wustl.edu |
|---|---|
| Subject: | Question re DF for testing covariate in lme() (W2k) |
| From: | KINLEY_ROBERT@LILLY.COM |
| Date: | Mon, 01 Sep 2003 10:20:39 +0100 |
|
I'm puzzled by the fact that lme() reports 1 DF for the denominator in the t-test , and also for the t-value used in calculating the confidence interval, for the slope term in the model described below . ... if anyone can point me in the right direction I'd be very grateful. Details :- I have the dataset called 'TF' :- Strength type AnalDay Assay StorageDays AssayOcc 1 4 FinalSol 0 40.4355 0 1 2 4 FinalSol 4 41.0857 4 2 3 4 FinalSol 7 40.8318 7 3 4 10 FinalSol 0 103.6812 0 1 5 10 FinalSol 4 104.5013 4 2 6 10 FinalSol 7 104.1406 7 3 7 20 FinalSol 0 195.6095 0 1 8 20 FinalSol 4 196.6439 4 2 9 20 FinalSol 7 196.1706 7 3 and I'm trying to fit the model Assay = Strength + AssayOcc + StorageDays where Strength is a fixed effect factor, StorageDays is a covariate ... a fixed slope effect , assumed common to all 3 strengths AssayOcc is a random effect for the variability between assay occasions [ Perhaps the model itself is inappropriate ..? ] > TF.lme <- lme(fixed = Assay ~ Strength + StorageDays, data = "" random = ~ 1 | AssayOcc) > summary(TF.lme) Linear mixed-effects model fit by REML Data: TF AIC BIC logLik 18.22203 15.87866 -3.111015 Random effects: Formula: ~ 1 | AssayOcc (Intercept) Residual StdDev: 0.4563897 0.09662428 Fixed effects: Assay ~ Strength + StorageDays Value Std.Error DF t-value p-value (Intercept) 113.4024 0.4309204 4 263.163 <.0001 Strength1 31.6617 0.0394467 4 802.645 <.0001 Strength2 41.2318 0.0227746 4 1810.431 <.0001 StorageDays 0.0751 0.0925766 1 0.811 0.5661 -- Standardized Within-Group Residuals: Min Q1 Med Q3 Max -1.019119 -0.107893 -0.0007243514 0.08046544 1.127012 Number of Observations: 9 Number of Groups: 3 > intervals(TF.lme) Approximate 95% confidence intervals Fixed effects: lower est. upper (Intercept) 112.205987 113.40241396 114.598841 Strength1 31.552162 31.66168333 31.771205 Strength2 41.168540 41.23177222 41.295005 StorageDays -1.101194 0.07510225 1.251399 Random Effects: Level: AssayOcc lower est. upper sd((Intercept)) 0.1115825 0.4563897 1.866705 Within-group standard error: lower est. upper 0.04832077 0.09662428 0.193214 If anyone can help I'd be most grateful thanks Bob |
| <Prev in Thread] | Current Thread | [Next in Thread> |
|---|---|---|
| ||
| Previous by Date: | Re: How S chooses the .Random.seed, Patrick Burns |
|---|---|
| Next by Date: | Re: Question re DF for testing covariate in lme() (W2k), Spencer Graves |
| Previous by Thread: | Mixed-Model fit - DF mismatch v. SAS/JMP, KINLEY_ROBERT |
| Next by Thread: | Re: Question re DF for testing covariate in lme() (W2k), Spencer Graves |
| Indexes: | [Date] [Thread] [Top] [All Lists] |