首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >使用fda软件包phaseplanePlot命令在R中绘制相平面

使用fda软件包phaseplanePlot命令在R中绘制相平面
EN

Stack Overflow用户
提问于 2016-03-25 20:32:26
回答 1查看 520关注 0票数 2

所以我试图计算一些函数数据的相平面图,它是相对于时间测量的,所以我想让图上的标签每12小时标记一次。

代码语言:javascript
复制
attach(alldata)

#identify the time measure form the data, and make a vector
tdata<-data.matrix(alldata[,1:1])

timev<-as.vector(tdata)

#identify the control group
control<-data.matrix(alldata[,2:25])

#create basis and smooth and create functional object
basis17<-create.fourier.basis(rangeval=c(0,5.24677083333333), nbasis=17)

csbasis17<-smooth.basis(argvals=tdata,control,basis17)

controlfd<-fd(coef(csbasis17),basis17) 

#creaste list of points whre i was labels to go and name them 

th<-list("12 hours"=timev[1:1],"24 hours"=timev[15:15],"36 hours"=timev[29:29],"48 hours"=timev[43:43],
         "60 hours"=timev[50:50], "72 hours"=timev[57:57], "84 hours"=timev[71:71],
         "96 hours"=timev[85:85] ,"108 hours"=timev[100:100],"120 hours"=timev[114:114], "136 hours"=timev[128:128])
#crerate list of 2 containing the points i want plotted and the names
ws<- list(th,labels(th))

#phase plane plot command 
phaseplanePlot(timev,controlfd,Lfdobj1=1, Lfdobj2=2,
               lty=c("longdash", "solid"),
               labels=ws)

这将生成以下不带标签的图?

这就是错误栏中出现的内容

代码语言:javascript
复制
Error in text.default(D1., D2., labels$labels) : 
  no coordinates were supplied

我想不出该怎么修??

EN

回答 1

Stack Overflow用户

发布于 2020-10-14 01:37:57

不需要将th放入列表中。由于您的示例不可重现,因此假设长度为1: 100,曲线的前100个评估点,我们希望在数据的原始"O"、中间"M"和末端"E"中绘制一个标签。然后,您应该按如下方式填充label参数:

代码语言:javascript
复制
arg <- 1:100

phaseplanePlot(timev,controlfd, Lfdobj1 = 1, Lfdobj2 = 2,
               lty = c("longdash", "solid"),
               label = list(evalarg = seq(arg[1], max(arg), length=3), 
                            labels = c("O","M","E"))

这里,

代码语言:javascript
复制
seq(arg[1], max(arg), length=3) 
[1]   1.0  50.5 100.0

然后,

票数 0
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/36220011

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档