[edit]
1 蟆曙誤 #
- R蠍磯 る 覿/螳覈,蟾覈蠏/蟲一
一企 一危
- Windows R [] -> [襴覲蟆] 一危郁 襴襦 覲蟆]
> cost.d <-read.table("cost.d", header=T)
> cov.wt(cost.d, cor=T)
$cov
fuel repair capital
fuel 23.013361 12.366395 2.906609
repair 12.366395 17.544111 4.773082
capital 2.906609 4.773082 13.963334
$center
fuel repair capital
12.218611 8.112500 9.590278
$n.obs
[1] 36
$cor
fuel repair capital
fuel 1.0000000 0.6154424 0.1621444
repair 0.6154424 1.0000000 0.3049570
capital 0.1621444 0.3049570 1.0000000
> pairs(cost.d)
#譟郁唄覿 一 #覲瑚 譟郁唄朱 譯殊伎 磯觜(襴暑) 襴觜(覦覲)伎 蟯螻 coplot(repair~fuel|capital, data=cost.d)
> set.seed(2)
> library(mvtnorm)
> x <- rmvnorm(20, sigma=matrix(c(1, 0.98, 0.98, 1), 2))
> x
[,1] [,2]
[1,] 0.62872170 1.05125106
[2,] -0.61625134 -0.81208810
[3,] 2.23543105 2.23568460
[4,] 0.36178518 0.79807392
[5,] -0.05901189 -0.04696427
[6,] -1.44903706 -1.81448776
[7,] 0.85017069 0.81754232
[8,] -0.56312456 -0.61359221
[9,] 2.03787185 1.86925931
[10,] 0.07583804 0.13642630
[11,] 0.79101426 0.83645122
[12,] 0.96199934 0.86826634
[13,] 0.37699900 0.58472713
[14,] -0.98481735 -0.87797193
[15,] 0.88841158 0.52652788
[16,] -2.16636324 -1.92376783
[17,] -0.41200937 -0.78035322
[18,] -0.54347793 -0.67618649
[19,] 0.43040268 0.20810378
[20,] 0.17868536 0.08269170
> plot(x[,1], x[,2], xlab="X1", ylab="X2")
mcc.f<-function(x, var1, var2)
{
s<-var(x)
s0<-sqrt(s[var1,var1])
s0q<-s[var1,var2]
Sq<-s[var2,var2]
sqrt(s0q%*%solve(Sq)%*%s0q)/s0
}
cost.d <-read.table("cost.d", header=T)
source("mcc.f")
mcc.f(cost.d, 1, c(2,3))
[,1]
[1,] 0.6160264
pcc.f<-function(x,z1,z2,x2)
{
x<-as.matrix(x)
s<-var(x)
s22<-s[x2,x2]
y1<-x[,z1]-x[,x2]%*%solve(s22)%*%s[z1,x2]
y2<-x[,z2]-x[,x2]%*%solve(s22)%*%s[z2,x2]
cor(y1,y2)
}
soil.d <-read.table("soil.d", header=T)
source("pcc.f")
pcc.f(soil.d, 1, 2, c(3,4))
[,1]
[1,] 0.7475448
#る襯れ 蠏碁(碁Μ 蠍, 襭螳 蟇磯Μ襯 蟇磯Μ襯 讌 伎)
andrews.f<-function(x,v)
{
n<-nrow(x)
p<-ncol(x)
nv<-length(v)
x<-as.matrix(x)
m<-cov.wt(x)$center
M<-matrix(m*rep(1,n*p),n,p,byrow=T)
d<-diag(1/sqrt(diag(cov.wt(x)$cov)))
x<-(x-M)%*%d
l<-100
tt<-seq(from=-pi,to=pi,length=l);
andr.v<-c(rep(1/sqrt(2),l))
for(i in 1:((p-1)/2))
andr.v<-c(andr.v,sin(i*tt),cos(i*tt))
if( p%%2 == 0 )
andr.v<-andr.v[-( ((p-1)*l+1): p*l )]
andr.m<-matrix(andr.v,nrow=l,ncol=p)
y <- andr.m%*%t(x)
z <- matrix(0,nv*l,2)
lab<-vector("numeric",nv*l)
plot(rep(tt,n),as.vector(y),type="n",xlab="t",ylab="Andrews plot")
for(i in 1:nv) {
lines(tt,y[,v[i]])
for(j in 1:l ) {
z[(i-1)*l+j,1]<-tt[j]
z[(i-1)*l+j,2]<-y[j,v[i]]
lab[(i-1)*l+j]<-v[i]
}
}
abline(h=0)
identify(z,labels=lab)
}
andrews.f(cost.d, c(1:36))
andrews.f(cost.d, c(9, 21, 23))
