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 -1 to 1. So, I
wrote the following.
f1-function(theta,j){
aa-seq(-1,1,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