我正在尝试创建一个函数来在R中执行以下操作:
输入参数
n
,s
,x-bar
所有整数。
- 使用
n - 1
生成rnorm( n - 1 , 0 , 1)
个随机数。- 取出每个数字,将其平方,然后将它们相加并将总和存储到
x
中。将整数s
除以x
,并取结果的平方根 存入sigma
- 从N(0,1)生成单个值。
rnorm( 1 , 0 , 1)
并存储到Z
。- 然后将
mu
的值存储到Z x ( sigma / sqrt(n) ) + x-bar
。- 打印
醇>mu
。
任何帮助将不胜感激!
我的尝试(这是对的吗?)
g <- function( n , s , xb ) {
vals <- rnorm( n - 1 , 0 , 1)
sum <- 0
for( i in vals )
{
sum <- sum + (i * i )
print( sum )
}
sigma <- sqrt ( s / sum )
z <- rnorm(1, 0 , 1)
mu <- ( z * ( sigma / sqrt(n) ) + xb )
print ( mu)
}
答案 0 :(得分:0)
g <- function( n , s , xb ) {
vals <- rnorm( n - 1 , 0 , 1)
sum<- sum(vals^2)
sigma <- sqrt ( s / sum )
z <- rnorm(1, 0 , 1)
mu <- ( z * ( sigma / sqrt(n) ) + xb )
print ( mu)
}
g(10,6,3)
[1] 2.89561