[R] Integration of vector syntax unknown

Benjamin Sabatini art_archaeology at yahoo.com
Fri Mar 22 18:36:11 CET 2013


Hello,

I'm very new to using R, but I was told it could do what I want. I'm not sure how best to enter the information but here goes...

I'm trying to transfer the following integral into R to solve for ln(gamma_1), on the left, for multiple instances of gamma_i and variable N_i.

gamma_i is, for example, (0, 0.03012048, 0.05000000, 0.19200000, 0.44000000, 0.62566845)

N_i (N_1 or N_2) is between 0 and 1 so that N_1+N_2=1, so if N_1=(0,.166,.180,.250,.325,.374), then N_2=(1.000, 0.834, 0.820, 0.750, 0.675, 0.626)

a_i (a_1 or a_2)


So, for gamma_i (in this case gamma_2), N_i (N_2), and a_i (a_2) first the following

a_i = ln(gamma_i)/(1-N_i)^2

then,

ln(gamma_1) = -a_2*N_1*N*2 - integration (from N_1=1, to N_1) a_2 dN_1


I hope that makes sense...



More information about the R-help mailing list