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]]
__
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?
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.
__
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
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:
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
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?
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
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
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
10 matches
Mail list logo