分类：r语言
2021
05-20
2021
05-20
数据分组，根据数据分析对象的特征，按照一定的数值指标，把数据分析对象划分为不同的区间部分来研究，以揭示内在的联系和规律性；在R中，我们常用ifelse函数来进行数据的分组，跟excel中的if函数是同一种用法。ifelse(condition,TRUE,FALSE)>data<-read.table('1.csv',sep='|',header=TRUE);>>level<-ifelse(+data$cost<=20,"(0,20]",+ifelse(+data$cost<=40,"(20,40]",+ifelse(+data$cost<=60...
继续阅读 >
可以先自定义函数也可以用的时候再定义。>mat<-matrix(c(1:3,7:9,4:6),byrow=T,nc=3)>mat[,1][,2][,3][1,]123[2,]789[3,]456>apply(mat,2,function(x){order(x,decreasing=T)[1]})#查找每一列[1]222>apply(mat,1,function(x){order(x,decreasing=T)[1]})#查找每一行[1]333>apply(mat,1,function(x){which.max(x)})#查找每一行...
继续阅读 >
2021
05-19
看代码吧~par(mar=c(5,5,3,4)+0.1)#似乎是设置图片位置bar<-barplot(gu[1:22,6],xlim=c(0.5,26),ylim=c(0,200000),ylab="交易量",col="blue",col.axis="blue",col.lab="blue")mtext(c(1:22),side=1,line=1,at=bar,cex=0.8,col="black")mtext("time",side=1,line=3,col="black")par(new=T)plot(bar,gu[1:22,2],axes=F,xlim=c(0.5,26),ylim=c(17,23),xlab="",ylab="",col="green",type="o")#ylim设置不好的话就会看不见折线ax...
继续阅读 >
R语言中，for循环运行比较慢for(iin1:1000){print(i^2)}补充：R语言：for循环使用小结基本结构展示：vals=c(5,6,7)for(vinvals){print(v)}#即把大括号里的内容对vals里的每一个值都循环run一遍实例展示：1.paste()命令是把几个字符连接起来如paste("A","B","C",sep="")得到的就是“ABC”，在次基础上写如下forloop:partnumber=c(1,2,5,78)for(iinpartnumber){print(paste("participantnumber",i,sep=""))}...
继续阅读 >
2021
05-15
在导入大批量数据时，如果没有显式地指定“stringsAsFactors=FALSE”，默认会将所有的字符串转换为因子，导致数据处理速度较慢。示例数据如下：name,math,english,sex,year"yiifaa",65,68,"M",2018"yiifee",95,98,"F",2018"guagua",75,78,"M",2018"MM",85,88,"F",2018查看数据概要，发现默认将字符串转换为因子，并进行了分组计数(这也是处理速度较慢的原因之一)概要如下：namemathenglishsexye...
继续阅读 >
PHP 代码分享
MongoDB 教程
正则表达式
jQuery 教程
JavaScript 代码片段
PHP7 教程
首先，让我们来进入常规步骤我安装的是recharts包，正常的写法呢，就是以下这个样子：install.packages("devtools")#devtools::install_github("madlogos/recharts")第一个问题：然而对于今天的我来说，那就太天真了，首先踏入的第一个坑：无法打开URL'http://yihui.name/xran/src/contrib/PACKAGES'Warningininstall.packages:从zip文件中抽取1时出了错类似于无法访问的错误，第一步的解决办法，就是改变包的镜像。在Tools-...
继续阅读 >
AR(p)模型与MA(q)实际上是ARMA(p,q)模型的特例。它们都统称为ARMA模型，而ARMA（p,q）模型的统计性质也是AR（p）与MA（q）模型的统计性质的有机组合。平稳系列建模假如某个观察值序列通过序列预处理可以判定为平稳非白噪声序列，就可以利用ARMA模型对序列建模。1.求出该观察值序列的样本自相关系数（ACF）与偏相关系数（PACF的值。2.根据根样本自相关系数和偏自相关系数的性质，选择阶数适当的ARMA（p,q）模型进行拟合。3.估计模型...
继续阅读 >
fivenum()函数：返回五个数据：最小值、下四分位数数、中位数、上四分位数、最大值对于奇数个数字=5，fivenum()先排序，依次返回最小值、下四分位数、中位数、上四分位数、最大值>fivenum(c(1,12,40,23,13))[1]112132340对于奇数个数字>5，fivenum()先排序，我们可以求取最小值，最大值，中位数。在排序中，最小值与中位数中间，若为奇数，取其中位数为下四分位数，若为偶数，取最中间两个数的平均值为下四分位数；在排序中...
继续阅读 >
卡方检验在数据统计中，卡方检验是一种很重要的方法。通常卡方检验的应用主要为：1、卡方拟合优度检验2、卡方独立性检验本文主要通过使用自己编程的方法实现相关检验。卡方拟合优度检验理论：1、我们先做出0假设：H0：总体服从假定的理论分布2、我们再构造一个统计量：3、当n充分大时4、我们得到该拒绝域代码#Chi_squareGoodnessOfFitTest#函数说明：#n为所得样本数据；p为理论概率#alpha为置信水平，df为自由度cgoft<-func...
继续阅读 >
看代码吧~#求单正态均值mu的置信区间#参数依次为置信水平alpha，正态样本x，已知总体方差（默认为未知）mu<-function(alpha,x,sigma=NA){n<-length(x)meanx<-mean(x)if(is.na(sigma)){t1<-qt(1-alpha/2,n-1)t2<-qt(1-alpha,n-1)mu11<-meanx-t1*sqrt(sum((x-meanx)^2)/(n-1))/sqrt(n)mu12<-meanx+t1*sqrt(sum((x-meanx)^2)/(n-1))/sqrt(n)mu21<-meanx+t2*sqrt(sum((x-meanx)^2)/(n...
继续阅读 >
查询百度之后，发现在R上一次显示多张图的函数很多，比如layout()或者分屏函数，但是这些都不是我想要的结果。之后，发现了line()函数可以保留原来图片继续作图，在括号中填入所需画图的部分即可（我是将一个矩阵作图）。先运行plot()函数，再注释掉plot()函数运行line()函数即可。效果如下：补充：R语言：在同一张图作不同曲线R语言将两条曲线作在同一张图的方法是library(ggplot2)year<-c(1993,1998,2003,2008)Res<-c(0.0227,0.0...
继续阅读 >
看代码吧~hist(data$MEDV,col="grey")#可以看到目标数据的总体分布legend(35,170,c("min:5","median:21.2","mean:22.51","max:50"),fill=c("lightblue","wheat"))补充：R语言作图——histogram(直方图)最近小仙同学很是烦恼，本以为自己已经掌握了ggplot2作图的语法，用read.csv(),ggplot()+geom_point()/boxplot()/violinplot()…就可以画遍天下图表，结果却发现到真正画图的时候，还是会出现不少的小问题。比如小仙最近要画一个...
继续阅读 >
程序实在是调不出来了，我决定破釜沉舟，直接把所有表格都打印出来，看看数据到底哪儿有问题。然后就开始了闹心的矩阵导出。。。首先，百度了一下，数据导出的代码为：write.table(x,file="",sep="",row.names=TRUE,col.names=TRUE,quote=TRUE)其中：x：需要导出的数据file：导出的文件路径sep：分隔符，默认为空格（""），也就是以空格为分割列row.names：是否导出行序号，默认为TRUE，也就是导出行序号col.names：是...
继续阅读 >