- 博客(7)
- 收藏
- 关注
RSS订阅原创 抽样（蒙特卡洛法）到底是在干啥？
最近在学马尔可夫蒙特卡洛法，期间学习了各种各样的抽样方法，也在《概率图模型基于R》这本书中看到了一些代码实现，但是说实话我都还没搞清楚抽样到底是在干什么。昨天在B站关注到一位UP主感觉讲的很不错，这里按照上图的分布我解释一下抽样到底是在干什么。因为我的统计基础薄弱，所以很多基础概念都不扎实。之前我看f(x)，就觉得这就是一个函数，抽样是不是就是从这个函数上取点？随机抽样是不是就是随机在这个函数上取点？那比如说上面这个函数的定义域是（-6，6），那我直接在（-6，...
2022-05-05 08:11:51
4534
1
原创 统计学习方法 EM算法的一个例子
参考《统计学习方法》第九章三硬币模型 假设有3枚硬币，分别记作A,B,C.这些硬币正面出现的概率分别是π，p和q·进行如下掷硬币试验：先掷硬币A,根据其结果选出硬币B或硬币C,正面选硬币B,反面选硬币C:然后掷选出的硬币，掷硬币的结果，出现正面记作1，出现反面记作0：独立地重复n次试验（这里，n=10),观测结果如下：1,1,0,1,0,0,1,0,1,1假设只能观测到掷硬币的结果，不能观测掷硬币的过程.问如何估计三硬币正面出现的概率，即三硬币模型的参数。题中一...
2022-04-20 14:54:51
1226
原创 R apply()、lapply()，sapply()
目录一、 apply()二、lapply() 和 sapply()首先我们给出三个函数的应用区别：apply() 作用在矩阵上 lapply() 作用在向量或列表上，返回列表 sapply() 作用在向量或列表上，返回向量一、 apply()> x<-array(c(1:24),c(3,4))> x [,1] [,2] [,3] [,4][1,] 1 4 7 10[2,] 2 5 8...
2022-04-19 12:58:48
304
原创 R ddply()
在学习《概率图模型基于R》的过程中遇到了ddply()，而书中讲解较少，本文内容作为理解书中代码做的补充。也就是服务于这篇博客。若有错误之处，希望大家在留言区指出。简单地说，ddply(data,variables,fun)就是对一个表格使用函数 在函数作用前先将表格继续分组首先我们建一个简单的表格student = data.frame(name=c("a","b","c","d"), years=c(12,13,14,12),
2022-04-18 21:08:10
3929
2
原创 R adj
无向图按R文档中的例子来理解这个函数。set.seed(123)gR3 <- randomGraph(LETTERS[1:4], M<-1:2, p=.5)plot(gR3)我们构造的图是：adj(gR3, "A") #就是输出在图gR3中与节点A相连的节点$A[1] "B" "D"adj(gR3, c(2,3)) #分别输出在图gR3中与第2个和第3个节点相连的节点$B[1] "A" "D"$C[1] "D"后来在有向图中我发现.
2022-04-17 19:11:30
249
原创 R graphNEL
graphNEL有nodes、edgeL和edgemod参数，其中nodes：节点名称edgeL：该节点的子节点edgemod：边的类型用两个例子说明这个函数例一data0 <- data.frame(x=c(1,2),y=c(3,4),z=c(5,6))#这里edges=2指的是第一个节点要连向第二个节点edges0 <- list(x=list(edges=2),y=list(edges=3),z=list()) g0 <- graphNEL(
2022-04-17 16:34:15
308
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人