Как я могу оценить 95% доверительные интервалы, используя профилирование для параметров, оцениваемых путем максимизации логарифмической функции правдоподобия с использованием optim в R?
Я знаю, что могу асимптотически оценить ковариационную матрицу, инвертировав гессиан , но я обеспокоен тем, что мои данные не соответствуют предположениям, необходимым для того, чтобы этот метод был действительным. Я бы предпочел оценить доверительные интервалы, используя другой метод.
Является ли приемлемым метод вероятности профиля , как обсуждалось в Stryhn и Christensen , а также в книге MASS Venables и Ripley, §8.4, с. 220-221?
Если так, есть ли пакеты, которые могут помочь мне сделать это в R? Если нет, то как бы выглядел псевдокод для такого метода?