啊啊啊啊啊吖
2018-10-24
#转成时间序列类型 x = rnorm(2) charvec = c(“2010-01-01”,”2010-02-01”) zoo(x,as.Date(charvec)) #包zoo xts(x, as.Date(charvec)) #包xts timeSeries(x,as.Date(charvec)) #包timeSeries #规则的时间序列，数据在规定的时间间隔内出现 t
啊啊啊啊啊吖
2018-10-24
R与金融时间序列分析常见问题集 library(zoo) #时间格式预处理 library(xts) #同上 library(timeSeires) #同上 library(urca) #进行单位根检验 library(tseries) #arma模型 library(fUnitRoots) #进行单
啊啊啊啊啊吖
2018-10-24
现在有一个list。list里每一个元素都是一个时间序列，元素名称为该股票代码，时间序列为该股票在一段时间的数据。现在想把list的元素合并为一个整的时间序列，并在每观测值添加一个变量，代表该观测值对应的股票。 a lapply() dplyr::bind_rows()
啊啊啊啊啊吖
2018-10-24
如图1，这个里面显示有三个internal函数，convexMin,setupLambda和loss.ncvreg，但是我想调用convexMin的时候它说找不到， 如图2，上面help可以查到帮助文件，但是想看源码就找不到了。 图2 a 图1 a 磕磕碰碰请求学长终于找到一个解决办法，那就是在R-CRAN上下载用tar.gz封装的ncvreg包，解压之后在R文件夹里有convexM
啊啊啊啊啊吖
2018-10-24
> str(query_GPS) 'data.frame': 13150 obs. of 9 variables: $ id : chr "1010852" "1010852" "1010852" "1010852" ... $ 车牌号码: chr "鄂A3QD28" "鄂A3QD28" "鄂A3QD28" "鄂A3QD28" ... $ 速度 : ch
啊啊啊啊啊吖
2018-10-24
library(tseries) #加载算法包tseries library(forecast) #加载算法包forecast YEAR <- 2012 #设置初始时间年份 Frequency<- 12 #设置数据类型（月度数据为12，季度数据为4） p<-1 q<-1 d<-0 P<-2 Q<-0 D<-1 step<-4 d
啊啊啊啊啊吖
2018-10-24
假设： 1. 妹子们一生中可以遇到100个追求者，追求者的优秀程度符合正态分布； 2. 每个妹子都具备判断并比较追求者优秀程度的能力； 3. 接受或拒绝一个追求者后永远无法后悔。 问题：当遇到追求者时，如何选择才能获得最优结果？ 下面介绍选择方法： 首先，为了不错过在未来可以接受更优秀的追求者，理性的妹子会拒绝最早的一批追求者，并且采用第一批追求者做样本量k，理性地判断出追求者中最
啊啊啊啊啊吖
2018-10-24
SAS中Proc logisitc过程提供了很完善的logistic回归的分析功能，学习R中完成此过程只是想比较一下两个软件在完成此过程的差别。虽然有很多帖子介绍如何采用R完成logistic回归过程，但是都相对过于简单，对于以下常用细节很少涉及。 1、模型筛选方法 2、如何简单设定哑变量 3、针对分类变量，如何选取特定水平作为参考水平 4、如何简单输出OR值及置信区间 5、如何构建条件logis
啊啊啊啊啊吖
2018-10-24
读取文件是最基本的能力同时也是事贼多的问题，例如前几天有同学的读取R文件如下状况： > d<-read.table("C:\\Rdata\\cuel.txt",header=T) Error in file(file, "rt") : 无法打开链结 此外: Warning message: In file(file, "rt") : 无法打开文件'C:\Rdata\cuel.txt':
啊啊啊啊啊吖
2018-10-24
randomForest 原创者 Leo Breiman原创性论文集，及R中randomForest变种列表 很多人，用完randomForest ,觉得它很多不可思议的地方，因此想着改善。但是，没有成功的，大量所谓改善的"randomforest"变种，至今没有一个能超越Leo Breiman 原创！！！ 要用好randomForest ，请读读Leo Breiman的原创性论文，任
啊啊啊啊啊吖
2018-10-22
pie(x, labels = names(x), edges = 200, radius = 0.8, clockwise = FALSE, init.angle = if(clockwise) 90 else 0, density = NULL, angle = 45, col = NULL, border = NULL, lty = NULL, main = NULL, .
啊啊啊啊啊吖
2018-10-22
barpot：适用于数据表格化后的情形。 barplot(height, width = 1, space = NULL, names.arg = NULL, legend.text = NULL, beside = FALSE, horiz = FALSE, density = NULL, angle = 45, col = NULL, border =
啊啊啊啊啊吖
2018-10-22
R软件实用例子，咱高校学生做题的好题库，注重统计应用：http://www.stathome.cn/html/S-plus_R/Rrumen/ 各种关于R的讲义，身为老师的你真的不来参考学习一下吗：http://research.stowers-institute.org/efg/index.html R可视化的网站，你能想到的图都在这啦：http://addictedtor.free.fr/
啊啊啊啊啊吖
2018-10-22
箱线图函数：boxplot(x, ...) ## S3 method for class 'formula' boxplot(formula, data = NULL, ..., subset, na.action = NULL) ## Default S3 method: boxplot(x, ..., range = 1.5, width = NULL, varwidth = FALSE,
啊啊啊啊啊吖
2018-10-22
1）之前在网上查看关于时间序列分析的文章，在分析之前首先需要去掉非平稳的部分。这个非平稳的部分，包含趋势和季节性。那这样的话，季节性是不是就是没有研究价值的部分？ 2）是不是对本身就有季节性变化的数据，才能在预测时，在R的forecast函数部分指定关于季节性的参数？ 对本身就没有季节性变化的数据，在做预测时，就不能在R的forecast函数部分指定关于季节性的参数？ 3）R的auto.
啊啊啊啊啊吖
2018-10-22
R语言怎么使用mgcv包做面板数据的非参数广义可加模型的估计？ library(mgcv) set.seed(0) n<-400 sig2<-4 x0 <- runif(n, 0, 1) x1 <- runif(n, 0, 1) x2 <- runif(n, 0, 1) x3 <- runif(n, 0, 1) pi <- asin(1) * 2 f <- 2 * sin(pi
啊啊啊啊啊吖
2018-10-22
a TXT中是factor类型，里面有文字和数字，如“混2”，我想让“混”赋值到STRU中，“2”赋值到FLOOR中，如何实现呢？数字不是只有一位，可能一位也可能两位，不定长度的。 library("dplyr") library("stringr") myData %>% mutate(TXT = as.character(TXT), FLOOR = str_ex
啊啊啊啊啊吖
2018-10-19
> head(sam) long lat bin pro 12464 115.2391 33.91198 a 0 10608 112.5319 27.75921 a 0 2848 112.8600 25.70828 a 0 10433 108.7582 24.06754 a
啊啊啊啊啊吖
2018-10-19
首先保证json数据文件都在当前文件夹下 读取json数据 library(jsonlite)#加载需要的包（安装包需要先安装好） 批量导入当前文件夹下的所有文件 filelist <- list.files(pattern=".*.json")#列出当前文件夹下所有的json数据文件名（每个文件名以.json结尾） length(filelist)#文件个数 time1<-Sys.time()
啊啊啊啊啊吖
2018-10-18
