Rで条件付きパーセンタイル解析
>>
wakekijyun_d<-
function (x, y, kugiri = 32, minmin = 1, last = F,points=F,pch=1,xlab="x",ylab="y"){
d <- exp(seq(min(log(x)), max(log(x)), length.out = kugiri))
st <- d[1:(length(d) - 1)]
fi <- d[2:length(d)]
if (minmin != 1) {
length2 <- 0
for (i in 1:length(st)) {
xb2 <- x[x >= st[i] & x < fi[i]]
length2[i] <- length(xb2)
}
f12345 <- max(which*1 {
xb2 <- x[x >= st[i] & x < fi[i]]
yb2 <- y[x >= st[i] & x < fi[i]]
median1[i] <- median(xb2)
median2[i] <- median(yb2)
logmean1[i] <- mean(log(xb2[xb2 != Inf]))
logmean2[i] <- mean(log(yb2[yb2 != Inf]))
l123[i] <- 0
mean1[i] <- mean(xb2, p = 0.5)
mean2[i] <- mean(yb2, p = 0.5)
q25[i] <- quantile(yb2, 0.25)
q75[i] <- quantile(yb2, 0.75)
q5[i] <- quantile(yb2, 0.05)
q95[i] <- quantile(yb2, 0.95)
q50[i] <- quantile(yb2, 0.5)
IQR[i] <- IQR(yb2)
try(d123 <- density(log(yb2)))
length1[i] <- length(yb2)
try(mode[i] <- exp(d123$x[which(d123$y == max(d123$y))]))
}
jj<- sqrt(st * fi)
if(points==F){
plot(jj,q5,ylim=c(min(q5[!is.na(q5) ]),max(q95[!is.na(q95)])),log="xy",col=0,xlab=xlab,ylab=ylab,cex.axis=1.5,cex.lab=1.5)
}
points(jj,q5,col=1,pch=pch)
points(jj,q25,col=2,pch=pch)
points(jj,q50,col=3,pch=pch)
points(jj,q75,col=4,pch=pch)
points(jj,q95,col=5,pch=pch)
data.frame(sqrt(st * fi), median2, mean1, mean2, mode, q25,
q75, IQR, logmean1, logmean2, l123, length1, median1,q5,q95)
}
*1:length2 >= minmin))) if (last == T) { st <- st[1:(f12345 - 1)] fi <- c(fi[2:(f12345 - 1)], fi[length(fi)]) } else { st <- st[1:(f12345 - 1)] fi <- c(fi[2:(f12345)]) } } median1 <- 0 median2 <- 0 logmean1 <- 0 logmean2 <- 0 mean1 <- 0 mean2 <- 0 mode <- 0 length1 <- 0 q25 <- 0 q75 <- 0 IQR <- 0 l123 <- 0 j123 <- y[x >= 30] q1<-0 q5<-0 q25<-0 q50<-0 q75<-0 q95<-0 for (i in 1:length(st
