ANOVA является аббревиатурой от «дисперсионного анализа». Это широко используемый статистический метод для анализа источника отклонений в наборе данных.
Я хочу рассчитать надежный ANOVA, как описано на https://rcompanion.org/rcompanion/d_08a.html а> . Я получаю сообщение об ошибке «Ошибка в list_to_dataframe (res, attr (.data, «split_labels»), .id, id_as_factor): результаты не имеют одинаковой длины».
Что означает ошибка? Что вызывает это?
Я исполь...
Хотя я дважды запускал один и тот же тест, но оказалось, что результаты разные.
Сначала я создал таблицу непредвиденных обстоятельств из своих данных:
> SAtab1 <- table(dt2$ACE_SA_BIN, dt2$IHD_other_healthy)
> SAtab1
0 1 2
0 345 577 29
1 10 25 2
А затем я запускаю тест Фишера...
16 Мар 2023 в 21:59
У меня есть датафрейм:
sub3 <- df1[, c('Attrition', "Age", "DistanceFromHome", "MonthlyIncome", "NumCompaniesWorked", "PercentSalaryHike", "TotalWorkingYears", "TrainingTimesLastYear", "YearsAtCompany", "YearsSinceLastPromotion", "YearsWithCurrManager")]
sub3
Где Истощение — это переменная о...
Я использую R для запуска ANOVA (обычный код aov(), который дает следующую информацию: Df, Sum Sq, Mean Sq, F value , Pr(>F)), и есть часть информации, которую я не смог получить из этих таблиц: Оценочные предельные средние для зависимой переменной.
Я знаю, что можно получить эту дополнительную инф...
Я пытаюсь запустить тест ANOVA для своих данных, используя этот фрейм данных.
df = structure(list(Total_Devices = c(2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2...
Я строю график, на котором показан средний ответ для уровней двух факторов, один на трех уровнях, а другой на семи. Однако с этим кодом:
library(dplyr)
library(ggplot2)
library(tidyr)
desfactorial %>%
mutate(across(Solvent:Water_added, as.factor)) %>%
pivot_longer(Solvent:Water_added) %>%
mut...
Я выполняю ортогональный план L9, и чтобы визуально увидеть, как различаются уровни факторов, мне интересно построить график средних значений. С кодом:
plot.design(PCR ~ Temperature+Time+Ratio+Pretreatment, data = scenedesmus)
Я получаю:
Я хотел бы улучшить этот график с помощью ggplot2, но мне э...
Я работаю над двухсторонним повторным измерением, используя пакет rstatix.
library(rstatix)
res.aov <- anova_test(data=kolding, dv=count, wid=id, within=treatment)
get_anova_table(res.aov)
А вот как выглядят данные:
id treatment species count
1 K1 biohut Goldsinny Wrasse 0
...
Это изображение имен моих столбцов в значении, а затем анова, которую я настроил для нее, которая работает.
Я хочу, чтобы у меня была возможность переключать имя каждого столбца в моей анове в цикле, поэтому мне не нужно переписывать код и просто менять имя столбца.
Пытался настроить цикл for и не с...
27 Фев 2023 в 17:39
У меня есть эти данные в R:
V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 experiment
1 87 50 37 50 50 50 75 75 87 75 1010v1
2 100 50 62 50 37 62 75 87 75 87 1010v1
3 100 100 100 100 100 100 100 ...
27 Фев 2023 в 17:10
Я пытался одновременно оценить несколько ANOVA с помощью цикла. Но я хочу перебрать как несколько предикторов, так и несколько результатов. Поэтому я пытался сделать вложенный цикл.
#data
test<-structure(list(Alcohol = c(1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L,
0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, ...
У меня возникли проблемы с поиском краткого способа повторения одностороннего анализа нескольких фреймов данных в R. Все фреймы данных, которые я использую, имеют одинаковое форматирование, но я пытался понять, как поместить их все в список, а затем выберите переменную в данных, чтобы выбрать для а...
20 Фев 2023 в 18:55
Я хочу нарисовать точно такой же график в R. Однако я хочу рассмотреть два варианта: (1) с одной осью x для каждого пола & (2) два разных xax для каждого пола. Вот также ссылка, где я нашел изображение: https://rpubs.com/WhataBurger/Anovatype3
Спасибо, что поделились знаниями.
Вот случайно сгене...
19 Фев 2023 в 20:50
Я использовал пример набора данных, который я загружаю в фрейм данных. Затем я использую statsmodels OLS, сравнивая текстуру в результате смешивания, а затем использую эту модель для таблицы ANOVA.
import numpy as np
import pandas as pd
import statsmodels.api as sm
from statsmodels.formula.api impo...
3 Фев 2023 в 01:06
Я хочу манипулировать набором данных, чтобы сделать его пригодным для тестирования ANOVA. Текущий способ структурирования df похож на df1, с множеством точек данных, нескольких типов и разделенных контекстуальными категориями. Насколько я понимаю (что может быть неправильно), мне нужно изменить ст...
Это вопрос R-статистики. У меня есть данные по многим предметам. Моя зависимая переменная - это некоторая мера крови, скажем, количество лейкоцитов (непрерывная переменная). bc = 5,6 Меня интересует независимая переменная группа Dx (3 уровня: контроль, депрессия, ремиссия). Я хочу «исправить» (д...
Как конкретизировать контрасты (точечные оценки, 95CI и p-значения) для межгрупповых различий внутригрупповых дельта-изменений?
В приведенном ниже примере меня бы интересовали межгрупповые (группа = 1 минус группа = 2) дельта-изменения (время = 3 минус время = 1).
ДФ и модель:
demo3 <- read.csv("htt...
17 Янв 2023 в 00:40
Я хочу использовать однократные повторные измерения в своем наборе данных, чтобы проверить, различаются ли значения 5 пациентов между измеренными 3 днями.
Я использую AnovaRM из statsmodels.stats.anova, и результатом является объект «AnovaResults». Я вижу p-значение с помощью функции print(), но не...
13 Янв 2023 в 23:57
Я пытаюсь провести тест Даннета. У меня есть два разных фактора: 1) группа, 2) условие Я могу запустить код с одним фактором, используя этот код
Dunnett <- glht(anova_growth,
linfct = mcp(group= "Dunnett"))
summary(Dunnett)
Однако мне нужно учитывать оба фактора, я пробую этот к...
Я пытаюсь запустить анова-тест для описательных переменных с 4 разными группами, эти 4 группы сгруппированы в зависимости от наличия или отсутствия 2 осложнений.
Мои данные
structure(list(values = c("F", "F", "M", "F", "F", "M", "F",
"F", "F", "F", "F", "F", "F", "M", "M", "F", "F", "F", "F", "M"
)...
11 Янв 2023 в 01:10
Я управлял этой моделью МЛМ
model1.a <- lmer(log10_number_of_sensations ~ Tone*Location + Tone*Intensity
+(1+Tone+Intensity|Subject),
neutralBlock_data,REML=TRUE)
И я использовал «lmertest», чтобы получить таблицу вариаций повторных измерений:
a1.a = anova(model1.a, ddf="K...
Программный выбор/запуск групповых сравнений (включая апостериорные)
Я хотел бы выполнить цикл в r, где, если значение p теста Шаприо-Уилка меньше 0,05, выполняется тест Краскела-Уоллиса, но если оно больше 0,05, выполняется односторонний ANOVA. Расширяя это, если результат Краскела-Уоллиса или одн...
26 Дек 2022 в 12:11
Я пытаюсь сделать коробочную диаграмму из набора данных WeightLoss из автомобильного пакета. Интересующие переменные разделены по месяцам, поэтому я сделал диаграммы по месяцам отдельно, показывая все группы (Контроль, Диета, Диета и Упражнения). Я только хочу получить блок-график групповой диеты ...
19 Дек 2022 в 03:47
Я хотел бы провести апостериорные тесты Tukey HSD для повторного измерения ANOVA. Введенная формула "TukeyHSD" возвращает мне ошибку. Не могу найти ответ на форуме. Могу ли я попросить о помощи?
«лечить» — фактор повторных измерений, «vo2» — зависимая переменная.
Ниже приведен скрипт, который выд...
Я использую tidymodels для соответствия нескольким моделям Random Forest. Затем я следовал этому руководству, чтобы сравнить модель Результаты. Проблема в том, что я получаю сообщение об ошибке: Ошибка в
UseMethod("anova") :
no applicable method for 'anova' applied to an object of class "ran...
12 Дек 2022 в 02:28