[R] t.test() with missing values

Petr PIKAL petr.pikal at precheza.cz
Fri Sep 14 14:12:42 CEST 2007


Hi

r-help-bounces at r-project.org napsal dne 14.09.2007 13:50:58:

> Birgit Lemcke wrote:
> >
> > Am 14.09.2007 um 12:05 schrieb Peter Dalgaard:
> >
> >> S Ellison wrote:
> >>>
> >>>>>> Peter Dalgaard <P.Dalgaard at biostat.ku.dk> 14/09/2007 09:26:16 >>>
> >>>>>>


<snip>


> >>>
> >> Yes, but....
> >>
> >> I don't think lme() will do better than what you can do by hand: Get 
two
> >> independent estimates of mu1-mu2 (one estimate from the pairs and one
> >> from the singlets), compute a weighted average using the s.e.'s and 
test
> >> that against zero (possibly after testing them for equality for good
> >> measure). This is easy if you use a plug-in approach: first assume 
that
> >> the s.e. are known, then plug in their empirical value. The tricky 
bit
> >> is to calculate the DF in the style of Welch's test.
> >
> >             I apologise but I really can not follow your explanations.
> > I am  R and Stastistics Beginner.
> >
> >             What do you mean with mu1-mu2 and what are s.e.´s?
> >
> That was a reply to S. Ellison. If you don't understand it, don't worry;
> you'll probably need to read a book chapter or more about weighted
> analyses to get up to speed for that. 
> 
> mu1, mu2 : (theoretical) mean for group 1, 2
> s.e.: standard error

But as Birgit actually does not have any paired values, according to the 
data she had sent, she can not do paired t.test at all. The only way is to 
compare averages from each vector by non paired t.test or to get some new 
values for which she have counterparts. 

Regards
Petr



> >             Once again thank you for your help.
> >
> >             Birgit
> >
> > 
> >
> > 
> >>
> >> --   O__  ---- Peter Dalgaard             Øster Farimagsgade 5, 
Entr.B
> >>   c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
> >>  (*) \(*) -- University of Copenhagen   Denmark          Ph:  (+45)
> >> 35327918
> >> ~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)                  FAX: (+45)
> >> 35327907
> >>
> >>
> >
> > Birgit Lemcke
> > Institut für Systematische Botanik
> > Zollikerstrasse 107
> > CH-8008 Zürich
> > Switzerland
> > Ph: +41 (0)44 634 8351
> > birgit.lemcke at systbot.uzh.ch
> >
> >
> >
> >
> >
> >
> 
> 
> -- 
>    O__  ---- Peter Dalgaard             Øster Farimagsgade 5, Entr.B
>   c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
>  (*) \(*) -- University of Copenhagen   Denmark          Ph:  (+45) 
35327918
> ~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)                  FAX: (+45) 
35327907
> 
> ______________________________________________
> 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.



More information about the R-help mailing list