<div>Dear Dr. L<span class="Apple-style-span" style="font-family: arial, sans-serif, &#39;Arial Unicode MS&#39;; font-size: 13px; border-collapse: collapse; white-space: pre-wrap; -webkit-border-horizontal-spacing: 2px; -webkit-border-vertical-spacing: 2px; ">obry and <span class="Apple-style-span" style="font-family: arial, helvetica, clean, sans-serif; white-space: normal; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; line-height: 15px; ">Colleagues,</span></span></div>
<div><font class="Apple-style-span" face="arial, helvetica, clean, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse; line-height: 15px;"><br></span></font></div><div><font class="Apple-style-span" face="arial, helvetica, clean, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse; line-height: 15px;">The code is running fine and I can reproduce the output as same as yours. It seems that the data is accessible. </span></font></div>
<div><font class="Apple-style-span" face="arial, helvetica, clean, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse; line-height: 15px;"><br></span></font></div><div><font class="Apple-style-span" face="arial, helvetica, clean, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse; line-height: 15px;">Yes, these are absolutely those genomes that I intend to work with. I want <span class="Apple-style-span" style="border-collapse: separate; font-family: arial; line-height: normal; ">to consider only all the CDS (or Open Reading Frame) for the statistical observation.</span></span></font></div>
<div><br></div><div>Thank you very much</div><div><br></div><div>Sincerely yours</div><div>Sourav Roy Choudhury</div><div><br></div><div><span class="Apple-style-span" style="font-family: arial, helvetica, clean, sans-serif; font-size: 13px; border-collapse: collapse; color: rgb(136, 136, 136); line-height: 15px; ">SOURAV ROY CHOUDHURY.<br style="line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">
Dept of Neuroscience<br style="line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">University of Calcutta <br style="line-height: 1.2em; outline-style: none; outline-width: initial; outline-color: initial; ">
West Bengal, INDIA</span></div><div><br></div><br><br><div><br></div><div><br><div class="gmail_quote">On Wed, Nov 4, 2009 at 12:12 AM, Jean lobry <span dir="ltr">&lt;<a href="mailto:lobry@biomserv.univ-lyon1.fr">lobry@biomserv.univ-lyon1.fr</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Dear Sourav Roy Choudhury,<br>
<br>
The first thing to do is to check that you have access to data.<br>
Complete virus genomes are in the refseqVirus bank. Try this:<br>
<br>
###<br>
 library(seqinr)<br>
 choosebank(&quot;refseqViruses&quot;)<br>
 query(&quot;herp&quot;, &quot;TID=548681&quot;)<br>
 mnemo &lt;- getName(herp)<br>
 mnemo<br>
 [1] &quot;NC_000898&quot; &quot;NC_001347&quot; &quot;NC_001348&quot; &quot;NC_001350&quot; &quot;NC_001491&quot;<br>
 [6] &quot;NC_001493&quot; &quot;NC_001650&quot; &quot;NC_001664&quot; &quot;NC_001716&quot; &quot;NC_001798&quot;<br>
[11] &quot;NC_001806&quot; &quot;NC_001826&quot; &quot;NC_001844&quot; &quot;NC_001847&quot; &quot;NC_001987&quot;<br>
[16] &quot;NC_002229&quot; &quot;NC_002512&quot; &quot;NC_002531&quot; &quot;NC_002577&quot; &quot;NC_002641&quot;<br>
[21] &quot;NC_002665&quot; &quot;NC_002686&quot; &quot;NC_002794&quot; &quot;NC_003401&quot; &quot;NC_003409&quot;<br>
[26] &quot;NC_003521&quot; &quot;NC_004065&quot; &quot;NC_004367&quot; &quot;NC_004812&quot; &quot;NC_005261&quot;<br>
[31] &quot;NC_005264&quot; &quot;NC_005881&quot; &quot;NC_006146&quot; &quot;NC_006150&quot; &quot;NC_006151&quot;<br>
[36] &quot;NC_006273&quot; &quot;NC_006560&quot; &quot;NC_006623&quot; &quot;NC_007016&quot; &quot;NC_007605&quot;<br>
[41] &quot;NC_007646&quot; &quot;NC_007653&quot; &quot;NC_008210&quot; &quot;NC_008211&quot; &quot;NC_009127&quot;<br>
[46] &quot;NC_009333&quot; &quot;NC_009334&quot; &quot;NC_011587&quot; &quot;NC_011644&quot;<br>
###<br>
<br>
Can You reproduce this?<br>
<br>
Are those the genomes You want to work with?<br>
<br>
Do You want to consider all CDS as a statististical observation<br>
or do You want to consider all genomes as a statistical observation?<br>
<br>
Best,<br><font color="#888888">
<br>
-- <br>
Jean R. Lobry            (<a href="mailto:lobry@biomserv.univ-lyon1.fr" target="_blank">lobry@biomserv.univ-lyon1.fr</a>)<br>
Laboratoire BBE-CNRS-UMR-5558, Univ. C. Bernard - LYON I,<br>
43 Bd 11/11/1918, F-69622 VILLEURBANNE CEDEX, FRANCE<br>
allo  : +33 472 43 27 56     fax    : +33 472 43 13 88<br>
<a href="http://pbil.univ-lyon1.fr/members/lobry/" target="_blank">http://pbil.univ-lyon1.fr/members/lobry/</a><br>
<br>
<br>
</font></blockquote></div><br></div>