R
paulsweet_M
你不知道的事!
展开
-
2015年09月22日
#ABG,P80> TEM > > names(TEM) [1] "Sample" "Date" "DateNr" "dDay1" "dDay2" [6] "dDay3" "Station" "Area" "X31UE_ED50" "X31UN_ED50" [11] "Year" "Mont原创 2016-07-23 18:11:25 · 252 阅读 · 0 评论 -
基础补缺:rep,dim,length
#加括号是可以直接让赋值语句计算出结果(a#rep中有each和无each的区别rep(a,each=2)> rep(a,each=2) [1] 0 0 5 5 10 10 15 15 20 20 25[12] 25 30 30 35 35 40 40 45 45 50 50[23] 55 55 60 60 65 65 70 70原创 2016-07-23 18:15:09 · 544 阅读 · 0 评论 -
自相关函数acf pacf
income "income_end.txt",header = T)y 7]x 1]acf(y,na.action = na.pass) #缺失值出现?acfpacf(y,na.action = na.pass)原创 2016-07-23 18:15:03 · 10630 阅读 · 0 评论 -
初识--循环
x for (i in seq(from=1,to=100,by=1)) print(x[i])原创 2016-07-23 18:14:58 · 144 阅读 · 0 评论 -
unbiase application plot
set.seed(1)x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13 x14 x15 x16 x17 x18 x19 x20 x21 x22 x23 x24 x25 x26 x27 x28 x29 x30 x31 x32 x33 x34 x35 x36 x37 x38 x39 x40 x41 x42 x43 x44 x45 x46 x47 x48 x49 x5原创 2016-07-23 18:14:24 · 1458 阅读 · 0 评论 -
如何运用循环赋值
#循环中,每次循环都打印出i的结果。#循环结束后,i不再更新,打印i输出的就是i最后一次的取值,#并且一直都是这个值,除非你再次给i赋值。set.seed(1)for (i in seq(from=1,to=500,by=1))i #i,运行得出500print(i) #print(i),运行得出1-500原创 2016-07-23 18:14:20 · 1399 阅读 · 0 评论 -
2015年10月21日
#补之前一条参考线之遗憾,添加参考线,水平h,垂直vabline(v=c(0),lty=2,lwd=2)原创 2016-07-23 18:12:55 · 148 阅读 · 0 评论 -
2015年10月21日
#改天写注释set.seed(1)x1x2x3x4x5opar par(mfrow=c(2,2))hist(x1,freq = F,breaks=10,main="n=10,t(4)",xlab = "x")lines(density(x1),col="red",lty=5,lwd=2)hist(x2,freq = F,breaks=20,main="n=20,t(4)",xlab = "x")l原创 2016-07-23 18:12:52 · 190 阅读 · 0 评论 -
2015年10月19日
#循环eg.ALLNests for (i in 1:27) { Nest.i Owls.i YourFileName jpeg(file=YourFileName) plot(x=Owls.i$ArrivalTime,y=Owls.i$NegPerChick, xlab="Arrival Time",ylab = "Negotiation behaviour",main =原创 2016-07-23 18:12:43 · 216 阅读 · 0 评论 -
2015年10月19日
#循环初谈,jpeg的输出!!!Owls names(Owls)str(Owls)unique(Owls$Nest)#setwd不是在编入的时候就不能改,而是调用入Rstudio的数据库bottomleft window储存起来;#只是在读取原数据时设定一下;#如果编过程中想选择别的出输出文件夹就可以换;#mac直接用bottom right的window来设置路径,win则用setwd语句;#注原创 2016-07-23 18:12:36 · 174 阅读 · 0 评论 -
2015年10月19日
> #3.比较计算: ==, >, =, isTRUE, identical> > a> a==a;a!=b;a>b;a[1] TRUE[1] TRUE[1] TRUE[1] FALSE> > #判断是否为TRUE,没太搞懂!!> isTRUE(a)[1] FALSE> isTRUE(!a)[1] FALSE> > #精确比较两个对象> identical(1,as.integer(1))[1]原创 2016-07-23 18:12:32 · 160 阅读 · 0 评论 -
2015年10月18日
> #1.基本计算(四则运算: 加减乘除, 余数, 整除, 绝对值, 判断正负)> > a> > #加减乘除> a+b;a-b;a*b;a/b[1] 15[1] 5[1] 50[1] 2> > #余数,整除> a%%b;a%/%b[1] 0[1] 2> > #绝对值> abs(-a)[1] 10> > #判断正负> sign(-2:3)[1] -1 -1 0 1 1 1> > #2.幂,原创 2016-07-23 18:12:27 · 324 阅读 · 0 评论 -
2015年10月15日
#生成随机数编号为1set.seed(1)x #生成序列,序列长度1000,可以理解为序列的精细度 y1 #生成N(O,1)的PDF#?xaxs"i";lwd是线宽,递增plot(x,y1,col="red",xlim=c(-5,5),ylim=c(0,0.45),type='l', xaxs="i", yaxs="i",ylab='Density',xlab='x',原创 2016-07-23 18:12:09 · 210 阅读 · 0 评论 -
2015年10月8日
#ABG中循环部份例题Owls names(Owls)str(Owls)unique(Owls$Nest)Owls.ATV plot(x=Owls.ATV$ArrivalTime,y =Owls.ATV$NegPerChick,xlab = "Arrival Time", ylab ="Negotiation behaviour", main="AutavauxTV")原创 2016-07-23 18:11:51 · 257 阅读 · 0 评论 -
2015年09月26日
Veg plot(x=Veg$BARESOIL,y=Veg$R,xlim = c(0,45),ylim = c(4,19))M.Loess #对数据进行计算,Loess平滑FIT #拟合Ord1 lines(Veg$BARESOIL[Ord1],FIT[Ord1],lwd=3,lty=2,col=14) #在图上添线原创 2016-07-23 18:11:36 · 249 阅读 · 0 评论 -
2015年9月25日
#plot,其中对颜色按数据分析Veg names(Veg)?parVeg$Time2 Veg$Time2 [Veg$Time Veg$Time2 [Veg$Time >1974] Veg$Col2 Veg$Col2 [Veg$Time Veg$Col2 [Veg$Time >1974] plot(x=Veg$BARESOIL,y=Veg$R,xlab = "Exposed soil",原创 2016-07-23 18:11:33 · 223 阅读 · 0 评论 -
文本分词,记词频,做词云
a b ?table lecture=read.csv("1.segment.txt",sep=",",header=TRUE,fileEncoding="UTF-8") # 查看前几行,看是否有字符编码问题 head(lecture); # 获取数据集长度 n=length(lecture[,1]); print(n) # == 文本预处理 res=lecture[lecture原创 2016-07-23 18:15:45 · 1004 阅读 · 0 评论