Year <- c(1888,1895,1901,1909,1917,1925,1933,1941,1949,1957,1965,1973,1981,1989,1997) Naru <-c(0.026349206,0.09689298,0.194725544,0.50310302,0.616606405,0.846628182,0.856009695,0.88542422,0.970925784,0.992300941,0.985845662,0.991990764,0.986659622,0.992662148,0.995891426) fm <- nls(Naru~a/(1+b*exp(c*1:15)),start=c(a=1,b=1,c=-1),trace=TRUE) plot(Year,Naru,cex=2) lines(Year,fitted(fm),col=2,lty=2,lwd=2) legend(locator(1),c("実測値","予測値"),col=1:2,lty=1:2,lwd=2)
library(MASS) data <- read.csv("y:\\20110203_Ba_cType2.csv",header=T) # 元データ res.mca <- mca(data,nf=3,abbrev=T) res.mca$d^2 / sum(res.mca$d^2) # 寄与率 res.mca$cs res.mca$cs biplot(res.mca$rs,res.mca$cs, var.axes=F)