

系統樹は自分で以下のように作成しています。  
どうやってもエラーが解消されなくて困っています。作成方法がいけないのでしょうか？？

[code]

    > newick_tree<-"(A:4,(B:2,(C:2.953469601,D:1.443532938):1.136357167):1.104667235);"
    > phylo_tree <- read.tree(text = newick_tree)
    > vcv(phylo_tree)
      A        B        C        D
    A 4 0.000000 0.000000 0.000000
    B 0 3.104667 1.104667 1.104667
    C 0 1.104667 5.194494 2.241024
    D 0 1.104667 2.241024 3.684557
    >  gls.ancova(log_lifespan ~ body_size, vcv(phylo_tree), Model, Model_I)
    Error in solve.default(Sigma) : 
      Lapack routine dgesv: system is exactly singular: U[1,1] = 0
    
[/code]


CC BY-SA 4.0
2023-11-19T07:52:31.457
96831
2023-11-19T08:31:33.020
2023-11-19T08:31:33.020
3054
60468
question
r
行列の逆行列が計算できない（Rの系統樹の処理）