2012-05-15 26 views
0

我试图在this数据中绘制变量“ipen”对“年份”的增长样条。我使用的代码是:R:ylim和xlab/ylab在plot()for grofit package not working

grofit <- read.csv('http://dl.dropbox.com/u/1791181/grofit.csv') 
# install.packages(c("grofit"), dependencies = TRUE) 
require(grofit) 
growth = gcFitSpline(grofit$year, grofit$ipen) 
plot(growth) 

这工作正常,并产生下面的情节。但问题是(a)我无法使用xlabylab选项更改默认标签,并且(b)plot()语句中的ylim=c(0,100)无法将x轴的比例更改为(0,100)。我在grofit()文档中找不到任何指针。 this

+1

我会使用'plot.gcFitSpline()'函数并对其进行自定义。 –

回答

2

不要直接使用plotgcFitSpline结果,而是使用str(growth)探索曲线拟合对象的结构,并确定要绘制什么,或者看plot.gcFitSpline代码(只需键入该函数的名称,而不括号中的控制台,然后按Enter)

例如:

growth = gcFitSpline(grofit$year, grofit$ipen) 

plot (grofit$year, grofit$ipen, pch=20, cex=0.8, col="gray", ylim=c(0, 100), 
    ylab="Growth", xlab="Year", las=1) 
points(growth$fit.time, growth$fit.data, t="l", lwd=2) 

这给:

result grofit plot

+1

绝对完美!非常感谢,尼科。 – user702432

1

其实,我觉得grofit文档确实有一个关于这个指针:

...其他图形参数也可以作为参数传递[原文]。这目前没有效果,只是为了满足通用函数的要求。

您可能需要重写具体到你想要的图形参数plot.gcFitSpline功能,但有可能是一个原因,功能是它的方式。

顺便说一句,我不确定给你的数据和你使用的软件包同名是个好主意。