Displaying 1 result from an estimated 1 matches for "hatx".
Did you mean:
hat
2012 Apr 03
3
regression for poisson distributed data
...er in the mail the code with the model and visualisations I use to
check my results. I also copied the test dataset from my txt file.
I am using R 2.15 and Rstudio to visualise it.
plot(FR~N0)
x <-
nls(FR~(exp(i+j*N0)/(1+exp(i+j*N0)))*(k*N0/(k+N0)),start=list(i=0.02,j=0.002,k=6))
summary(x)
hatx <- predict(x)
lines(spline(N0,hatx))
N0 FR
10 2
10 3
10 2
10 4
10 2
10 2
10 1
10 2
10 2
10 2
20 2
20 3
20 3
20 3
20 4
20 2
20 4
20 2
20 3
20 2
30 1
30 2
30 3
30 4
30...