![](https://img-blog.csdnimg.cn/20201014180756930.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Rcode
TOMOCAT
知乎可私信: TOMOCAT
展开
-
[Rcode]管道操作
原文地址:http://blog.csdn.net/fairewell/article/details/72878107例子:取10000个随机数符合,符合正态分布。求这个10000个数的绝对值,同时乘以50。把结果组成一个100*100列的方阵。计算方阵中每行的均值,并四舍五入保留到整数。把结果除以7求余数,并话出余数的直方图。法一:#转载 2017-11-28 20:05:54 · 19492 阅读 · 0 评论 -
[Rcode]R语言数据统计量处理
1、提取百分比数值,不含百分号(unname)2、构造一个任意长度的向量用于给循环赋值 x<-vector(mode="numeric",length=0) for(i in 1:100) { x[i]<-i } 构造已知长度的向量:x<-vector(length=n)...原创 2020-01-28 00:43:47 · 610 阅读 · 0 评论 -
[Rcode]三种判别分析比较及R代码
判别方法的比较:#距离判别:按照方差是否相等比较x与总体均值的距离。#Bayes判别:假定对研究对象已经有一定的认识,但这种认识常用先验概率来描述,取得样本后,就可以用样本修正已有的先验概率,得到后验概率。#Fisher判别:按照类内方差尽可能小,类间方差尽可能大来求判别函数。#距离判别: #二分类问题:discriminiant.distance <- function(T原创 2017-12-05 18:11:08 · 21320 阅读 · 0 评论 -
[Rcode]时间序列(整理自R in action)【待更新】
一 生成时序对象sales<-c(18,33,41,7,34)tsales###常规操作得到基础信息plot(tsales) #绘制时序图start(tsales);end(tsales) #起止点frequency(tsales) ###取子集tsales.subset二 时间原创 2017-12-03 23:28:11 · 16914 阅读 · 0 评论 -
[Rcode]线性回归
#简单线性回归: ##常用绘图: fit<-lm(weight~height,data=women) summary(fit) plot(women$height,women$weight,xlab="Height (in inches)",ylab="Weight (in pounds)") abline(fit) fit2<-lm(mpg~wt+I(wt^2),data=mtcar原创 2017-12-12 17:51:45 · 20964 阅读 · 0 评论 -
[Rcode]subset函数
subset函数,从某一个数据框中选择出符合某条件的数据或是相关的列(1)单条件查询> selectresult=subset(df1,name=="aa")> selectresult name age sex1 aa 20 f> df1 name age sex1 aa 20 f2 bb 29 m3 cc转载 2017-11-24 17:06:27 · 27917 阅读 · 1 评论 -
[Rcode]数据框的基础操作【不定期更新】
一 数据的筛选与处理1 添加新的一列,把连续型变量(如身高)转化为类别型变量data$height[data$height>=3]data heightcat heightcat[height>1.8] heightcat[height>=1.6&height heightcat[height原创 2017-12-02 16:13:06 · 15359 阅读 · 0 评论 -
[Rcode]实现文字包含某一字符的批量替换和删除
一 向量中包含空格的字符批量删除grep函数x <- c("A b","abcdfd","123 4","1234fd")x[-grep(" ",x)]二 对向量中包含某一字符(例删除包含空格的变量)的批量替换x<-c("a b","c d")y#当字符串中存在通配符,如"."时。按如下方式使用gsub(),强制取消通配。x <- "1.4"y <- gsub(".",原创 2017-11-16 19:56:45 · 23182 阅读 · 0 评论 -
[Rcode]构造函数把数据中某项状态的持续时间转化为新的向量
countDurationTime <- function(x){ # Examples: x <- c('A', 'A', 'A', 'B', 'B', 'C') #> countDurationTime(x) # items freq # 1 A 3 # 2 B 2 # 3 C 1 if (length(x)原创 2017-11-30 22:18:53 · 19530 阅读 · 0 评论 -
[Rcode]ifelse函数:用平均值替代缺失值,或对数据不同水平分组
一 处理缺失值dta$PM2.5 #把NA值用均值替代data$x 二 对数据分组dta$Level<- ifelse(dta$PM2.5<=35, '优良', ifelse(dta$PM2.5<=75, '轻度污染', ifelse(dta$PM2.5原创 2017-11-30 21:52:54 · 24771 阅读 · 0 评论 -
[Rcode]聚类分析
#聚类分析一般步骤[见R in action343]#计算距离,d<-dist(x,method=)#层次聚类分析(以平均联动聚类为例)if(!require('flexclust')){ install.packages('flexclust') require('flexclust')}data(nutrient,package="flexclust")row.name原创 2017-12-08 11:53:38 · 16020 阅读 · 0 评论