zoukankan      html  css  js  c++  java
  • 方差分析(R),偏eta平方 , 非中心F分布

    #计算方差分析
    #数据来源于data文件
    #wrtter: mike
    #time: 2020,10,24
    
    facA <- factor(c(rep(1,15),rep(2,15)))
    facB <- factor(c(rep(1,5),rep(2,5),rep(3,5),rep(1,5),rep(2,5),rep(3,5)))
    value <- c(2,1,4,4,2,
               3,4,4,5,1,
               5,7,7,3,6,
               8,10,12,11,6,
               6,4,3,3,2,
               3,4,4,5,3
               )
    data <- data.frame(facA,facB,value)
    
    
    #数据分析
    res <- aov(value~facA*facB)
    summary(res)
    
    #手动算出eta方,lamda
    #facA F=8.6 df=1, facB  F=5.9, df=2, residual  MS=2.6 df=24
    #计算偏eta方,也就是效应量
    eta_squ_A <- (8.6*1)/(8.6*1+24)
    eta_squ_B <- (5.9*2)/(5.9*2+24)
    
    #计算非中心F的参数lamda
    lamda_A <- 30 * (eta_squ_A)/(1-eta_squ_A)
    lamda_B <- 30*(eta_squ_B)/(1-eta_squ_B)

  • 相关阅读:
    第四次
    jsp
    2021.3.4
    第八次作业
    第七次作业
    第六周作业
    第五周作业
    第四周
    第四次jsp作业
    第二次作业
  • 原文地址:https://www.cnblogs.com/zijidefengge/p/13869917.html
Copyright © 2011-2022 走看看