[R] any function for monotone nonparametric regression?

2006-06-16 Thread Cunningham Kerry
I am wondering if there is any package in R that can fit a nonparametric regression model with monotone constraints on the fitted results. - [[alternative HTML version deleted]] __

[R] looking for ENAR partner

2006-01-16 Thread Cunningham Kerry
Please forgive my abuse of this list. One of my friend is going to ENAR spring meeting and has booked air tickets. As a student with only limited support, he is looking for a partner to share the hotel costs. If you happen to know any information, could you please contact me at this email?

[R] solving a complicated equation

2005-11-06 Thread Cunningham Kerry
I want to solve the following equation for x p=a*exp(-x^2/2)+b*P(Zx) where p,a,b are known, Z is a standard normal variable. Clearly there is no analytic form for P(Zx). I am wondering if any expert could direct one easy way on this. Thank you. __

[R] question about sm.density

2005-10-27 Thread Cunningham Kerry
How can I draw a 95% contour in sm.density? For example, y - cbind(rnorm(50), rnorm(50)) sm.density(y, display = slice) will give 25%, 50% and 75% contours automatically, but no reference on other values. __ R-help@stat.math.ethz.ch mailing

[R] Fwd: predictive interval in nlme

2005-10-20 Thread Cunningham Kerry
I sent out this message a couple of days ago. Yet received no reply. There is a possibility that my description is not very clear or example is not highly specific. I am just wondering if anybody could point out (or have seen) a similar function that can do this job for lme model. Thanks. Note:

[R] predictive interval in nlme

2005-10-18 Thread Cunningham Kerry
Suppose I have the following data: y x id 44 0 104 48 58 104 48 55 204 47 105 204 41 275 206 18 67 209 ... I fit the model fit=lme(y~x+I(x^2),random=~1|id) Now I want to make a prediction plot: time=seq(0,300,len=100) plot(predict(fit,data.frame(x=time),level=0)) Very fine. It gives me

[R] how to keep very small or large number?

2005-09-22 Thread Cunningham Kerry
When I was computing some joint probabilities, I found that R reported most of the results to to -Inf and thus didn't record the value. I guess it is b/c the joint log(probability) can be extremely small. Is there a way in R to keep the values even if they are small?

Re: [R] how to keep very small or large number?

2005-09-22 Thread Cunningham Kerry
for(i in 1:n) { p - p*exp((y[i]*log(plogis(pta0+pta1*x1[i]+pta2*x2[i]))+(1-y[i])*log(1-plogis(pta0+pta1*x1[i]+pta2*x2[i] } but the result is shown to be zero. --- Cunningham Kerry [EMAIL PROTECTED] wrote: When I was computing some joint probabilities, I found that R reported most

[R] matrix calculation?

2005-09-13 Thread Cunningham Kerry
Hi, I am wondering if anybody has a simpler solution to calculate the following: ma a b [1,] 1 4 [2,] 2 3 [3,] 3 2 [4,] 4 5 pa [,1] [,2] [,3] [,4] a1234 b4325 diag(ma%*%pa) [1] 17 13 13 41 I only want to calculate the product of the row of the first

[R] how to merge several data sets?

2005-08-27 Thread Cunningham Kerry
Dear R-helpers, I want to merge several data sets into one single data set. For example, there are three separate data sets like: Set 1: id age gender 01 12 M 03 15 F 04 19 M ... Set 2: id time x1 01 1 0.25 01 2 0.27 01 3 0.29 03 1 0.15 03 2 0.18 04 2 0.22 04 3 0.54 ... Set 3: id