Sorry, I replied to the wrong email.  Here it is again:

Try this where g is f summed over j and k for given scalars
theta and rho and gv is g vectorized over theta.  I have not
checked this carefully so be sure you do:

f <- function(theta = 0, rho = 0, j = 0, k = 0)
       dnorm(theta+2*pi*j,0,1)*pnorm(2*pi*(k+1)-rho*(theta+2*pi*j))
g <- function(theta = 0, j = 0, k = 0)
       sum(sapply(k, function(k) sum(f(theta, j = j, k = k))))
gv <- function(theta = 0, ...) sapply(theta, g, ...)

x <- seq(0, pi, length = 100)
plot(x, gv(x, j = -1:1, k = -1:1))

On 4/9/06, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> Hi, thanks for your reply.
> Here I would like to ask you again more directly.
> The following is what I had for now.
> The function to begin with is
> dnorm(theta+2*pi*j,0,1)*(pnorm(((2*pi*(k+1)-rho*
> (theta+2*pi*j)).
>
> Now, I wanted to sum it over k from -10000 to 10000. So, I
> wrote the following.
> f1<-function(theta,j){
> aa<-seq(-10000,10000,lenght=20001)
> sum(dnorm(theta+2*pi*j,0,1)*(pnorm(((2*pi*(aa+1)-rho*
> (theta+2*pi*j))
> }
>
> Next, I would like sum over j as
> (This is a way I thouhgt about computing a double sum)
> f2<-function(theta){
> aa<-c(-1,0,1)
> sum(f(theta,aa))
> }
> expecting f2 be a function of theta.
>
> Then, I wanted to plot f2 for theta from 0 to 2*pi. But, f2
> (seq(0,2*pi,length=10)) is a scalar, not a vector with the
> dimension of 10x1. I would like to ask you how I can use
> rowSums with outer or some other method to correct my coding.
>
> I would like to thank you in advance.
> Sungsu.
> UC riverside.
>
> ______________________________________________
> 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

Reply via email to