「統計学関連なんでもあり」の過去ログ--- 039

No.01809　Rのplot；個別タイトルと軸要素の順序　　【Rはじめました】　2006/12/08(Fri) 10:21

おかげさまで，だいぶ分かってきたのですが，また壁にぶつかっております。何とか自分で解決しようと試みておりましたが，どうしてもできません。初歩的なことで恐縮ですが，お力添えお願いします。

species <- c("A","B","A","C","D","E","A","B","E","A","C","B","D")
length <- c(13,20,14,18,18,13,14,17,18,14,19,20,15)
color <- c(1,1,1,2,2,1,1,2,1,2,2,2,1)
sex <- c(1,1,1,1,1,1,1,2,2,2,2,2,2)
data <- data.frame(species = species, length = length, color = color, sex = sex)
data
data <- data[data$species =="A"|data$species =="B"|data$species =="D",]
data

data $species <- factor(data$species)
par(mfrow=c(2,2), mar=c(2,2,1,1), oma=c(1,3,1,1))
for(i in 1:2) {
for (j in 1:2) {
selected <- is.element(data$color,i) & is.element(data$sex,j)
print(sum(selected))
boxplot(length[selected]~(species[selected]),data=data,ylim=c(10,30),mtext("length", side=2, line=1,outer=TRUE))
}}
としたときに，sexとcolorで抽出した4つのグラフにそれぞれ別のタイトルを入れたいのです。
ここで，

title <- c("male,red","male,bule","female,red","femele,bule")
for(k in 1:4){
mtext=(title[,k],side=3)

など，色々やってみましたが，うまくいきません。ご教授頂けないでしょうか。

また，A，B，Dの順に並んでいるx軸の要素をA，D，Bに並べ換えることは可能でしょうか。

No.01811　Re: Rのplot；個別タイトルと軸要素の順序　　【takahashi】　2006/12/08(Fri) 12:30

species <- c("A","B","A","C","D","E","A","B","E","A","C","B","D")
length <- c(13,20,14,18,18,13,14,17,18,14,19,20,15)
color <- c(1,1,1,2,2,1,1,2,1,2,2,2,1)
sex <- c(1,1,1,1,1,1,1,2,2,2,2,2,2)
data <- data.frame(species = species, length = length, color = color, sex = sex)
data <- data[data$species =="A"|data$species =="B"|data$species =="D",]

data $species <- factor(data$species)
tbl<-c(a=1,b=3,c=2)
data$species.index<-factor(tbl[data$species])

par(mfrow=c(2,2), mar=c(2,2,1,1), oma=c(1,3,1,1))
title <- c("male,red","male,bule","female,red","femele,bule")
for(i in 1:2) {
for (j in 1:2) {
selected <- is.element(data$color,i) & is.element(data$sex,j)
boxplot(length[selected]~(species.index[selected]),data=data,names=c("A","D","B"),ylim=c(10,30))
#boxplot(length[selected]~(species[selected]),data=data,ylim=c(10,30))
mtext(title[(i-1)*2+j],side=3)
}}
mtext("length", side=2, line=1,outer=TRUE)

もっと簡単な方法があるかもしれませn

No.01822　Re: Rのplot；個別タイトルと軸要素の順序　　【青木繁伸】　2006/12/08(Fri) 18:05

boxplot には main 引数を与えることができます
boxplot(length[selected]~(species.index[selected]), data=data, 
names=c("A","D","B"), ylim=c(10,30), main=title[(i-1)*2+j])
でよろしいかと。。。

No.01850　Re: Rのplot；個別タイトルと軸要素の順序　　【Rはじめました】　2006/12/10(Sun) 02:04

青木先生，takahasi様，いつも有難うございます。
今回も非常に助かりました。感謝感激です。