Tracy数据-CSDN博客

转载 RStudio创建R程序包

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

2017-02-17 16:32:31 6337

原创 R-Garch异方差

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

2014-08-18 18:39:36 1805

原创 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 18936 9

原创 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 56504 3

转载用R检验配对股票的协整性

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

2014-08-12 14:41:30 5642

转载 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 2110

转载揭开R语言中环境空间的神秘面纱

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

2014-07-03 22:48:25 4780 1

转载 [概率统计]—中心极限定理

什么是中心极限定理中心极限定理Central Limit Theorem：设从均值为μ、方差为σ^2;（有限）的任意一个总体中抽取样本量为n的样本，当n充分大时，样本均值的抽样分布近似服从均值为μ、方差为σ^2/n的正态分布。注意：原来的分布不一定要符合正态分布，可以是任何的分布，可以是离散也可以是连续，即无要求。n为sample size，每次取n个样本，每次样本的mean，分别

2014-03-10 19:02:42 18361

原创 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 2488

转载 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 1934 1

原创 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 3461

原创 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 3536

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 3385

原创 R-常见操作符

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

2013-12-21 23:24:23 2922

原创 R-$和@符号的意义

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

2013-12-21 22:35:01 27201

转载 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 3219

原创 R-数据类型及特殊值判断

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

2013-12-14 18:28:16 8909

原创 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 29178 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 13683

原创 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 26296

原创 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 9202

原创 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 71564 4

原创 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 51456

原创 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 11082 1

转载用source函数代替繁冗的R语言打包过程

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

2013-10-31 16:33:51 18455

原创 R-数据保存和输出

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

2013-10-31 16:15:44 2727

原创 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 5755

原创 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 14025

原创 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 12193

原创 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 10343

u012543538的专栏