read.ms.output()

[suiyuano]

when using gap to read an input file of msms, it reports that:

"Error in gap::read.ms.output(hap_file, verbose = verbose) :
不是所有的length(marker) == ndraws都是TRUE"

could u plz tell me how to solve the issue?

[jinghuazhao]

Thanks for letting me know .

I haven't used it for a long time and just experimented with it. It worked as follows,

1. Compiled from inst/ms/ with bash clms which gives ms.
2. Move ms into a folder/directory that is in the search path, in my case ~/bin
3. Start R, load gap and run the documentation example,

    system("ms 5 4 -s 5 > ms.out")
    msout1 <- read.ms.output("ms.out")
   
    system("ms 50 4 -s 5 > ms.out")
    msout2 <- read.ms.output("ms.out",outfile="out",outfileonly=TRUE)
   
    msout <- system("ms 5 4 -s 5 -L", intern=TRUE)
    msout3 <- read.ms.output(msout,FALSE)

Hope this helps

[suiyuano]

Thank you very much for your reply, in fact, I checked the code you published for the "read.ms.output()" section and I found that the problem occurred when read.ms.output() processed the headers in the ms output, because I was using msms. I have solved the problem. Thank you very much again! mango ***@***.***  

…

------------------&nbsp;原始邮件&nbsp;------------------ 发件人: "jinghuazhao/R" ***@***.***&gt;; 发送时间:&nbsp;2024年6月24日(星期一) 晚上6:38 ***@***.***&gt;; ***@***.******@***.***&gt;; 主题:&nbsp;Re: [jinghuazhao/R] read.ms.output() (Issue #7) Thanks for letting me know . I haven't used it for a long time and just experimented with it. It worked as follows, I compiled ms from inst/ms/ with bash clms which gives ms. Move ms into a folder/directory that is in the search path, in my case ~/bin Start R, load gap and run the documentation example, system("ms 5 4 -s 5 &gt; ms.out") msout1 <- read.ms.output("ms.out") system("ms 50 4 -s 5 &gt; ms.out") msout2 <- read.ms.output("ms.out",outfile="out",outfileonly=TRUE) msout <- system("ms 5 4 -s 5 -L", intern=TRUE) msout3 <- read.ms.output(msout,FALSE) Hope this helps — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you authored the thread.Message ID: ***@***.***&gt;

[jinghuazhao]

Thanks for the feedback. In fact read.ms.output() was based author of ms himself and it would indeed be nice to mirror a read.msms.output().