r - IRLS関数で「Model.frame.default(formula = y ~x)変数の長さが異なる(見つかった)のエラー」の対処方法

Question

反復的に再重み付けされた最小二乗法を使用して、方程式の係数を推定しようとしています:

y = a0 + a1 * cos(2*pi/365 * x) + b1 * sin(2*pi/365 * x) + c1 * x

IRLS を使用してこの式を当てはめると、エラーが発生します。

x <- c(69,77,117,149,165,197,213,221,229,237,261,269,301,309,40,104,128,144,168,176,192,208,224,240,256,352) 
y <- c(1812,1603,1689,1901,1858,1834,2138,1830,1819,1944,1895,1775,1700,1712,1876,1657,2094,2116,2121,2276,2198,2211,2003,1944,1795,1698) 
a0=1 
a1=2 
b1=3 
c1=4  
fit <- irls(y ~ a0 + a1 * cos(2*pi/365 * x) + b1 * sin(2*pi/365 * x) + c1 * x, data= heart)    
Error in model.frame.default(formula, data) : 
  variable lengths differ (found for 'a0')   

どうすればこの問題に対処できますか? IRLS を使用して係数を推定する他の方法はありますか?