我正在尝试运行混合效果模型,F2_difference
该模型预测其余列作为预测变量,但我收到一条错误消息
固定效应模型矩阵排名不足,因此下降7列/系数.
从这个链接,固定效果模型是排名不足,我想我应该findLinearCombos
在R包中使用caret
.但是,当我尝试时findLinearCombos(data.df)
,它给了我错误信息
qr.default(object)中的错误:外部函数调用中的NA/NaN/Inf(arg 1)另外:警告消息:在qr.default(object)中:强制引入的NAs
我的数据没有任何NA - 可能导致这种情况的原因是什么?(对不起,如果答案很明显 - 我是R的新手).
我的所有数据都是除了我试图预测的数值之外的因素.这是我的数据的一小部分样本.
sex <- c("f", "m", "f", "m")
nasal <- c("TRUE", "TRUE", "FALSE", "FALSE")
vowelLabel <- c("a", "e", "i", "o")
speaker <- c("Jim", "John", "Ben", "Sally")
word_1 <- c("going", "back", "bag", "back")
type <- c("coronal", "coronal", "labial", "velar")
F2_difference <- c(345.6, -765.8, 800, 900.5)
data.df <- data.frame(sex, nasal, vowelLabel, speaker,
word_1, type, F2_difference
stringsAsFactors = TRUE)
Run Code Online (Sandbox Code Playgroud)
编辑:这是一些更多的代码,如果它有帮助.
formula <- F2_difference ~ …
Run Code Online (Sandbox Code Playgroud) 我正在尝试对我的数据进行子集化,并且在我的一生中我无法弄清楚为什么这不起作用。这是我的数据的链接,它是一个 csv 文件Data。我已经尝试过这些方法来进行子集化,但它们都只是创建了一个没有观察的标题。
alveolar.df <- combine.df[combine.df$articulation == "velar_coronal", ]
alveolar.df <- subset(combine.df, articulation == "velar_coronal")
alveolar.df <- combine.df[combine.df$articulation %in% c("velar_coronal"), ]
Run Code Online (Sandbox Code Playgroud)
这是数据的示例 -
combine.df <- structure(list(start = c(215.76, 278.78, 862.35, 1050.22, 512.59,
552.97), end = c(217.18, 280.68, 863.85, 1051.72, 515.02, 555.67
), articulation = c(" velar_coronal", " velar_labial", " velar_coronal",
" velar_coronal", " velar_coronal", " velar_coronal"), hometown = structure(c(1L,
NA, 2L, 3L, 4L, 4L), .Label = c(" Montana ", " ", " San Leandr ",
" Southern California …
Run Code Online (Sandbox Code Playgroud) 我有以下代码绘制点并在它们之间画一条线.
ggplot (data = subset(df, vowel == "O" & gender == "f"), aes (x = time, y = val, color = formant)) +
geom_point()+
geom_line(aes(group=interaction(formant, number)))
Run Code Online (Sandbox Code Playgroud)
它产生了这个:
有没有办法按颜色/线型对负斜率与这些线的正斜率进行分组?
编辑:这是我的数据:
number <- c(1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3)
formant <- c("F2", "F2", "F2", "F2", "F2", "F2", "F3", "F3", "F3", "F3", "F3", "F3")
time <- c(50, 50, 50, 99, 99, 99, 50, 50, 50, 99, 99, 99)
val <- c(400, 500, 600, 450, 550, 650, 300, …
Run Code Online (Sandbox Code Playgroud)