thr3ads.net - R help - [R] r programming help [Jun 2005]

If this information is useful, please help other people find it:
Share via:

Mohammad Ehsanul Karim

2005-Jun-22 13:36 UTC

[R] r programming help

Dear list,

Is there anyway i can make the following formula short
by r-programming?

CYCLE.n<-c(NA,
WET[1]*DRY[1],
WET[1]*DRY[2]+WET[2]*DRY[1],
WET[1]*DRY[3]+WET[2]*DRY[2]+WET[3]*DRY[1],
WET[1]*DRY[4]+WET[2]*DRY[3]+WET[3]*DRY[2]+WET[4]*DRY[1],
WET[1]*DRY[5]+WET[2]*DRY[4]+WET[3]*DRY[3]+WET[4]*DRY[2]+WET[5]*DRY[1],
WET[1]*DRY[6]+WET[2]*DRY[5]+WET[3]*DRY[4]+WET[4]*DRY[3]+WET[5]*DRY[2]+WET[6]*DRY[1],
WET[1]*DRY[7]+WET[2]*DRY[6]+WET[3]*DRY[5]+WET[4]*DRY[4]+WET[5]*DRY[3]+WET[6]*DRY[2]+WET[7]*DRY[1],
WET[1]*DRY[8]+WET[2]*DRY[7]+WET[3]*DRY[6]+WET[4]*DRY[5]+WET[5]*DRY[4]+WET[6]*DRY[3]+WET[7]*DRY[2]+WET[8]*DRY[1],
WET[1]*DRY[9]+WET[2]*DRY[8]+WET[3]*DRY[7]+WET[4]*DRY[6]+WET[5]*DRY[5]+WET[6]*DRY[4]+WET[7]*DRY[3]+WET[8]*DRY[2]+WET[9]*DRY[1],
WET[1]*DRY[10]+WET[2]*DRY[9]+WET[3]*DRY[8]+WET[4]*DRY[7]+WET[5]*DRY[6]+WET[6]*DRY[5]+WET[7]*DRY[4]+WET[8]*DRY[3]+WET[9]*DRY[2]+WET[10]*DRY[1],
WET[1]*DRY[11]+WET[2]*DRY[10]+WET[3]*DRY[9]+WET[4]*DRY[8]+WET[5]*DRY[7]+WET[6]*DRY[6]+WET[7]*DRY[5]+WET[8]*DRY[4]+WET[9]*DRY[3]+WET[10]*DRY[2]+WET[11]*DRY[1],
WET[1]*DRY[12]+WET[2]*DRY[11]+WET[3]*DRY[10]+WET[4]*DRY[9]+WET[5]*DRY[8]+WET[6]*DRY[7]+WET[7]*DRY[6]+WET[8]*DRY[5]+WET[9]*DRY[4]+WET[10]*DRY[3]+WET[11]*DRY[2]+WET[12]*DRY[1],
WET[1]*DRY[13]+WET[2]*DRY[12]+WET[3]*DRY[11]+WET[4]*DRY[10]+WET[5]*DRY[9]+WET[6]*DRY[8]+WET[7]*DRY[7]+WET[8]*DRY[6]+WET[9]*DRY[5]+WET[10]*DRY[4]+WET[11]*DRY[3]+WET[12]*DRY[2]+WET[13]*DRY[1],
WET[1]*DRY[14]+WET[2]*DRY[13]+WET[3]*DRY[12]+WET[4]*DRY[11]+WET[5]*DRY[10]+WET[6]*DRY[9]+WET[7]*DRY[8]+WET[8]*DRY[7]+WET[9]*DRY[6]+WET[10]*DRY[5]+WET[11]*DRY[4]+WET[12]*DRY[3]+WET[13]*DRY[2]+WET[14]*DRY[1],
WET[1]*DRY[15]+WET[2]*DRY[14]+WET[3]*DRY[13]+WET[4]*DRY[12]+WET[5]*DRY[11]+WET[6]*DRY[10]+WET[7]*DRY[9]+WET[8]*DRY[8]+WET[9]*DRY[7]+WET[10]*DRY[6]+WET[11]*DRY[5]+WET[12]*DRY[4]+WET[13]*DRY[3]+WET[14]*DRY[2]+WET[15]*DRY[1],
WET[1]*DRY[16]+WET[2]*DRY[15]+WET[3]*DRY[14]+WET[4]*DRY[13]+WET[5]*DRY[12]+WET[6]*DRY[11]+WET[7]*DRY[10]+WET[8]*DRY[9]+WET[9]*DRY[8]+WET[10]*DRY[7]+WET[11]*DRY[6]+WET[12]*DRY[5]+WET[13]*DRY[4]+WET[14]*DRY[3]+WET[15]*DRY[2]+WET[16]*DRY[1],
WET[1]*DRY[17]+WET[2]*DRY[16]+WET[3]*DRY[15]+WET[4]*DRY[14]+WET[5]*DRY[13]+WET[6]*DRY[12]+WET[7]*DRY[11]+WET[8]*DRY[10]+WET[9]*DRY[9]+WET[10]*DRY[8]+WET[11]*DRY[7]+WET[12]*DRY[6]+WET[13]*DRY[5]+WET[14]*DRY[4]+WET[15]*DRY[3]+WET[16]*DRY[2]+WET[17]*DRY[1],
WET[1]*DRY[18]+WET[2]*DRY[17]+WET[3]*DRY[16]+WET[4]*DRY[15]+WET[5]*DRY[14]+WET[6]*DRY[13]+WET[7]*DRY[12]+WET[8]*DRY[11]+WET[9]*DRY[10]+WET[10]*DRY[9]+WET[11]*DRY[8]+WET[12]*DRY[7]+WET[13]*DRY[6]+WET[14]*DRY[5]+WET[15]*DRY[4]+WET[16]*DRY[3]+WET[17]*DRY[2]+WET[18]*DRY[1],
WET[1]*DRY[19]+WET[2]*DRY[18]+WET[3]*DRY[17]+WET[4]*DRY[16]+WET[5]*DRY[15]+WET[6]*DRY[15]+WET[7]*DRY[13]+WET[8]*DRY[12]+WET[9]*DRY[11]+WET[10]*DRY[10]+WET[11]*DRY[9]+WET[12]*DRY[8]+WET[13]*DRY[7]+WET[14]*DRY[6]+WET[15]*DRY[5]+WET[16]*DRY[4]+WET[17]*DRY[3]+WET[18]*DRY[2]+WET[19]*DRY[1],
WET[1]*DRY[20]+WET[2]*DRY[19]+WET[3]*DRY[18]+WET[4]*DRY[17]+WET[5]*DRY[16]+WET[6]*DRY[15]+WET[7]*DRY[14]+WET[8]*DRY[13]+WET[9]*DRY[12]+WET[10]*DRY[11]+WET[11]*DRY[10]+WET[12]*DRY[9]+WET[13]*DRY[8]+WET[14]*DRY[7]+WET[15]*DRY[6]+WET[16]*DRY[5]+WET[17]*DRY[4]+WET[18]*DRY[3]+WET[19]*DRY[2]+WET[20]*DRY[1],
WET[1]*DRY[21]+WET[2]*DRY[20]+WET[3]*DRY[19]+WET[4]*DRY[18]+WET[5]*DRY[17]+WET[6]*DRY[16]+WET[7]*DRY[15]+WET[8]*DRY[14]+WET[9]*DRY[13]+WET[10]*DRY[12]+WET[11]*DRY[11]+WET[12]*DRY[10]+WET[13]*DRY[9]+WET[14]*DRY[8]+WET[15]*DRY[7]+WET[16]*DRY[6]+WET[17]*DRY[5]+WET[18]*DRY[4]+WET[19]*DRY[3]+WET[20]*DRY[2]+WET[21]*DRY[1],
WET[1]*DRY[22]+WET[2]*DRY[21]+WET[3]*DRY[20]+WET[4]*DRY[19]+WET[5]*DRY[18]+WET[6]*DRY[17]+WET[7]*DRY[16]+WET[8]*DRY[15]+WET[9]*DRY[14]+WET[10]*DRY[13]+WET[11]*DRY[12]+WET[12]*DRY[11]+WET[13]*DRY[10]+WET[14]*DRY[9]+WET[15]*DRY[8]+WET[16]*DRY[7]+WET[17]*DRY[6]+WET[18]*DRY[5]+WET[19]*DRY[4]+WET[20]*DRY[3]+WET[21]*DRY[2]+WET[22]*DRY[1],
WET[1]*DRY[23]+WET[2]*DRY[22]+WET[3]*DRY[21]+WET[4]*DRY[20]+WET[5]*DRY[19]+WET[6]*DRY[18]+WET[7]*DRY[17]+WET[8]*DRY[16]+WET[9]*DRY[15]+WET[10]*DRY[14]+WET[11]*DRY[13]+WET[12]*DRY[12]+WET[13]*DRY[11]+WET[14]*DRY[10]+WET[15]*DRY[9]+WET[16]*DRY[8]+WET[17]*DRY[7]+WET[18]*DRY[6]+WET[19]*DRY[5]+WET[20]*DRY[4]+WET[21]*DRY[3]+WET[22]*DRY[2]+WET[23]*DRY[1],
WET[1]*DRY[24]+WET[2]*DRY[23]+WET[3]*DRY[22]+WET[4]*DRY[21]+WET[5]*DRY[20]+WET[6]*DRY[19]+WET[7]*DRY[18]+WET[8]*DRY[17]+WET[9]*DRY[16]+WET[10]*DRY[15]+WET[11]*DRY[14]+WET[12]*DRY[13]+WET[13]*DRY[12]+WET[14]*DRY[11]+WET[15]*DRY[10]+WET[16]*DRY[9]+WET[17]*DRY[8]+WET[18]*DRY[7]+WET[19]*DRY[6]+WET[20]*DRY[5]+WET[21]*DRY[4]+WET[22]*DRY[3]+WET[23]*DRY[2]+WET[24]*DRY[1],
WET[1]*DRY[25]+WET[2]*DRY[24]+WET[3]*DRY[23]+WET[4]*DRY[22]+WET[5]*DRY[21]+WET[6]*DRY[20]+WET[7]*DRY[19]+WET[8]*DRY[18]+WET[9]*DRY[17]+WET[10]*DRY[16]+WET[11]*DRY[15]+WET[12]*DRY[14]+WET[13]*DRY[13]+WET[14]*DRY[12]+WET[15]*DRY[11]+WET[16]*DRY[10]+WET[17]*DRY[9]+WET[18]*DRY[8]+WET[19]*DRY[7]+WET[20]*DRY[6]+WET[21]*DRY[5]+WET[22]*DRY[4]+WET[23]*DRY[3]+WET[24]*DRY[2]+WET[25]*DRY[1],
WET[1]*DRY[26]+WET[2]*DRY[25]+WET[3]*DRY[24]+WET[4]*DRY[23]+WET[5]*DRY[22]+WET[6]*DRY[21]+WET[7]*DRY[20]+WET[8]*DRY[19]+WET[9]*DRY[18]+WET[10]*DRY[17]+WET[11]*DRY[16]+WET[12]*DRY[15]+WET[13]*DRY[14]+WET[14]*DRY[13]+WET[15]*DRY[12]+WET[16]*DRY[11]+WET[17]*DRY[10]+WET[18]*DRY[9]+WET[19]*DRY[8]+WET[20]*DRY[7]+WET[21]*DRY[6]+WET[22]*DRY[5]+WET[23]*DRY[4]+WET[24]*DRY[3]+WET[25]*DRY[2]+WET[26]*DRY[1],
WET[1]*DRY[27]+WET[2]*DRY[26]+WET[3]*DRY[25]+WET[4]*DRY[24]+WET[5]*DRY[23]+WET[6]*DRY[22]+WET[7]*DRY[21]+WET[8]*DRY[20]+WET[9]*DRY[19]+WET[10]*DRY[18]+WET[11]*DRY[17]+WET[12]*DRY[16]+WET[13]*DRY[15]+WET[14]*DRY[14]+WET[15]*DRY[13]+WET[16]*DRY[12]+WET[17]*DRY[11]+WET[18]*DRY[10]+WET[19]*DRY[9]+WET[20]*DRY[8]+WET[21]*DRY[7]+WET[22]*DRY[6]+WET[23]*DRY[5]+WET[24]*DRY[4]+WET[25]*DRY[3]+WET[26]*DRY[2]+WET[27]*DRY[1],
WET[1]*DRY[28]+WET[2]*DRY[27]+WET[3]*DRY[26]+WET[4]*DRY[25]+WET[5]*DRY[24]+WET[6]*DRY[23]+WET[7]*DRY[22]+WET[8]*DRY[21]+WET[9]*DRY[20]+WET[10]*DRY[19]+WET[11]*DRY[18]+WET[12]*DRY[17]+WET[13]*DRY[16]+WET[14]*DRY[15]+WET[15]*DRY[14]+WET[16]*DRY[13]+WET[17]*DRY[12]+WET[18]*DRY[11]+WET[19]*DRY[10]+WET[20]*DRY[9]+WET[21]*DRY[8]+WET[22]*DRY[7]+WET[23]*DRY[6]+WET[24]*DRY[5]+WET[25]*DRY[4]+WET[26]*DRY[3]+WET[27]*DRY[2]+WET[28]*DRY[1],
WET[1]*DRY[29]+WET[2]*DRY[28]+WET[3]*DRY[27]+WET[4]*DRY[26]+WET[5]*DRY[25]+WET[6]*DRY[24]+WET[7]*DRY[23]+WET[8]*DRY[22]+WET[9]*DRY[21]+WET[10]*DRY[20]+WET[11]*DRY[19]+WET[12]*DRY[18]+WET[13]*DRY[17]+WET[14]*DRY[16]+WET[15]*DRY[15]+WET[16]*DRY[14]+WET[17]*DRY[13]+WET[18]*DRY[12]+WET[19]*DRY[11]+WET[20]*DRY[10]+WET[21]*DRY[9]+WET[22]*DRY[8]+WET[23]*DRY[7]+WET[24]*DRY[6]+WET[25]*DRY[5]+WET[26]*DRY[4]+WET[27]*DRY[3]+WET[28]*DRY[2]+WET[29]*DRY[1]
)




Thank you for your time. 
PS: if the above formula can not be seen correctly,
one can see details at
http://www.angelfire.com/ab5/get5/alt.ren.txt

----------------------------------

Mohammad Ehsanul Karim 

Web: http://snipurl.com/ehsan 
ISRT, University of Dhaka, BD

Peter Dalgaard

2005-Jun-22 14:03 UTC

head link

[R] r programming help

Mohammad Ehsanul Karim <wildscop at yahoo.com> writes:
> Dear list,
> 
> Is there anyway i can make the following formula short
> by r-programming?
> 
> CYCLE.n<-c(NA,
> WET[1]*DRY[1],
> WET[1]*DRY[2]+WET[2]*DRY[1],
> WET[1]*DRY[3]+WET[2]*DRY[2]+WET[3]*DRY[1],....

As far as I can see: 

z <- toeplitz(DRY)
z[upper.tri(z)] <- 0
c(NA, z %*% WET)

or convolve() with suitable options, padding, and/or cutting (but
beware, there could be devils in the details). convolve(WET,DRY,
type="o") gives you about twice what you need, I believe.


-- 
   O__  ---- Peter Dalgaard             ??ster Farimagsgade 5, Entr.B
  c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
 (*) \(*) -- University of Copenhagen   Denmark          Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)                  FAX: (+45) 35327907

Peter Wolf

2005-Jun-22 15:53 UTC

head link

[R] r programming help

Try:

 > DRY<-c(2,5,3,7,11)
 > WET<-(1:5)*10
 > print(filter(c(rep(0,length(WET)),DRY),WET))

Time Series:
Start = 1
End = 10
Frequency = 1
[1]  NA  NA   0  20  90 190 360 640  NA  NA

 > CYCLE.n<-c(NA,
 > WET[1]*DRY[1],
 > WET[1]*DRY[2]+WET[2]*DRY[1],
 > WET[1]*DRY[3]+WET[2]*DRY[2]+WET[3]*DRY[1],
 > WET[1]*DRY[4]+WET[2]*DRY[3]+WET[3]*DRY[2]+WET[4]*DRY[1],
 > WET[1]*DRY[5]+WET[2]*DRY[4]+WET[3]*DRY[3]+WET[4]*DRY[2]+WET[5]*DRY[1])
 > print(CYCLE.n)
[1]  NA  20  90 190 360 640

Peter Wolf

Mohammad Ehsanul Karim wrote:
>Dear list,
>
>Is there anyway i can make the following formula short
>by r-programming?
>
>CYCLE.n<-c(NA,
>WET[1]*DRY[1],
>WET[1]*DRY[2]+WET[2]*DRY[1],
>WET[1]*DRY[3]+WET[2]*DRY[2]+WET[3]*DRY[1],
>WET[1]*DRY[4]+WET[2]*DRY[3]+WET[3]*DRY[2]+WET[4]*DRY[1],
>WET[1]*DRY[5]+WET[2]*DRY[4]+WET[3]*DRY[3]+WET[4]*DRY[2]+WET[5]*DRY[1],
>WET[1]*DRY[6]+WET[2]*DRY[5]+WET[3]*DRY[4]+WET[4]*DRY[3]+WET[5]*DRY[2]+WET[6]*DRY[1],
>WET[1]*DRY[7]+WET[2]*DRY[6]+WET[3]*DRY[5]+WET[4]*DRY[4]+WET[5]*DRY[3]+WET[6]*DRY[2]+WET[7]*DRY[1],
>WET[1]*DRY[8]+WET[2]*DRY[7]+WET[3]*DRY[6]+WET[4]*DRY[5]+WET[5]*DRY[4]+WET[6]*DRY[3]+WET[7]*DRY[2]+WET[8]*DRY[1],
>WET[1]*DRY[9]+WET[2]*DRY[8]+WET[3]*DRY[7]+WET[4]*DRY[6]+WET[5]*DRY[5]+WET[6]*DRY[4]+WET[7]*DRY[3]+WET[8]*DRY[2]+WET[9]*DRY[1],
>WET[1]*DRY[10]+WET[2]*DRY[9]+WET[3]*DRY[8]+WET[4]*DRY[7]+WET[5]*DRY[6]+WET[6]*DRY[5]+WET[7]*DRY[4]+WET[8]*DRY[3]+WET[9]*DRY[2]+WET[10]*DRY[1],
>WET[1]*DRY[11]+WET[2]*DRY[10]+WET[3]*DRY[9]+WET[4]*DRY[8]+WET[5]*DRY[7]+WET[6]*DRY[6]+WET[7]*DRY[5]+WET[8]*DRY[4]+WET[9]*DRY[3]+WET[10]*DRY[2]+WET[11]*DRY[1],
>WET[1]*DRY[12]+WET[2]*DRY[11]+WET[3]*DRY[10]+WET[4]*DRY[9]+WET[5]*DRY[8]+WET[6]*DRY[7]+WET[7]*DRY[6]+WET[8]*DRY[5]+WET[9]*DRY[4]+WET[10]*DRY[3]+WET[11]*DRY[2]+WET[12]*DRY[1],
>WET[1]*DRY[13]+WET[2]*DRY[12]+WET[3]*DRY[11]+WET[4]*DRY[10]+WET[5]*DRY[9]+WET[6]*DRY[8]+WET[7]*DRY[7]+WET[8]*DRY[6]+WET[9]*DRY[5]+WET[10]*DRY[4]+WET[11]*DRY[3]+WET[12]*DRY[2]+WET[13]*DRY[1],
>WET[1]*DRY[14]+WET[2]*DRY[13]+WET[3]*DRY[12]+WET[4]*DRY[11]+WET[5]*DRY[10]+WET[6]*DRY[9]+WET[7]*DRY[8]+WET[8]*DRY[7]+WET[9]*DRY[6]+WET[10]*DRY[5]+WET[11]*DRY[4]+WET[12]*DRY[3]+WET[13]*DRY[2]+WET[14]*DRY[1],
>WET[1]*DRY[15]+WET[2]*DRY[14]+WET[3]*DRY[13]+WET[4]*DRY[12]+WET[5]*DRY[11]+WET[6]*DRY[10]+WET[7]*DRY[9]+WET[8]*DRY[8]+WET[9]*DRY[7]+WET[10]*DRY[6]+WET[11]*DRY[5]+WET[12]*DRY[4]+WET[13]*DRY[3]+WET[14]*DRY[2]+WET[15]*DRY[1],
>WET[1]*DRY[16]+WET[2]*DRY[15]+WET[3]*DRY[14]+WET[4]*DRY[13]+WET[5]*DRY[12]+WET[6]*DRY[11]+WET[7]*DRY[10]+WET[8]*DRY[9]+WET[9]*DRY[8]+WET[10]*DRY[7]+WET[11]*DRY[6]+WET[12]*DRY[5]+WET[13]*DRY[4]+WET[14]*DRY[3]+WET[15]*DRY[2]+WET[16]*DRY[1],
>WET[1]*DRY[17]+WET[2]*DRY[16]+WET[3]*DRY[15]+WET[4]*DRY[14]+WET[5]*DRY[13]+WET[6]*DRY[12]+WET[7]*DRY[11]+WET[8]*DRY[10]+WET[9]*DRY[9]+WET[10]*DRY[8]+WET[11]*DRY[7]+WET[12]*DRY[6]+WET[13]*DRY[5]+WET[14]*DRY[4]+WET[15]*DRY[3]+WET[16]*DRY[2]+WET[17]*DRY[1],
>WET[1]*DRY[18]+WET[2]*DRY[17]+WET[3]*DRY[16]+WET[4]*DRY[15]+WET[5]*DRY[14]+WET[6]*DRY[13]+WET[7]*DRY[12]+WET[8]*DRY[11]+WET[9]*DRY[10]+WET[10]*DRY[9]+WET[11]*DRY[8]+WET[12]*DRY[7]+WET[13]*DRY[6]+WET[14]*DRY[5]+WET[15]*DRY[4]+WET[16]*DRY[3]+WET[17]*DRY[2]+WET[18]*DRY[1],
>WET[1]*DRY[19]+WET[2]*DRY[18]+WET[3]*DRY[17]+WET[4]*DRY[16]+WET[5]*DRY[15]+WET[6]*DRY[15]+WET[7]*DRY[13]+WET[8]*DRY[12]+WET[9]*DRY[11]+WET[10]*DRY[10]+WET[11]*DRY[9]+WET[12]*DRY[8]+WET[13]*DRY[7]+WET[14]*DRY[6]+WET[15]*DRY[5]+WET[16]*DRY[4]+WET[17]*DRY[3]+WET[18]*DRY[2]+WET[19]*DRY[1],
>WET[1]*DRY[20]+WET[2]*DRY[19]+WET[3]*DRY[18]+WET[4]*DRY[17]+WET[5]*DRY[16]+WET[6]*DRY[15]+WET[7]*DRY[14]+WET[8]*DRY[13]+WET[9]*DRY[12]+WET[10]*DRY[11]+WET[11]*DRY[10]+WET[12]*DRY[9]+WET[13]*DRY[8]+WET[14]*DRY[7]+WET[15]*DRY[6]+WET[16]*DRY[5]+WET[17]*DRY[4]+WET[18]*DRY[3]+WET[19]*DRY[2]+WET[20]*DRY[1],
>WET[1]*DRY[21]+WET[2]*DRY[20]+WET[3]*DRY[19]+WET[4]*DRY[18]+WET[5]*DRY[17]+WET[6]*DRY[16]+WET[7]*DRY[15]+WET[8]*DRY[14]+WET[9]*DRY[13]+WET[10]*DRY[12]+WET[11]*DRY[11]+WET[12]*DRY[10]+WET[13]*DRY[9]+WET[14]*DRY[8]+WET[15]*DRY[7]+WET[16]*DRY[6]+WET[17]*DRY[5]+WET[18]*DRY[4]+WET[19]*DRY[3]+WET[20]*DRY[2]+WET[21]*DRY[1],
>WET[1]*DRY[22]+WET[2]*DRY[21]+WET[3]*DRY[20]+WET[4]*DRY[19]+WET[5]*DRY[18]+WET[6]*DRY[17]+WET[7]*DRY[16]+WET[8]*DRY[15]+WET[9]*DRY[14]+WET[10]*DRY[13]+WET[11]*DRY[12]+WET[12]*DRY[11]+WET[13]*DRY[10]+WET[14]*DRY[9]+WET[15]*DRY[8]+WET[16]*DRY[7]+WET[17]*DRY[6]+WET[18]*DRY[5]+WET[19]*DRY[4]+WET[20]*DRY[3]+WET[21]*DRY[2]+WET[22]*DRY[1],
>WET[1]*DRY[23]+WET[2]*DRY[22]+WET[3]*DRY[21]+WET[4]*DRY[20]+WET[5]*DRY[19]+WET[6]*DRY[18]+WET[7]*DRY[17]+WET[8]*DRY[16]+WET[9]*DRY[15]+WET[10]*DRY[14]+WET[11]*DRY[13]+WET[12]*DRY[12]+WET[13]*DRY[11]+WET[14]*DRY[10]+WET[15]*DRY[9]+WET[16]*DRY[8]+WET[17]*DRY[7]+WET[18]*DRY[6]+WET[19]*DRY[5]+WET[20]*DRY[4]+WET[21]*DRY[3]+WET[22]*DRY[2]+WET[23]*DRY[1],
>WET[1]*DRY[24]+WET[2]*DRY[23]+WET[3]*DRY[22]+WET[4]*DRY[21]+WET[5]*DRY[20]+WET[6]*DRY[19]+WET[7]*DRY[18]+WET[8]*DRY[17]+WET[9]*DRY[16]+WET[10]*DRY[15]+WET[11]*DRY[14]+WET[12]*DRY[13]+WET[13]*DRY[12]+WET[14]*DRY[11]+WET[15]*DRY[10]+WET[16]*DRY[9]+WET[17]*DRY[8]+WET[18]*DRY[7]+WET[19]*DRY[6]+WET[20]*DRY[5]+WET[21]*DRY[4]+WET[22]*DRY[3]+WET[23]*DRY[2]+WET[24]*DRY[1],
>WET[1]*DRY[25]+WET[2]*DRY[24]+WET[3]*DRY[23]+WET[4]*DRY[22]+WET[5]*DRY[21]+WET[6]*DRY[20]+WET[7]*DRY[19]+WET[8]*DRY[18]+WET[9]*DRY[17]+WET[10]*DRY[16]+WET[11]*DRY[15]+WET[12]*DRY[14]+WET[13]*DRY[13]+WET[14]*DRY[12]+WET[15]*DRY[11]+WET[16]*DRY[10]+WET[17]*DRY[9]+WET[18]*DRY[8]+WET[19]*DRY[7]+WET[20]*DRY[6]+WET[21]*DRY[5]+WET[22]*DRY[4]+WET[23]*DRY[3]+WET[24]*DRY[2]+WET[25]*DRY[1],
>WET[1]*DRY[26]+WET[2]*DRY[25]+WET[3]*DRY[24]+WET[4]*DRY[23]+WET[5]*DRY[22]+WET[6]*DRY[21]+WET[7]*DRY[20]+WET[8]*DRY[19]+WET[9]*DRY[18]+WET[10]*DRY[17]+WET[11]*DRY[16]+WET[12]*DRY[15]+WET[13]*DRY[14]+WET[14]*DRY[13]+WET[15]*DRY[12]+WET[16]*DRY[11]+WET[17]*DRY[10]+WET[18]*DRY[9]+WET[19]*DRY[8]+WET[20]*DRY[7]+WET[21]*DRY[6]+WET[22]*DRY[5]+WET[23]*DRY[4]+WET[24]*DRY[3]+WET[25]*DRY[2]+WET[26]*DRY[1],
>WET[1]*DRY[27]+WET[2]*DRY[26]+WET[3]*DRY[25]+WET[4]*DRY[24]+WET[5]*DRY[23]+WET[6]*DRY[22]+WET[7]*DRY[21]+WET[8]*DRY[20]+WET[9]*DRY[19]+WET[10]*DRY[18]+WET[11]*DRY[17]+WET[12]*DRY[16]+WET[13]*DRY[15]+WET[14]*DRY[14]+WET[15]*DRY[13]+WET[16]*DRY[12]+WET[17]*DRY[11]+WET[18]*DRY[10]+WET[19]*DRY[9]+WET[20]*DRY[8]+WET[21]*DRY[7]+WET[22]*DRY[6]+WET[23]*DRY[5]+WET[24]*DRY[4]+WET[25]*DRY[3]+WET[26]*DRY[2]+WET[27]*DRY[1],
>WET[1]*DRY[28]+WET[2]*DRY[27]+WET[3]*DRY[26]+WET[4]*DRY[25]+WET[5]*DRY[24]+WET[6]*DRY[23]+WET[7]*DRY[22]+WET[8]*DRY[21]+WET[9]*DRY[20]+WET[10]*DRY[19]+WET[11]*DRY[18]+WET[12]*DRY[17]+WET[13]*DRY[16]+WET[14]*DRY[15]+WET[15]*DRY[14]+WET[16]*DRY[13]+WET[17]*DRY[12]+WET[18]*DRY[11]+WET[19]*DRY[10]+WET[20]*DRY[9]+WET[21]*DRY[8]+WET[22]*DRY[7]+WET[23]*DRY[6]+WET[24]*DRY[5]+WET[25]*DRY[4]+WET[26]*DRY[3]+WET[27]*DRY[2]+WET[28]*DRY[1],
>WET[1]*DRY[29]+WET[2]*DRY[28]+WET[3]*DRY[27]+WET[4]*DRY[26]+WET[5]*DRY[25]+WET[6]*DRY[24]+WET[7]*DRY[23]+WET[8]*DRY[22]+WET[9]*DRY[21]+WET[10]*DRY[20]+WET[11]*DRY[19]+WET[12]*DRY[18]+WET[13]*DRY[17]+WET[14]*DRY[16]+WET[15]*DRY[15]+WET[16]*DRY[14]+WET[17]*DRY[13]+WET[18]*DRY[12]+WET[19]*DRY[11]+WET[20]*DRY[10]+WET[21]*DRY[9]+WET[22]*DRY[8]+WET[23]*DRY[7]+WET[24]*DRY[6]+WET[25]*DRY[5]+WET[26]*DRY[4]+WET[27]*DRY[3]+WET[28]*DRY[2]+WET[29]*DRY[1]
>)
>
>
>
>
>Thank you for your time. 
>PS: if the above formula can not be seen correctly,
>one can see details at
>http://www.angelfire.com/ab5/get5/alt.ren.txt
>
>----------------------------------
>
>Mohammad Ehsanul Karim 
>
>Web: http://snipurl.com/ehsan 
>ISRT, University of Dhaka, BD
>
>______________________________________________
>R-help at stat.math.ethz.ch mailing list
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
>  
>

Seemingly Similar Threads

Search for more reasonably related threads

R help - Jun 2005 - r programming help

[R] r programming help

[R] r programming help

[R] r programming help

Seemingly Similar Threads