R语言实现多元线性回归

k=seq(0.1,0.9,0.1) #[1] 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
par(mfrow=c(3,3),mai=c(0.6,0.5,0.2,0.1)) #实现一页多图的功能
for(i in 1:9)
  + barplot(dbinom(0:5,5,k[i]),
            xlab = "x",ylab = "p",ylim = c(0,0.6),
            main = substitute(B(5,b),list(b=k[i])),
            col = "lightblue")
#没有不合格的概率
dbinom(0,5,0.06)


#恰有一个不合格的概率
dbinom(1,5,0.06)


#有三个及三个一下的不合格品概率
pbinom(3,5,0.06)
par(mfrow=c(1,1),mai=c(0.7,0.7,0.2,0.1),cex=0.8)
curve(dnorm(x,-2,1),
      from=-6,to=2,
      xlim = c(-6,6),ylab = "Density")#绘制μ=-2，σ=1正态曲线
abline(h=0)#绘制y=0的一条线
curve(dnorm(x,-2,1.5),
      from=-7,to=3,
      add=T,xlim = c(-6,6),
      col = "red")
curve(dnorm(x,2,1.5),
      from=-3,to=7,add=T,
      xlim = c(-6,6),col = "blue")
segments(-2,0,-2,
         dnorm(2,mean=2,sd=1),
         lty=2,col = "red")#绘制中线
#计算正态分布概率
#X~N(50,10^2),P(X<=40)
pnorm(40,mean = 50,sd=10)

#X~N(50,10^2),P(30<=X<=40)
pnorm(40,mean = 50,sd=10)-pnorm(30,mean = 50,sd=10)

#Z~N(0,1),P(Z<=2.5)
pnorm(2.5,mean = 0,sd=1)

#Z~N(0,1),P(-1.5<=Z<=2)
pnorm(2,mean = 0,sd=1)-pnorm(-1.5,mean = 0,sd=1)


#计算正态分布分位数
#累积概率为0.025时的反函数值z
qnorm(0.025,mean=0,sd=1) #左侧概率值

par(mfrow=c(2,3),mai=c(0.6,0.6,0.2,0.1))#分块
n=5000
df=c(2,5,10,15,20,30)
for(i in 1:6){
	x<-rchisq(n,df[i])
    hist(x,xlim=c(0,60),#设置不同自由度n
	prob=T,col='lightblue',
	xlab=expression(chi^2),ylab="Density",#xy图例
	main=paste("df=",df[i]))#标题
   curve(dchisq(x,df[i]),col=2,add=T)}
#计算概率
#自由度为15，卡方值小于10的概率
pchisq(10,df=15)

#自由度为25，卡方值大于15的概率
1-pchisq(15,df=25)


#计算分位数
#自由度为10，卡方分布右尾概率为0.05时的反函数值
#左值
qchisq(0.95,df=10)
par(mfrow=c(1,1),mai=c(0.8,0.8,0.1,0.2))
curve(dnorm(x,0,1),from=-3,to=3,xlim=c(-4,4),ylab="f(x)",lty=1,col=1)#黑色正态分布曲线
abline(h=0)#y=0水平线
segments(0,0,0,dnorm(0),col="blue",lty=2)#对称轴蓝虚线
curve(dt(x,2),from=-4,to=4,add=T,lty=3,col=2)#自由度为2的红线
curve(dt(x,5),from=-4,to=4,add=T,lty=4,col=3)#自由度为5的绿线
#计算概率
#自由度为10，t值小于-2的概率
 pt(-2,df=10)

#自由度为10，t值大于3的概率
 1-pt(3,df=15)


#计算分位数
#自由度为10，t分布双尾概率为0.025时的t值
 qt(0.975,df=25) #左值
curve(df(x,10,20),from=0,to=5,xlim=c(0,5),xlab="F",ylab="f(x)",lty=1,col=1)
curve(df(x,5,10),from=0,to=5,add=T,lty=2,col=2)
abline(h=0,v=0)
#计算概率
#分子自由度为10，分母自由度为8，F值小于3的概率
 pf(3,df1=10,df2=8)

#分子自由度为18，分母自由度为15，F值大于2.5的概率
 1- pf(2.5,df1=18,df2=15)


#计算分位数 #左值
#分子自由度为25，分母自由度为20，F分布累积概率为0.95时的F值
 qf(0.95,df1=25,df2=20)

Workspace restored from D:/sltj/.RDataFiles/1527078135.RData 
> k=seq(0.1,0.9,0.1)
> par(mfrow=c(3,3),mai=c(0.6,0.5,0.2,0.1))
> for(i in 1:9)
+   + barplot(dbinom(0:5,5,k[i]),
+             xlab = "x",ylab = "p",ylim = c(0,0.6),
+             main = substitute(B(5,b),list(b=k[i])),
+             col = "lightblue")
> dbinom(0,5,0.06)
[1] 0.733904
> dbinom(1,5,0.06)
[1] 0.2342247
> pbinom(3,5,0.06)
[1] 0.9999383
> par(mfrow=c(1,1),mai=c(0.7,0.7,0.2,0.1),cex=0.8)
> curve(dnorm(x,-2,1),
+       from=-6,to=2,
+       xlim = c(-6,6),ylab = "Density")
> abline(h=0)
> curve(dnorm(x,-2,1.5),
+       from=-7,to=3,
+       add=T,xlim = c(-6,6),
+       col = "red")
> curve(dnorm(x,2,1.5),
+       from=-3,to=7,add=T,
+       xlim = c(-6,6),col = "blue")
> segments(-2,0,-2,
+          dnorm(2,mean=2,sd=1),
+          lty=2,col = "red")
> pnorm(40,mean = 50,sd=10)
[1] 0.1586553
> pnorm(40,mean = 50,sd=10)-pnorm(30,mean = 50,sd=10)
[1] 0.1359051
> pnorm(2.5,mean = 0,sd=1)
[1] 0.9937903
> pnorm(2,mean = 0,sd=1)-pnorm(-1.5,mean = 0,sd=1)
[1] 0.9104427
> qnorm(0.025,mean=0,sd=1)
[1] -1.959964
> par(mfrow=c(2,3),mai=c(0.6,0.6,0.2,0.1))
> n=5000
> df=c(2,5,10,15,20,30)
> for(i in 1:6){
+ 	x<-rchisq(n,df[i])
+     hist(x,xlim=c(0,60),#设置不同自由度n
+ 	prob=T,col='lightblue',
+ 	xlab=expression(chi^2),ylab="Density",#xy图例
+ 	main=paste("df=",df[i]))#标题
+    curve(dchisq(x,df[i]),col=2,add=T)}
> pchisq(10,df=15)
[1] 0.1802601
> 1-pchisq(15,df=25)
[1] 0.9413826
> qchisq(0.95,df=10)
[1] 18.30704
> par(mfrow=c(1,1),mai=c(0.8,0.8,0.1,0.2))
> curve(dnorm(x,0,1),from=-3,to=3,xlim=c(-4,4),ylab="f(x)",lty=1,col=1)
> abline(h=0)
> segments(0,0,0,dnorm(0),col="blue",lty=2)
> curve(dt(x,2),from=-4,to=4,add=T,lty=3,col=2)
> curve(dt(x,5),from=-4,to=4,add=T,lty=4,col=3)
> pt(-2,df=10)
[1] 0.03669402
> 1-pt(3,df=15)
[1] 0.004486369
> qt(0.975,df=25)
[1] 2.059539
> curve(df(x,10,20),from=0,to=5,xlim=c(0,5),xlab="F",ylab="f(x)",lty=1,col=1)
> curve(df(x,5,10),from=0,to=5,add=T,lty=2,col=2)
> abline(h=0,v=0)
> pf(3,df1=10,df2=8)
[1] 0.9335491
> 1- pf(2.5,df1=18,df2=15)
[1] 0.03944963
> qf(0.95,df1=25,df2=20)
[1] 2.07392
> qf(0.95,df1=25,df2=20)
[1] 2.07392