threeAFC(sensR)
threeAFC()所属R语言包:sensR
Create 3-AFC binomial family
创建3-AFC二项式族
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Creates af copy of the binomial family with the inverse link function changed to equal the 3-AFC psychometric function and correspondingly changed link function and derivative of the inverse link function.
创建AF副本的二项式家庭的改变等于3-AFC的心理功能,并相应地改变了链接功能及衍生的逆联接函数的反向链接功能。
用法----------Usage----------
threeAFC()
值----------Value----------
A binomial family object for models. Among other things it inludes the psychometric function as linkinv and the inverse psychometric function (for direct dprime computation) as linkfun.
一个的二项式家庭对象的模型。除其他事项外里包括linkinv和的逆心理测量的功能(的直接dprime计算)linkfun的心理功能。
注意----------Note----------
Several functions in this package makes use of the function, but it may also be used on its own—see the example below.
在此包中的几个功能,使用该功能,但它也可以用来在其自己的例子下面。
(作者)----------Author(s)----------
Rune Haubo B Christensen and Per Bruun Brockhoff
参考文献----------References----------
models for sensory discrimination tests as generalized linear models.
参见----------See Also----------
triangle, twoAFC, duotrio, discrim, discrimPwr, discrimSim, AnotA, discrimSS,
triangle,twoAFC,duotrio,discrim,discrimPwr,discrimSim,AnotA,discrimSS,
实例----------Examples----------
xt <- matrix(c(10, 5), ncol = 2) # data: 10 correct answers, 5 incorrect[数据:10个正确答案,不正确]
res <- glm(xt ~ 1, family=threeAFC)
summary(res)
## Extended example plotting the profile likelihood[#扩展的例子策划的档案可能性]
## data: 10 correct answers, 5 incorrect[#数据:10个正确答案,5不正确]
xt <- matrix(c(10, 2), ncol = 2)
summary(res <- glm(xt ~ 1, family = threeAFC))#, etastart = etastart))[etastart = etastart))]
N <- 100
dev <- double(N)
level <- c(0.95, 0.99)
delta <- seq(1e-4, 5, length = N)
for(i in 1:N)
dev[i] <- glm(xt ~ -1 + offset(delta[i]),
family = threeAFC)$deviance
plot(delta, exp(-dev/2), type = "l",
xlab = expression(delta),
ylab = "Normalized Profile Likelihood")
## Add Normal approximation:[#添加师范大学近似:]
lines(delta, exp(-(delta - coef(res))^2 /
(2 * vcov(res))), lty = 2)
## Add confidence limits:[#添加的置信区间:]
lim <- sapply(level, function(x)
exp(-qchisq(x, df=1)/2) )
abline(h = lim, col = "grey")
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|