方差分析,二元的,可以分为无交互效应和有交互效应。
代码如下:

import numpy as np
import pandas as pd
import scipy.stats as stats
import matplotlib.pyplot as plt
plt.style.use("bmh")import statsmodels.api as sm
from statsmodels.formula.api import olsybar = 87.5
alpha = 2.5
beta = 5A = [1]*10 + [-1]*10
B = [1]*5 + [-1] * 5 +  [1]*5 + [-1] * 5y = np.array([ybar]*20) +np.array(A) * 2.5 +np.array(B) * 5 + np.random.normal(loc=0, scale=2, size=20)df = pd.DataFrame({"group1":A, "group2":B, "y":y})
mod = ols("y~group1+group2+group1*group2", data=df).fit()ano_table = sm.stats.anova_lm(mod)
ano_tabletable = pd.pivot_table(df, index=['group1', 'group2'])
table
table2= table.reset_index()fig, ax = plt.subplots()
y = table2['y']
xall = [1, 2]
ax.plot(xall[0], y[0], 'bo')
ax.plot(xall[0], y[1], 'bo')
ax.plot(xall[1], y[2], 'bo')
ax.plot(xall[1], y[3], 'bo')ax.plot([1, 2], [y[0], y[2]], 'r')
ax.plot([1, 2], [y[1], y[3]], 'g')table2
fig, ax = plt.subplots()
ax.plot(1, y[0], 'bo')
ax.plot(1, y[2], 'bo')
ax.plot(2, y[1], 'bo')
ax.plot(2, y[3], 'bo')ax.plot([1, 2], [y[0], y[1]], 'r')
ax.plot([1, 2],[y[2], y[3]] , 'g')##########################################
#################有交互作用
##########################################y = np.array([ybar]*20) +np.array(A) * 2.5 +np.array(B) * 5 + np.array(A) * np.array(B) * 3 + np.random.normal(loc=0, scale=2, size=20)df = pd.DataFrame({"group1":A, "group2":B, "y":y})
mod = ols("y~group1+group2+group1*group2", data=df).fit()ano_table = sm.stats.anova_lm(mod)
ano_tabletable = pd.pivot_table(df, index=['group1', 'group2'])
table
table2= table.reset_index()fig, ax = plt.subplots()
y = table2['y']
xall = [1, 2]
ax.plot(xall[0], y[0], 'bo')
ax.plot(xall[0], y[1], 'bo')
ax.plot(xall[1], y[2], 'bo')
ax.plot(xall[1], y[3], 'bo')ax.plot([1, 2], [y[0], y[2]], 'r')
ax.plot([1, 2], [y[1], y[3]], 'g')table2
fig, ax = plt.subplots()
ax.plot(1, y[0], 'bo')
ax.plot(1, y[2], 'bo')
ax.plot(2, y[1], 'bo')
ax.plot(2, y[3], 'bo')ax.plot([1, 2], [y[0], y[1]], 'r')
ax.plot([1, 2],[y[2], y[3]] , 'g')

方差分析 交互效应和无交互效应相关推荐

  1. SPSS 方差分析 两因素 无交互项

    操作 Analyze→General Linear Model→Univariate 选择因变量 自变量 在Model中选择Build terms 选择自变量1的Main effect 自变量2的Ma ...

  2. 《数学建模:基于R》一一2.2 方差分析

    本节书摘来自华章计算机<数学建模:基于R>一书中的第2章,第2.2节,作者:薛 毅 更多章节内容可以访问云栖社区"华章计算机"公众号查看. 2.2 方差分析 方差分析是 ...

  3. R语言学习笔记(七)方差分析

    文章目录 写在前面 方差分析(ANOVA) 一些术语 单因素方差分析 多重比较 评估检验的假设条件 双因素方差分析 双因素方差分析的实现 可视化方法实例 正交试验设计与方差分析 正交表 考虑样本间的交 ...

  4. spss多因素方差分析

    多因素方差分析 多因素方差分析是对一个独立变量是否受一个或多个因素或变量影响而进行的方差分析.SPSS调用"Univariate"过程,检验不同水平组合之间因变量均数,由于受不同因 ...

  5. 【数据分析】双因素方差分析

    0.双因素方差分析的分类 无交互作用的方差分析 假定因素AAA和因素BBB的效应之间是相互独立的,不存在相互关系. 有交互作用的方差分析 假定因素AAA和因素BBB的结合会产生出一种新的效应. 无交互 ...

  6. python 方差齐性检验_Python数据科学:方差分析

    之前已经介绍的变量分析: ①相关分析:一个连续变量与一个连续变量间的关系. ②双样本t检验:一个二分分类变量与一个连续变量间的关系. 本次介绍: 方差分析:一个多分类分类变量与一个连续变量间的关系. ...

  7. 【统计模型】ToothGrowth数据集双因素方差分析

    目录 ToothGrowth数据集双因素方差分析 一.研究目的 二.数据来源和相关说明 三.描述性分析 3.1 样本描述 3.2 样本均值 3.3 箱线图 四.数学建模 五.结论与建议 5.1 结论 ...

  8. 第四节 单因素、多因素方差分析

    第四节 单因素.多因素方差分析 单因素方差分析 适用条件:单因素方a差分析用来检验3组以上的样本数据是否来自均值相等的总体. 原理:单因素方差分析是独立样本t检验的拓展性分析内容,独立样本t检验只能检 ...

  9. 【应用统计学】方差分析

    一.问题的提出 [例7-1] 三台设备平均灌装时间分别是15.82秒.16.67秒和14.97秒.试用样本数据检验这3台机器灌装过程的时间是否存在显著不同,以便对设备的购买做出决策.( α=0.05 ...

  10. 【R】快速实现统计推断

        如君愿,开门见山,直入主题吧! 1 t检验 1.1 单样本t检验 对总体均值的假设检验     单样本 t 检验是最基础的假设检验,利用来自总体的样本数据,推断总体均值于假设的检验值之间是否存 ...

最新文章

  1. 你哪来这么多事(四):职工信息排序
  2. javascript判断输入的值是不是正整数
  3. tensorflow-Inception-v3模型训练自己的数据代码示例
  4. SqlServer 2014还原数据库时提示:无法在已有的上还原文件,请重新发出RESTORE语句,用WITH REPLACE来覆盖原先存在的文件
  5. 日记20190415
  6. c语言实战1200例 pdf6,C语言程序设计6.6.2.pdf
  7. 【Hibernate】Hibernate系列6之HQL查询
  8. 服务器怎么控制忽略样式_使用JavaScript来编写你的CSS样式代码——JSS
  9. java static method_java 中static的几种用法
  10. Ubuntu文本编辑器适合linux平台的其它版本
  11. linux下的shell操作mysql
  12. 没有电脑可以学习python_有没有适合初学者学习Python的方法?
  13. LoRaWan 硬件和信道特点 TDMA的MAC协议优势
  14. 如何“谨慎”使用“数据驱动”的风控模型(三)——监控篇
  15. 胖哥和几个群友写了个好用的OAuth2授权服务器
  16. 教你玩转商业字体设计
  17. 数据库和数据库软件的安装
  18. springboot mybatis 分页
  19. 十个著名的思想实验-黑客帝国思想原来是(Brain in a Vat)
  20. K8S中iptables和ipvs区别

热门文章

  1. 惠普计算机图标不在桌面,惠普笔记本电脑桌面图标显示图标不正常如何还原 惠普笔记本电脑桌面图标显示图标不正常还原的方法...
  2. 七夕送女友什么礼物有意义、送女朋友实用的七夕礼物清单
  3. C++——以任意数为底的对数
  4. axure 自学入门教程
  5. Win7远程桌面 提示 登录没有成功解决方法
  6. 常见计算机网络协议汇总
  7. 微信小程序-腾讯地图报错:鉴权失败,请传入正确的key
  8. BZOJ 1213 HNOI2004 高精度开根 二分+高(Py)精(thon)度
  9. 局域网中的通信子网和资源子网
  10. 泰坦尼克号沉没之谜,用数据还原真相——Titanic获救率分析(用pyecharts)