Hello R users,
I am having difficulting getting multcomp to run.
I have a dataframe attached with a numeric variable q12a and a numeric variable
quota (which is really a classification variable).
quota has 10 levels and unequal sample sizes.
a12a has some missing data.
I am interested in doing pairwise testing across the 10 quota groups on q12a.
Using the ctest package the following code ran generating the pvalue matrix part
of which I list below. When I use the multcomp package and attempt to replicate
this analysis, I cannot get it to work. Below I show three attempts that
failed.
Any help would be much appreciated.> pairwise.t.test(q12a, quota, p.adj = "fdr")
$method
[1] "t tests with pooled SD"
$data.name
[1] "q12a and quota"
$p.value
1 2 3 4 5 6 7 8 9
2 4.805732e-09 NA NA NA NA NA NA NA NA
> simtest(q12a ~ quota)
Error in parseformula(formula, data, subset, na.action, whichf, ...) :
at least one factor required
> simtest(q12a ~ factor(quota))
Error in parse(file, n, text, prompt) : parse error
> simtest(q12a ~ factor(quota),na.action=na.exclude)
Error in parse(file, n, text, prompt) : parse error
Greg Blevins
Partner, The Market Solutions Group
Windows XP, version 1.9.
[[alternative HTML version deleted]]