Displaying 1 result from an estimated 1 matches for "ple3".
Did you mean:
ple
2005 Feb 18
1
Contingency tables profiles
...like my solution about legend in profiles graphics: I inserted empty
extra columns in order to avoid tue superimposed of legend.
#Data
N <- matrix(0,3,6)
N[1,] <- c(7,7,5,0,4,4)
N[2,] <- c(0,0,0,5,5,5)
N[3,] <- c(4,4,0,0,3,0)
rownames(N) <- c("Ple1","Ple2","Ple3")
colnames(N) <- c("A","B","C","D","E","F")
N
# Row profiles
PF <- N/rowSums(N)
# Columns profiles
PC <- t(N)/rowSums(t(N))
# Graphics
par(mfrow=c(2,2))
barplot(cbind(cbind(t(PF),0),0),legend=colnames(F),density=100,cex.names...