我正在R中使用quickpsy软件包。我想有一个参数,该参数可缩放quickpsy为累积正态分布函数cum_normal_fun
计算的标准偏差。有没有办法做到这一点?我无法传递初始参数,因为quickpsy软件包显然正在使用概率分析来计算参数(均值和std dev),因此参数随每个数据集而变化。这是我第一次使用R,所以我可能会遗漏一些明显的东西。
fit <- quickpsy(avg, Gain, nYes, n, grouping = (ID), prob = levels, xmin = 0, xmax = 2, fun = cum_normal_fun)
我查看了源代码中的quickpsy,发现cum_normal_function只是调用pnorm()
,所以我尝试定义自己的函数,例如:
scaled_dev_fun <- function(x,p) { pnorm(x, p[1],p[2]) }
但我收到错误消息:
Error in fitpsy(d, x, k, n, random, within, between, grouping, xmin, xmax, :
parini (initial parameters) must be specified.
据我了解,我无法为parini
指定初始参数,因为这将使我无法完成通过概率分析来计算它们的快速检查。任何建议都非常感谢。