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?