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]]


