I am trying to figure out what is the best way to estimate beta why accounting for the uncertainty in x and y. For example, I have
x = rnorm(10,0,2) x.se = rnorm(10,0,0.7) y = 20*x y.se = rnorm(10,0,1) fit <- lm(y~x+0)
However, I would like to account for the standard errors of x and y. Ideally, I would like to estimate beta which takes into account this uncertainty. Any suggestions on how best to do this?