[R] Strange problem.

[Ko-Kang Kevin Wang]

On Mon, 1 Sep 2003, Yannong.Dong-1 wrote:

> >std<-1000
> >mean<-8000
> >prior<-function(n){1/(sqrt(2*pi)*std)*exp(-1.0*(n-mean)^2/(2*std^2))}
> >plot(prior,1,15000)
> >post<-function(n){
> + if(n < 9543) 
> +  0
> + else
> +  prior(n)/n
> + }
> >plot(post,1,15000)  # This command didn't work. The error message is " x and y lengths differ "

Try replacing the post() function with something like:
post<-function(n){
 ifelse(n < 9543, 0, prior(n) / n)
}

-- 
Cheers,

Kevin

------------------------------------------------------------------------------
"On two occasions, I have been asked [by members of Parliament],
'Pray, Mr. Babbage, if you put into the machine wrong figures, will
the right answers come out?' I am not able to rightly apprehend the
kind of confusion of ideas that could provoke such a question."

-- Charles Babbage (1791-1871) 
---- From Computer Stupidities: http://rinkworks.com/stupid/

--
Ko-Kang Kevin Wang
Master of Science (MSc) Student
SLC Tutor and Lab Demonstrator
Department of Statistics
University of Auckland
New Zealand
Homepage: http://www.stat.auckland.ac.nz/~kwan022
Ph: 373-7599
    x88475 (City)
    x88480 (Tamaki)