我已经使用R中的permuco对某些时间序列数据进行了聚类分析。(对控制/治疗条件的标签进行置换,并计算F统计量,这些显着差异的时间聚类是偶然发生的可能性。) 到目前为止一切顺利。
我使用此软件包随附的内置函数plot.clusterlm生成了许多图。但是,数据来自不同的组,并且在每个图中重新缩放y轴上的F值,即,根据效果的强弱重置值和刻度。
这是有问题的,因为基于不同聚类分析的不同图在视觉上不具有可比性。 我想重新缩放y轴,以使所有聚类都沿相同的F值可视化(例如0-10)。
我还没有做到这一点,我想知道是否有办法将任何其他函数传递到plot.clusterlm中。 这是函数的用法,但是我看不到重新调整y轴的方法。 (尽管可以通过操纵nbbaselinepts和nbptsperunit来缩放x轴,但这不是我想要的...)
plot(x, effect = "all", type = "statistic",
multcomp = "clustermass", alternative = "two.sided",
enhanced_stat = FALSE, nbbaselinepts = 0, nbptsperunit = 1, ...)
如果对此有任何想法,请告诉我。 谢谢!
答案 0 :(得分:1)
感谢您使用permuco
!我在GitHub上发布了一个问题,以提供实现这些功能的解决方案。您可以预期permuco
的后续发行版中会发生更改。
但是,plot()
方法显示出F统计量,这不是衡量效果大小的好方法。效果大小的更好衡量方法是在afex package
答案 1 :(得分:0)
在基本R中,绘图设备的轴是这样改变的:
x<-1:10; y=x*x
# Simple graph
plot(x, y)
# Enlarge the scale
plot(x, y, xlim=c(1,15), ylim=c(1,150))
# Log scale
plot(x, y, log="y")
这是STHDA中的一个示例,您可以在其中找到许多有用的教程。