Displaying 1 result from an estimated 1 matches for "r_manuf".
Did you mean:
manuf
2007 Mar 28
0
nlsystemfit: Errors with reproducing the manual example
...hg
The following object(s) are masked from data ( position 3 ) :
ba cr dbh dg elev hg smi tht
Error in array(x, c(length(x), 1), if (!is.null(names(x))) list(names(x), :
dim<- : invalid first argument
/
2) When I try to do something similar on my own data (called R_manuf and
containing 3 columns named va98, ts98 and gom98) with the code
library(systemfit)
dati <- read.table("R_manuf",header=T,na.strings="nan")
attach(dati)
y1 <- va98/ts98 - 1;
y2 <- gom98 - va98;
eq1.formula <- y1.98 ~ a1*ts98^(b1-1) ;
eq2.formula &...