Julia中的二项式GLM-如何指定命中和未命中

Julia中的二项式GLM-如何指定命中和未命中,julia,glm,Julia,Glm,在Julia中,我想用族二项式()和LogitLink()计算GLM。我的数据是三个线性关系:xvalues,点击次数,以及未命中次数。我想通过它们在x轴上的位置来解释二元分布的命中和未命中。我有多个具有相同x坐标的样本(因为数据最初来自被展平的2D数组) 在R中,我必须在两列矩阵中提供命中和未命中。类似于以下工作: glm1虽然将数据集膨胀到一个~100k行的数据帧并不好看,但它确实让它工作起来了。要使用以下代码,请首先将数据集加载到xVAL、点击次数和未点击次数(如问题中的链接所示),然后:

在Julia中,我想用族
二项式()
LogitLink()
计算GLM。我的数据是三个线性关系:
xvalues
点击次数
,以及
未命中次数
。我想通过它们在x轴上的位置来解释二元分布的命中和未命中。我有多个具有相同x坐标的样本(因为数据最初来自被展平的2D数组)

在R中,我必须在两列矩阵中提供命中和未命中。类似于以下工作:


glm1虽然将数据集膨胀到一个~100k行的数据帧并不好看,但它确实让它工作起来了。要使用以下代码,请首先将数据集加载到
xVAL
点击次数
未点击次数
(如问题中的链接所示),然后:

# spreading dataset to one row per trial...   
data = DataFrame(
    xvals = vcat(rep(xvals,hits),rep(xvals,misses)), 
    outcome = vcat(rep(1,sum(hits)),rep(0,sum(misses))))

glm1 = glm(outcome ~ xvals, data, Binomial(),LogitLink())

我粗略地看了一下,结果似乎与数据相符。另请注意,
二项式
逻辑链接
是位置参数,而不是命名参数。

酷,有效。不过我还是希望能有更优雅的东西。但也许优雅还没有实现。。。