Displaying 2 results from an estimated 2 matches for "pol1".
Did you mean:
col1
2010 Jun 07
1
fit data with y = x^-1
Dear list,
I am getting weired with fitting data with a 1/x-polynomial. Suggest I have
the following data:
x <- c(1,2,3,4,5,6,7)
y <- c(100,20,4,2,1,.3,.1)
I may fit this with a linear model
fit1 = lm(y ~ I(x))
Getting plot out of this model I applied
library(polynom)
pol1 = polynomial(fit1$coefficients)
f1 = as.function(pol1)
plot(x,y)
lines(x, f1(x), col = 2)
Clearly, this model appears insufficient. Therefore I intended to aplly
fit2 = lm(y ~ I(1/x))
pol2 = polynomial(fit2$coefficients)
f2 = as.function(pol2)
lines(x, f2(x), col = 3)
as this would be much more...
2010 Oct 15
1
Problem using BRugs
...me figure out the problem and
hopefully get the code to work?
Many thanks for your help. --- Maomao
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
data<-list("y","cap2","pol2","cap1","pol1","g","wo","wd","ww","mu","tau")
inits<-function() {list(beta=beta0, rho_o=rho_o_0, rho_d=rho_d_0,
rho_w=rho_w_0)}
parameters<-c("beta", "rho_o", "rho_d", "rho_w")
probit.sim<-BRugsFi...