Tôi thực hiện đo của hai biến x và y . Cả hai đều có sự không chắc chắn đã biết σ x và σ y liên quan đến chúng. Tôi muốn tìm mối quan hệ giữa x và y . Tôi làm nó như thế nào?
EDIT : mỗi có khác nhau σ x , i liên kết với nó, và cùng với các y i .
Ví dụ R sinh sản:
## pick some real x and y values
true_x <- 1:100
true_y <- 2*true_x+1
## pick the uncertainty on them
sigma_x <- runif(length(true_x), 1, 10) # 10
sigma_y <- runif(length(true_y), 1, 15) # 15
## perturb both x and y with noise
noisy_x <- rnorm(length(true_x), true_x, sigma_x)
noisy_y <- rnorm(length(true_y), true_y, sigma_y)
## make a plot
plot(NA, xlab="x", ylab="y",
xlim=range(noisy_x-sigma_x, noisy_x+sigma_x),
ylim=range(noisy_y-sigma_y, noisy_y+sigma_y))
arrows(noisy_x, noisy_y-sigma_y,
noisy_x, noisy_y+sigma_y,
length=0, angle=90, code=3, col="darkgray")
arrows(noisy_x-sigma_x, noisy_y,
noisy_x+sigma_x, noisy_y,
length=0, angle=90, code=3, col="darkgray")
points(noisy_y ~ noisy_x)
## fit a line
mdl <- lm(noisy_y ~ noisy_x)
abline(mdl)
## show confidence interval around line
newXs <- seq(-100, 200, 1)
prd <- predict(mdl, newdata=data.frame(noisy_x=newXs),
interval=c('confidence'), level=0.99, type='response')
lines(newXs, prd[,2], col='black', lty=3)
lines(newXs, prd[,3], col='black', lty=3)
Vấn đề với ví dụ này là tôi nghĩ rằng nó giả định rằng không có sự không chắc chắn trong . Làm thế nào tôi có thể sửa lỗi này?
Deming
hàm trong gói R MethComp .
lm
phù hợp với mô hình hồi quy tuyến tính, đó là: mô hình kỳ vọng của đối với P ( Y | X ) , trong đó rõ ràng Y là ngẫu nhiên và X được coi là đã biết. Để đối phó với sự không chắc chắn trong X, bạn sẽ cần một mô hình khác.