在R中重现先前的密度图

Question

我仍然习惯于将公式转换为R中的代码。我正在尝试通过Computer Age Statistical Inference中的经验贝叶斯章节（第6章）。我想制作图6.4，但为了做到这一点，我需要找到数据函数的边际分布。

为了获得数据并绘制它，我做了以下工作。

myscript.R

然后他尝试以下模型

通过以下计算

要获得我想要重现的情节。

我习惯使用标准R函数进行推理，例如nodes <- read.table("https://web.stanford.edu/~hastie/CASI_files/DATA/nodes.txt", header = T) nodes %>% ggplot(aes(x=x/n))+ geom_histogram(bins = 30)+ theme_bw()+ labs(x = "nodes", n = "p=x/n") 和。{ lm()，但现在我需要最大限度地找到一个边缘分布，我迷失了。