作者buski (......)
看板Statistics
标题[问题] 在R中读入外部资料库的问题
时间Wed Feb 22 05:53:34 2006
我目前用R进行时间序列的分析
我都是用下列语法读入
a1=read.table ('d:xxxx.txt)
data=t(c(a1)) ----->我一直不太清楚这行的目的何在?
work=ts(data)
转成时间序列的形式之後进行分析....
我手上的资料库是每日的雨量 (daily precipitation)
他资料格式做成每个月都有31个栏位,如此一来小月跟二月当然就有misssing
他是用-9999标示
我一开始是将这些栏位都delete掉,不过读入R的时候问题就来了
他会显示第二列没有31个栏位的错误 (看来栏位是被固定住了)
如果我在第一行加入 spe'/' 这个指令 他是能读,不过转成TS的时候
他会用NA取代那些空格。 问题是这样我分析资料会有问题
因为ACF或PACF的指令都不容许有字元形式的资料存在。
请问有指令可以命令他不要去读(-9999)的栏位 或者可以让他连续往下读值而
不受栏位影响吗?
还有一个问题是如果我想计算每周或每月的资料
在R下面该怎麽写呢吗?(每七天相加或按照每月不同天数相加取均值)
谢谢!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 68.43.58.122