Hi all,
I would like to fit a gamma pdf to my data using the method of RSS (Residual
Sum-of-Squares). Here are the data:
x <- c(86, 90, 94, 98, 102, 106, 110, 114, 118, 122, 126, 130, 134, 138,
142, 146, 150, 154, 158, 162, 166, 170, 174)
y <- c(2, 5, 10, 17, 26, 60, 94, 128, 137, 128, 77, 68, 65, 60, 51, 26, 17, 9,
5, 2, 3, 7, 3)
I have typed the following code, using nls method:
fit <- nls(y ~ (1/((s^a)*gamma(a))*x^(a-1)*exp(-x/s)), start = c(s=3, a=75,
x=86))
But I have the following message error (sorry, this is in German):
Fehler in qr(.swts * attr(rhs, "gradient")) :
Dimensionen [Produkt 3] passen nicht zur Länge des Objektes [23]
Zusätzlich: Warnmeldung:
In .swts * attr(rhs, "gradient") : Länge des längeren Objektes
ist kein Vielfaches der Länge des kürzeren Objektes
Could anyone help me with the code?
I would greatly appreciate it.
Sincerely yours,
Vincent Laperrière.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
