PopArt中介连接网络R脚本分支长度异常修复求助

针对你在绘制线粒体种群突变分析的中介连接（Median Joining, MJ）网络时，遇到的长分支过短、短分支过长，且需保留PopArt拓扑结构并体现突变步数的问题，可通过以下几种方案解决：

方案1：基于pegas+igraph生成并调整网络 #

直接在R中用与PopArt同源的MJ算法生成网络，手动绑定突变步数到分支长度，同时保留拓扑：

# 加载依赖包
library(pegas)
library(igraph)

# 1. 读取线粒体序列（FASTA格式）
dna <- read.dna("mt_sequences.fasta", format = "fasta")

# 2. 生成中介连接网络（拓扑与PopArt一致）
mj_net <- median.joining(dna)

# 3. 转换为igraph对象，提取每条边的突变步数（权重）
g <- as.igraph(mj_net)
edge_mutations <- E(g)$weight

# 4. 自定义分支长度缩放规则，平衡短/长分支显示
# 示例：突变步数为1时固定长度，步数>1时按比例放大，避免短分支冗余
scaled_length <- ifelse(edge_mutations == 1, 2, (edge_mutations - 1)*1.8 + 2)

# 5. 绘图，用权重约束布局以保留拓扑
plot(g,
     edge.length = scaled_length,
     edge.width = edge_mutations * 0.7,  # 分支宽度也对应突变步数
     vertex.size = 14,
     vertex.label = V(g)$name,
     layout = layout_with_fr(g, weights = edge_mutations)  # 力导向布局结合突变权重
)

方案2：导入PopArt导出的网络文件调整 #

如果已经在PopArt中得到满意的拓扑，可导出为NEXUS格式，在R中导入后调整分支长度：

library(phangorn)

# 1. 读取PopArt导出的NEXUS网络文件
popart_net <- read.networx("popart_mj_network.nex")

# 2. 提取突变步数作为分支长度基准
edge_mutations <- popart_net$edge.length

# 3. 调整缩放系数绘图
plot(popart_net,
     edge.length = edge_mutations * 3,  # 全局缩放系数按需调整
     edge.width = edge_mutations,
     vertex.size = 12,
     show.node.label = TRUE
)

方案3：用poppr的单倍型网络工具快速实现 #

poppr包内置的haploNet支持直接将突变步数映射为分支长度，操作更简洁：

library(poppr)

# 1. 生成单倍型网络
hap_net <- haploNet(dna)

# 2. 绘图时直接绑定突变步数到分支属性
plot(hap_net,
     edge.length = "mutations",  # 分支长度严格对应突变步数
     edge.width = "mutations",
     vertex.size = 10,
     layout = layout_with_fr(hap_net, weights = haploNet.stats(hap_net)$mutations)
)

关键注意事项 #

• 拓扑一致性：确保R中使用的MJ算法与PopArt一致（pegas和poppr的MJ实现均遵循标准算法），或直接导入PopArt导出的网络文件。
• 分支缩放逻辑：避免使用默认的自动缩放，通过自定义规则或直接绑定突变步数，平衡短分支（如1步突变）的冗余显示和长分支的细节展示。
• 布局约束：使用带权重的力导向布局（如layout_with_fr），让节点位置随突变步数调整，保留PopArt的拓扑结构。

内容的提问来源于stack exchange，提问作者n rank