回复 第2楼 的 nuomin：找到了，谢谢
要用乘法模型调整季节性数据，怎么计算季节指数呀。请教
回复 第3楼 的 nuomin：应该是不行，只能用AIC准则
回复 第2楼 的 diyiziran：嗯，就得再做一步去掉了
用step逐步回归，结果选择了三个变量，其中一个变量的p值=0.06，怎么设置参数使结果只含有两个p值小于0.05的变量，谢谢！
回复 第2楼 的 nuomin：好的，谢谢
rt
用lattice包的parallelplot函数画平行坐标图
library(lattice)
parallelplot(iris[1:10,],horizontal.axis = FALSE)结果中有一条直线是在最上方的一条水平线，这应该是这个观测的各个变量值都是最大值导致的。要怎么调整参数才能使它不是水平线，而是具体的折线呢？？
回复 第2楼 的 ahhhh：好像是这样。这样直接用读入的数据求mean()函数报错。。
用read.csv读数据时，总是在左侧有序号一列，就是1,2,3,...这样的，要怎么才能只读取数据而不独进这些呢
rt
回复 第4楼 的 doctorjxd：是在一个半圆内画的那种
有没有现成的命令，还是要自己编写函数？
回复 第5楼 的 shelly_dan：
library(ismev)
mrl.plot() #MEF
gdp.fitrange() #根据稳定性确定阈值的图回复 第3楼 的 shelly_dan：先载入evir包
回复 第1楼 的 shelly_dan：
library(evir)
help(hill)回复 第1楼 的 David TJ：grep()
回复 第5楼 的 superdesolator：
比较笨的方法。。
>value=numeric(26)
>for(i in 1:26){if(b[i,]>= a[1,'startdate'] & b[i,]< a[1,'enddate']) value=a[1,'value'] else value=a[2,'value']}
>cbind(b,value)
或者
>library(xts)
>b=xts(numeric(26),seq(as.POSIXct("1999-01-01"),length=26,by="months"))
>b["1999-01-01/2001-12-01"]=1
>b["2001-01-01/"]=2
同求好方法出现。。
回复 第1楼 的 warrenzhang：时间可以直接比较吧
> "2001-01-01" < "2001-02-01"
[1] TRUE
> "2001-01-01" > "2001-02-01"
[1] FALSE
是这意思吗回复 第8楼 的 meilinvhai：应该根据具体问题分析吧。如果某一观测确实的变量比较多就删去，因为填充会有较大的误差；对于不能处理具有缺失值的如线性回归等要进行填充，可以用众数或者通过分析相关性相似性等填充；对于能处理具有确实值的像回归树等应该不用管也行吧。。。个人见解。。。