PTT
Submit
Submit
選擇語言
正體中文
简体中文
PTT
R_Language
[問題] bioinformatics, XStringset格式
作者:
cb1040
(Absurde.Y)
2015-11-09 17:12:54
[程式諮詢] Bioinformatics的XStringset格式轉換
[入門]
2010修課玩過VBA,不過之後就沒了直到今日,也算新手吧
此版首po,請前輩們指教
[問題敘述]:
小弟我剛踏入bioinformatics,用Bioinfomatics with R cookbook 自學中
在MSA, Multiple Sequence Alignment部分,
使用下面這個package時卡關了
muscle, (multiple sequence comparison by log-expectation)
我的問題是:
如何轉換XStringset格式,且是多筆資料轉換
(有查詢Biostrings 中AAString功能,但仍不知道如何使用)
[程式範例]
這是原書步驟
http://imgur.com/eFWyDjk
http://imgur.com/94zPqQQ
這是會用到的fasta檔案
http://tinyurl.com/q6kjpru
以下是我的內容
pastie版 (第一次使用,還不熟,請見諒)
http://pastie.org/10540288
自己截圖版
http://imgur.com/ecLpoF1
http://imgur.com/4BqdEng
http://imgur.com/rT9PEqr (後面是其他物種的序列,共十種,我就不貼了)
http://imgur.com/0WvRNXJ
因為教材有點年份,code有些出入
例如install.packages("muscle")會找不到東西
或是沒有muscle::read.fasta這個功能
這裏再說一次我的問題
有查詢Biostrings 中 XStringset、AAString功能,但仍不會使用
想請問: 如何將多筆資料轉換成XStringset格式
[環境敘述]:
直接截圖
http://imgur.com/u42f95X
[關鍵字]: R, muscle, Biostrings, multiple sequence alignment, XStringset
作者:
Godkin
(山裡的人)
2015-11-09 23:57:00
用readDNAStringSet("fastaMSA.fasta", format="fasta")讀這個package請用bioconductor上的新版本多筆資料可以直接通通塞到同一個fasta檔裡頭跑
作者:
cb1040
(Absurde.Y)
2015-11-10 10:25:00
我用source(
http://bioconductor.org....)
然後biocLite ("muscle") 安裝,這樣應該是最新版對吧我用readXStringSet成功了!!謝謝!!!關於第三點,是我必須額外存一個新的fasta檔案對嗎
作者:
Godkin
(山裡的人)
2015-11-10 17:16:00
對, 可以用cat去串聯每個fasta檔
繼續閱讀
[分享] Rcpp 中如何儲存和讀取C++物件
Wush978
[情報] GSOC 2016
Wush978
Re: [問題] 比對兩個list
CKPILOT
[問題] postscript 產生的 eps 圖檔
Chris7462
Re: [問題] 比對兩個list
celestialgod
[問題] 比對兩個list
angela79979
Re: [問題] 依照條件,計算欄位,合併欄位
celestialgod
[問題] 依照條件,計算欄位,合併欄位
JJabrams
Re: [問題] 質因數分解
showfeb
Re: [問題] 資料讀入
ror
Links
booklink
Contact Us: admin [ a t ] ucptt.com