Próbuję narysować linię regresji najmniejszych kwadratów przy użyciu abline(lm(...))
, która również musi przejść przez określony punkt. Widzę, że this question jest powiązany, ale nie do końca tego, co chcę. Oto przykład:R draw (abline + lm) linia najlepszego dopasowania przez dowolny punkt
test <- structure(list(x = c(0, 9, 27, 40, 52, 59, 76), y = c(50, 68,
79, 186, 175, 271, 281)), .Names = c("x", "y"))
# set up an example plot
plot(test,pch=19,ylim=c(0,300),
panel.first=abline(h=c(0,50),v=c(0,10),lty=3,col="gray"))
# standard line of best fit - black line
abline(lm(y ~ x, data=test))
# force through [0,0] - blue line
abline(lm(y ~ x + 0, data=test), col="blue")
ta wygląda następująco:
Teraz jak bym go o zmuszanie linię przez oznaczonego dowolnego punktu (x=10,y=50)
jednocześnie minimalizując dystans do innych punktów?
# force through [10,50] - red line
??
Nice. Wydaje się, że nie jest to szorstkie rozwiązanie. – thelatemail
Szorstki? Nie z komputerowego punktu widzenia, ale raczej statystycznego. Znajdziesz kilka rants na ten temat przez Billa Venablesa (jak z Venables/Ripley MASS) –
To właśnie miałem na myśli. – mnel