- 博客(30)
- 收藏
- 关注
RSS订阅转载 RStudio创建R程序包
RStudio Server的安装及配置，请参考：多人在线协作R开发RStudio ServerRStudio已经为我们准备了一个有界面的环境，方便我们制作R包。创建一个工程sayHello编辑DESCRIPTION创建sayHello.R的脚本编辑sayHello.Rd的文档1). 创建一个工程sayHello2). 编辑rstudio/sayHell
2017-02-17 16:32:31
6430
原创 R-ARMA(p,q)如何选找最小AIC的p,q值
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
2014-08-18 17:41:56
19151
7
原创 R-时间序列自相关acf，偏自相关pacf
关于自相关、偏自相关：一、自协方差和自相关系数 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
2014-08-18 17:31:12
56857
3
转载 用R检验配对股票的协整性
基于统计套利的配对交易策略是一种市场中性策略。具体地说，是指从市场上找出历史股价走势相近的股票进行配对，当配对股票价格差（Spread）偏离历史均值时，则做空股价偏高的股票，同时做多股价偏低的股票，等待它们回归到长期均衡关系，由此赚取两股票价格收敛的报酬。 进行配对交易，第一步也是最关键的一步是寻找符合配对条件的股票，即两支历史价格走势相近，具有长期稳定关系的股票。本文解释如何用R来实现协
2014-08-12 14:41:30
5788
转载 R语言自定义启动环境
平时做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
2014-07-04 11:14:40
2170
转载 揭开R语言中环境空间的神秘面纱
前言环境空间(environment)对于大部分的R使用者来说，都是比较陌生的。虽然我们不了解它的运行原理，但也不影响我们使用R语言。环境空间是R语言中关于计算机方面的底层设计，主要用于R语言是环境加载器。通过环境空间，封装了加载器的运行过程，让使用者在不知道底层细节的情况下，可以任意加载使用到的第三方的R语言程序包。本文将揭开R语言中环境空间的神秘面纱。目录R
2014-07-03 22:48:25
5102
1
转载 [概率统计]—中心极限定理
什么是中心极限定理中心极限定理Central Limit Theorem：设从均值为μ、方差为σ^2;（有限）的任意一个总体中抽取样本量为n的样本，当n充分大时，样本均值的抽样分布近似服从均值为μ、方差为σ^2/n的正态分布。注意：原来的分布不一定要符合正态分布，可以是任何的分布，可以是离散也可以是连续，即无要求。n为sample size，每次取n个样本，每次样本的mean，分别
2014-03-10 19:02:42
18942
原创 R-web爬取一般步骤
R网页爬取步骤library(XML)#第一步：读取url并解析成R结构url1=“http:。。。。"url=htmlParse(url1,encoding="UTF-8")#第二步：找网页节点test#第三步：读取节点内容test_text_listtest_text#第四步：读取节点属性content1content1
2014-01-20 16:20:10
2544
转载 R语言—豆瓣搜索电影
豆瓣搜索想要的电影名并返回电影评分，走起library(RCurl)library(XML)movieScore <- function(x) { stopifnot(is.character(x)) # 提交搜索豆瓣表单 search <- getForm("http://movie.douban.com/subject_search", search_text
2014-01-20 15:51:52
2000
原创 R-自画蜡烛图R CandleChart
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)
2014-01-06 17:15:31
3680
原创 R-自画蜡烛图CandleChart
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)
2014-01-06 17:13:09
3860
原创 R-什么是slot插槽？
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",
2013-12-23 00:31:47
3576
原创 R-常见操作符
R语言常见操作符 ‘:: :::’ 访问命名空间变量 ‘$ @’ 访问组分或者插槽 ‘[ [[’ 索引 ‘^’ 指数操作符，从左至右 ‘- +’ 减和加 ‘:’
2013-12-21 23:24:23
3135
原创 R-$和@符号的意义
R语言中，$表示从一个dataframe中取出某一列数据@是从R的类实例里面读取数据，bg=x@colors$bg.col就是从对象实例x中取出colors，而这个colors本身又是个dataframe，所以需要进一步用$读取bg.col列。
2013-12-21 22:35:01
27951
转载 R-各种分布函数
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,
2013-12-19 13:19:48
3346
原创 R-数据类型及特殊值判断
R中获取数据类型信息的一些有用函数R语言的对象“类”很多，虽然我们不可能一一去详细学习，但接触到一类新数据时我们需要了解一些基本信息才能进行进一步的操作。R提供了一些非常有用的方法（函数）。getClass( )函数我们前面已经见过了，它的参数是表示类的字符串。class( )可获取一个数据对象所属的类，它的参数是对象名称。str( )可获取数据对象的结构组成，这很有用。mod
2013-12-14 18:28:16
8999
原创 R-找出两个数列中的相同元素
R-找出两列数中的相同数通常在用交叉判断，比如满足A条件的点以及同时满足B条件的点a b c 要找出a,b,c中的相同数字解决办法一：intersect(intersect(a,b),c)解决办法二：Reduce(intersect, list(a,b,c))
2013-12-09 15:39:57
29927
3
原创 R-xts()时间序列
xts()1. 所属包：xts> install.packages("xts")> library(xts)载入需要的程辑包：zoo载入程辑包：‘zoo’下列对象被屏蔽了from ‘package:base’: as.Date, as.Date.numericWarning message:程辑包‘xts’是用R版本3.0.2 来建造的
2013-12-02 22:45:48
13862
原创 R-quantile()
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))
2013-11-29 15:51:15
26483
原创 R-list类型直接转data.frame()之后是factor类型，直接按行列位置取数会出错
实践举例：> recordBuyPosition_ru1405_1021[1]$`20625` TIMENO DATETIME BUY1AMT cumBUY1AMT [1,] "201310210915001" "2013-10-21 09:15:00" "3" "3" [2,] "20131021
2013-11-26 11:02:52
9277
原创 R-往一个list某列中写入新的数据append()，查找列名 names(sapply()及list相关操作
之前写过一篇文章是往文本文件中不断写入数据时用到了append=T,今天用到append()问题：如何往list某列末尾添加数据？直接用append(某list，要写入的数据）append(x, values, after = length(x))Argumentsxthe vector to be modified.val
2013-11-18 23:50:01
72557
1
原创 R-排序order(),sort(),rank()
R中有三个排序函数：order(),sort(),rank() sort(x)是对向量x进行排序，返回值排序后的数值向量。 rank()是求秩的函数，它的返回值是这个向量中对应元素的“排名”。 而order()的返回值是对应“排名”的元素所在向量中的位置。注意三者返回值的差异。Example：> x> order(x)#x中排第2的数字在原数列中的位
2013-11-15 15:37:48
52645
原创 R-timeDate() 根据日期查星期几
#目标：输入一个日期，需要判断是否周末/或者周几package(timeDate)#报错： x x[1] "2013-11-15" library(timeDate) dayOfWeek(x)错误: 不是所有的inherits(x, "timeDate")都是TRUE修正：timeDate()要求的数据格式必须是GMT格式的，所以要先对.Date
2013-11-15 14:29:45
11531
1
转载 用source函数代替繁冗的R语言打包过程
用source函数代替繁冗的R语言打包过程 熊荣川 xiong rongchuan 六盘水师范学院生物信息学实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz 经过初级的学习和使用R语言之后我们渐渐的开始动手写自己的R语言小程序，这些小程序因为和自己的工作非常契合而变得通用性
2013-10-31 16:33:51
18876
原创 R-数据保存和输出
1. 如果用write.table(....,.csv,....)中要加上sep=",",否则excel表中数据会不分列，writ.csv则不用加此2. 保存数据时记得+ row.names=F,否则会自动添加一列序号列
2013-10-31 16:15:44
2867
原创 R-write.table族存出数据多一列
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即可
2013-10-28 22:28:27
5958
原创 R-数据写入，不覆盖 append=T，以及数据存出后多出一列
遇到问题：(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
2013-10-27 23:21:47
14276
原创 R-让函数返回多个对象
R函数中如何返回多个对象的问题：You can return a vector, or a list. 你可以返回一个变量或者一个list 点击打开链接举例：func ... result return(result) # you can omit this } if you omit the "return(result)"
2013-10-23 17:21:38
12529
原创 R-workspace相关操作
R空间映像相关操作#保存save.image() #存为.RData格式，文件夹为getwd()对应的文件夹save(d, file = "c:/data/foo.Rdata") #存到指定文件夹和指定空间名下#加载load(".RData") load("dumData.Rdata") #加载#删除方法一：file.remove('.RData')方法
2013-10-23 10:37:48
10703
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人