学习R点滴记录（基础）

学习过程中常用的语句，记录下来时常添加和复习

数据库操作

library(DBI)//载入DBI
library(RMySQL)//载入RMySQL
conn<-dbConnect(MySQL(),dbname="dbname",username="root",password="1qaz2wsx",client.flag=CLIENT_MULTI_STATEMENTS)
s<-dbGetQuery(conn,"SELECT distinct S_id,S_country FROM qogismonitor.a_timenew;")
view(mm)
dbDisconnect(conn)

数据处理

排序

newdata<-mm[order(S_id,M_id),]

选择

new<-mm[which(S_id==31),]
a1<-s[a,]  

合并

b<-cbind(b,coef(fit))
total<-merge(a,b,by.x<-"",by.y<-"")//表的连接

重命名

names(a)[3]<-"Ctime"

删除变量

rm()

转置

a<-t(a)

读入数据

write.table(D,"F:/D.txt",row.names=FALSE,col.names=FALSE)
write.csv(D,"F:/D.csv")

数据类型转换

Dis<-as.matrix(dis)

回归

fit<-lm(A$"2"~A$distance)
anova(fit1,fit2)
coef(fit1)//系数
summary(fit1)

描述性数据统计

sapply(Bt[2],mean)
sapply(Bt[2],sd)
Bcountry<-aggregate(Bt$new.distance,by=list(am=Bt$country),mean)//分组统计变量

分组数据是否有明显的差异，如果P<0.001,就认为分组有明显的差异。

> class<-Bt$country
> var<-Bt$new.distance
> mydata<-as.data.frame(cbind(class,var))
> kruskal.test(mydata.var ~ mydata.class, data=mydata)