ミグストラノート

ラノベ読みの日記です。その他、地質。

R compositions 三角ダイアグラム

cran.r-project.org/web/packages/compositions/compositions.pdfのP237とP145,その他諸々のページ

datファイルは.以下の内容をdat形式で保存.標準化して拡張子をdatへ変更.

10.0 80.0 10.0 A
20.0 20.0 60.0 B
50.0 20.0 30.0 C 

値のどれかが0だとバグる???

以下は,P236のコピペを一部改変.英語の説明を読んでも正直分からない部分がある.値をいじって試してみる必要がある.

# require("compositions")
# table <- read.table("san.dat", row.names = 4)
# plot(acomp(table), pca=TRUE,col.pca="red", axes=TRUE)

#require("compositions")
#table <- read.table("san.dat", row.names = 4)
#plot(acomp(table),axes=TRUE)
#ternaryAxis(side=1:3,pos=0,col.axis="black",col.lab="red")
#角を赤に 

#require("compositions")
#table <- read.table("san.dat", row.names = 4)
#plot(acomp(table),axes=TRUE)
#ternaryAxis(side=1:3,tick=TRUE,pos=0,lty="solid",lwd=10,col.axis="black",col.lab="red",Xlab="A",Ylab="B",Zlab="C",small=TRUE,)
#角を赤に lwdは枠の太さ P236の説明を読め

require("compositions")
table <- read.table("san.dat", row.names = 4)
plot(acomp(table),axes=TRUE)
ternaryAxis(side=1:3,at=seq(0.2,0.8,by=0.2),
tick=TRUE,pos=0,
font.axis=par("font.axis"),
font.lab=par("font.lab"),
lty="solid",lwd=1,
len.tck=0.025,dist.lab=0.03,
dist.axis=0.03,
lty.tck="solid",
col.axis=par("col.axis"),
col.lab=par("col.lab"),
cex.axis=par("cex.axis"),
cex.lab=par("cex.lab"),
Xlab="A",Ylab="B",Zlab="C",small=FALSE,
xpd=NA,aspanel=FALSE)
#Xlabが極の名前の変更 smallをFALSEにした方がバランスが良い