[R] "taking the log then later exponentiate the result" query

2009-04-12 Thread Mary Winter
Hi, I am trying to figure out the observed acceptance rate and M, using generalised rejection sampling to generate a sample from the posterior distribution for p. I have been told my code doesn't work because I need to "take the log of the expression for M, evaluate it and then exponentia

[R] Generalised Rejection Sampling

2009-04-12 Thread Mary Winter
Hi, I am trying to figure out the observed acceptance rate and M, using generalised rejection sampling to generate a sample from the posterior distribution for p. I have been told my code doesn't work because I need to "take the log of the expression for M, evaluate it and then exponent

[R] (no subject)

2009-04-11 Thread Mary Winter
Hi, Is it possible to increase just the text size of the main title on my graph? If so, how?! Along the same lines..along the x axis at a lot of the labels overlap each other, is it possible to make this text smaller of write it on two lines so it doesn't take up as much space along

Re: [R] Frequency table to histogram

2009-04-07 Thread Mary Winter
de a copy of the first few lines of your frequency > table (with only relevant columns)? I think this could help many > potential contributors understand in more details what you're trying to do. > > Cheers, > > Luc > > Mary Winter wrote: > > > > > &

[R] Frequency table to histogram

2009-04-07 Thread Mary Winter
I have read a frequency table in to R called "temp." I now want to create a histogram table from it, but I obviously first have to expand the data - to the sample size of 100. I want to use the command rep(), but I'm not sure how to go about it..I tried using the code: temp1<-rep(temp$C

Re: [R] Residual resampling for non linear reg model

2008-05-04 Thread Mary Winter
Since sending this message I have now solved the problem - needed to alter the initial values of alpha and beta! > From: [EMAIL PROTECTED]> To: r-help@r-project.org> Date: Mon, 5 May 2008 > 00:03:53 +0100> Subject: [R] Residual resampling for non linear reg model> > > > > > I was attempting t

[R] Residual resampling for non linear reg model

2008-05-04 Thread Mary Winter
I was attempting to use the residual resampling approach to generate 999 bootstrap samples of alpha and beta and find their confidence intervals. However, I keep getting the error message:Error in nls(resample.mp ~ cases/(alpha + (beta * cases)), start = init.values, : singular gra