数据 (林鼠, 包装= “猿” ) 猿 :: 作为字符。 DNAbin公司 (木老鼠[ 1 : 5 , 1 : 5 ]) #> [,1] [,2] [,3] [,4] [,5]
#>编号305
#>304号a“t”t“c”g
#>编号306“a”“t”“c”“g”
#>编号0906S“a”“t”“c”“g”
#>编号0908S“a”“t”“c”“g”
山鼠<- 林鼠[, 应用 (林鼠, 2 , 功能 (v) ! 任何 (v) == 0xf0(0xf0) ))]
###计算完整距离矩阵并打印前几行和前几列 图书馆 (公里) 伍德穆斯.kdist<- k距离 (林鼠, k个= 6 ) 打印 ( as矩阵 (木鼠.kdist)[ 1 : 7 , 1 : 7 ], 数字= 2 ) #>编号305编号304编号306编号0906S编号0908S编号0909S编号0910S
#>编号305 0.000 0.0322 0.0295 0.033 0.036 0.037 0.037
#>编号304 0.032 0.0000 0.0051 0.020 0.022 0.032 0.023
#>编号306 0.030 0.0051 0.0000 0.016 0.017 0.026 0.018
#>编号0906S 0.033 0.0202 0.0162 0.000 0.024 0.033 0.014
#>编号0908S 0.036 0.0224 0.0171 0.024 0.000 0.033 0.025
#>编号0909S 0.037 0.0322 0.0264 0.033 0.033 0.000 0.034
#>编号0910S 0.037 0.0233 0.0176 0.014 0.025 0.034 0.000 ###计算并打印嵌入的距离矩阵 禁止显示警告 ( RNG版本 ( "3.5.0" )) 设置种子 ( 999 ) 种子<- 样品 ( 1 : 15 , 尺寸= 三 ) 木老鼠.mbed<- 被动式爆炸装置 (林鼠, 种子= 种子, k个= 6 ) #>已注册的S3方法被“openssl”覆盖:
#>方法来自
#>打印字节Rcpp
打印 (woodmouse.mbed[,]), 数字= 2 ) #>编号0909S编号0913S编号304
#>编号305 0.0368 0.0391 0.0322
#>编号304 0.0322 0.0102 0.0000
#>编号306 0.0264 0.0098 0.0051
#>编号0906S 0.0332 0.0215 0.0202
#>编号0908S 0.0332 0.0273 0.0224
#>编号0909S 0.0000 0.0368 0.0322
#>编号0910S 0.0341 0.0176 0.0233
#>编号0912S 0.0242 0.0322 0.0273
#>编号0913S 0.0368 0.0000 0.0102
#>编号1103S 0.0171 0.0251 0.0202
#>编号1007S 0.0046 0.0368 0.0322
#>编号1114S 0.0451 0.0428 0.0373
#>编号1202S 0.0345 0.0176 0.0233
#>编号1206S 0.0304 0.0251 0.0202
#>编号1208S 0.0046 0.0409 0.0359
##计算两两距离矩阵 蒸馏器1<- 猿 :: 距离.dna (林鼠, 模型= “K80” ) 蒸馏器2<- k距离 (林鼠, k个= 7 ) ##构建邻接树 物理1<- 猿 :: 新泽西州 (距离1) 物理2<- 猿 :: 新泽西州 (距离2) ##以阶梯式的方式重新排列树木 物理1<- 猿 :: 用梯子把 (物理1) 物理2<- 猿 :: 用梯子把 (物理2) ##将phylo对象转换为树状图 dnd1号机组<- as.dend程序 (物理1) dnd2型<- as.dend程序 (物理2) ##绘制缠结图 dndlist(数据清单)<- 右旋糖 :: 拒绝列表 (dnd1、dnd2) 右旋糖 :: 缠结图 (dndlist, 快速= 真的 , 内部边距= 5 )
禁止显示警告 ( RNG版本 ( "3.5.0" )) 设置种子 ( 999 ) 山鼠。 OTU(OTU)<- 非处方药 (林鼠, k个= 5 , 阈值= 0.97 , 方法= “最远” , 斯特哈特= 20 ) 山鼠。 OTU(OTU) #>编号305*编号304编号306*编号0906S编号0908S编号0909S*编号0910S编号0912S
#> 3 1 1 1 1 2 1 2
#>编号0913S编号1103S编号1007S编号1114S编号1202S编号1206S编号1208S
#> 1 2 2 3 1 1 2