Dear list,
I tried to use the nparcomp to run some post hoc non-parametric comparison
and got and error.
Error in uniroot(pfct, interval = interval) :
f() values at end points not of opposite sign
Appreciate any comments.
the command line:
>nparcomp(Ulceration~Group,data=test,type='Dunnett',control='Non-treated')
Jun
===================================================================
data as follows
structure(list(Group = c("Duoderm", "Duoderm", "Duoderm", "Duoderm",
"Duoderm", "Duoderm", "Duoderm", "Duoderm", "Duoderm", "Duoderm",
"Duoderm", "Duoderm", "Duoderm", "Duoderm", "Duoderm", "Duoderm",
"Duoderm", "Duoderm", "Duoderm", "Duoderm", "Duoderm", "Duoderm",
"Duoderm", "Duoderm", "Fibrase", "Fibrase", "Fibrase", "Fibrase",
"Fibrase", "Fibrase", "Fibrase", "Fibrase", "Fibrase", "Fibrase",
"Fibrase", "Fibrase", "Fibrase", "Fibrase", "Fibrase", "Fibrase",
"Fibrase", "Fibrase", "Fibrase", "Fibrase", "Fibrase", "Fibrase",
"Fibrase", "Fibrase", "Kollagenase", "Kollagenase", "Kollagenase",
"Kollagenase", "Kollagenase", "Kollagenase", "Kollagenase", "Kollagenase",
"Kollagenase", "Kollagenase", "Kollagenase", "Kollagenase", "Kollagenase",
"Kollagenase", "Kollagenase", "Kollagenase", "Kollagenase", "Kollagenase",
"Kollagenase", "Kollagenase", "Kollagenase", "Kollagenase", "Kollagenase",
"Kollagenase", "Non-treated", "Non-treated", "Non-treated", "Non-treated",
"Non-treated", "Non-treated", "Non-treated", "Non-treated", "Non-treated",
"Non-treated", "Non-treated", "Non-treated", "Non-treated", "Non-treated",
"Non-treated", "Non-treated", "Non-treated", "Non-treated", "Non-treated",
"Non-treated", "Non-treated", "Non-treated", "Non-treated", "Non-treated",
"Stimulen", "Stimulen", "Stimulen", "Stimulen", "Stimulen", "Stimulen",
"Stimulen", "Stimulen", "Stimulen", "Stimulen", "Stimulen", "Stimulen",
"Stimulen", "Stimulen", "Stimulen", "Stimulen", "Stimulen", "Stimulen",
"Stimulen", "Stimulen", "Stimulen", "Stimulen", "Stimulen", "Stimulen",
"Vehicle", "Vehicle", "Vehicle", "Vehicle", "Vehicle", "Vehicle",
"Vehicle", "Vehicle", "Vehicle", "Vehicle", "Vehicle", "Vehicle",
"Vehicle", "Vehicle", "Vehicle", "Vehicle", "Vehicle", "Vehicle",
"Vehicle", "Vehicle", "Vehicle", "Vehicle", "Vehicle", "Vehicle"
), Ulceration = c(5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5), Inflamation = c(3, 4, 3, 2, 3, 3, 4, 4, 2, 2, 3, 3,
3, 3, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 2, 3, 3, 4, 3, 3, 2, 4, 4,
4, 4, 4, 4, 4, 3, 3, 4, 3, 5, 3, 3, 4, 4, 3, 3, 2, 4, 2, 3, 3,
4, 3, 4, 3, 3, 4, 3, 4, 2, 3, 3, 4, 2, 3, 4, 3, 2, 3, 3, 3, 2,
3, 2, 2, 2, 2, 4, 3, 2, 3, 3, 4, 3, 3, 4, 3, 4, 2, 4, 3, 4, 2,
4, 3, 4, 3, 2, 2, 2, 2, 3, 2, 3, 2, 4, 3, 2, 4, 4, 4, 2, 2, 3,
3, 2, 4, 3, 2, 3, 2, 2, 2, 4, 2, 3, 2, 3, 2, 3, 3, 3, 4, 3, 3,
4, 4, 2, 3, 2, 3), Fibroplasia = c(4, 4, 4, 4, 4, 3, 4, 4, 4,
3, 4, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 2, 4, 4, 3,
2, 4, 4, 4, 4, 4, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 3, 3, 3, 4, 4, 3, 4, 4, 4, 3, 4,
4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4,
4, 4, 3, 4, 4, 3, 4, 3, 2, 3, 3, 4, 3, 3, 4, 4, 3, 3, 3, 4, 3,
3, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 3, 3, 4, 4, 4, 4, 4, 3, 4,
3, 4, 4, 4, 4, 4, 4, 4, 4), Fibrosis.and.Adexnal.Atrophy = c(4,
4, 4, 3, 4, 4, 4, 4, 4, 3, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 3, 4, 3, 4, 4, 3, 3, 4, 4, 4, 4, 4, 4, 4, 3, 3, 4, 3, 4,
4, 4, 3, 4, 3, 3, 3, 3, 3, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 3, 3,
4, 4, 3, 3, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4,
4, 4, 3, 4, 4, 4, 3, 4, 3, 4, 3, 4, 4, 3, 4, 3, 2, 3, 3, 4, 4,
3, 4, 4, 3, 3, 3, 4, 3, 3, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 4,
3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4), Inflammation = c(2,
2, 2, 1, 1, 1, 2, 3, 1, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1,
1, 1, 1, 1, 1, 1, 1, 2, NA, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 2,
1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1,
1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, NA, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1,
2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 2), Fibroplasia.1 =
c(4,
4, 4, 4, 4, 3, 4, 4, 4, 3, 4, 3, 3, 3, 4, 4, 4, 4, 3, 3, 3, 3,
4, 4, 4, 4, 3, 3, 4, 3, NA, 4, 4, 4, 4, 4, 4, 4, 3, 3, 3, 3,
3, 3, 4, 4, 4, 4, 3, 3, 3, 3, 4, 4, 4, 3, 4, 3, 4, 4, 4, 4, 3,
3, 4, 4, 4, 3, 4, 3, 4, 4, 3, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4,
4, 3, 3, 3, NA, 4, 4, 4, 4, 3, 4, 3, 3, 3, 3, 3, 4, 2, 4, 3,
4, 4, 3, 4, 4, 2, 3, 2, 3, 3, 3, 3, 4, 4, 4, 4, 3, 3, 3, 3, 3,
4, 3, 3, 4, 4, 4, 4, 3, 3, 4, 3, 3, 4, 4, 4, 4, 3, 4, 4), Fibrosis = c(3,
3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,
3, 3, 3, 3, 3, 3, 3, 3, NA, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,
3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,
3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2, 3, 3, 3, 3, 3,
3, 3, 3, 3, NA, 3, 3, 3, 3, 3, 2, 2, 3, 3, 3, 3, 3, 2, 3, 3,
3, 3, 3, 3, 3, 2, 3, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2, 3,
3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3), Total = c(25,
26, 25, 22, 24, 22, 26, 27, 23, 21, 24, 21, 22, 23, 25, 24, 25,
26, 23, 23, 24, 24, 24, 24, 21, 24, 20, 24, 24, 22, 12, 25, 25,
25, 26, 26, 26, 24, 22, 21, 24, 22, 24, 22, 25, 23, 25, 24, 22,
21, 24, 21, 24, 24, 25, 23, 25, 23, 24, 24, 24, 24, 20, 21, 24,
26, 21, 22, 25, 23, 21, 24, 23, 24, 23, 22, 23, 23, 22, 23, 25,
23, 23, 24, 24, 24, 23, 22, 17, 25, 25, 21, 25, 22, 25, 19, 24,
23, 22, 23, 21, 16, 22, 20, 25, 22, 22, 24, 25, 19, 21, 20, 25,
23, 20, 22, 25, 24, 23, 25, 24, 22, 22, 21, 23, 23, 23, 20, 24,
23, 24, 23, 23, 22, 24, 23, 23, 24, 26, 25, 23, 23, 24, 25)), .Names =
c("Group",
"Ulceration", "Inflamation", "Fibroplasia", "Fibrosis.and.Adexnal.Atrophy",
"Inflammation", "Fibroplasia.1", "Fibrosis", "Total"), row.names = c(NA,
144L), class = "data.frame")
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
