Hi, > I first thought that this p=103 was the error, however I did the following: > > wnew <-W[-1,] You need to remove one row AND one column to have a 102x102 matrix. In the code above you are just removing one row. I believe that the error is there. But you should check why you W is 103x103 if you only have 102 data points… Best, Virgilio