对于bootstrap的一些粗浅认识

彭友松 2011-5-30 22:11

什么是bootstrap？ bootstrap就是从一个原始样本中进行有放回的重复采样，采样大小和原始样本大小相同，采样次数根据计算量而定。从每个重新样的样本中可以计算某个统计量的bootstrap 分布，比如说均值，多个重采样样本的均值构成了原始样本均值的bootstrap分布。在采样完后需要检查待研究统计量的bootstrap分布是不是符合 ...

个人分类: 统计学与R语言学习|23522 次阅读|没有评论

R语言中的机器学习-分类树和回归树

热度 1 彭友松 2011-5-18 16:31

在R中包tree专门用于分类树和回归树的构建，其中的函数tree()和线性回归的函数lm()类似，只需要指定函数形式和数据。如： data - read.table("test") model - tree(V1~.-1,data) 这里需要注意的是，如果公式中的Y是数值变量，那么tree()默认为回归树，否则是分类树。 包中的函数prune.tree()用于修建树，predic ...

个人分类: 统计学与R语言学习|18996 次阅读|1 个评论 热度 1

对于流行病学建模的粗浅认识-1

热度 2 彭友松 2011-5-12 23:28

之前没有关心过流行病学方面，一直以为就是折腾几个微分方程，没有什么意思。现在有一个想法需要用到流行病学建模，才想去关注一下这个领域。不看不知道，一看吓一跳，竟然有这么多人在做流行病学，而且各种各样的模型都有，好晕好晕，不禁感叹这水真深，有点不敢趟了。不过水深归水深，还是想多了解一下在我这个外行看来处 ...

个人分类: 流行病学建模|6550 次阅读|1 个评论 热度 2

R语言中的最大似然估计

彭友松 2011-5-12 16:37

对于最大似然估计，都是依赖于似然函数的，因此关键在于写出似然函数，然后对该似然函数进行优化。似然函数依个人的问题而定，比如说下面的是正态分布的似然函数： normal - function(theta,x){ mu - theta sigma2 - theta n - length(x) logL - -0.5*n*log(2*pi)- ...

个人分类: 统计学与R语言学习|31277 次阅读|没有评论

关于回归的一些问题

彭友松 2011-4-22 17:52

1 如何理解回归得到的R-squared和残差标准误？ 下面是一个简单的回归模型，给出了模型的上面三个量，如粗体所示。 Call: lm(formula = data ~ data + data ) Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) -124.626 92 ...

个人分类: 统计学与R语言学习|4806 次阅读|没有评论

R语言的图像输出和结果输出

彭友松 2011-4-22 15:19

如果需要保存图像文件，我现在一直在用的图像输出方式是打开一个图像设备，如png("文件名")类似的方式，其他图像格式类似，如tiff格式的可以用tiff("文件名"),画完图之后需要把它关掉，用dev.off()命令。下面是一个例子： x - 1:10 y - 1:10 png("test.png") plot(x,y) dev.off() 如果需要保存R程序的输出结 ...

个人分类: 统计学与R语言学习|21702 次阅读|没有评论

R中package的安装

热度 1 彭友松 2011-3-28 14:58

之前一直用install.packages("packageName")安装R包的，最近文件服务器扩容之后发现这样安装出现了问题，窗口出现的警告信息： 无法在貯藏處/tjjiang/pengyousong/R/rattle_2.6.5.tar.gz/src/contrib中读写索引 我也不知道该怎么解决，于是只有寻找新的方法。在参考了网上的一些博客之后，现在总结下面两种方法： 1 从 ...

个人分类: 统计学与R语言学习|41448 次阅读|1 个评论 热度 1

Windows中的copy命令

彭友松 2010-11-18 22:25

在Linux中Cat命令可以把两个文件file1和file2连接起来，通过管道输出到另外一个文件file3，格式如下： file3 但是在Windows中没有对应的Cat命令，在网上找了好久，终于发现原来copy命令可以完成这个功能。copy命令可以把一个文件复制到另外一个文件上，这是大家都知道的，但是它也能把多个文件连接起来，只要在 ...

个人分类: Linux|6088 次阅读|没有评论

在windows中如何用perl调用可执行文件

彭友松 2010-11-18 22:04

假设test是可执行文件，在目录D：/work/里面，下面是其中可能的三种方式： 1exec('D：/work/test');用绝对路径.exec调用程序会使当前的perl程序终止执行，但是它好像它不能从可执行文件跳出来 2直接用``运行，即`test` 3用system()命令，即system('test'); 4fork(linux环境)和pr ...

个人分类: Linux|6488 次阅读|没有评论

本页有 1 篇博文因作者的隐私设置或未通过审核而隐藏