csaps()
в Matlab выполняет кубический сплайн в соответствии с конкретным определением параметра сглаживания p
. Вот некоторый код Matlab и его результат:
% x variable
age = 75:99
% y variable
diffs = [-39 -2 -167 -21 -13 32 -37 -132 -143 -91 -93 -88 -62 -112 -95 -28 -90 -40 -27 -23 -28 -11 -8 -6 1]
% 0.0005 is the parameter p, and the later specification of
% age are the desired x for prediction
csaps(age,diffs,0.0005,age)
% result (column headers removed):
-63.4604 -64.0474 -64.6171 -65.1397 -65.6111 -66.0165 -66.3114
-66.4123 -66.2229 -65.6726 -64.7244 -63.3582 -61.5676 -59.3568
-56.7364 -53.7382 -50.4086 -46.7922 -42.9439 -38.9183 -34.7629
-30.5180 -26.2186 -21.8912 -17.5532
Я хотел бы получить тот же результат в R. Я пробовал base::smooth.spline()
, но параметр сглаживания spar
указан по-другому, что я не могу связать с p
Matlab (можете вы ?). Ближайший результат, который мне удалось получить, был с функцией smooth.Pspline()
пакета pspline
. Вот некоторый код, чтобы все закрутилось в R:
age <- 75:99
diffs <- c(-39L, -2L, -167L, -21L, -13L, 32L, -37L, -132L, -143L, -91L,
-93L, -88L, -62L, -112L, -95L, -28L, -90L, -40L, -27L, -23L,
-28L, -11L, -8L, -6L, 1L)
predict(pspline::smooth.Pspline(
x = age,
y = diffs,
norder = 2,
method = 1,
spar = 1 / 0.0005 # p given in MP and matlab as 0.0005
),age)
# which gives something close, but not exactly the same:
[1] -63.46487 -64.05103 -64.61978 -65.14158 -65.61214 -66.01662 -66.31079
[8] -66.41092 -66.22081 -65.67009 -64.72153 -63.35514 -61.56447 -59.35372
[15] -56.73367 -53.73584 -50.40680 -46.79098 -42.94333 -38.91850 -34.76393
[22] -30.51985 -26.22131 -21.89474 -17.55757
Страница помощи csaps() находится здесь
smooth.spline()
справку можно найти здесь (код не указан, потому что я думаю, что связь между spar
и p
довольно туманна, поэтому, возможно, не стоит идти по этому пути)
pspline::smooth.Pspline()
помощь находится здесь
Этот запрос другого человека из 2008 года, похоже, остался без ответа. , заставляя меня чувствовать себя этим парнем.
R переполнен сплайнами, поэтому, если сообразительные среди вас могут указать мне на тот, который делает то же самое, что и csaps()
Matlab (или трюк в этом направлении), я был бы очень признателен.
[РЕДАКТИРОВАТЬ 19-8-2013] spar
необходимо указать как (1-p)/p
(а не 1/p
), и тогда результаты будут согласовываться с точностью, насколько это возможно. Смотрите ответ ниже.
edit csaps
и посмотрите на код. - person horchler   schedule 17.08.2013