R，如何在数据帧上进行循环中的anova？
在R中，可以使用anova()函数在数据帧上进行循环中的方差分析（ANOVA）。
ANOVA是一种统计方法，用于比较两个或多个组之间的均值是否存在显著差异。在R中，可以使用anova()函数来执行方差分析。
首先，确保你的数据帧中包含了需要进行方差分析的变量。假设你的数据帧名为df，其中包含了一个因变量y和一个自变量x。你可以使用以下代码进行方差分析：
model <- lm(y ~ x, data = df) # 创建线性回归模型
anova_result <- anova(model) # 执行方差分析在上述代码中，首先使用lm()函数创建了一个线性回归模型，其中y是因变量，x是自变量。然后，使用anova()函数执行方差分析，并将结果保存在anova_result变量中。
执行完上述代码后，你可以通过打印anova_result来查看方差分析的结果。anova_result将包含各个组之间的方差分析统计量，如F值、p值等。
在R中，还有其他一些函数可以用于执行不同类型的方差分析，如单因素方差分析（one-way ANOVA）、多因素方差分析（two-way ANOVA）等。你可以根据具体需求选择适合的函数进行分析。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云官网：https://cloud.tencent.com/
- 云服务器（CVM）：https://cloud.tencent.com/product/cvm
- 云数据库 MySQL 版：https://cloud.tencent.com/product/cdb_mysql
- 人工智能平台（AI Lab）：https://cloud.tencent.com/product/ailab
- 云存储（COS）：https://cloud.tencent.com/product/cos
- 区块链服务（Tencent Blockchain）：https://cloud.tencent.com/product/tencent_blockchain
相关·内容
我在R和编码方面都是新手。我计算了多列的多个方差分析(总共16列)。为此，"Purr“方法帮助我：
anova_results_5sector <- purrr::map(df_anova_ch[,3:18], ~aov(.x ~ df_anova_ch$Own_5sector))
summary(anova_results_5sector[[1]])
因此，检索输出(p值等)的最愚蠢的方法是以下方法
summary(anova_results_5sector$Env_Pillar)
summary(anova_results_5sector$Gov_Pillar)
summ
浏览 6提问于2022-01-28得票数 0
我有一个定义的函数，它计算我需要的所有必要的统计数据(例如，双向anova和多重比较)。
def stats(F1_para1,F2_para2,para3,para4,value):
#MEAN,SEM,COUNT
msc = df.groupby([F1_para1,F2_para2,para3,para4])[value].agg(['mean','sem','count'])
msc.reset_index(inplace=True) #converts any columns in index as column
浏览 0提问于2021-11-24得票数 0
2回答
我正在对我当前的数据集执行Anova测试，这些数据集有多个列，我试图循环这些列以使事情变得更容易，但在我看来，我似乎总是面临着相同的错误，称为“可变长度不同”。 下面是我的循环代码： for(i in 5:125){
WL<- colnames(NB[i])
model <- lm(WL ~ Treatment , data = NB)
if(!exists("aovNB")){
aovNB<-anova(model)
}
if(exists("aovNB")){
aovNB <- rbind(aovNB,a
浏览 13提问于2021-10-05得票数 0
回答已采纳
2回答
我该如何打印每次迭代？它只打印第一次迭代多次。或者这是不是正确地在每一列上循环？ 编辑:我刚刚意识到我一直在anova中指定相同的变量。但是现在我不知道如何在每一列上运行这个anova。 CR <- data2[data2$cv == "CR",]
colu <- CR[,5:36]
for (i in colnames(colu)){
ANOVA <- aov(i ~ block + loc, data=CR)
ANOVA_SUMMARY <- summary(ANOVA)
test <- HSD.test(ANOVA ,
浏览 19提问于2021-08-24得票数 0
回答已采纳
2回答
我正在尝试在R中运行anova模型。我有一个包含3行12列的数据文件。每一行都是特定级别的解释变量的数据。单元格i，j是第i级响应。文件扩展名为".dat“。我正在运行以下R代码，以尝试获得一个36 x 2的数据帧来运行anova模型，而不是3 x 12的原始数据帧：
data <- read.table("usedcar.dat", row.names = 1)
young <- data[1,]
med <- data[2,]
old <- data[3,]
Price <- c(young, med, old)
Age <- a
浏览 0提问于2017-11-06得票数 12
回答已采纳
1回答
我在记事本中保存了一个名为ANOVA的数据帧：
143 141 150 146 148
152 149 137 143 0
134 136 132 127 0
129 127 132 129 130
当我在R控制台中使用read.table()命令时，
> read.table("ANOVA.txt")
V1 V2 V3 V4 V5
1 143 141 150 146 148
2 152 149 137 143 0
3 134 136 132 127 0
4 129 127 132 129 130
Warning message:
In read.tab
浏览 5提问于2013-05-21得票数 0
回答已采纳
我有三个不同的数据框架，我需要将它们组合成一个来运行ANOVA，创建一些图表等等。我还需要创建一个新的变量来标识每个数据来自哪个数据帧。他们被称为NNA2，TRB和AP。每个人都有不同的长度，因为每一个参与者的数量不同。我在R Studio上跑。
通过使用子集()，我已经将这些数据帧从较大的数据帧减少了。现在我要把它们结合起来。我尝试创建一个新的数据框架作为测试，它给出了以下错误，即有关长度的错误：
prueba <- data.frame(NNA2 2$Sexo，TRB$Sexo，AP$Sexo) data.frame中的错误(NNA2 2$Sexo，TRB$Sexo，AP$Sexo
浏览 1提问于2019-06-20得票数 2
1回答
使用R，我使用下面的脚本创建了一个包含39个anovas的向量
w3<-lapply(split(q1, q1$taxa), aov, formula=trig ~ d)
我可以通过摘要(w3$“type”)分别从每个值中提取p值，但是如何一次对所有39个值进行提取，最好是将所有的p值写入一个单独的数据框架中？
此外，还可以为每个单独的anova创建一个包含摘要的列表。
e7<-lapply(w3, FUN=summary)
但是有什么方法可以从这个列表中用p值来生成数据帧吗？
浏览 2提问于2020-06-05得票数 1
回答已采纳
我正在解决这个问题，不确定如何继续。 加载R数据集mtcar作为pandas数据帧。考虑自变量wt的log和因变量mpg的log，建立线性回归模型。将模型与数据进行拟合。 对上一步得到的线性模型进行方差分析。(提示:使用anova.anova_lm) 显示F统计值。 我在下面的另一篇文章中看到了解决方案。但这似乎行不通。 import statsmodels.api as sm
import numpy as np
mtcars = sm.datasets.get_rdataset('mtcars')
mtcars_data = mtcars.data
liner_model
浏览 70提问于2020-01-17得票数 0
在sklearn 的文档中有一个例子。我想对支持向量机的超视距、i.d.、C和gamma做进一步的GridSearchCV，对于示例中使用的每一个百分位数的特征，如下所示：
transform = feature_selection.SelectPercentile(feature_selection.f_classif)
clf = Pipeline([('anova', transform),
('normal',preprocessing.StandardScaler()),
('s
浏览 4提问于2015-04-19得票数 1
回答已采纳
1回答
我有一个很大的数据帧，其中一些列由于取0的对数而具有NA。
我一直在对数据进行各种测试(ANOVA，Tukey，Kruskal Wallis，Mann Whitney)，但我无法弄清楚NA值发生了什么。
R是否完全排除了这些值？
浏览 3提问于2016-08-05得票数 1
从线性模型中运行的变量pointA，pointB，pointC和anova中，我想收集数据帧中的p值。在我的真实数据集中，大约有30个变量，所以我在寻找一个运行这些变量的函数，并收集数据框中的p值。这是一种比遍历每个数据并手动将它们放在一个数据帧中更有效的方法。
set.seed(1)
id <- rep(1:3,each=4)
trt <- rep(c("A","OA", "B", "OB"),3)
pointA <- sample(1:10,12, replace=TRUE)
pointB<- samp
浏览 13提问于2020-02-28得票数 1
回答已采纳
1回答
我知道通过使用anova()可以得到一个anova表。但是，如果您已经有r中的数据(即您有一个文件并将其读取到r中)。但我有一张印刷的桌子如下所示：
如何使用r计算一个anova表？
谢谢!
浏览 6提问于2022-09-09得票数 0
我认为这是一个愚蠢的问题，但我有点卡住了。我想把我的posthoc Tukey分析放在一个表格中，但是，我不能这样做。我遵循的步骤是:首先，我将posthoc结果转换为数据帧，然后使用gt包创建表。 问题是我无法在表中显示我的组，在本例中是岛屿的名称 这是我正在使用的代码 Tuk_anova1 <- TukeyHSD(AnovaModel.1)
Tuk_anova1 <- TukeyHSD(AnovaModel.1)
Tuk_anova1
TK <-as.data.frame(Tuk_anova1$Island)
gt_tbl <- gt(TK)
gt_tbl The
浏览 12提问于2021-11-01得票数 0
我正在寻找一种方法来从Python中的数据帧执行ANOVA和HSD测试。我试着在论坛和教程上阅读一些例子，但我没有实现将其应用到我的工作中。
这是一个简单的Pandas数据帧：
Date Density Hour Repetition Glucose
A HD AM 1 6.7
A HD AM 2 6.8
A HD PM 2 9.6
A HD PM 3 11.9
B HD AM 1 23
B
浏览 1提问于2014-09-11得票数 1
我正在分析混合设计中的数据，使用来自aov_car()包的afex，具有1中间和1在因子内的数据。
我尝试使用来自anova_stats()包的sjstats，它在设计内部或设计之间非常好，但是抛出了错误：
Fehler in UseMethod("anova") :
nicht anwendbare Methode für 'anova' auf Objekt der Klasse "Anova.mlm" angewendet
英文错误(谷歌翻译)：
Error in UseMethod ("anova"):
non-appl
浏览 0提问于2019-05-21得票数 1
回答已采纳
到目前为止，我的代码如下所示：
Points = readOGR(dsn = "./Data/filename.shp",layer = "layername",stringsAsFactors = FALSE)
Points$LDI = extract(LDI, Points)
LDI = raster("./Data/filename2.tif")
Points$LDI = extract(LDI, Points)
PointsDF = Points@data
for(i in PointsDF) {
Mod1 = lm(LD
浏览 5提问于2016-01-26得票数 2
回答已采纳
1回答
我使用statsmodel库运行了ANOVA分析，如下所示：
import numpy as np
import pandas as pd
from scipy import stats
from scipy.stats import norm
from scipy.signal import savgol_filter
from matplotlib import rcParams
import scipy.stats as stats
#ANOVA ONE WAY
mc = MultiComparison(df['NDVI'], df['line'])
浏览 1提问于2020-02-18得票数 0
1回答
常见的任务是必须对数据集的不同子集执行特定的统计分析(如anova、glm或混合模型)，并将输出表与汇总系数和p值组合在单个数据帧中。我正在寻找一个通用函数，它将接受模型的类型(例如aov(...)或lm(...)或glm(...)或glmer(...) )和特定的输出项，对于这些项，必须根据数据集中的某个分组变量为每个重复分析返回系数和p值。
假设我有一个数据帧，我想在其中对一个数据帧data中不同级别的因子“复制”进行某种分析：
data(iris)
library(car)
data=data.frame()
for (i in 1:10) {data=rbind(data,cbind(
浏览 1提问于2015-09-08得票数 4
2回答
我有一个数据帧，它看起来像下面这样，但是有很多额外的行。"Content“和"Process”都可以是"a“或"n”。
dv1 dv2 dv3 content process
1 3 4 a a
2 1 2 a n
5 3 5 n a
4 3 2 n n
我知道如何对以这种方式设置的数据运行几个方差分析。我将运行以下命令：
dv <- df[,1:3]
output <- lm(as.matrix(dv) ~ con
浏览 0提问于2018-09-06得票数 2
我喜欢在一个数据帧的多个列上运行多个单向ANOVA。我这样做的方法是使用for循环。我的数据框架的第一列包含组。为了提供一个可重复的例子，我在这里取虹膜数据集。我想使用rstatix::anova_test()而不是f.ex。aov()，因为rstatix::anova_test()是管道友好的，似乎是一个更好的选择，不平衡的数据(就像我有它)，也允许定义的平方和类型的方差。
当我用aov()编写for循环时，它可以工作。不幸的是，到目前为止，我在rstatix::anova_test()方面做类似的工作都失败了。有人能帮我吗？
data <- iris %>% relocate(
浏览 2提问于2022-07-09得票数 0
回答已采纳
1回答
我试图在一个数据帧的多个子集上运行一个双向的方差分析，而不必对数据进行实际的子集，因为这是有效的。
示例数据：
DF<-structure(list(Sample = c(666L, 676L, 686L, 667L, 677L, 687L,
822L, 832L, 842L, 824L, 834L, 844L), Time = c(300L, 300L, 300L,
300L, 300L, 300L, 400L, 400L, 400L, 400L, 400L, 400L), Ploidy = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 2L, 2L,
浏览 1提问于2014-05-30得票数 0
回答已采纳
1回答
数据帧示例：
df
ID B C D
1 A 1 1 3
2 B 2 3 1
3 C 1 1 1
4 D 3 1 1
5 E 1 0 0
我用这段代码(在网上找到)对数据帧上的各种变量进行了anova循环。
数据:df
library(car)
LLA <- rep(NA, ncol (df))
sink("dfresults.doc")
for (i in 4:ncol(df)) {
column <- names(df[i])
contrasts(df$Group)<-cont
浏览 2提问于2020-03-16得票数 2
1回答
我使用似然比检验(在R中)在我的模型中寻找三个固定因素(地点，年份，栖息地)的主要影响，如下所示：
model1<-glm(tot.mass~hab, data=biom, family = Gamma(link = "log"))
anova(model1, test="Chisq")
model2<-glm(tot.mass~year, data=biom, family = Gamma(link = "log"))
anova(model2, test="Chisq")
model3<-glm(tot.ma
浏览 4提问于2017-07-09得票数 0
1回答
如何使用python比较两个回归模型。在R，能够比较使用anova(md1，md2)。但在试图使用sm.stats.anova_lm(md1，md2)的python中，显示AttributeError：'MixedLMResults‘对象没有属性’has‘。
R例子：
md1 = lmer(A ~ B + (1|Z), data = df, REML = FALSE)
md2 = lmer(A ~ B + C + (1|Z), data = df, REML = FALSE)
anova(md1,md2)
python示例：
md1= smf.mixedlm("A ~ B
浏览 1提问于2019-08-14得票数 0
回答已采纳
2回答
我有一个dataframe，它包含一个包含lm公式的列。当我对特定行[2]运行此列时，我得到了该LM的汇总输出。这可以很好地工作，但由于我在该列中有959行，所以我想编写一个for循环，以便对这些回归进行方差分析。如何指定在for循环中对该列表中的所有对象进行寻址？
为了让你更好地理解，这里有一个MWE：
数据帧：
structure(list(Week = 7:17, Category = c("2", "2", "2", "2",
"2", "2", "2", "2&#
浏览 1提问于2017-12-15得票数 0
回答已采纳
我有一份数据帧列表。我想在for循环中执行一系列操作，但在此之前，我需要提取每个数据集的字符串名称，以用作变量/数据帧名称后缀。
for(i in dflist) {
suffix<- deparse(substitute(i))
print(suffix)
}
但是，我的输出如下所示：
[1] "i"
[1] "i"
[1] "i"
我知道这是因为R的懒惰评估框架。但是，我如何绕过这个限制，并让R将dflist中数据帧的名称分配给后缀变量？
浏览 0提问于2021-11-27得票数 0
运行重复测量后，ANOVA和命名输出
RM_test <- ezANOVA(data=test_data, dv=var_test, wid = .(subject),
within = .(water_year), type = 3)
我使用str(RM_test)查看了命名对象的内部结构，并收到以下信息：
List of 3
$ ANOVA :List of 3
..$ ANOVA :'data.frame': 1 obs. of 7 variables:
.. ..$ Ef
浏览 0提问于2016-06-22得票数 0
回答已采纳
1回答
我正在测试几个组之间的线性趋势，但是，我的数据违反了组间方差相等的假设(由Levene的方差同质性检验)。
在SPSS中，随着线性趋势假设等方差的显着性，在不假定等方差的情况下，显着性有自动输出。正在做什么‘测试’或‘调整’？我能在R里做这个吗?怎么做？
图片SPSS输出：()
我很难找到SPSS到底在做什么，但这可能是某种形式的韦尔奇修正吗？
# TEST homogeneity of variance
leveneTest(ICECAP_A ~ SFMental_f, data = SCI)
P< 0.001，因此我们不接受方差同质性的零值。
# Use built-in contr
浏览 0提问于2019-06-12得票数 0
1回答
我想运行lm模型，保存模型比较结果并提取p值.我想把所有的信息保存在一个数据。以钻石数据集为例：
diamonds %>%
group_by(cut) %>%
do(model1 = lm(price~carat, data=.),
model2 = lm(price~carat+depth, data=.)) %>%
mutate(anova = anova(model2,model1)) %>%
mutate(pval= anova$'Pr(>F'[2])
我收到了下面的错误消息：
浏览 0提问于2017-08-04得票数 3
回答已采纳
这是我第一次在这里，所以放轻松！老实说，以我花在故障排除上的时间来说，我本可以手工完成每一项工作，但我下定决心要做到这一点……
我想运行71种不同的单因素方差分析，比较实验组(组)和实验组(组) 71个大脑感兴趣区(ROI)的体积。分为三组:对照组(1)，恢复组(2)和易感组(3)。
我的tibble数据帧如下所示：
> head(df)
# A tibble: 6 x 74
ID group whole_brain amygdala arbor_vitae auditory_cortex bed_nuclei
<chr> <dbl> &
浏览 0提问于2020-04-09得票数 1
2回答
我希望在R中运行多个ANOVA，所以我希望编写一个函数。
df = iris
run_anova <- function(var1,var2,df) {
fit = aov(var1 ~ var1 , df)
return(fit)
}
在虹膜数据集中，列名为"Sepal.Length“"Sepal.Width”"Petal.Length“"Petal.Width”“Petal.Width”。
假设我想在方程中使用这些列，如何将它们传递到run_anova函数中？我试着把它们当作字符串传递
run_anova("Sepal.Width&#
浏览 4提问于2020-02-24得票数 0
回答已采纳
我有一个带有丢失数据的数据帧，我将其归因于mice。
我不能完全理解我得到的输出。SSQ =平方和。但是平方和是多少呢？剩余平方和？或者总平方和-残差平方和？我假设是后者，但不确定。
df1是有意义的-它是每个变量中的组数-1。df2，我不明白。我有473个变量，在老鼠身上有20个推算。但即使是473x20也不等于下面第一列中的500万！
另外，eta2是指单向anova的结果，而partial-eta2是指我的多向anova的结果吗？
底部残差是什么？
非常感谢您的所有建议！！
imput <- mice(bdd, seed=1, pred = pred1, meth = me
浏览 0提问于2017-09-18得票数 1
我是新的R和目前我工作的参与者数据，我想运行的每一个类别的方差分析。
我使用以下代码从CSV读取我的数据：
CSV <- read.csv("data.csv", header=TRUE)
矢量CSV现在有很多。12个变量中。最后一个变量描述组，其中前11个变量是每个类别。我希望通过根据第12个变量的值将数据分成组来运行ANOVA，并对每个变量2-11运行一个ANOVA。
如何根据第12个变量将数据分成N个组，并对每个变量2-11进行方差分析？
浏览 5提问于2014-10-21得票数 0
回答已采纳
1回答
我正在尝试设置一个函数，该函数对我的主数据帧的特定子集中的选择变量(data.var1，data.var2)运行回归分析，但无论我是在函数内部还是外部定义子集，该函数都是在整个数据帧上运行，而不是仅在我希望它的子集上运行。
#Function subsetting data by temp and running regression
varReg21C <- function(data.var1,data.var2) {
data21C <- subset(allPursuit,allPursuit$temp == 21)
fitData <- lm(data.va
浏览 16提问于2020-05-30得票数 1
回答已采纳
我想要定义一个包含基本信息和统计测试结果的简单数据结构。每个测试都包含一个名称和一个结果。有多种选择可供选择，我无法思考什么选择是好的。一些简介代码：
type Calculation a = Either Failed a
data Failed = NoVariance | NotEnoughData | ...
type Name = String
type P = Double
以及正反两方面的例子：
data Test = Test Name (Calculation P)
优点:它可能是最简单的数据结构--每个测试包含一个名称和一个P值
缺点:没有类型的安全。ANOVA将被视为与
浏览 0提问于2019-02-07得票数 2
2回答
我正在对我的数据集中的多个组执行简单的线性回归。但是，我希望从每个回归中提取摘要，并将它们按组组织到主表中。我可以这样运行它(并且它是有效的)：
fit_basic <- rs2_anova %>% #Run multiple simple linear regressions
group_by(quant_method) %>%
nest() %>%
mutate(model = map(data, ~lm(recoveries ~ treatment, data = .)))
fit_basic_A <- fit_basic
浏览 0提问于2020-03-13得票数 0
1回答
我需要一些帮助，在R中执行N-way ANOVA，以捕获不同因素之间的相互依赖关系。在我的数据中，大约有100个不同的因素，我使用以下代码执行方差分析。
model.lm<-lm(y~., data=data)
anova(model.lm)
据我所知(可能我错了)，这只在每个因素上执行单因素方差分析。由于某些原因，我需要在所有100个组之间进行N-way ANOVA，即从x1到x100。我是否需要像下面这样指定每个因子，或者有一个速记符号来表示？
model.lm<-lm(y~x1*x2*x3....,x100, data=data)
anova(model.lm)
浏览 0提问于2012-11-26得票数 0
1回答
我创建了一个类似下面的脚本来做一些我称之为“加权”回归的事情：
library(plyr)
set.seed(100)
temp.df <- data.frame(uid=1:200,
bp=sample(x=c(100:200),size=200,replace=TRUE),
age=sample(x=c(30:65),size=200,replace=TRUE),
weight=sample(c(1:10),size=200,replace=TRUE
浏览 0提问于2012-04-22得票数 13
回答已采纳
我希望能够在使用SelectKBest时，通过自己执行网格搜索简历来再现学习GridSearchCV的结果。但是，我发现我的代码可以产生不同的结果。下面是一个可重复的例子：
import numpy as np
from sklearn.datasets import make_classification
from sklearn.feature_selection import SelectKBest
from sklearn.model_selection import GridSearchCV, StratifiedKFold
from sklearn.metrics import r
浏览 4提问于2021-03-05得票数 2
回答已采纳
我对R语言比较陌生，并且使用它来分析我的数据。我正在使用一个名为'agricolae‘的包和它的HSD.test组件来输出基于处理它的anova的tukeyHSD结果。我简化了下面的代码，但本质上讲，tukey代码是在一个for循环中运行的，该循环遍历了一系列因素/交互。它工作得很好，但是，我只希望循环在相应因子的anova结果显着的地方运行。(即<= 0.05)。
f1 <- as.formula(paste(m,'~ Treatment + Genotype + GENOTYPExTREATMENT', sep = ''))
anova.
浏览 3提问于2020-07-31得票数 2
我有三个因素:单词、类型和语域。在SPSS中，很容易在SPSS中进行成对比较(或简单比较)，语法是：
/EMMEANS=TABLES(word*register*type) COMPARE(type) ADJ (BONFERRONI)
它会给我一个这样的结果：
但是我如何在R中使用Multcomp包实现这一点呢？假设我有一个这样的anova结果：
library (multcomp)
anova <- aov (data, min~word*register*type)
summary (anova)
如何进行类似于SPSS的成对比较？
浏览 3提问于2016-06-22得票数 0
2回答
我在r中有这个数据，我需要用它来做一个变量。
我的问题是，我只需要Faktor 0和2中的数据，而不需要1。
FSQ$total_Sum[FSQ$Faktor=="0"]
但我不知道如何用Faktor 0,2创建一个anova。
anova_FSQ_Vor_Nach <- aov(FSQ$Faktor ~FSQ$total_Sum[FSQ$Faktor=="0"])
谢谢你的帮助。
浏览 2提问于2022-01-16得票数 1
回答已采纳
我正在尝试循环通过数据帧的每一列的TukeyHSD测试，并比较处理级别。下面是一些模拟数据，它们是我拥有的数据的简化版本(我的数据大约有350列)： df1 <- data.frame(cmpd1 = c(500,436,1,1,1,1),
cmpd2 = c(1,1,1,1,1,253),
cmpd3 = c(1,1,300,57,150,260),
treatment=c("W","W","A","A","
浏览 21提问于2021-09-13得票数 1
回答已采纳
在R中，我试图将一个数据帧的多个列折叠为两列，将第一个数据帧的列名复制到结果数据框架中的自己的列中。例如，我有以下数据框架df：
A B C D
1 2 3 4
5 6 7 8
我试图得到这个输出，这在进行ANOVA测试时很有帮助：
DV IV
A 1
A 5
B 2
B 6
C 3
C 7
D 4
D 8
我一直在手动声明新的数据框架，如下所示：
df2 <- data.frame("DV"
浏览 4提问于2015-03-30得票数 5
回答已采纳
1回答
我在R中有几个数据帧，其结构如下
> df1
messy_col_name1 messy_group_name1
numeric data "group1"
... ...
numeric data "group1"
> df2
messy_col_name2 messy_group_name2
numeric data "group2"
... ...
numeric data "group2"
.
.
.
浏览 0提问于2019-03-11得票数 0
回答已采纳
本学期我一直在使用R exams包为我的统计学入门课程创建考试。它真的是一个很棒的工具！我可以从头开始创建几个问题，并将它们导入到canvas中，而不会出现问题。然而，当我尝试导入它们时，有一些问题会给我带来问题(例如，包中包含的anova和boxplot示例)。如果我使用以下命令，可以成功导入： R> library("exams")
R> set.seed(1)
R> exams2canvas("anova.Rmd") 然而，当我尝试创建同一问题的多个版本时，有时会遇到问题： R> library("exams")
浏览 33提问于2020-11-17得票数 1
回答已采纳
我有一个非常大的数据帧(12,000行x 95列)，我想遍历每一列并运行ANOVA，将结果附加到列表中。
results = []
for i in df.columns[1:]:
mod = ols('var ~ C({})'.format(i), data = df).fit()
aov = sm.stats.anova_lm(mod, typ = 2)
results.append(aov)
我分别在几列上运行了相同的代码，并且在这个循环中也运行了前5列，这两个列都可以工作。但是，当我在整个数据集上运行时，我得到这个错误：
ValueError
浏览 15提问于2021-02-19得票数 0
3回答
到目前为止，我知道如何在R上运行ANOVA，但是我总是要复制代码来运行另一个变量的de ANOVA，我想知道是否可以将它以循环的方式传递给aov()变量的名称，并将ANOVA的结果存储在变量中，这样我就不必通过复制代码块来手动更改它们。
例如：
我想测试的变量: Z，Y，X
分类变量:治疗
VectorVariables = c(Z，Y，X)
for (i in Vector) {
AnovaZ <- aov(Z ~ Treatment) #then
AnovaY <- aov(Y ~ Treatment) # and so on..
}
这在某种程度上是可能的？
浏览 1提问于2018-09-01得票数 2
回答已采纳
3回答
这是我的小函数和数据。请注意，我想设计一个功能，而不是一般用途的个人使用。
dataf <- data.frame (A= 1:10, B= 21:30, C= 51:60, D = 71:80)
myfun <- function (dataframe, varA, varB) {
daf2 <- data.frame (A = dataframe$A*dataframe$B,
B= dataframe$C*dataframe$D)
anv1 <- lm(varA ~ varB,
浏览 0提问于2011-11-14得票数 4
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
