0

apply 一連の関数のいずれかを使用して記述できるはずのいくつかの R コードを使用していますが、その方法がわかりません。複数の列を持つデータフレームがあり、関数を呼び出したいのですが、関数の入力はデータフレームの複数の列を使用しています。このデータと関数 f があるとします。

data<- data.frame(T=c(1,2,3,4), S=c(3,7,8,4), K=c(5,6,11,9))
data
V<-c(0.1,0.2,0.3,0.4,0.5,0.6)

f<-function(para_h,S,T,a,t,b){
  r<- V
  steps<-T
  # Recursive form: Terminal condition for the A and B at time T
  A_T=0
  B_T=0
  A=c()
  B=c()
  #  A and B a time T-1
  A[1]= r[steps]*a 
  B[1]= a*para_h[5]+ ((para_h[4])^(-2))         
  # Recursion back to time t      
  for (i in 2:steps){
    A[i]= A[i-1]+ r[steps-i+1]*a + para_h[1]*B[i-1]
    B[i]= para_h[2]*B[i-1]+a*para_h[5]+ (para_h[4]^(-2))
}
  f = exp(log(S)*a + A[t] + B[t]*b )

  return(f)
}

この関数は、いくつかの特定の値に対してうまく機能します:

> para_h<-c(0.1,0.2,0.3,0.4,0.5,0.7)
> f(para_h,S=3,T=2,a=0.4,t=1,b=0.1)
[1] 3.204144

データ フレームの各列 S と T に関数を適用したいと考えています。したがって、私のコードは次のようになります。

mapply(function(para_h,S,T,a,t,b) f(para_h,S,T,a,t,b) ,para_h,S=data$S,T=data$T,a=0.4,t=1,b=0.1)

これによりエラーが発生します。

> mapply(function(para_h,S,T,a,t,b) f(para_h,S,T,a,t,b) ,para_h,S=data$S,T=data$T,a=0.4,t=1,b=0.1)
Error in A[i] = A[i - 1] + r[steps - i + 1] * a + para_h[1] * B[i - 1] : 
  replacement has length zero

問題は次のとおりだと確信しています。「ステップ」はベクトルです。エレガントなソリューションに本当に感謝します。これが何らかの意味を成していることを願っています。アドバイスをいただければ幸いです。

4

1 に答える 1

1

いくつかのこと：

para_h1）関数の各呼び出しは完全なベクトルを期待しますが、mapplyコードでは一度に1つの値しか受け取ることができないため、おそらく次のようなものが必要です：

 mapply(function(S,T) f(para_h,S,T,a=0.4,t=1,b=0.1), data$S, data$T)

またはこれ：

 apply(data,1,function(d) f(para_h,d['S'],d['T'],a=0.4,t=1,b=0.1))

2) 関数は次の場合にエラーをスローしますT==1(これは の最初の行の場合ですdata)。そのため、このコードを実行できるようにサンプル データ セットを変更する必要がある場合があります。

于 2015-02-19T06:49:01.100 に答える