Блог

Мне было интересно узнать о некоторых вещах, в dummy.coef() которых преобразуются оцененные параметры (контрасты) в моделях ANOVA в исходные. Это работает только для одномерных моделей, но изменения, необходимые для того, чтобы заставить его работать и для многомерных моделей, кажутся незначительными. В dummy.coef.lm() :

• строка 52 coef <- object$coefficients должна быть предназначена coef <- as.matrix(object$coefficients) для размещения одномерных и многомерных моделей ( coef(object) является вектором в 1-м случае и матрицей во 2-м)
• строка 60 ans <- drop(mm[rn == tl[j], keep, drop = FALSE] %*% coef[keep]) должна ans <- drop(mm[rn == tl[j], keep, drop = FALSE] %*% coef[keep, ]) содержать все столбцы в coef
• строка 61 names(ans) <- rnn[rn == tl[j]] может names(ans) <- rep(rnn[rn == tl[j]], ncol(coef)) содержать имена строк всех столбцов

Метод печати потребует некоторых изменений, но, похоже, это все. Кто-нибудь знает, почему dummy.coef() не был разработан для обработки многомерных моделей?

Еще одна вещь, на которую я наткнулся: строки 20-22

 for (i in vars) args[[i]] <- if (nxl[[i]] == 1)
    rep.int(1, nl)
else factor(rep.int(xl[[i]][1L], nl), levels = xl[[i]])
  

Это безопасно? Т.Е., Если if() предложение есть TRUE , не будет ли неожиданного else ? Я бы ожидал чего-то вроде

 for (i in vars) args[[i]] <- if (nxl[[i]] == 1) {
    rep.int(1, nl)
} else { factor(rep.int(xl[[i]][1L], nl), levels = xl[[i]]) }
  

Это касается только вашего второго вопроса (вопроса об if() x else y операторе в строках 20-22 в коде).

Для начала попробуйте вырезать и вставить эти два блока в сеанс R.:

 test <- TRUE

# Block 1 -- Doesn't work
if(test) 
    cat("test is TRUEn")
    else
    cat("test is FALSEn")

# Block 2 -- Works
{
if(test) 
    cat("test is TRUEn")
    else
    cat("test is FALSEn")
}
  

Что происходит? В {} этом вся разница. Блок 1 «чтение-анализ-оценка» кода построчно, вызывая именно ту проблему, которую вы ожидаете. Блок 2, с другой стороны, считывается и полностью анализируется перед выполнением какой-либо оценки. Это часть того, что {} предписывает R делать. Когда он получает блок кода в целом, анализатор явно анализирует if() x else y блок как одно выражение.

Код, который вы процитировали, был взят изнутри тела функции (т. Е. Внутри {} пары). В этом контексте он обрабатывается правильно (т. Е. Как Блок 2).

HTH

У нас есть официальный отчет об ошибке dummy.coef с «своего рода» патчем, на который я сейчас смотрю — он кажется хорошим, хотя и меняет немного больше вещей сразу, чем я ожидал.

Когда я искал другие проблемы в Google dummy.coef , я получил этот старый поток SO. Позвольте мне ответить (как член основной команды R), что мы заинтересованы в расширении до многомерных моделей, если изменения в коде (и документации) невелики. Пожалуйста, ответьте, если вы заинтересованы в помощи.