u012543538的专栏

R语言中，$表示从一个dataframe中取出某一列数据@是从R的类实例里面读取数据，bg=x@colors$bg.col就是从对象实例x中取出colors，而这个colors本身又是个dataframe，所以需要进一步用$读取bg.col列。

Slots are linked to S4 objects. A slot can be seen as a part, element or a "property" of an object. Say you have a car object, then you can have the slots "price", "number of doors", "type of engine",

R语言常见操作符             ‘:: :::’           访问命名空间变量       ‘$ @’              访问组分或者插槽       ‘[ [[’             索引       ‘^’                指数操作符，从左至右       ‘- +’              减和加       ‘:’

getCandleChart<-function(x,theme='white'){ x<-na.omit(x) xax<-x$DATETIME print(paste("xax[1]=",xax[1])) Opens <- as.numeric(x$Open) Highs <- as.numeric(x$High)

getCandleChart x<-na.omit(x) xax<-x$DATETIME print(paste("xax[1]=",xax[1])) Opens <- as.numeric(x$Open) Highs <- as.numeric(x$High) Lows <- as.numeric(x$Low)

豆瓣搜索想要的电影名并返回电影评分，走起library(RCurl)library(XML)movieScore <- function(x) { stopifnot(is.character(x)) # 提交搜索豆瓣表单 search <- getForm("http://movie.douban.com/subject_search", search_text

R网页爬取步骤library(XML)#第一步：读取url并解析成R结构url1=“http:。。。。"url=htmlParse(url1,encoding="UTF-8")#第二步：找网页节点test#第三步：读取节点内容test_text_listtest_text#第四步：读取节点属性content1content1

平时做R开发时，每次启动R环境都会要加载很多package，甚是繁琐。最近学习《R in action》时，介绍了自定义R启动环境的方法。windows环境下，R启动时会到R_Home\etc目录下找Rprofile.site文件，其中“R_Home”指的是R安装目录，例如c:\R。可以用notepad等文本编辑器打开c:\R\etc\Rprofile.site，进行修改 # Thin

前言环境空间(environment)对于大部分的R使用者来说，都是比较陌生的。虽然我们不了解它的运行原理，但也不影响我们使用R语言。环境空间是R语言中关于计算机方面的底层设计，主要用于R语言是环境加载器。通过环境空间，封装了加载器的运行过程，让使用者在不知道底层细节的情况下，可以任意加载使用到的第三方的R语言程序包。本文将揭开R语言中环境空间的神秘面纱。目录R

基于统计套利的配对交易策略是一种市场中性策略。具体地说，是指从市场上找出历史股价走势相近的股票进行配对，当配对股票价格差（Spread）偏离历史均值时，则做空股价偏高的股票，同时做多股价偏低的股票，等待它们回归到长期均衡关系，由此赚取两股票价格收敛的报酬。 进行配对交易，第一步也是最关键的一步是寻找符合配对条件的股票，即两支历史价格走势相近，具有长期稳定关系的股票。本文解释如何用R来实现协

当p，q确定（参见上篇文章）以后，即可进行

关于自相关、偏自相关：一、自协方差和自相关系数      p阶自回归AR(p)      自协方差 r(t,s)=E[X(t)-EX(t)][X(s)-EX(s)]      自相关系数ACF=r(s,t)/[(DX(t).DX(s))^0.5]二、平稳时间序列自协方差与自相关系数      1、平稳时间序列可以定义r(k)为时间序列的延迟k自协方差函数：      r

RCode;#ARMA Modeling寻找AIC值最小的p,qarmaSearch     armacoef     for (p in 0:5){         for (q in 0:5) {             #data.arma = arima(diff(data), order = c(p, 0, q))             data.arm

R语言中各种分布函数举例：curve(dnorm(x,m=10,sd=2),from=0,to=20,main="Normal distribution")curve(dgamma(x, scale=1.5, shape=2),from=0, to=15, main="Gamma distribution")curve(dweibull(x,

R中获取数据类型信息的一些有用函数R语言的对象“类”很多，虽然我们不可能一一去详细学习，但接触到一类新数据时我们需要了解一些基本信息才能进行进一步的操作。R提供了一些非常有用的方法（函数）。getClass( )函数我们前面已经见过了，它的参数是表示类的字符串。class( )可获取一个数据对象所属的类，它的参数是对象名称。str( )可获取数据对象的结构组成，这很有用。mod

遇到问题：(pending)write.csv(a,"a.csv",append=T)提示：Warning message:In write.csv(a,"a.csv",  :  不能修改'append'于是改存txt格式，write.csv(a,"a.txt",append=T) #a.txt已经存在并有数据Warning message:In wri

write.table族(a,file="data.txt",sep=",",quote=F,col.names=T)，但导出来的txt文件第一列总是R标记的行数，1，2，3，4.....而且，当用excel或文本打开时，由于多了一列，会造成列名对应错误。解决办法：这是R自动添加的行号，在保存数据write.table()中加入row.names=F即可

R函数中如何返回多个对象的问题：You can return a vector, or a list. 你可以返回一个变量或者一个list 点击打开链接举例：func   ...   result   return(result)  # you can omit this }   if you omit the "return(result)"

1. 如果用write.table(....,.csv,....)中要加上sep=",",否则excel表中数据会不分列，writ.csv则不用加此2. 保存数据时记得+ row.names=F,否则会自动添加一列序号列

用source函数代替繁冗的R语言打包过程 熊荣川 xiong rongchuan 六盘水师范学院生物信息学实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz 经过初级的学习和使用R语言之后我们渐渐的开始动手写自己的R语言小程序，这些小程序因为和自己的工作非常契合而变得通用性

#目标：输入一个日期，需要判断是否周末/或者周几package(timeDate)#报错： x x[1] "2013-11-15" library(timeDate) dayOfWeek(x)错误: 不是所有的inherits(x, "timeDate")都是TRUE修正：timeDate()要求的数据格式必须是GMT格式的，所以要先对.Date

R中有三个排序函数：order(),sort(),rank()   sort(x)是对向量x进行排序，返回值排序后的数值向量。   rank()是求秩的函数，它的返回值是这个向量中对应元素的“排名”。   而order()的返回值是对应“排名”的元素所在向量中的位置。注意三者返回值的差异。Example：> x> order(x)#x中排第2的数字在原数列中的位

实践举例：> recordBuyPosition_ru1405_1021[1]$`20625`      TIMENO            DATETIME              BUY1AMT cumBUY1AMT [1,] "201310210915001" "2013-10-21 09:15:00" "3"     "3"        [2,] "20131021

R中取百分位比用quantile()1. 求某个百分位比> x> quantile(x,0.9)     90% 1.377796 > quantile(x,c(0.1,0.9))      10%       90% -1.317574  1.377796 2. 产生一序列百分位比值> quantile(x,seq(0.1,1,0.1))   

R-找出两列数中的相同数通常在用交叉判断，比如满足A条件的点以及同时满足B条件的点a b c 要找出a,b,c中的相同数字解决办法一：intersect(intersect(a,b),c)解决办法二：Reduce(intersect, list(a,b,c))

xts()1. 所属包：xts> install.packages("xts")> library(xts)载入需要的程辑包：zoo载入程辑包：‘zoo’下列对象被屏蔽了from ‘package:base’:    as.Date, as.Date.numericWarning message:程辑包‘xts’是用R版本3.0.2 来建造的 

之前写过一篇文章是往文本文件中不断写入数据时用到了append=T,今天用到append()问题：如何往list某列末尾添加数据？直接用append(某list，要写入的数据）append(x, values, after = length(x))Argumentsxthe vector to be modified.val

R空间映像相关操作#保存save.image() #存为.RData格式，文件夹为getwd()对应的文件夹save(d, file = "c:/data/foo.Rdata") #存到指定文件夹和指定空间名下#加载load(".RData") load("dumData.Rdata") #加载#删除方法一：file.remove('.RData')方法

RStudio Server的安装及配置，请参考：多人在线协作R开发RStudio ServerRStudio已经为我们准备了一个有界面的环境，方便我们制作R包。创建一个工程sayHello编辑DESCRIPTION创建sayHello.R的脚本编辑sayHello.Rd的文档1). 创建一个工程sayHello2). 编辑rstudio/sayHell