我正在处理原始测序数据,我想使用平滑参数为 10^-16 的三次样条对其进行平滑处理。我想使用 predict(smooth.pspline( parameters )) 创建一个包含平滑数据的数据框,但我不确定要在 pspline 函数中放入哪些参数来执行我想要的操作,因为 x 必须升序。如果我正在使用遵循以下格式的非常大的数据集(存储在名为“ipavg”的 df 中):
Chromosome Location Timing
1 33985 -4.402
1 11749 -0.4502
...
13 189184 0.2031
...
问题是我不知道如何将 x 和 y 值放入参数中,因为它需要按升序排列,并且此处“Timing”的值会上下波动(有许多高峰和低谷)。如何使用具有平滑参数 10^-16 的 pspline 平滑 ipavg 并将预测的平滑值放入相同格式的新数据帧(稍后可以绘制)中?任何帮助表示赞赏!谢谢你。