0

多くのdata.frameがあり、それぞれに多くの列が含まれています。私の最初のdata.framecol1= a、col2 = b、col3=cと言ってください

x-axis = b/aとy-axis=aをプロットしたいと思います。なんとかプロットしました（散布図）

plot (dataframe$b/dataframe$a, dataframe$a, xlim=...,ylim=..)

ここで、分散データのパターンを取得する必要があります（xとyの両方が変化しているため、線形回帰は必要ありません）。コマンドloess(..)を使用して、パターンを表示することができました。

lo_smooth<-loess(x,y, f=number, iter=number)

信頼区間（CI）をグラフに追加するにはどうすればよいですか？私の目標は、2つのdata.frameが互いにCI内にあるかどうかを確認することです。

4

1 に答える 1

0

あなたの試み（よくやった！）とあなたの説明を使用する解決策

いくつかのダミーデータ

dppm <- data.frame(a = runif(100, 1, 100), b = runif(100,1, 100))
dppm_2 <- data.frame(a = runif(100, 1, 75), b = runif(100,1,75))
dppm_3 <- data.frame(a = runif(100, 1,50), b = runif(100,1,50))

reshape2を使用してこれらのデータを単一のデータフレームライブラリにマージします（reshape2）

data_list <- list(dppm1 = dppm, dppm2 = dppm_2, dppm3 = dppm_3)
all_data <- melt(data_list, id.vars = c('a','b'))

この単一のデータフレームにL1は、識別子（のリストコンポーネントの名前）である列がありますdata_list。

head(all_data)
          a        b    L1
## 1 83.202896 36.94026 dppm1
## 2 42.618987 11.23863 dppm1
## 3 29.505029 11.91742 dppm1
## 4 63.569487 59.07395 dppm1
## 5 94.499772 47.32779 dppm1
## 6  4.535389 64.11570 dppm1

次に、この識別子を使用して、この結合されたデータセットと色をプロットできます。またfill、CIが同じように色付けされるように、smoothのを同じ識別子に設定します。

ggplot(all_data,aes(x = b/a, y = a, colour = L1)) + 
  geom_point() + 
  stat_smooth(method = "loess", se = TRUE,level = 0.90, aes(fill = L1))+
  coord_cartesian(ylim = c(0, 100))

ここに画像の説明を入力してください

于 2012-07-24T01:25:39.507 に答える