Hello!
I am trying to use the function "plot" in agricolae package. I used it a while ago and it gave me good results. Now it gaves the warning "In plot.group(out, xlab = "x... y...)", :
NAs introduced by coercion"
The result is a plot with no actual bars or boxes, just letters. I tried it with the example code and gave me the same result.
Basically I try to analyze my data with anova, then apply post-hoc test and simplyl visualize the results e.g. in a boxplot or barchart, without apllying the letters manually on top of each treatment. Something shown here, but it worked with simpler code earlier (plus that code does not work for me).
Let's say I try to use the LSD.test function with code given in the help.
library(agricolae)
data(sweetpotato)
model<-aov(yield~virus, data=sweetpotato)
out <- LSD.test(model,"virus", p.adj="bonferroni")
#stargraph
# Variation range: max and min
plot(out)
#endgraph
# Old version LSD.test()
df<-df.residual(model)
MSerror<-deviance(model)/df
out <- with(sweetpotato,LSD.test(yield,virus,df,MSerror))
#stargraph
# Variation interquartil range: Q75 and Q25
plot(out,variation="IQR")
#endgraph
out<-LSD.test(model,"virus",p.adj="hommel",console=TRUE)
plot(out,variation="SD") # variation standard deviation
And it gives no actual plots either. I feel puzzled and a bit helpless at the moment.
If anyone could help me out it would be much appreciated!
Dome