R中的极大似然预计 1. 数据与模子 我们要利用的数据来自于“MASS”包中的geyser数据。先把数据调出来，看看它长什么样子。 > geyser     waiting  duration 1        80 4.0166667 2        71 2.1500000 3        57 4.0000000 4        80 4.0000000 5        ... 继续阅读R中的极大似然预计

R语言利器之ddply ddply()函数位于plyr包，用于对data.frame举办分组统计，与tapply有些雷同 # 利用stringsAsFactors=F来防备data.frame把向量转为factor apache = data.frame(httpCode=c(200,200,200,404,404,500), time=c(100,... 继续阅读R语言利器之ddply

use R for fun系列之小应用建造篇 系列以use R for fun为主题，以COS论坛上的英华帖、相关的package以及本身的一些code为素材，团结自身的一些编程体会，从而整合成文。本文是第二篇小应用建造篇。*本文素材出处均已在正文注明本文继承承接上一篇的话题(小游戏开拓篇)，继承在交互操纵上做文章，差异的是这里引入了更富厚的操纵和idea，仅仅做些小游戏还远远达不到我们的胃口... 继续阅读use R for fun系列之小应用建造篇

R语言：抓交通闯祸犯 一辆卡车违反交通法则，撞人后逃跑。现场有三人眼见事件，但都没有记着车号，只记下车号的一些特征。甲说：牌照的前两位数字是沟通的；乙说：牌照的后两位数字是沟通的，但与前两位差异； 丙是数学家，他说：四位的车号恰好是一个整数的平方。请按照以上线索求出车号。 x=1:9 y=0:9 foo = function(x,y) if(x!=y) x*1000+x*100+y*10+y... 继续阅读R语言：抓交通闯祸犯

R语言：新娘和新郞 三对情侣介入婚礼，三个新郞为A、B、C，三个新娘为X、Y、Z。有人不知道谁和谁成婚，于是询问了六位新人中的三位，但听到的答复是这样的：A说他将和X成婚；X说她的未婚夫是C；C说他将和Z成婚。这人听后知道他们在恶作剧，全是谎言。请编程找出谁将和谁成婚。 a <- expand.grid(rep(list(1:3),3)) names(a) <- c(‘x’,... 继续阅读R语言：新娘和新郞

Windows下R的编辑器（RWinEdt） 作为一个建模人员，并且平时的事情情况就是Windows，预计不大大概迁移到Linux平台，也不想记那么多的Emacs下的快捷键，因此，RWinEdt才是我的归宿。RWinEdt和R Console协同的方法很像当年的GAUSS，是我以为最舒服的一种建模调算法的方法，在一个文本编辑器中写好大段代码，看那边不爽就用鼠标点已往修改，对付差异的实现方法还可... 继续阅读Windows下R的编辑器（RWinEdt）

主身分阐明的步调及R代码 拿到原始数据X=[x1,x2,…xp]； 个中xi 是列向量（1..N）;这就是说X是N×P的矩阵，有N个样本，P个指标 接下来，要做的是尺度化处理惩罚，因为差异的指标大概4量纲可能数据区间差异等等 （1）xij*=xij-mean(xj)/var(xj)(尺度差); 个中xij*是尺度化的数据值，xj暗示第j列数据； （2）计较相关矩阵 （3）计较特征值和特征矩阵 ... 继续阅读主身分阐明的步调及R代码

R语言绘制维恩图 Gplots包 > library(gplots) Loading required package: gtools Loading required package: gdata Attaching package: ‘gdata’ The following object(s) are masked from ‘package:utils’: object.siz... 继续阅读R语言绘制维恩图

R 常用统计检讨要领 R语言：常用统计检讨要领<一> 正态总体均值的假设检讨t检讨单个总体例一某种元件的寿命X（小时），听从正态漫衍，N（mu,sigma^2），个中mu,sigma^2均未知，16只元件的寿命如下：问是否有来由认为元件的平均寿命大于255小时。呼吁：X<-c(159, 280, 101, 212, 224, 379, 179, 264,222, 362, 1... 继续阅读R 常用统计检讨要领

R语言-将数值的矩阵映射为一个颜色方格的矩阵 data(mtcars)　　fit = lm(mpg~., mtcars)　　cor = summary(fit, correlation = TRUE)$correlation　　#相关阵上下倒转再转置，是为了让画出的图和相关阵偏向一致　　cor2 = t(cor[11:1,])　　#下面颜色取自ellipse包中的plotcorr函数的示例　　... 继续阅读R语言-将数值的矩阵映射为一个颜色方格的矩阵