Re: [R] r programming help

From: Peter Wolf <s-plus_at_wiwi.uni-bielefeld.de>
Date: Thu 23 Jun 2005 - 01:53:03 EST

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@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
>
>



R-help@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 Received on Thu Jun 23 01:56:37 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:57 EST