1

センタイル回帰で特定のセンタイルを超える値をカウントする方法を理解しようとしています。パッケージ lmsqreg を使用して基準値を使用して百分位回帰を作成し、その上に実験値をプロットしています。例えば。


#reference data

male.weight <- lmsqreg.fit(log(males$heart),
                           log(males$weight),
                           pvec = (centiles <- c(0.03, 0.1, 0.5, 0.9, 0.97)),
                           lam.fixed=1)


plot(log(males$weight),log(males$heart),
     axes=FALSE,
     frame.plot=TRUE,
     xlab="Weight (kg)",
     ylab="Heart weight (g)",
     main="Men", 
     col = "white",
     xlim = c(3.2,5.2),
     ylim = c(4.7, 7))

axis(1,at=log(xat <- c(3,5,7,9,11,13,15,17,19)*10),
     labels=xat)

axis(2,at=log(yat <- (1.5:10)*100),labels=yat)

lpars <- list(lty=rep(c(2,1,2),
                      c(2,1,2)),
              col=c("#444444","#999999","#000000")
              [c(1,2,3,2,1)])

for(l in 1:nrow(male.weight$qsys$outmat) ){
  lines(male.weight$qsys$targetx,
        male.weight$qsys$outmat[l,],
        lty=lpars$lty[l],col=lpars$col[l])

}



 #experimental data


points(log(m.heart$weight[m.heart$cod.g == "nCA"]),
       log(m.heart$heart.w[m.heart$cod.g == "nCA"]),
       pch = 21, col = "black")

points(log(m.heart$weight[m.heart$cod.g == "CA"]),
       log(m.heart$heart.w[m.heart$cod.g == "CA"]),
       pch = "♥", col = "red")

私が知りたいのは、m.heart$weight[m.heart$cod.g == "nCA"](実験データ) の何パーセントが 0.9 および 0.97 センタイルより上にあるかということです。似たような回答が他にもあるかもしれませんが、私はこの分野の命名法に詳しくないので、何を検索すればよいかわかりません。

前もって感謝します、

ローラ

4

0 に答える 0