编辑
这个答案解决的问题是找到对他们不喜欢的学生给予较少分数的年级生。
原始帖子
我认为很容易实现的方法如下:
μk,ikiyk,i
1个
假设模型
yk,i=μk,i+α+τek,iααiα
2
Giiy~k,i
yk,i−μk,i−α=y~k,i=Gi+σie~k,i
并对进行11个单独估计 Gσ
3
现在不寻常的观察是这样的数量
T=|y~−Giσi|
注意
eT
R码
下面是R中的代码。请注意,在您的情况下,将同时给出mu和y,因此应忽略为它们分配了rnorm-numbers时生成的行。我包括它们,以便能够在没有数据的情况下评估脚本。
mu_0 <- 50;
alpha <- 5;
tau<- 10;
# 0 Generate data for first assignment
mu <- matrix(rnorm(605, mu_0, tau), 11)
# 1 Generate data for second assignment and estimate alpha
G <- rnorm(11, 0)*10;
for(i in 1:11){
y[i,] <- rnorm(55, 0, sigma) + mu[i,] + alpha + G[i];
}
alpha_hat <- mean(y-mu)
alpha_hat
# 2 Form \tilde{y} and find unsual observations
ytilde <- y - mu - alpha_hat
T <- matrix(0, 11, 55);
for(i in 1:11){
G_hat <- mean(ytilde[i,]);
sigma_hat <- sd(ytilde[i,]);
T[i,] <- order(abs(ytilde[i,] - G_hat)/sigma_hat)
}
# 3 Look at grader number 2 by
T[2,]