R语言初识

Posted 2020-06-13

# 创建数据集&基本数据管理
1、向量 创建函数 c() a <- c(1,2,3,4)
a[c(i,j)] :[]给定元素所处位置的数值,即向量a中第i和第j个元素，a[2]第二个元素即2
length(a):返回向量a中元素的个数
2、矩阵 创建函数 X<-matrix（数据，nrow=n,ncol=m）
+使用下标及方括号[]来选择矩阵中的行、列或元素，X[i,],第i行
+将矩阵转化为数据框 dataframe <-as.data.frame(matircname)
3、数组 创建函数 <-array（）
4、数据框：将不同类型数据集放入一个矩阵中， mydata<-data.frame(变量1，变量2，。。)
+指明数据框中的变量方式：mydata$变量名
+使用下标及方括号[],mydata[1：i],1到i列
+指定列明，,mydata[c("column1","column2")]
----
+利用函数attach()将数据框添加到R的搜索中，也可以利用with()
+names(dataframe):生成一个包含所有变量名的字符型向量
5、因子：名义型变量及有序型变量称为因子，<-factor(因子),（把因子转译成数字）
6、列表：可以整合向量、矩阵、数据框及其他列表的组合，创建函数 <-list()
7、创建新变量
+ dataframe$varietyname <-dataframe$variety1 +/-/*/ dataframe$varerity2;;
+ attach(dataframe) dataframe$varietyname <-variety1+variety2
detach(dataframe)
+dataframe <-transform(dataframe,variety=variety1+/-/ variety2)
7、变量的重编吗 数据框$变量[dataframe$varity condition] <-expression:
当变量满足什么条件时，将expression赋值给新变量，或者用函数with（）
8、变量的重命名 rename() 或name()
9.函数is.na（）：检测缺失值是否存在
10.na.rm=true,可在计算之前移除缺失值并使用剩余值进行计算
11、na.omit()可以删除所有含有缺失数据的行
12、as.date(X,"input_format")用于将日期以字符串形式输入R中转化为以数值形式存储的日期变量
13、Sys.date()返回当天的日期，date()则返回当前的日期和时间
14、format(x,format="%B %d %Y"):输出指定格式的日期值
15、函数difftime(newday,oldday,units="week/day/year")计算时间间隔，并以星期，天，时等表示
16、函数as.character(),将日期转换为字符型变量
17、order（）函数对数据框变量进行排序
18、数据集的横向合并即添加列，使用函数merge（）,
+ 利用相同变量内连接 新数据框名 <- merge(dataframeA,dataframeB,by="变量"/by=c("bianliang1","bianliang2"))
+ 无相同变量的连接 cbline( A,B)
19、数据集的纵向合并函数添加观测 rbline(A,B),条件A与和变量数要相同
-----
#数据管理
1、scale()对矩阵或数据框的指定列进行均值为0、标准差为1的标准化
+ newdata <-scale(dataframe)
+ newdata <-transform(dataframe,variety=scale(variety)*SD+M)
2、概率函数
密度函数dorm,
分布函数pnorm,
分位数函数qnorm,
随机数生成函数rnorm
3、函数runif（）用来生成0~1区间上服从均匀分布的伪随机数。可以通过set.seed()设定随机数种子
4、根据指定的均值和协方差生成多元正态数据：MASS包中的mvrnorm()函数
mvrnorm(n,mean,sigma)，然后使用函数as.data.frame()矩阵转化成数据框
5、函数apply(矩阵/数组/数据框，维度下标，分析函数)，对指定的维度（行）进行分析
lapply/sapply将函数应用到list列表中
6、转置：函数t()对一个矩阵或者数据框进行转置
7、数据整合：将多组观测替换为根据这些观测计算的描述性统计量
* 函数aggregate（x待折叠的数据对象，by=list(变量1，变量2)，计算描述性统计函数）
* reshape包，融合melt（）函数及重铸cast（）函数
----
#基本统计分析
一、描述性统计量
1、summary()计算描述性统计量；
2、*Hmisc包中的describe（）函数可返回变量和观测的数量、缺失值和唯一值得数目、平均值、分位数以及五个最大的值和最小的值
*pastecs包中stat.desc()函数或者describe()函数计算描述性统计量
二、分组计算描述性统计量
1、分组计算描述性统计量，利用函数aggregate（）
+ aggregate(数据集中的[变量组]，by=list(变量名=数据集$变量)，function（比如mean.sd等）)，一次执行一个函数
+一次执行多个函数，使用自定义函数及by()，见P134
2、doBy包中的summaryBY()
3、psych包中的describe.by()
4、reshape包
三、频数表