inst/examples/R/stat/SchaetzungRegressionsparameter/Punktwolke.R

In sigbertklinke/shinyExample: Create and run simple (interactive) programs as Shiny apps

par(mar=c(2,2,2,1))
par(bg="white")
c = 1.4

X = c(1.47212467127845, 1.43803816908291, 1.39816119472708, 1.35265431871237, 1.30170078116712, 1.24550575400328, 1.18429551476157, 1.11831653547164, 1.04783449019622, 0.973133185255426, 0.894513416438906, 0.812291757807323, 0.726799286960364, 0.638380251904096, 0.547390684885741, 0.454196968777484, 0.359174361781949, 0.262705486399871, 0.165178788744325, 0.06698697440533, -0.03147457283696, -0.12980938308319, -0.227621496759157, -0.324517059005263, -0.42010590559103, -0.514003133968743, -0.605830653140166, -0.695218706095549, -0.781807358694595, -0.865247948994213, -0.945204491187104, -1.02135502849802, -1.09339292959007, -1.16102812326088, -1.22398826645696, -1.28201984090304, -1.33488917393073, -1.38238337939581, -1.42431121489567, -1.46050385183491, -1.49081555523854, -1.51512427057529, -1.53333211522811, -1.54536577263291, -1.55117678749854, -1.5507417609191, -1.54406244459308, -1.5311657337699, -1.51210355895223, -1.48695267679019, -1.4558143610094, -1.41881399461743, -1.37610056503069, -1.3278460641547, -1.27424479583335, -1.21551259345594, -1.15188595087228, -1.08362107011542, -1.01099282976651, -0.934293678115755, -0.853832455576358, -0.7699331510932, -0.682933597553698, -0.593184111454005, -0.501046082298121, -0.406890517409916, -0.311096548017573, -0.214049902625953, -0.116141353824037, -0.0177651447816395, 0.0806825982286982, 0.178805460892388, 0.276208337073337, 0.372499019766904, 0.467289780379078, 0.560198929972681, 0.650852356194005, 0.73888502969122, 0.823942473958736, 0.905682192688981, 0.983775048884155, 1.05790659017475, 1.12777831500824, 1.19310887460953, 1.2536352058731, 1.30911359062536, 1.35932063699166, 1.40405417891654, 1.44313409021509, 1.47640300987754, 1.50372697570645, 1.52499596373525, 1.54012433125592, 1.54905116167191, 1.55174050978779, 1.54818154654791, 1.53838860264111, 1.52240111079614, 1.50028344699998, 1.47212467127845)

Y = c(1.47212467127845, 1.50028344699998, 1.52240111079614, 1.53838860264111, 1.54818154654791, 1.55174050978779, 1.54905116167191, 1.54012433125592, 1.52499596373525, 1.50372697570645, 1.47640300987754, 1.44313409021509, 1.40405417891654, 1.35932063699166, 1.30911359062536, 1.2536352058731, 1.19310887460953, 1.12777831500824, 1.05790659017475, 0.983775048884155, 0.905682192688982, 0.823942473958736, 0.73888502969122, 0.650852356194005, 0.560198929972681, 0.467289780379077, 0.372499019766904, 0.276208337073337, 0.178805460892388, 0.0806825982286976, -0.0177651447816391, -0.116141353824037, -0.214049902625953, -0.311096548017574, -0.406890517409916, -0.501046082298121, -0.593184111454004, -0.682933597553698, -0.7699331510932, -0.853832455576358, -0.934293678115755, -1.01099282976651, -1.08362107011542, -1.15188595087228, -1.21551259345594, -1.27424479583335, -1.3278460641547, -1.37610056503069, -1.41881399461743, -1.4558143610094, -1.48695267679019, -1.51210355895223, -1.5311657337699, -1.54406244459308, -1.5507417609191, -1.55117678749854, -1.54536577263291, -1.53333211522811, -1.51512427057529, -1.49081555523854, -1.46050385183491, -1.42431121489567, -1.38238337939581, -1.33488917393073, -1.28201984090304, -1.22398826645696, -1.16102812326088, -1.09339292959007, -1.02135502849802, -0.945204491187104, -0.865247948994213, -0.781807358694595, -0.695218706095549, -0.605830653140167, -0.514003133968742, -0.42010590559103, -0.324517059005264, -0.227621496759157, -0.12980938308319, -0.03147457283696, 0.0669869744053296, 0.165178788744326, 0.262705486399871, 0.359174361781949, 0.454196968777483, 0.547390684885742, 0.638380251904096, 0.726799286960364, 0.812291757807323, 0.894513416438907, 0.973133185255426, 1.04783449019622, 1.11831653547164, 1.18429551476157, 1.24550575400328, 1.30170078116712, 1.35265431871237, 1.39816119472708, 1.43803816908291, 1.47212467127845)

plot(X, Y, type='n', lwd=2, xaxt="n", yaxt="n", xlim=c(-3,3), ylim=c(-3,3), axes=F, xlab="", ylab="")
polygon(X, Y, col="lightblue")

arrows(-3, -3, -3, 3, angle=20, length=0.2)
arrows(-3, -3, 3, -3, angle=20, length=0.2)

mtext(expression(bar("x")), side=1, line=0, at=0, cex=c)
mtext(expression(bar("y")), side=2, line=0, at=0, las=2, cex=c)
mtext("X", side=1, line=0, at=3, cex=c)
mtext("Y", side=2, line=0, at=3, las=2, cex=c)

set.seed(1)
points(rnorm(20,0,0.25), rnorm(20,0,0.25), cex=0.5, pch=19)
points(rnorm(10,-1,0.2), rnorm(10,-1,0.2), cex=0.5, pch=19)
points(rnorm(20,1,0.2), rnorm(20,1,0.2), cex=0.5, pch=19)
points(rnorm(20,-0.23,0.2), rnorm(20,-0.23,0.2), cex=0.5, pch=19)
points(rnorm(20,0.23,0.2), rnorm(20,0.23,0.2), cex=0.5, pch=19)
points(rnorm(20,0.6,0.1), rnorm(20,0.6,0.1), cex=0.5, pch=19)
points(rnorm(20,-0.4,0.1), rnorm(20,-0.5,0.1), cex=0.5, pch=19)

lines(c(0,0), c(-3,3), lwd=2)
lines(c(-3,3), c(0,0), lwd=2)
lines(c(-1,1), c(-3,3))
lines(c(-2,2), c(-3,3))
lines(c(-3,3), c(-3,3))
lines(c(-3,4), c(-2.2,3))