> Задача 4
> ...
> 2) Для параметров m=10, σ = 1, n=10 сгенерировать 100 выборок x1,...,xn. По каждой выборке вычислить оценку максимального правдоподобия параметров m и σ . Нанести полученные оценки в виде точек на плоскости, отметив их знаком ‘o’. Точку с координатами m=0, σ = 1 отметить на той же плоскости знаком Х.

R> X ← replicate(100, rnorm(10, 10, 1))
R> sdb ← function(x) sqrt(sum((x-mean(x))^2/length(x)))
R> d10 ← data.frame(m=apply(X, c(2), mean), s=apply(X, c(2), sdb))
R> X2 ← replicate(100, rnorm(300, 10, 1))
R> d300 ← data.frame(m=apply(X2, c(2), mean), s=apply(X2, c(2), sdb))
R> library(ggplot2)
R> plot1 ← ggplot(d10, aes(x=m, y=s)) + geom_point(shape="o") + geom_point(data=data.frame(m=10, s=1), shape="x", size=4, color="blue") + opts(title="n=10")
R> plot1
R> plot2 ← ggplot(d300, aes(x=m, y=s)) + geom_point(shape="o") + geom_point(data=data.frame(m=10, s=1), shape="x", size=4, color="blue") + opts(title="n=300")
R> plot2
R> plot1
R> plot1 + xlim(9,11) + ylim(0.5,1.5)
R> plot2 + xlim(9,11) + ylim(0.5,1.5)
R> postscript("Документи/ШАД/червоненкис_задание1.plot1.ps", paper="special", horizontal=F, width=5.417, height=5.417)
R> plot1 + xlim(9,11) + ylim(0.5,1.5)
R> dev.off()
R> postscript("Документи/ШАД/червоненкис_задание1.plot2.ps", paper="special", horizontal=F, width=5.417, height=5.417)
R> plot2 + xlim(9,11) + ylim(0.5,1.5)
R> dev.off()

Обожаю R и ggplot, и благодарен чуваку, который мне их показал и принудительно им научил. По своей воле я бы это говно никогда не полез даже пробовать. А так ок, няшнота.

#ozfefe

add comment recommend bookmark subscribe