Displaying 1 result from an estimated 1 matches for "koncentr".
Did you mean:
concentr
2010 Mar 22
1
maxNR - Error in p(a, b) : element 1 is empty; the part of the args list of '*' being evaluated was: (b, t)
...e version of R being used is "RGui for
Windows" if it helps to know this.
R-code below:
library(maxLik)
> require(maxLik)
>
> x <-
> read.table('http://www.math.ku.dk/kurser/2008-09/blok4/stat2/doku/data/Eksempel_6_3.txt',
> header = TRUE);
> t <- log(x$Koncentration);
> X <- x$Status;
>
> p <- function(a,b) exp(a+b*t)/(1+exp(a+b*t));
> S <- sum(X);
> SP <- sum(t*X);
>
> F <- function(a,b) {
+ c(sum(p(a,b)) - S,
+ sum(t*p(a,b)) - SP)
+ }
>
>
> z <- maxNR(F, start=1, print.level=2)
Error in p(a,...