组保留的标签 余下为需要删除的标签
unset($oldtag[$key]);
}
}
}
if (!empty($oldtag)) {
$tagids = array();
foreach ($oldtag as $tagid => $tagname) {
$tagids[] = $tagid;
}
well_oldtag_delete($tagids, $tid);
}
$r = well_tag_process($tid, $fid, $create_tag, $tagarr);
return $r;
}
// 删除标签和绑定的主题
function well_oldtag_delete($tagids, $tid)
{
$pagesize = count($tagids);
$arrlist = well_tag_find_by_tagids($tagids, 1, $pagesize);
$delete_tagids = array(); // 删除
$tagids = array();
$n = 0;
foreach ($arrlist as $val) {
++$n;
if (1 == $val['count']) {
// 只有一个主题
$delete_tagids[] = $val['tagid'];
} else {
$tagids[] = $val['tagid'];
}
}
!empty($delete_tagids) and well_tag_delete($delete_tagids);
$arlist = well_tag_thread_find_by_tid($tid, 1, $n);
if ($arlist) {
$ids = array();
foreach ($arlist as $val) $ids[] = $val['id'];
well_tag_thread_delete($ids);
}
!empty($tagids) and well_tag_update($tagids, array('count-' => 1));
}
// 标签数据处理 $arr=新提交的数组 $tagarr=保留的旧标签
function well_tag_process($tid, $fid, $new_tags = array(), $tagarr = array())
{
if (empty($tid)) return '';
// 新标签处理入库
if ($new_tags) {
$threadarr = array();
$tagids = array();
$i = 0;
$size = 5;
$n = count($tagarr);
$n = $n > $size ? $size : $size - $n;
foreach ($new_tags as $name) {
++$i;
$name = trim($name);
$name = stripslashes($name);
$name = strip_tags($name);
$name = str_replace(array(' ', '#', "@", "$", "%", "^", '&', '·', '<', '>', ';', '`', '~', '!', '¥', '……', ';', '?', '?', '-', '—', '_', '=', '+', '.', '{', '}', '|', ':', ':', '、', '/', '。', '[', ']', '【', '】', '‘', ' ', ' ', ' ', ' ', ' '), '', $name);
$name = htmlspecialchars($name, ENT_QUOTES);
if ($name && $i <= $n) {
// 查询标签
$read = well_tag_read_name($name);
if ($read) {
// 存在 count+1
$tagids[] = $read['tagid'];
} else {
// 入库
$arr = array('name' => $name, 'count' => 1);
$tagid = well_tag_create($arr);
FALSE === $tagid and message(-1, lang('create_failed'));
$read = array('tagid' => $tagid, 'name' => $name);
}
$tag_thread = array('tagid' => $read['tagid'], 'tid' => $tid);
$threadarr[] = $tag_thread;
$tagarr[$read['tagid']] = $read['name'];
}
}
!empty($threadarr) and tag_thread_big_insert($threadarr);
!empty($tagids) and well_tag_update($tagids, array('count+' => 1));
}
$json = empty($tagarr) ? '' : xn_json_encode($tagarr);
return $json;
}
?> $v = implode(",", $v);
$temp[] = $v;
}
// 去掉重复的字符串,也就是重复的一维数组
$temp = array_unique($temp);
// 再将拆开的数组重新组装
$output = array();
foreach ($temp as $k => $v) {
if ($stkeep) $k = $starr[$k];
if ($ndformat) {
$temparr = explode(",", $v);
foreach ($temparr as $ndkey => $ndval) $output[$k][$ndarr[$ndkey]] = $ndval;
} else $output[$k] = explode(",", $v);
}
return $output;
}
// 合并二维数组 如重复 值以第一个数组值为准
function array2_merge($array1, $array2, $key = '')
{
if (empty($array1) || empty($array2)) return NULL;
$arr = array();
foreach ($array1 as $k => $v) {
isset($v[$key]) ? $arr[$v[$key]] = array_merge($v, $array2[$k]) : $arr[] = array_merge($v, $array2[$k]);
}
return $arr;
}
/*
* 对二维数组排序 两个数组必须有一个相同的键值
* $array1 需要排序数组
* $array2 按照该数组key排序
* */
function array2_sort_key($array1, $array2, $key = '')
{
if (empty($array1) || empty($array2)) return NULL;
$arr = array();
foreach ($array2 as $k => $v) {
if (isset($v[$key]) && $v[$key] == $array1[$v[$key]][$key]) {
$arr[$v[$key]] = $array1[$v[$key]];
} else {
$arr[] = $v;
}
}
return $arr;
}
?>
群体分层软件STRUCTURE使用编程频道|福州电脑网
群体分层软件STRUCTURE使用
编程之家630
更新时间:2026-04-03 20:25:24 STRUCTURE软件使用马尔科夫链蒙特卡罗(MCMC)方法进行群体结构分析,在群体遗传学上得到了广泛应用。特别是在全基因组关联研究(GWAS)中,有时需要使用STRUCTURE软件计算获得群体结构数据。由于算法较为复杂,并且伴随几款其它软件使用,使得分析过程有些困难。这里以一个实例介绍软件安装、使用和分析计算过程,最后给出最优分群的Q矩阵。
一、STRUCTURE软件安装及应用
点击STRUCTURE网址,下载Windows版本Structure.msi安装。首先将电脑屏幕分辨率调整为较低水平(如800x600),然后运行STRUCTUE,结果如下图所示:
下面以STRUCTURE自带的sample-data数据为例进行计算分析,该数据位于安装目录下的Examples子目录。
(1) 在D盘建立一个Examples目录,将sample_data文件拷贝到此目录下。
(2) 点击File -> New Project,出现如下对话框:
以上三栏中依次填上StrucTest、D:\Examples和D:\Examples\sample_data,然后点击Next按钮。
(3) 第2步根据数据特征填写如下内容:
(4) 本例第3步不添加任何信息:
(5) 第4步勾选如下内容:
(6) 点击Finish,出现如下所填写的信息:
点击Proceed, 至此完成项目建立。
(7) 在菜单中点击Parameter Set -> New, 对话框填写如下内容:
(8) 点击OK后,命名参数集为ParameterSet01:
(9) 点击确定后,窗口情形如下:
(10) 点击Project -> Start a job, 选择参数集ParameterSet01, 填写群体数目K从1到10,重复运行(run)10次:
(11) 点击Start, 开始计算:
(12) 经过约1小时30分钟完成计算,结果为从ParameterSet01_run_1_f到ParameterSet01_run_100_f共100个文件,存储在D:\Examples\StrucTest\Paramete
本文发布于:2025-01-31,感谢您对本站的认可!
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。
本文标签:群体软件Structure
发布评论