<?xml version="1.0" encoding="UTF-8"?><rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>小生这厢有礼了(BioFaceBook Personal Blog) &#187; 兴趣杂项</title>
	<atom:link href="http://www.biofacebook.com/?cat=14&#038;feed=rss2" rel="self" type="application/rss+xml" />
	<link>http://www.biofacebook.com</link>
	<description>记录生物信息学点滴足迹（NGS,Genome,Meta,Linux)</description>
	<lastBuildDate>Sun, 23 Aug 2020 03:28:53 +0000</lastBuildDate>
	<language>en-US</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>https://wordpress.org/?v=4.1.41</generator>
	<item>
		<title>ssh, scp and rsync</title>
		<link>http://www.biofacebook.com/?p=1451</link>
		<comments>http://www.biofacebook.com/?p=1451#comments</comments>
		<pubDate>Sun, 23 Feb 2020 08:38:06 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[Linux相关]]></category>
		<category><![CDATA[兴趣杂项]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1451</guid>
		<description><![CDATA[<p>[Admin.DESKTOP-7JT504C] ➤ rsync -P &#8211;rsh=ssh /drives/d/Kraken_12.tar.gz cityu_jhli_1@172.16.22.11:/BIGDATA1/cityu_jhli_1/mhyleung/database/findfungi Warning: Permanently added &#8216;172.16.22.11&#8217; (RSA) to the list of known hosts.</p> <p>rsync -P -avz -e &#8220;ssh -p5566 -i /drives/C/Users/Admin/Desktop/cityu_jhli_1.id&#8221; /drives/d/Kraken_12.tar.gz cityu_jhli_1@172.16.22.11:/BIGDATA1/cityu_jhli_1/mhyleung/database/findfungi</p> <p>[Admin.DESKTOP-7JT504C] ➤ scp -P 5566 -i /drives/C/Users/Admin/Desktop/cityu_jhli_1.id -r /drives/d/Kraken_12.tar.gz cityu_jhli_1@172.16.22.11:/BIGDATA1/cityu_jhli_1/mhyleung/database/findfungi</p> <p>&#160;</p> <p>ssh -p5566 -i /drives/C/Users/Admin/Desktop/cityu_jhli_1.id cityu_jhli_1@172.16.22.11</p> ]]></description>
				<content:encoded><![CDATA[<p>[Admin.DESKTOP-7JT504C] ➤ rsync -P &#8211;rsh=ssh /drives/d/Kraken_12.tar.gz cityu_jhli_1@172.16.22.11:/BIGDATA1/cityu_jhli_1/mhyleung/database/findfungi Warning: Permanently added &#8216;172.16.22.11&#8217; (RSA) to the list of known hosts.</p>
<p>rsync -P -avz -e &#8220;ssh -p5566 -i /drives/C/Users/Admin/Desktop/cityu_jhli_1.id&#8221; /drives/d/Kraken_12.tar.gz cityu_jhli_1@172.16.22.11:/BIGDATA1/cityu_jhli_1/mhyleung/database/findfungi</p>
<p>[Admin.DESKTOP-7JT504C] ➤ scp -P 5566 -i /drives/C/Users/Admin/Desktop/cityu_jhli_1.id -r /drives/d/Kraken_12.tar.gz cityu_jhli_1@172.16.22.11:/BIGDATA1/cityu_jhli_1/mhyleung/database/findfungi</p>
<p>&nbsp;</p>
<p>ssh -p5566 -i /drives/C/Users/Admin/Desktop/cityu_jhli_1.id   cityu_jhli_1@172.16.22.11</p>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1451</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>ClinVAP: A reporting strategy from variants to therapeutic options</title>
		<link>http://www.biofacebook.com/?p=1424</link>
		<comments>http://www.biofacebook.com/?p=1424#comments</comments>
		<pubDate>Mon, 16 Dec 2019 09:37:26 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[兴趣杂项]]></category>
		<category><![CDATA[生物信息]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1424</guid>
		<description><![CDATA[<p>https://academic.oup.com/bioinformatics/advance-article/doi/10.1093/bioinformatics/btz924/5674039</p> <p>&#160;</p> ClinVAP: A reporting strategy from variants to therapeutic options <p>&#160;</p> Abstract Motivation <p>Next-generation sequencing (NGS) has become routine in oncology and opens up new avenues of therapies, particularly in personalized oncology setting. An increasing number of cases also implies a need for a more robust, automated, and reproducible processing of long lists of [...]]]></description>
				<content:encoded><![CDATA[<p><a href="https://academic.oup.com/bioinformatics/advance-article/doi/10.1093/bioinformatics/btz924/5674039">https://academic.oup.com/bioinformatics/advance-article/doi/10.1093/bioinformatics/btz924/5674039</a></p>
<p>&nbsp;</p>
<h1 class="wi-article-title article-title-main">ClinVAP: A reporting strategy from variants to therapeutic options</h1>
<p>&nbsp;</p>
<h2 id="190031274" class="abstract-title">Abstract</h2>
<section class="abstract"></section>
<section class="sec">
<div class="title">Motivation</div>
<p>Next-generation sequencing (NGS) has become routine in oncology and opens up new avenues of therapies, particularly in personalized oncology setting. An increasing number of cases also implies a need for a more robust, automated, and reproducible processing of long lists of variants for cancer diagnosis and therapy. While solutions for the large-scale analysis of somatic variants have been implemented, existing solutions often have issues with reproducibility, scalability, and interoperability.</p>
</section>
<section class="sec">
<div class="title">Results</div>
<p>ClinVAP is an automated pipeline which annotates, filters, and prioritizes somatic single nucleotide variants (SNVs) provided in variant call format. It augments the variant information with documented or predicted clinical effect. These annotated variants are prioritized based on driver gene status and druggability. ClinVAP is available as a fully containerized, self-contained pipeline maximizing reproducibility and scalability allowing the analysis of larger scale data. The resulting JSON-based report is suited for automated downstream processing, but ClinVAP can also automatically render the information into a user-defined template to yield a human-readable report.</p>
</section>
<section class="sec">
<div class="title">Availability and Implementation</div>
<p>ClinVAP is available at <a class="link link-uri" href="https://github.com/PersonalizedOncology/ClinVAP" target="">https://github.com/PersonalizedOncology/ClinVAP</a></p>
</section>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1424</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>dictionary to csv</title>
		<link>http://www.biofacebook.com/?p=1416</link>
		<comments>http://www.biofacebook.com/?p=1416#comments</comments>
		<pubDate>Thu, 07 Nov 2019 08:58:47 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[兴趣杂项]]></category>
		<category><![CDATA[脚本语言]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1416</guid>
		<description><![CDATA[<p>#!/usr/bin/env python import os,re,sys,string,commands,getopt,subprocess,glob,csv import prettytable as pt from os import path #SL335752_kneaddata_paired_1_kneaddata_paired_1.fastq.gz_rep.mpa.txt_bracken.txt #SL311013_1_kneaddata_paired_1_kneaddata_paired_1.fastq.gz_rep_mpa.txt_bracken.txt def main(): dic = {} unique_speciesname = [] speciesname = [] samplenames = [] for d in os.listdir(&#8216;.&#8217;): print(d) a = [] a = re.split(&#8220;_kneaddata_paired_1_&#8221;,d) if len(a) &#62;=2: samplename = a[0] dic[a[0]] = {} #os.rename(d, newname) fh = open(d, &#8216;r&#8217;) fhlines [...]]]></description>
				<content:encoded><![CDATA[<p>#!/usr/bin/env python<br />
import os,re,sys,string,commands,getopt,subprocess,glob,csv<br />
import prettytable as pt<br />
from os import path<br />
#SL335752_kneaddata_paired_1_kneaddata_paired_1.fastq.gz_rep.mpa.txt_bracken.txt<br />
#SL311013_1_kneaddata_paired_1_kneaddata_paired_1.fastq.gz_rep_mpa.txt_bracken.txt<br />
def main():<br />
dic = {}<br />
unique_speciesname = []<br />
speciesname = []<br />
samplenames = []<br />
for d in os.listdir(&#8216;.&#8217;):<br />
print(d)<br />
a = []<br />
a = re.split(&#8220;_kneaddata_paired_1_&#8221;,d)<br />
if len(a) &gt;=2:<br />
samplename = a[0]<br />
dic[a[0]] = {}<br />
#os.rename(d, newname)<br />
fh = open(d, &#8216;r&#8217;)<br />
fhlines = fh.readlines()<br />
fh.close()</p>
<p>for line in fhlines:<br />
line = line.strip()<br />
if re.search(&#8220;name&#8221;, line):<br />
continue<br />
else:<br />
b = []<br />
b = re.split(&#8220;\t&#8221;,line)<br />
speciesname.append(b[0])<br />
length = len(b)<br />
dic[a[0]][b[0]] = b[length-1]<br />
unique = [unique_speciesname.append(x) for x in speciesname if x not in unique_speciesname]<br />
for sample_name in dic.keys():<br />
samplenames.append(sample_name)<br />
for name in unique_speciesname:<br />
#print name<br />
if dic[sample_name].has_key(name):<br />
print_line = name + &#8220;\t&#8221; + dic[sample_name][name]<br />
else:<br />
dic[sample_name][name] = &#8220;0&#8221;<br />
print_line = name + &#8220;\t&#8221; + dic[sample_name][name]</p>
<p>#print dic<br />
csv_columns = [&#8216;Species&#8217;]<br />
for ele in samplenames:<br />
csv_columns.append(ele)<br />
print csv_columns<br />
#csv_columns = [&#8216;No&#8217;,&#8217;Name&#8217;,&#8217;Country&#8217;]<br />
csv_file = &#8220;sample_species.csv&#8221;<br />
#{&#8216;SL311014&#8242;: {&#8216;Acinetobacter sp. WCHA55&#8242;: &#8216;0.00029&#8217;, &#8216;Streptococcus sp. oral taxon 431&#8242;: &#8216;0.00007&#8217;, &#8216;Bacillus velezensis': &#8216;0.00001&#8217;, &#8216;Ahniella affigens': &#8216;0.00003&#8217;, &#8216;Arsenicicoccus sp. oral taxon 190&#8242;: &#8216;0.00077&#8217;, &#8216;Aureimonas sp. AU20&#8242;: &#8216;0&#8217;, &#8216;Mycobacterium sp. MS1601&#8242;: &#8216;0&#8217;, &#8216;Acinetobacter sp. WCHAc010052&#8242;: &#8216;0.00003&#8217;, &#8216;Candidatus Micrarchaeota archaeon Mia14&#8242;: &#8216;0.00000&#8217;, &#8216;Pseudomonas sp. MT-1&#8242;: &#8216;0.00005&#8217;, &#8216;Halothece sp. PCC 7418&#8242;: &#8216;0.00001&#8217;}, &#8216;SL311013&#8242;: {&#8216;Acinetobacter sp. WCHA55&#8242;: &#8216;0&#8217;, &#8216;Streptococcus sp. oral taxon 431&#8242;: &#8216;0.00009&#8217;, &#8216;Bacillus velezensis': &#8216;0.00004&#8217;, &#8216;Candidatus Micrarchaeota archaeon Mia14&#8242;: &#8216;0&#8217;, &#8216;Ahniella affigens': &#8216;0.00002&#8217;, &#8216;Arsenicicoccus sp. oral taxon 190&#8242;: &#8216;0.00059&#8217;, &#8216;Mycobacterium sp. MS1601&#8242;: &#8216;0.00029&#8217;, &#8216;Acinetobacter sp. WCHAc010052&#8242;: &#8216;0.00001&#8217;, &#8216;Aureimonas sp. AU20&#8242;: &#8216;0.00026&#8217;, &#8216;Pseudomonas sp. MT-1&#8242;: &#8216;0.00003&#8217;, &#8216;Halothece sp. PCC 7418&#8242;: &#8216;0.00001&#8217;}, &#8216;SL311012&#8242;: {&#8216;Acinetobacter sp. WCHA55&#8242;: &#8216;0&#8217;, &#8216;Streptococcus sp. oral taxon 431&#8242;: &#8216;0.00008&#8217;, &#8216;Bacillus velezensis': &#8216;0.00001&#8217;, &#8216;Candidatus Micrarchaeota archaeon Mia14&#8242;: &#8216;0&#8217;, &#8216;Ahniella affigens': &#8216;0.00001&#8217;, &#8216;Arsenicicoccus sp. oral taxon 190&#8242;: &#8216;0.00054&#8217;, &#8216;Mycobacterium sp. MS1601&#8242;: &#8216;0.00026&#8217;, &#8216;Acinetobacter sp. WCHAc010052&#8242;: &#8216;0.00002&#8217;, &#8216;Aureimonas sp. AU20&#8242;: &#8216;0.00024&#8217;, &#8216;Pseudomonas sp. MT-1&#8242;: &#8216;0.00003&#8217;, &#8216;Halothece sp. PCC 7418&#8242;: &#8216;0.00000&#8217;}}<br />
dict_data_mine = []<br />
#dic_ele = {}<br />
for name in unique_speciesname:<br />
dic_ele = {}<br />
dic_ele[&#8220;Species&#8221;]=name<br />
for key in dic.keys():<br />
#dic_ele[&#8220;speciesname&#8221;]=&gt;name<br />
dic_ele[key]=dic[key][name]<br />
dict_data_mine.append(dic_ele)<br />
dict_data = [<br />
{&#8216;No': 1, &#8216;Name': &#8216;Alex&#8217;, &#8216;Country': &#8216;India&#8217;},<br />
{&#8216;No': 2, &#8216;Name': &#8216;Ben&#8217;, &#8216;Country': &#8216;USA&#8217;},<br />
{&#8216;No': 3, &#8216;Name': &#8216;Shri Ram&#8217;, &#8216;Country': &#8216;India&#8217;},<br />
{&#8216;No': 4, &#8216;Name': &#8216;Smith&#8217;, &#8216;Country': &#8216;USA&#8217;},<br />
{&#8216;No': 5, &#8216;Name': &#8216;Yuva Raj&#8217;, &#8216;Country': &#8216;India&#8217;},<br />
]</p>
<p>#print dict_data_mine<br />
try:<br />
with open(csv_file, &#8216;w&#8217;) as csvfile:<br />
writer = csv.DictWriter(csvfile, fieldnames=csv_columns)<br />
writer.writeheader()<br />
for data in dict_data_mine:<br />
writer.writerow(data)<br />
except IOError:<br />
print(&#8220;I/O error&#8221;)<br />
if __name__== &#8216;__main__':<br />
main()</p>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1416</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Ubuntu中v2ray客户端配置实例</title>
		<link>http://www.biofacebook.com/?p=1411</link>
		<comments>http://www.biofacebook.com/?p=1411#comments</comments>
		<pubDate>Sat, 12 Oct 2019 07:51:14 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[兴趣杂项]]></category>
		<category><![CDATA[服务器管理]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1411</guid>
		<description><![CDATA[<p>首先使用bash &#60;(curl -L -s https://install.direct/go.sh)来快捷安装v2ray，如下：</p> root@vm:~# bash &#60;(curl -L -s https://install.direct/go.sh) Installing V2Ray v4.18.0 on x86_64 Downloading V2Ray: https://github.com/v2ray/v2ray-core/releases/download/v4.18.0/v2ray-linux-64.zip % Total % Received % Xferd Average Speed Time Time Time Current Dload Upload Total Spent Left Speed 100 608 0 608 0 0 581 0 --:--:-- 0:00:01 --:--:-- 581 100 10.5M 100 10.5M 0 0 [...]]]></description>
				<content:encoded><![CDATA[<p>首先使用<code>bash &lt;(curl -L -s https://install.direct/go.sh)</code>来快捷安装v2ray，如下：</p>
<pre><code>root@vm:~# bash &lt;(curl -L -s https://install.direct/go.sh)
Installing V2Ray v4.18.0 on x86_64
Downloading V2Ray: https://github.com/v2ray/v2ray-core/releases/download/v4.18.0/v2ray-linux-64.zip
  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100   608    0   608    0     0    581      0 --:--:--  0:00:01 --:--:--   581
100 10.5M  100 10.5M    0     0   172k      0  0:01:02  0:01:02 --:--:--  194k
Extracting V2Ray package to /tmp/v2ray.
Archive:  /tmp/v2ray/v2ray.zip
  inflating: /tmp/v2ray/config.json  
   creating: /tmp/v2ray/doc/
  inflating: /tmp/v2ray/doc/readme.md  
  inflating: /tmp/v2ray/geoip.dat    
  inflating: /tmp/v2ray/geosite.dat  
   creating: /tmp/v2ray/systemd/
  inflating: /tmp/v2ray/systemd/v2ray.service  
   creating: /tmp/v2ray/systemv/
  inflating: /tmp/v2ray/systemv/v2ray  
  inflating: /tmp/v2ray/v2ctl        
 extracting: /tmp/v2ray/v2ctl.sig    
  inflating: /tmp/v2ray/v2ray        
 extracting: /tmp/v2ray/v2ray.sig    
  inflating: /tmp/v2ray/vpoint_socks_vmess.json  
  inflating: /tmp/v2ray/vpoint_vmess_freedom.json  
PORT:51332
UUID:7378f6a4-790a-11e9-8f9e-2a86e4085a59
Created symlink /etc/systemd/system/multi-user.target.wants/v2ray.service → /etc/systemd/system/v2ray.service.
V2Ray v4.18.0 is installed.
</code></pre>
<p>然后编辑<code>/etc/v2ray/config.json</code>文件，如下设置：</p>
<pre><code>{
  "inbounds": [{
    "port": 10808,  // SOCKS 代理端口，在浏览器中需配置代理并指向这个端口
    "listen": "127.0.0.1",
    "protocol": "socks",
    "settings": {
      "udp": true
    }
  }],
  "outbounds": [{
    "protocol": "vmess",
    "settings": {
      "vnext": [{
        "address": "server", // 服务器地址，请修改为你自己的服务器 ip 或域名
        "port": 10086,  // 服务器端口
        "users": [{ "id": "b831381d-6324-4d53-ad4f-8cda48b30811" }]
      }]
    }
  },{
    "protocol": "freedom",
    "tag": "direct",
    "settings": {}
  }],
  "routing": {
    "domainStrategy": "IPOnDemand",
    "rules": [{
      "type": "field",
      "ip": ["geoip:private"],
      "outboundTag": "direct"
    }]
  }
}
</code></pre>
<p>编辑完成后保存，重新启动v2ray</p>
<pre><code>root@vm:~# service v2ray stop
root@vm:~# service v2ray start
root@vm:~# service v2ray status
● v2ray.service - V2Ray Service
   Loaded: loaded (/etc/systemd/system/v2ray.service; enabled; vendor preset: en
   Active: active (running) since Sat 2019-05-18 08:58:43 CST; 5s ago
 Main PID: 8025 (v2ray)
    Tasks: 7 (limit: 2311)
   CGroup: /system.slice/v2ray.service
           └─8025 /usr/bin/v2ray/v2ray -config /etc/v2ray/config.json

5月 18 08:58:43 vm systemd[1]: Started V2Ray Service.
5月 18 08:58:43 vm v2ray[8025]: V2Ray 4.18.0 (Po) 20190228
5月 18 08:58:43 vm v2ray[8025]: A unified platform for anti-censorship.
5月 18 08:58:44 vm v2ray[8025]: 2019/05/18 08:58:44 [Warning] v2ray.com/core: V2
</code></pre>
<p>然后Firefox设置代理如下：</p>
<p>设置-常规-网络设置 勾选<strong>手动代理配置</strong>，在<strong>SOCKS主机</strong>中填入<code>127.0.0.1</code>本地IP和<strong>端口</strong>，协议勾选<strong>SOCKS_v5</strong> 建议勾选<strong>使用SOCKSv5时代理DNS</strong></p>
<p><img src="https://unixetc.com/imgs/3178501606.webp" alt="20190518090412.webp" /></p>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1411</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Good software</title>
		<link>http://www.biofacebook.com/?p=1385</link>
		<comments>http://www.biofacebook.com/?p=1385#comments</comments>
		<pubDate>Wed, 24 Jul 2019 06:46:13 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[Linux相关]]></category>
		<category><![CDATA[二代测序]]></category>
		<category><![CDATA[兴趣杂项]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1385</guid>
		<description><![CDATA[<p>multiqc ranger ployly</p> <p>https://github.com/MultiQC</p> <p>https://github.com/ranger/ranger</p> <p>https://zhuanlan.zhihu.com/p/34369349</p> ]]></description>
				<content:encoded><![CDATA[<p>multiqc   ranger    ployly</p>
<p>https://github.com/MultiQC</p>
<p>https://github.com/ranger/ranger</p>
<p>https://zhuanlan.zhihu.com/p/34369349</p>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1385</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Histograms and Density Plots</title>
		<link>http://www.biofacebook.com/?p=1359</link>
		<comments>http://www.biofacebook.com/?p=1359#comments</comments>
		<pubDate>Tue, 26 Feb 2019 01:54:25 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[兴趣杂项]]></category>
		<category><![CDATA[生物信息]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1359</guid>
		<description><![CDATA[Histograms and Density Plots Histograms <p>You can create histograms with the function hist(x) where x is a numeric vector of values to be plotted. The option freq=FALSE plots probability densities instead of frequencies. The option breaks= controls the number of bins.</p> <p># Simple Histogram hist(mtcars$mpg)</p> <p> click to view</p> <p># Colored Histogram with Different Number [...]]]></description>
				<content:encoded><![CDATA[<h1>Histograms and Density Plots</h1>
<h2>Histograms</h2>
<p>You can create histograms with the function <strong>hist(</strong><em>x</em><strong>)</strong> where <em>x</em> is a numeric vector of values to be plotted. The option <strong>freq=FALSE</strong> plots probability densities instead of frequencies. The option <strong>breaks=</strong> controls the number of bins.</p>
<p><code># Simple Histogram<br />
hist(mtcars$mpg)</code></p>
<p><a href="https://www.statmethods.net/graphs/images/histogram1.jpg"><img src="https://www.statmethods.net/graphs/images/smhistogram1.jpg" alt="simple histogram" width="103" height="103" /></a> click to view</p>
<p><code># Colored Histogram with Different Number of Bins<br />
hist(mtcars$mpg, breaks=12, col="red")</code></p>
<p><a href="https://www.statmethods.net/graphs/images/histogram2.jpg"><img src="https://www.statmethods.net/graphs/images/smhistogram2.jpg" alt="colored histogram" width="103" height="103" /></a> click to view</p>
<p><code># Add a Normal Curve (Thanks to Peter Dalgaard)<br />
x &lt;- mtcars$mpg<br />
h&lt;-hist(x, breaks=10, col="red", xlab="Miles Per Gallon",<br />
main="Histogram with Normal Curve")<br />
xfit&lt;-seq(min(x),max(x),length=40)<br />
yfit&lt;-dnorm(xfit,mean=mean(x),sd=sd(x))<br />
yfit &lt;- yfit*diff(h$mids[1:2])*length(x)<br />
lines(xfit, yfit, col="blue", lwd=2)</code></p>
<p><a href="https://www.statmethods.net/graphs/images/histogram3.jpg"><img src="https://www.statmethods.net/graphs/images/smhistogram3.jpg" alt="histogram with normal curve" width="103" height="103" /></a> click to view</p>
<p>Histograms can be a poor method for determining the shape of a distribution because it is so strongly affected by the number of bins used.</p>
<p>To practice making a density plot with the hist() function, try <a href="https://campus.datacamp.com/courses/cleaning-data-in-r/chapter-1-introduction-and-exploring-raw-data?ex=10">this exercise.</a></p>
<h2>Kernel Density Plots</h2>
<p>Kernal density plots are usually a much more effective way to view the distribution of a variable. Create the plot using <strong>plot(density(</strong><em>x</em><strong>)) </strong>where<em> x </em>is a numeric vector.</p>
<p><code># Kernel Density Plot<br />
d &lt;- density(mtcars$mpg) # returns the density data<br />
plot(d) # plots the results</code></p>
<p><a href="https://www.statmethods.net/graphs/images/density1.jpg"><img src="https://www.statmethods.net/graphs/images/smdensity1.jpg" alt="simple density plot" width="103" height="103" /></a> click to view</p>
<p><code># Filled Density Plot<br />
d &lt;- density(mtcars$mpg)<br />
plot(d, main="Kernel Density of Miles Per Gallon")<br />
polygon(d, col="red", border="blue")</code></p>
<p><a href="https://www.statmethods.net/graphs/images/density2.jpg"><img src="https://www.statmethods.net/graphs/images/smdensity2.jpg" alt="colored density plot" width="103" height="103" /></a> click to view</p>
<h2>Comparing Groups VIA Kernal Density</h2>
<p>The <strong>sm.density.compare( ) </strong>function in the <strong><a href="http://cran.r-project.org/web/packages/sm/index.html">sm</a> </strong>package allows you to superimpose the kernal density plots of two or more groups. The format is <strong>sm.density.compare(</strong><em>x</em><strong>, </strong><em>factor</em>) where <em>x</em> is a numeric vector and <em>factor</em> is the grouping variable.</p>
<p><code># Compare MPG distributions for cars with<br />
# 4,6, or 8 cylinders<br />
library(sm)<br />
attach(mtcars)</p>
<p># create value labels<br />
cyl.f &lt;- factor(cyl, levels= c(4,6,8),<br />
labels = c("4 cylinder", "6 cylinder", "8 cylinder"))</p>
<p># plot densities<br />
sm.density.compare(mpg, cyl, xlab="Miles Per Gallon")<br />
title(main="MPG Distribution by Car Cylinders")</p>
<p># add legend via mouse click<br />
colfill&lt;-c(2:(2+length(levels(cyl.f))))<br />
legend(locator(1), levels(cyl.f), fill=colfill)</code></p>
<p><a href="https://www.statmethods.net/graphs/images/density3.png"><img src="https://www.statmethods.net/graphs/images/smdensity3.jpg" alt="comparing densities" width="103" height="103" /></a> click to view</p>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1359</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>ARGs-OAP: 抗性基因在线分析工具(转）</title>
		<link>http://www.biofacebook.com/?p=1348</link>
		<comments>http://www.biofacebook.com/?p=1348#comments</comments>
		<pubDate>Fri, 28 Dec 2018 05:53:21 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[兴趣杂项]]></category>
		<category><![CDATA[生物信息]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1348</guid>
		<description><![CDATA[ <p>ARGs_OAP_v2.0（步骤1）：https://github.com/biofuture/Ublastx_stageone ARGs-OAP在线分析网站（步骤2）： http://smile.hku.hk/SARGs</p> 无处不在的抗性基因 图片.png <p>环境中抗生素抗性基因（ARGs）的来源： 随机突变或表达潜在抗性基因等方式使细菌体内基因组上存在的抗性基因原型、准抗性基因或潜在抗性基因被表达出来，从而使细菌获得的抗生素抗性。 抗生素在人和动物肠道内诱导产生耐药菌，这些编码ARGs的耐药菌经由粪便排出并进入环境中，是环境中ARGs的重要来源。 抗性基因的水平转移是抗性基因在环境中传播的的主要机制，通过将包含抗性基因的质粒、转座子、整合子作为载体，通过细菌之间细胞与细胞的接触，将抗性基因从载体细胞转移到受体细胞。</p> 如何检测环境中抗生素抗性基因（ARGs）： PCR技术&#8212;定性。 qPCR技术&#8212;定量。 宏基因组测序：以环境样品中的整个微生物群体基因组为研究对象，检测环境样本微生物中的物种组成、丰度，基因预测、基因丰度，利用数据库进行注释，得到样本中ARGs的种类和丰度与样本的相关性。 ARDB数据库：主要包含细菌病原菌的多种抗性基因数据，不能为环境样本宏基因组数据提供详细的ARG概况（即对每个检测到的ARG提供type/subtype的ARG分类信息和丰度信息）。 CARD数据库：以Antibiotic Resistance Ontology（ARO）为分类单位的形式所构建，ARO用于关联抗生素模块及其目标、抗性机制、基因变异等信息。 ResFinder：需要较长的查询reads。对于在ResFinder中被检测为ARG的序列，其必须至少覆盖数据库中匹配ARG长度的五分之二，具有不小于50％的相似性。 ARGO：侧重于万古霉素和β-内酰胺抗性基因。 ARG-ANNOT：设计用于检测细菌基因组中的ARG而不是环境样品。 构建ARG综合数据库SARG v1.0 整合CARD和ARDB数据库 CARD数据库2,513条序列； ARDB数据库7,828条序列； 去除586条共享序列； SARG包含4246条ARGs参考序列。 去除非ARG序列 去除冗余序列（完整蛋白质序列具有100％同一性） 去除与SNP相关的ARG序列 去除描述为“假定蛋白质”或“未命名蛋白质”的序列 构建结构化ARG数据库SARG 构建ARG综合数据库SARG v2.0 图片.png 使用SARG v1.0作为从NCBI-NR获取潜在ARG序列的种子。 NCBI-NR序列BLASTP比对SARG v1.0数据库（e-value:1e-7, identity: 90％、80％、70％）； levels: Accurate, Moderate and Loose 。 基于序列相似度或关键字匹配将ARG序列分配给不同的Subtype。 合并时，删除有多个分类结果的序列，只保留具有匹配分类（type和subtype）的序列。 <p>Number of ARGs reference genes [...]]]></description>
				<content:encoded><![CDATA[<div>
<div>
<p>ARGs_OAP_v2.0（步骤1）：<a href="https://github.com/biofuture/Ublastx_stageone" target="_blank" rel="nofollow">https://github.com/biofuture/Ublastx_stageone</a><br />
ARGs-OAP在线分析网站（步骤2）： <a href="http://smile.hku.hk/SARGs" target="_blank" rel="nofollow">http://smile.hku.hk/SARGs</a></p>
<h6>无处不在的抗性基因</h6>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="960" data-height="548"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-8916eb4c325acfd1.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/960/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-8916eb4c325acfd1.png" data-original-width="960" data-original-height="548" data-original-format="image/png" data-original-filesize="103004" /></div>
</div>
<div class="image-caption">图片.png</div>
</div>
<p>环境中抗生素抗性基因（ARGs）的来源：<br />
随机突变或表达潜在抗性基因等方式使细菌体内基因组上存在的抗性基因原型、准抗性基因或潜在抗性基因被表达出来，从而使细菌获得的抗生素抗性。<br />
抗生素在人和动物肠道内诱导产生耐药菌，这些编码ARGs的耐药菌经由粪便排出并进入环境中，是环境中ARGs的重要来源。<br />
抗性基因的水平转移是抗性基因在环境中传播的的主要机制，通过将包含抗性基因的质粒、转座子、整合子作为载体，通过细菌之间细胞与细胞的接触，将抗性基因从载体细胞转移到受体细胞。</p>
<h6>如何检测环境中抗生素抗性基因（ARGs）：</h6>
<ul>
<li>PCR技术&#8212;定性。</li>
<li>qPCR技术&#8212;定量。</li>
<li>宏基因组测序：以环境样品中的整个微生物群体基因组为研究对象，检测环境样本微生物中的物种组成、丰度，基因预测、基因丰度，利用数据库进行注释，得到样本中ARGs的种类和丰度与样本的相关性。</li>
</ul>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="448" data-height="91"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-d29bae0fb6b71e22.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/448/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-d29bae0fb6b71e22.png" data-original-width="448" data-original-height="91" data-original-format="image/png" data-original-filesize="5386" /></div>
</div>
<div class="image-caption"></div>
</div>
<ul>
<li>ARDB数据库：主要包含细菌病原菌的多种抗性基因数据，不能为环境样本宏基因组数据提供详细的ARG概况（即对每个检测到的ARG提供type/subtype的ARG分类信息和丰度信息）。</li>
<li>CARD数据库：以Antibiotic Resistance Ontology（ARO）为分类单位的形式所构建，ARO用于关联抗生素模块及其目标、抗性机制、基因变异等信息。<br />
ResFinder：需要较长的查询reads。对于在ResFinder中被检测为ARG的序列，其必须至少覆盖数据库中匹配ARG长度的五分之二，具有不小于50％的相似性。</li>
<li>ARGO：侧重于万古霉素和β-内酰胺抗性基因。<br />
ARG-ANNOT：设计用于检测细菌基因组中的ARG而不是环境样品。</li>
</ul>
<h6>构建ARG综合数据库SARG v1.0</h6>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="547" data-height="576"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-6463988bf8d38e6e.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/547/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-6463988bf8d38e6e.png" data-original-width="547" data-original-height="576" data-original-format="image/png" data-original-filesize="45949" /></div>
</div>
<div class="image-caption"></div>
</div>
<ol>
<li>整合CARD和ARDB数据库<br />
CARD数据库2,513条序列；<br />
ARDB数据库7,828条序列；<br />
去除586条共享序列；<br />
SARG包含4246条ARGs参考序列。</li>
<li>去除非ARG序列</li>
<li>去除冗余序列（完整蛋白质序列具有100％同一性）</li>
<li>去除与SNP相关的ARG序列</li>
<li>去除描述为“假定蛋白质”或“未命名蛋白质”的序列</li>
<li>构建结构化ARG数据库SARG</li>
</ol>
<h6>构建ARG综合数据库SARG v2.0</h6>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="606" data-height="584"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-d1e1742d757f33bd.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/606/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-d1e1742d757f33bd.png" data-original-width="606" data-original-height="584" data-original-format="image/png" data-original-filesize="55799" /></div>
</div>
<div class="image-caption">图片.png</div>
</div>
<ol>
<li>使用SARG v1.0作为从NCBI-NR获取潜在ARG序列的种子。</li>
<li>NCBI-NR序列BLASTP比对SARG v1.0数据库（e-value:1e-7, identity: 90％、80％、70％）； levels: Accurate, Moderate and Loose 。</li>
<li>基于序列相似度或关键字匹配将ARG序列分配给不同的Subtype。</li>
<li>合并时，删除有多个分类结果的序列，只保留具有匹配分类（type和subtype）的序列。</li>
</ol>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="616" data-height="470"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-6459880dfed2e0fd.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/616/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-6459880dfed2e0fd.png" data-original-width="616" data-original-height="470" data-original-format="image/png" data-original-filesize="39448" /></div>
</div>
<div class="image-caption"></div>
</div>
<p>Number of ARGs reference genes in core SARG database (column ‘core SARG’) and updated SARG database using different cut off of identity (90%, 80% and 70%) for retrieving. A is the profile before using parallel classification to seat each sequence into hierarchical structure. B is the results of sequences amount after being classified into specific ARGs types and subtypes.</p>
<h6>ARGs-OAP概述</h6>
<p>ARGs-OAP是一个抗生素抗性基因分析平台、在线分析工具。</p>
<p>ARGs-OAP可以从宏基因组数据集中快速鉴定并定量分析抗生素抗性基因。<br />
ARGs-OAP中包含一个结构化ARG数据库SARG（type&#8211;subtype&#8211;reference sequence）。</p>
<p>ARGs-OAP 1.0版包括CARD及ARDB数据库的序列， 2.0版新纳入了NCBI-NR数据库中的ARG序列。</p>
<p>使用ARGs-OAP 进行注释后，对获得的ARGs：可以通过总reads数、16S rRNA基因拷贝数和细胞数量进行ARGs丰度标准化；2.0版优化了细胞数量定量分析过程。</p>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="480" data-height="360"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-6f172cd68aab2836.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/480/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-6f172cd68aab2836.png" data-original-width="480" data-original-height="360" data-original-format="image/png" data-original-filesize="42812" /></div>
</div>
<div class="image-caption"></div>
</div>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="760" data-height="414"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-0ecd4cb8d81295f7.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/760/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-0ecd4cb8d81295f7.png" data-original-width="760" data-original-height="414" data-original-format="image/png" data-original-filesize="158336" /></div>
</div>
<div class="image-caption"></div>
</div>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="762" data-height="390"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-b19ce4d819ea74fe.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/762/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-b19ce4d819ea74fe.png" data-original-width="762" data-original-height="390" data-original-format="image/png" data-original-filesize="124575" /></div>
</div>
<div class="image-caption"></div>
</div>
<h6>ARGs-OAP在线工具使用步骤</h6>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="757" data-height="385"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-4b3ab90b6e9c61c4.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/757/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-4b3ab90b6e9c61c4.png" data-original-width="757" data-original-height="385" data-original-format="image/png" data-original-filesize="119385" /></div>
</div>
<div class="image-caption"></div>
</div>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="806" data-height="504"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-26f5f9621a87cb80.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/806/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-26f5f9621a87cb80.png" data-original-width="806" data-original-height="504" data-original-format="image/png" data-original-filesize="117685" /></div>
</div>
<div class="image-caption"></div>
</div>
<p>1.本地计算机预先筛选潜在的ARG序列，以减少上传序列文件的大小；<br />
2.使用在线平台注释/分类ARG序列。</p>
<p>对于宏基因组数据，快速预筛选可去除总序 列&gt; 99.3％的不相关序列，显着减少上传文件的大小并加速在线BLASTX分析。</p>
<p>步骤2：上传预筛选后的ARG序列数据至online pipeline。<br />
ARGs_OAP_v2.0（步骤1）：<a href="https://github.com/biofuture/Ublastx_stageone" target="_blank" rel="nofollow">https://github.com/biofuture/Ublastx_stageone</a><br />
ARGs-OAP在线分析网站（步骤2）： <a href="http://smile.hku.hk/SARGs" target="_blank" rel="nofollow">http://smile.hku.hk/SARGs</a></p>
<p>The output files can be downloaded as tables listing the abundances of ARGs types/subtypes in different units:<br />
“ppm” (number of ARGs sequences in one million sequences) ;<br />
“copies of ARG per copy of 16S rRNA” ;<br />
“copies of ARG per prokaryote’s cell” .</p>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="759" data-height="517"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-45ecd61e1e620142.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/759/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-45ecd61e1e620142.png" data-original-width="759" data-original-height="517" data-original-format="image/png" data-original-filesize="205821" /></div>
</div>
<div class="image-caption"></div>
</div>
<p>当数据集包含新ARG时（即数据集2）：identity cutoff 设置为高于60％，则MCC值显著下降（图4a和4b），此水平下灵敏度也显著降低（图3d和4e），数据库的不完整性对注释精度影响不大（图4g和4h）。<br />
E-value 对这三个评估指标的影响：MCC值和精度随着E-value的减小而增加，但灵敏度没有太大变化。<br />
评估序列长度的影响：较长的读长导致较高的MCC和灵敏度（图3b和3c ）。<br />
最佳E-value 和 identity cutoff 值：与E值相比， identity值显示出更大的影响。蓝色箭头表示在以前ARGs注释（ E-value为1e-5， identity为90％）中对短读数宏基因组数据进行分析时，MCC值和灵敏度较低假阴性率很高，并且错过了许多ARG样序列。为了揭示更全面的ARG概况，基于使用模拟数据集2所示的MCC结果，如红色箭头所示，建议的最佳identity cutoff 为60％，E-value为1e-7。</p>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="733" data-height="322"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-6e336ceb68e7e97e.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/733/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-6e336ceb68e7e97e.png" data-original-width="733" data-original-height="322" data-original-format="image/png" data-original-filesize="41547" /></div>
</div>
<div class="image-caption"></div>
</div>
<p>序列覆盖度小于85％时，灵敏度和MCC值几乎没有影响。<br />
序列覆盖度从85％增加到100％时，灵敏度和MCC值急剧下降。<br />
更严格的序列覆盖度会错过更多类似ARG的序列。</p>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="1123" data-height="698"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-9d1ba3b494157f12.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/1000/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-9d1ba3b494157f12.png" data-original-width="1123" data-original-height="698" data-original-format="image/png" data-original-filesize="184851" /></div>
</div>
<div class="image-caption"></div>
</div>
<div class="image-package">
<div class="image-container">
<div class="image-container-fill"></div>
<div class="image-view" data-width="1200" data-height="790"><img class="" src="//upload-images.jianshu.io/upload_images/7600498-8f61f297ba929e5c.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/1000/format/webp" alt="" data-original-src="//upload-images.jianshu.io/upload_images/7600498-8f61f297ba929e5c.png" data-original-width="1200" data-original-height="790" data-original-format="image/png" data-original-filesize="369366" /></div>
</div>
<div class="image-caption"></div>
</div>
<p>参考文献：<br />
Yang Y, Jiang X, Chai B, et al. ARGs-OAP: online analysis pipeline for antibiotic resistance genes detection from metagenomic data using an integrated structured ARG-database[J]. Bioinformatics, 2016, 32(15):2346.<br />
Yin X, Jiang X T, Chai B, et al. ARGs-OAP v2.0 with an Expanded SARG Database and Hidden Markov Models for Enhancement Characterization and Quantification of Antibiotic Resistance Genes in Environmental Metagenomes[J]. Bioinformatics, 2018</p>
</div>
<p>作者：周运来就是我<br />
链接：https://www.jianshu.com/p/feb181e7888e<br />
來源：简书<br />
简书著作权归作者所有，任何形式的转载都请联系作者获得授权并注明出处。</p></div>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1348</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>pronunciation</title>
		<link>http://www.biofacebook.com/?p=1345</link>
		<comments>http://www.biofacebook.com/?p=1345#comments</comments>
		<pubDate>Mon, 15 Oct 2018 09:00:46 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[兴趣杂项]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1345</guid>
		<description><![CDATA[<p>发音第一个 n 带鼻音</p> <p>analysis (noun single )</p> <p>analyses (noun) (multiplex number of analysis)</p> <p>动词</p> <p>analyse a&#8217; n(ə)lʌɪz 发音a 不带鼻音</p> <p>th 在尾巴 发 f</p> <p>thanks</p> <p>中间ch发k的时候要弱化发音</p> ]]></description>
				<content:encoded><![CDATA[<p>发音第一个 n 带鼻音</p>
<p>analysis    (noun single  )</p>
<p>analyses (noun)  (multiplex number of analysis)</p>
<p>动词</p>
<p>analyse  a&#8217; n(ə)lʌɪz          发音a 不带鼻音</p>
<p>th 在尾巴 发 f</p>
<p>thanks</p>
<p>中间ch发k的时候要弱化发音</p>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1345</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>VOSEQ server start</title>
		<link>http://www.biofacebook.com/?p=1314</link>
		<comments>http://www.biofacebook.com/?p=1314#comments</comments>
		<pubDate>Thu, 12 Jul 2018 03:24:28 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[兴趣杂项]]></category>
		<category><![CDATA[服务器管理]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1314</guid>
		<description><![CDATA[ #### <p>[shenzy@LFE0530 VoSeq-2.1.1]$ source /usr/bin/virtualenvwrapper.sh [shenzy@LFE0530 VoSeq-2.1.1]$ workon voseq_environment (voseq_environment) [shenzy@LFE0530 VoSeq-2.1.1]$</p> <p>python voseq/manage.py runserver &#8211;settings=voseq.settings.local 143.89.29.80:8000</p> <p>&#160;</p> setup environmental variables, virtual environments Wai-Yin Kwan edited this page on Jul 5, 2015 · 12 revisions Pages 7 Home Common git commands Computer setup Github workflow Helpful links misc commands setup environmental variables, virtual environments [...]]]></description>
				<content:encoded><![CDATA[<div class="gh-header">
<div class="gh-header-show">
<h1 class="gh-header-title instapaper_title">####</h1>
<p>[shenzy@LFE0530 VoSeq-2.1.1]$ source /usr/bin/virtualenvwrapper.sh<br />
[shenzy@LFE0530 VoSeq-2.1.1]$ workon voseq_environment<br />
(voseq_environment) [shenzy@LFE0530 VoSeq-2.1.1]$</p>
<p>python voseq/manage.py runserver &#8211;settings=voseq.settings.local 143.89.29.80:8000</p>
<p>&nbsp;</p>
<h1 class="gh-header-title instapaper_title">setup environmental variables, virtual environments</h1>
<div class="gh-header-meta">Wai-Yin Kwan edited this page on Jul 5, 2015 · <a class="history" href="https://github.com/LearnTeachCode/marsrocks/wiki/setup-environmental-variables,--virtual-environments/_history">12 revisions</a></div>
</div>
</div>
<div id="wiki-content" class="wiki-content">
<div class="has-rightbar">
<div id="wiki-rightbar" class="wiki-rightbar">
<div class="wiki-pages-box readability-sidebar boxed-group flush js-wiki-pages-box ">
<h3 class="js-wiki-toggle-collapse wiki-auxiliary-content">Pages <span class="Counter">7</span></h3>
<div class="boxed-group-inner wiki-auxiliary-content wiki-auxiliary-content-no-bg">
<ul class="wiki-pages" data-filterable-for="wiki-pages-filter" data-filterable-type="substring">
<li><strong><a class="wiki-page-link" href="https://github.com/LearnTeachCode/marsrocks/wiki">Home</a></strong></li>
<li><strong><a class="wiki-page-link" href="https://github.com/LearnTeachCode/marsrocks/wiki/Common-git-commands">Common git commands</a></strong></li>
<li><strong><a class="wiki-page-link" href="https://github.com/LearnTeachCode/marsrocks/wiki/Computer-setup">Computer setup</a></strong></li>
<li><strong><a class="wiki-page-link" href="https://github.com/LearnTeachCode/marsrocks/wiki/Github-workflow">Github workflow</a></strong></li>
<li><strong><a class="wiki-page-link" href="https://github.com/LearnTeachCode/marsrocks/wiki/Helpful-links">Helpful links</a></strong></li>
<li><strong><a class="wiki-page-link" href="https://github.com/LearnTeachCode/marsrocks/wiki/misc-commands">misc commands</a></strong></li>
<li><strong><a class="wiki-page-link" href="https://github.com/LearnTeachCode/marsrocks/wiki/setup-environmental-variables,--virtual-environments">setup environmental variables, virtual environments</a></strong></li>
</ul>
</div>
</div>
<h5 class="mt-0 mb-2">Clone this wiki locally</h5>
<div class="width-full input-group"><input id="wiki-clone-url" class="form-control input-sm text-small text-gray input-monospace js-url-field" readonly="readonly" type="text" value="https://github.com/LearnTeachCode/marsrocks.wiki.git" /></div>
</div>
<div id="wiki-body" class="wiki-body gollum-markdown-content instapaper_body">
<div class="markdown-body">
<p>When developing and deploying a web app, different environments (local machine, live site) need different configurations (passwords, database names, etc). We can use environmental variables to setup the different environments. Here are two options to set up environmental variables: easy way with autoenv or harder way with virtualenvwrapper.</p>
<p>You need to use a text editor to edit these files. In this demo, I&#8217;m using atom, but you can use any text editor.</p>
<h2><a id="user-content-easy-way-using-autoenv" class="anchor" href="https://github.com/LearnTeachCode/marsrocks/wiki/setup-environmental-variables,--virtual-environments#easy-way-using-autoenv"></a>Easy way using autoenv</h2>
<h4><a id="user-content-part-1-only-do-this-once" class="anchor" href="https://github.com/LearnTeachCode/marsrocks/wiki/setup-environmental-variables,--virtual-environments#part-1-only-do-this-once"></a>Part 1. Only do this once.</h4>
<ol>
<li>Download <a href="https://github.com/kennethreitz/autoenv">autoenv</a>.</li>
</ol>
<pre><code>$ git clone git://github.com/kennethreitz/autoenv.git ~/.autoenv

$ echo 'source ~/.autoenv/activate.sh' &gt;&gt; ~/.bashrc
</code></pre>
<h4><a id="user-content-part-2-do-this-for-every-project" class="anchor" href="https://github.com/LearnTeachCode/marsrocks/wiki/setup-environmental-variables,--virtual-environments#part-2-do-this-for-every-project"></a>Part 2. Do this for every project.</h4>
<ol>
<li>create .env file in the root directory of the project.</li>
</ol>
<pre><code>$ cd &lt;path/to/project&gt;

$ touch .env
</code></pre>
<p>open the .env file.</p>
<pre><code>$ atom .env
</code></pre>
<p>Put environmental variables into the .env file, then save the file.</p>
<pre><code>export VARIABLE_NAME="value"

</code></pre>
<ol start="2">
<li>reload shell</li>
</ol>
<pre><code> source ~/.bashrc
</code></pre>
<ol start="3">
<li>type <code>y</code> when you will see a message like:</li>
</ol>
<pre><code>autoenv: This is the first time you are about to source /path/.env:

autoenv: Are you sure you want to allow this? (y/N)</code></pre>
</div>
</div>
</div>
</div>
<p>##########################################</p>
<pre><code class="language-text" data-lang="text">sudo pip install virtualenvwrapper
</code></pre>
<pre><code class="language-text" data-lang="text">export WORKON_HOME=~/venvs</code></pre>
<pre><code class="language-text" data-lang="text">source /usr/bin/virtualenvwrapper.sh
</code></pre>
<pre>mkvirtualenv -p /usr/bin/python3 voseq_environment

[shenzy@LFE0530 VoSeq]$ mkvirtualenv -p /usr/bin/python3 voseq_environment
Running virtualenv with interpreter /usr/bin/python3
Using base prefix '/usr'
New python executable in /home/shenzy/envs/voseq_environment/bin/python3
Also creating executable in /home/shenzy/envs/voseq_environment/bin/python
Installing setuptools, pip, wheel...done.
virtualenvwrapper.user_scripts creating /home/shenzy/envs/voseq_environment/bin/predeactivate
virtualenvwrapper.user_scripts creating /home/shenzy/envs/voseq_environment/bin/postdeactivate
virtualenvwrapper.user_scripts creating /home/shenzy/envs/voseq_environment/bin/preactivate
virtualenvwrapper.user_scripts creating /home/shenzy/envs/voseq_environment/bin/postactivate
virtualenvwrapper.user_scripts creating /home/shenzy/envs/voseq_environment/bin/get_env_details
(voseq_environment) [shenzy@LFE0530 VoSeq]$ workon voseq_environment</pre>
<p><strong>pip install django</strong></p>
<pre>pip install django-suit</pre>
<pre class="default prettyprint prettyprinted"><code><span class="pln">pip install </span><span class="pun">-</span><span class="pln">r requirements</span><span class="pun">.</span><span class="pln">txt</span></code></pre>
<pre><code class="language-text" data-lang="text">#########################33</code></pre>
<p>cd /home/shenzy/software/VoSeq</p>
<p>workon voseq_environment</p>
<p>source ~/.autoenv/activate.sh      #  使得能够调用 blastn 等 通过  export path</p>
<p><strong>pip install django</strong></p>
<p>pip install django-suit</p>
<pre class="default prettyprint prettyprinted"><code><span class="pln">pip install </span><span class="pun">-</span><span class="pln">r requirements</span><span class="pun">.</span><span class="pln">txt</span></code></pre>
<p>make serve</p>
<p>(voseq_environment) [shenzy@LFE0530 VoSeq]$ make serve<br />
python voseq/manage.py create_stats &#8211;settings=voseq.settings.local<br />
python voseq/manage.py runserver &#8211;settings=voseq.settings.local<br />
Performing system checks&#8230;</p>
<p>System check identified no issues (0 silenced).<br />
July 16, 2018 &#8211; 06:10:13<br />
Django version 1.10.4, using settings &#8216;voseq.settings.local&#8217;<br />
Starting development server at http://127.0.0.1:8000/<br />
Quit the server with CONTROL-C.</p>
<p>&nbsp;</p>
<p>(voseq_environment) [shenzy@LFE0530 VoSeq]$ make serve<br />
python voseq/manage.py create_stats &#8211;settings=voseq.settings.local<br />
python voseq/manage.py runserver &#8211;settings=voseq.settings.local<br />
Performing system checks&#8230;</p>
<p>System check identified no issues (0 silenced).<br />
July 16, 2018 &#8211; 07:02:22<br />
Django version 1.10.4, using settings &#8216;voseq.settings.local&#8217;<br />
Starting development server at http://127.0.0.1:8000/<br />
Quit the server with CONTROL-C.<br />
^C<br />
(voseq_environment) [shenzy@LFE0530 VoSeq]$ vi setup.py<br />
(voseq_environment) [shenzy@LFE0530 VoSeq]$ vi runserver.py<br />
(voseq_environment) [shenzy@LFE0530 VoSeq]$ python voseq/manage.py runserver &#8211;settings=voseq.settings.local 143.89.29.80:8000<br />
Performing system checks&#8230;</p>
<p>System check identified no issues (0 silenced).<br />
July 16, 2018 &#8211; 07:03:30<br />
Django version 1.10.4, using settings &#8216;voseq.settings.local&#8217;<br />
Starting development server at http://143.89.29.80:8000/</p>
<p>&nbsp;</p>
<p>########################</p>
<p>[shenzy@LFE0530 VoSeq]$ sudo -u postgres -i<br />
[postgres@LFE0530 ~]$ psql<br />
psql (9.2.23)<br />
Type &#8220;help&#8221; for help.</p>
<p>postgres=#</p>
<p>sudo -u postgres -i</p>
<p>&nbsp;</p>
<div class="votecell post-layout--left">
<div class="vote">
<h1 class="grid--cell fs-headline1 fl1"><a class="question-hyperlink" href="https://stackoverflow.com/questions/26846093/postgresql-9-3-on-centos-7-with-custom-pgdata">Postgresql 9.3 on Centos 7 with custom PGDATA</a></h1>
</div>
</div>
<div class="answercell post-layout--right">
<div class="post-text">
<p>920 sudo yum install postgresql postgresql-contrib postgresql-server-dev-9.3<br />
922 sudo yum uninstall postgresql<br />
923 sudo yum remove postgresql<br />
924 sudo yum install postgresql<br />
925 sudo yum install postgresql*<br />
927 sudo yum reinstall postgresql*<br />
931 sudo -u postgres createuser owning_user<br />
932 sudo -u postgres createuser shenzy<br />
933 sudo -u postgres createuser postgres<br />
935 sudo -u postgres -i<br />
945 sudo -u postgres -i</p>
<p>&nbsp;</p>
<p>try this:</p>
<pre class="lang-sql prettyprint prettyprinted"><code> <span class="pun">##</span><span class="pln"> Login </span><span class="kwd">with</span><span class="pln"> postgres </span><span class="kwd">user 
</span><span class="pln">sudo -u postgres -i</span><span class="pln">
 export PGDATA</span><span class="pun">=/</span><span class="pln">your_path</span><span class="pun">/</span><span class="pln">data
 pg_ctl </span><span class="pun">-</span><span class="pln">D </span><span class="pun">$</span><span class="pln">PGDATA </span><span class="kwd">start</span> <span class="pun">&amp;</span></code></pre>
</div>
</div>
<p>&nbsp;</p>
<pre class="lang-sql prettyprint prettyprinted"><code><span class="pln">service postgresql </span><span class="kwd">start/status</span></code></pre>
<p>SHOW data_directory;</p>
<p>postgres=# \q<br />
[postgres@LFE0530 ~]$ pwd<br />
/var/lib/pgsql<br />
[postgres@LFE0530 ~]$ psql<br />
psql (9.2.23)<br />
Type &#8220;help&#8221; for help.</p>
<p>postgres=# SHOW data_directory;<br />
data_directory<br />
&#8212;&#8212;&#8212;&#8212;&#8212;-<br />
/home/pgsql<br />
(1 row)</p>
<p>&nbsp;</p>
<p>虽然还是显示下面错误，但其实上面的已经可以让 pgsql启动，并服务！！！</p>
<p>Redirecting to /bin/systemctl status postgresql.service<br />
● postgresql.service &#8211; PostgreSQL database server<br />
Loaded: loaded (/usr/lib/systemd/system/postgresql.service; disabled; vendor preset: disabled)<br />
Active: failed (Result: exit-code) since 三 2018-07-18 22:31:30 EDT; 5min ago<br />
Process: 8035 ExecStart=/usr/bin/pg_ctl start -D ${PGDATA} -s -o -p ${PGPORT} -w -t 300 (code=exited, status=1/FAILURE)<br />
Process: 8029 ExecStartPre=/usr/bin/postgresql-check-db-dir ${PGDATA} (code=exited, status=0/SUCCESS)</p>
<p>7月 18 22:31:29 LFE0530 pg_ctl[8035]: HINT: Is another postmaster already running on port 5432? If not, wait a few seconds and retry.<br />
7月 18 22:31:29 LFE0530 pg_ctl[8035]: LOG: could not bind IPv4 socket: Address already in use<br />
7月 18 22:31:29 LFE0530 pg_ctl[8035]: HINT: Is another postmaster already running on port 5432? If not, wait a few seconds and retry.<br />
7月 18 22:31:29 LFE0530 pg_ctl[8035]: WARNING: could not create listen socket for &#8220;localhost&#8221;<br />
7月 18 22:31:29 LFE0530 pg_ctl[8035]: FATAL: could not create any TCP/IP sockets<br />
7月 18 22:31:30 LFE0530 pg_ctl[8035]: pg_ctl: could not start server<br />
7月 18 22:31:30 LFE0530 systemd[1]: postgresql.service: control process exited, code=exited status=1<br />
7月 18 22:31:30 LFE0530 systemd[1]: Failed to start PostgreSQL database server.<br />
7月 18 22:31:30 LFE0530 systemd[1]: Unit postgresql.service entered failed state.<br />
7月 18 22:31:30 LFE0530 systemd[1]: postgresql.service failed.<br />
[root@LFE0530 data]#</p>
<p>&nbsp;</p>
<p>voseq=# \COPY public_interface_vouchers(code,notes) FROM &#8216;/home/shenzy/software/VoSeq/U-RVDBv13.0.voucher_10_5000top.csv&#8217; DELIMITER &#8216;,&#8217; CSV HEADER;<br />
voseq=#</p>
<p>voseq=#<br />
\COPY public_interface_sequences FROM &#8216;/home/shenzy/software/VoSeq/all.gene_fasta_10test_import.csv&#8217; DELIMITER &#8216;,&#8217; CSV;<br />
voseq=#<br />
\COPY public_interface_sequences FROM &#8216;/home/shenzy/software/VoSeq/all.gene_fasta_10test_import.csv&#8217; DELIMITER &#8216;,&#8217; CSV;<br />
voseq=#<br />
\COPY public_interface_vouchers FROM &#8216;/home/shenzy/software/VoSeq/U-RVDBv13.0.voucher_import.csv&#8217; DELIMITER &#8216;,&#8217; CSV;</p>
<p>&nbsp;</p>
<p>###########################<br />
export</p>
<p>voseq=# COPY public_interface_sequences TO &#8216;/home/shenzy/software/VoSeq/testseq.csv&#8217; WITH CSV;</p>
<p>&nbsp;</p>
<p>############33</p>
<p>#empty table<br />
voseq=# truncate table public_interface_sequences  CASCADE;<br />
NOTICE:  truncate cascades to table &#8220;public_interface_primers&#8221;<br />
TRUNCATE TABLE<br />
voseq=# truncate table public_interface_vouchers CASCADE;<br />
NOTICE:  truncate cascades to table &#8220;public_interface_flickrimages&#8221;<br />
NOTICE:  truncate cascades to table &#8220;public_interface_localimages&#8221;<br />
NOTICE:  truncate cascades to table &#8220;public_interface_sequences&#8221;<br />
NOTICE:  truncate cascades to table &#8220;public_interface_primers&#8221;</p>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1314</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Sci-Hub 最新可用网址</title>
		<link>http://www.biofacebook.com/?p=1311</link>
		<comments>http://www.biofacebook.com/?p=1311#comments</comments>
		<pubDate>Mon, 15 Jan 2018 13:27:01 +0000</pubDate>
		<dc:creator><![CDATA[szypanther]]></dc:creator>
				<category><![CDATA[兴趣杂项]]></category>

		<guid isPermaLink="false">http://www.biofacebook.com/?p=1311</guid>
		<description><![CDATA[<p>Pubmedy：显示影响因子+引用数、Sci-hub全文下载的浏览器扩展</p> <p>很多人都知道一个下载文献的好地方 Sci-Hub，但是现在 Sci-Hub 的站长 Alexandra Elbakyan 女士面临着爱思唯尔的诉讼，让其面临数百万美元的赔偿。因此现在 Sci-Hub 在某些时候可能不可用了。不过 Alexandra Elbakyan 女士并未妥协，希望大家能多支持一下这位女士（貌似只能接受比特币捐款，我也不知道怎么弄，要是paypal或者payoneer就好了）。</p> <p>相信很多人都知道下面这两个网址：</p> <p>Sci-Hub 官网：http://www.sci-hub.org/</p> <p>Sci-Hub 替代网址：http://www.sci-hub.io/</p> <p>http://sci-hub.ac/和http://sci-hub.cc/也都已经阵亡（2017月12月6日测试）</p> <p>以上几个网址都挂了，现在可用的是</p> <p>http://sci-hub.bz/（12月8日测试 已经不可用了）</p> <p>http://sci-hub.la/</p> <p>http://sci-hub.hk/</p> <p>http://sci-hub.tw/</p> <p>http://sci-hub.tv/</p> <p>http://80.82.77.83/</p> <p>http://80.82.77.84/</p> ]]></description>
				<content:encoded><![CDATA[<p><a href="http://blog.sciencenet.cn/blog-3196388-1049657.html" target="_blank">Pubmedy：显示影响因子+引用数、Sci-hub全文下载的浏览器扩展</a></p>
<p>很多人都知道一个下载文献的好地方 Sci-Hub，但是现在 Sci-Hub 的站长 Alexandra Elbakyan 女士面临着爱思唯尔的诉讼，让其面临数百万美元的赔偿。因此现在 Sci-Hub 在某些时候可能不可用了。不过 Alexandra Elbakyan 女士并未妥协，希望大家能多支持一下这位女士（貌似只能接受比特币捐款，我也不知道怎么弄，要是paypal或者payoneer就好了）。</p>
<p>相信很多人都知道下面这两个网址：</p>
<p>Sci-Hub 官网：http://www.sci-hub.org/</p>
<p>Sci-Hub 替代网址：http://www.sci-hub.io/</p>
<p>http://sci-hub.ac/和http://sci-hub.cc/也都已经阵亡（2017月12月6日测试）</p>
<p>以上几个网址都挂了，现在可用的是</p>
<p><strong>http://sci-hub.bz/（12月8日测试 已经不可用了）</strong></p>
<p><strong>http://sci-hub.la/</strong></p>
<p><strong>http://sci-hub.hk/</strong></p>
<p><strong>http://sci-hub.tw/</strong></p>
<p><strong>http://sci-hub.tv/</strong></p>
<p><strong>http://80.82.77.83/</strong></p>
<p><strong>http://80.82.77.84/</strong></p>
]]></content:encoded>
			<wfw:commentRss>http://www.biofacebook.com/?feed=rss2&#038;p=1311</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
	</channel>
</rss>
