写入R代码以从N(0,1)和N(5,1)各自绘制大小为100的样本.在同一图表上绘制两个图形.将图形保存Sample.pdf在工作目录中I.
我的尝试:
pdf("SampleGraph.pdf",width=7,height=5)
x=rnorm(100)
y=rnorm(100,5,1)
plot(x,lty=2,lwd=2,col="red")
lines(y,lty=3,col="green")
dev.off()
Run Code Online (Sandbox Code Playgroud)
它不起作用.
MATLAB中是否有任何内置函数可供计算
二项分布的第0.025个分位数,参数n = 60,p = 0.4
0.975分位数的标准正态分布?
数据
x=[33 44 29; 16 25 45; 33 19 54; 22 21 49; 11 24 56];
Run Code Online (Sandbox Code Playgroud)
包含在三个不同的商店销售产品5天.因此该列表示商店1,2和3,行表示周一至周五.
现在我计算了每一天的汇总统计数据.
meanday=mean(x');
medianday=median(x');
minday=min(x');
maxday=max(x');
varianceday=var(x');
summaryday=[meanday;medianday;minday;maxday;varianceday];
Run Code Online (Sandbox Code Playgroud)
现在我想添加这样的标题:
summaryday =
Mean Median Minimum Maximum variance
Monday 32.3333 28.6667 35.3333 27.6667 30.3333
Tuesday 33.0000 25.0000 33.0000 22.0000 24.0000
Wednesday 20.0000 16.0000 19.0000 21.0000 11.0000
Thursday 44.0000 45.0000 54.0000 40.0000 56.0000
Friday 144.3333 220.3333 310.3333 114.3333 536.3333
Run Code Online (Sandbox Code Playgroud)
我该如何添加标题?
我想计算方差= 36的正态分布均值μ的最大似然估计。
但是以下过程显示错误:
set.seed(36)
x <- rnorm(50,mean=20,sd=6)
loglikelikelihood<- function(mu,x){
L = (-50)*log(sqrt(2*pi)*6)-((1/72)*(sum(x^2)-2*mu*sum(x)+50*mu^2))
}
optimize(f=loglik,X=x,interval=c(0,100),maximum=T)
Run Code Online (Sandbox Code Playgroud)