Jun Shen
2011-Apr-29 21:49 UTC
[R] Use nparcomp function from nparcomp library to run post hoc
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]]
Dennis Murphy
2011-Apr-29 22:49 UTC
[R] Use nparcomp function from nparcomp library to run post hoc
Hi: Is this the function nparcomp() in the nparcomp package or the one from the mutoss package? When using functions from packages, it is useful to indicate the package name. I'm assuming you're using the nparcomp package, because your code worked for me when that package was loaded:> library(nparcomp)Loading required package: multcomp Loading required package: mvtnorm Loading required package: survival Loading required package: splines> nparcomp(Ulceration~Group,data=test,type='Dunnett',control='Non-treated')Nonparametric Multiple Comparison Procedure based on relative contrast effects , Type of Contrast : Dunnett NOTE: *-------------------Weight Matrix------------------* - Weight matrix for choosen contrast based on all-pairs comparisons *-----------Analysis of relative effects-----------* - Simultaneous Confidence Intervals for relative effects p(i,j) with confidence level 0.95 - Method = Multivariate Delta-Method (Logit) - p-Values for H_0: p(i,j)=1/2 *----------------Interpretation--------------------* p(a,b) > 1/2 : b tends to be larger than a *--------------Mult.Distribution-------------------* - Equicoordinate Quantile - Global p-Value *--------------------------------------------------* $weight.matrix < snipped for brevity - all zeros > $Data.Info Sample Size 1 Duoderm 24 2 Fibrase 24 3 Kollagenase 24 4 Non-treated 24 5 Stimulen 24 6 Vehicle 24 $Analysis.of.relative.effects Comparison rel.effect confidence.interval t.value 1 p(Non-treated,Duoderm) 0.5 [ 0.499 ; 0.501 ] 0 2 p(Non-treated,Fibrase) 0.5 [ 0.499 ; 0.501 ] 0 3 p(Non-treated,Kollagenase) 0.5 [ 0.499 ; 0.501 ] 0 4 p(Non-treated,Stimulen) 0.5 [ 0.499 ; 0.501 ] 0 5 p(Non-treated,Vehicle) 0.5 [ 0.499 ; 0.501 ] 0 p.value.adjusted p.value.unadjusted 1 1 1 2 1 1 3 1 1 4 1 1 5 1 1 $Mult.Distribution Quantile p.Value.global 1 2.568766 1 $Correlation [1] NA A graphic also appears indicating zero effect, which is what one would expect since Ulceration = 5 for every observation in the data frame.> sessionInfo()R version 2.13.0 (2011-04-13) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=en_US.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] splines stats graphics grDevices utils datasets methods [8] base other attached packages: [1] nparcomp_1.0-1 multcomp_1.2-5 survival_2.36-9 mvtnorm_0.9-999 [5] sos_1.3-0 brew_1.0-6 plyr_1.5.2 loaded via a namespace (and not attached): [1] tcltk_2.13.0 tools_2.13.0 Check your version of R and the nparcomp package against this. If you have an older version of R or nparcomp, perhaps an upgrade is sufficient to fix the problem. HTH, Dennis On Fri, Apr 29, 2011 at 2:49 PM, Jun Shen <jun.shen.ut at gmail.com> wrote:> 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 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. >