Teykovo There seems to be a bug in a R package named pgam:periodogram(). The function intensity happened to have two cos. One of them should be sin. Here is the correct version. function (w, x) { n