# Observations on flowers. Middle Russia, 1998-2005
# Possible way of analysis:
if(0){
ss <- read.table("open/stamens.txt", h=TRUE, sep="\t")
str(ss)
rr <- as.character(as.roman(3:5))
rr
for (i in 3:5) ss <- cbind(ss, assign(rr[i], ss$STAMENS %% i == 0))
names(ss)[3:5] <- rr
ss$Other <- rowSums(ss[,3:5]) == 0
str(ss)
(ss.t <- xtabs(cbind(III, IV, V, Other) ~ SPECIES, data=ss))
dotchart(ss.t, lcolor=1, pch=19)
VTcoeffs(ss.t)[2,] # asmisc.r
chisq.test(ss.t, simulate.p.value=TRUE)
pairwise.Table2.test(ss.t, exact=TRUE) # asmisc.r
}
#
SPECIES	Species name
STAMENS	Number of stamens