3

質問を少し修正しました。

     integg <- function(t,a,b,con,s){   
  u <- ifelse(t - a < 0 ,0, t - a) 
  l <- ifelse(t - (a+b) < 0,0, t - (a+b))  
  s * integrate(Vectorize(function(foo,x){x}),lower=l,upper=u,x=con)      
}

この方程式は私に積分値を与えます、そして私は3つの配列を持っています：引数'a'、'b'、および's'を表すAs、Bs、およびSs。

配列が次のようになっているとしましょう。

As <- seq(from=50,to=60,by=0.01)
Bs <- seq(from=130,to=140,by=0.01)
Ss <- seq(from=0.0001,to=0.01,by=0.0001)
# con is a constant
con <- 55
# I have 7 values for t and I want to do one at a time, 
# so for this example I have t=360
t <- 360
# although I'll want to do also for my other values c(0,20,40,60,120,240)

私の最終的な目標は、これらの配列As、Bs、およびSsのすべての組み合わせを相互にテストすることです。私はアウターを使おうとしていて、その後ループが失敗しました。

# first make one array w/ all possible combinations
all_poss <- outer(As,Bs,paste)
# now include the third array
all_poss <- outer(all_poss,Ss,paste)

head(all_poss)
> [1] "50 130 1e-04"    "50.01 130 1e-04" "50.02 130 1e-04" 
  [4] "50.03 130 1e-04" "50.04 130 1e-04" "50.05 130 1e-04"


    ### I would have to change my integg function a little bit, to deal w/ the pasted items
      integg2 <- function(t,con,all){  
    a <- strsplit(h,split=' ')[[1]][1]
    b <- strsplit(h,split=' ')[[1]][2]
    s <- strsplit(h,split=' ')[[1]][3]   

    u <- ifelse(t - a < 0 ,0, t - a) 
    l <- ifelse(t - (a+b) < 0,0, t - (a+b))  
    s * integrate(Vectorize(function(foo,x){x}),lower=l,upper=u,x=con)   
}

     ### I would then need to loop integg2() somehow through my list of all possibilities

      all_vals <- sapply(all_poss,integg2)
      # I haven't gotten this to work, but i'm not sure this is 
      # even an efficient way to do what I want

最後に、ある種のループが必要です。これらの配列のすべての可能性を組み合わせる方法について、より良いアイデアがあれば、ループのより効率的な方法を教えてください。

どんな助けでも素晴らしいでしょう！

4

1 に答える 1

1

この関数expand.gridは、入力したベクトル/係数のすべての組み合わせを含むデータフレームを作成します。もう1つの方法は、`data.tableで結合を使用することです。

integrateリストを返します。おそらく、を使用して値が必要です$value。

また、個人的な好みのために変数名を変更しました（そして、のような組み込み関数と名前が重複するのは嫌いですt）。

av <- seq(from=50,to=60,by=0.01)
bv <- seq(from=130,to=140,by=0.01)
sv <- seq(from=0.0001,to=0.01,by=0.0001)
tv <- c(seq(from=0,to=60,by=20),seq(from=120,to=360,by=120))
con <- 55

##method 1: using built-in functions (warning: can be slower and memory-intensive)
cmb <- expand.grid(list(av=av,bv=bv,sv=sv,tv=tv))
cmb <- within(cmb,{
    u <- ifelse(tv - av < 0 ,0, tv - av)
    l <- ifelse(tv - (av+bv) < 0,0, tv - (av+bv))
    value <- sv * mapply(function(...){integrate(...)$value},
        lower=l,upper=u,
        MoreArgs=list(f=Vectorize(function(x,constant){constant}),constant=con))
})

##method 2: using package data.table (for speed and efficient memory use)
dt.av <- data.table(av,k=1,key="k")
dt.bv <- data.table(bv,k=1,key="k")
dt.sv <- data.table(sv,k=1,key="k")
dt.tv <- data.table(tv,k=1,key="k")
cmb <- dt.av[dt.bv[dt.sv[dt.tv]]] #joins together
cmb[,u := ifelse(tv - av < 0 ,0, tv - av)]
cmb[,l := ifelse(tv - (av+bv) < 0,0, tv - (av+bv))]
cmb[,value:=mapply(function(...){integrate(...)$value},
    lower=l,upper=u,
    MoreArgs=list(f=Vectorize(function(x,constant){constant}),constant=con)
)]
于 2012-08-31T19:24:06.320 に答える