R语言自定义极大似然函数函数、假设检验函数
目录
一、自定义极大似然函数二、自定义假设检验函数
目录
一、自定义极大似然函数
1、求出似然函数,以正太分布、指数分布为例 正太分布似然函数: 指数分布似然函数(假设x>0): 2、函数代码: 正太分布求极大似然函数的函数代码
norm.fun
n
N=qnorm(a/2) # 分位数
area1=xbar-(sigma/n)**0.5*N
area2=xbar+(sigma/n)**0.5*N
area=c(area1,area2)
# 上面是由公式写出置信区间
NT=(xbar-mu)/(sigma/n)**0.5 # 计算N检验值
TF='reject' # 先默认TF 为拒绝
if ( abs(NT)
N=qnorm(a)
area1=xbar+(sigma/n)**0.5*N
area2=Inf # 无穷
area=c(area1,area2)
NT=(xbar-mu)/(sigma/n)**0.5
TF='reject'
if (NT>N)
{
TF='accept'
}
result=list('accept_or_reject'=TF,'interval estimation'=area)
}
if(alternative=='greater')
{
N=-qnorm(a) #N应该是a的上分位,对下分位求负数是相等的
area1=-Inf
area2=xbar-(sigma/n)**0.5*N
area=c(area1,area2)
NT=(xbar-mu)/(sigma/n)**0.5
TF='reject'
if (NT |