2

次のデータセット (CEU) があります。

group  x      y
1     -23     100
1     -0.90   69.62
1     -0.90   72.03
2     -23     100
2      0.69   48.01
2      0.69   45.63

グループの値ごとに、以下に示す関数を x 値と y 値の各サブセットに適用したいと考えています。次に、すべての結果を結合し、それらをテーブルに書き込んでエクスポートしたいと思います。

これを行うためにplyr関数を正確に適用する方法がわかりません...それが実際に正しい行動方針である場合。

x<-c(-23.0000,-0.9031,-0.9031)
y<-c(100,85.72,86.65)

par<-c(16.88,100.28,-.75,4.129)

dcrit<-function(d) { 
    sumsq<-0
    for (i in 1:length(x)){
      sumsq<-sumsq+ (y[i]-(par[1]+(par[2]-par[1])/(1+10^((x[i]-par[3])*d))))^2      
    }
    sumsq
}

S<-function(par) { 
    a<-par[1]
    b<-par[2]
    c<-par[3]
    d<-par[4]
    sumsq<-0
    for (i in 1:length(x)){
      sumsq<-sumsq+ (y[i]-(a+(b-a)/(1+10^((x[i]-c)*d))))^2      
    }
    sumsq
}
optim(par,S)

CEU <- read.csv(file="C:/files/CEU.csv",head=TRUE,sep=",")
CEU

data <- ddply(CEU,.(group),function(xy) 
{
par[1]<-min(y)
par[2]<-100
par[3]<-x[[which.min(abs(y-50))]]
par[4]<-optimize(dcrit,interval=c(-100,100))$minimum

o<-optim(par,S)
par<-o$par

a<-par[1];
b<-par[2];
c<-par[3];
d<-par[4];

k<-(b-a)/(20-a)-1
if (k>0) ec20<-c+1/d*log10(k) else ec20<-NA
ec20

z<-(b-a)/(50-a)-1
 if (z>0) ec50<-c+1/d*log10(z) else ec50<-NA
ec50

j<-(b-a)/(80-a)-1
if (j>0) ec80<-c+1/d*log10(j) else ec80<-NA
ec80

data.frame(ec20, ec50, ec80)

})

data

コードはエラーなしで実行されますが、次のように設定された元の x 値と y 値に対してのみ実行されます。

 x<-c(-23.0000,-0.9031,-0.9031)
 y<-c(100,85.72,86.65)

データセット CEU の x と y の値は、ddply では使用されません。グループ値の場合のように、元の x と y を反復的に置き換えることはありません。データには適切な数のグループがあり、ec20/ec50/ec80 の値は正確ですが、元の x と y についてのみです。 

> data
   group       ec20       ec50       ec80
1      1 -0.3652977 -0.6843279 -0.8530892
2      2 -0.3652977 -0.6843279 -0.8530892
3      3 -0.3652977 -0.6843279 -0.8530892
4      4 -0.3652977 -0.6843279 -0.8530892
5      5 -0.3652977 -0.6843279 -0.8530892
4

1 に答える 1

3

出力を生成する必要があるだけです。

これがあなたの出力がある場所だと思いますか？

k<-(b-a)/(20-a)-1
if (k>0) ec20<-c+1/d*log10(k) else ec20<-NA
ec20

z<-(b-a)/(50-a)-1
 if (z>0) ec50<-c+1/d*log10(z) else ec50<-NA
ec50

j<-(b-a)/(80-a)-1
if (j>0) ec80<-c+1/d*log10(j) else ec80<-NA
ec80

それらをdata.frame関数の最後に入れます：

    ...
    data.frame(ec20, ec50, ec80)
}

これで、data.frameそれらすべてを含むが得られます。ec20ec50ec80

optim に関する問題について: 問題は次の部分にあると思います

par[3]<-x[which.min(abs(y-50))]

R のSingle[は、通常の添字ではありません。スライスを取得します。この場合、data.frame列です。その行はpar、数値ベクトルからlist. さらにブラケットを追加します。

par[3]<-x[[which.min(abs(y-50))]]
于 2011-08-29T23:02:42.013 に答える