I am trying to write a table of posterior probabilities from a Centipede function and having some trouble. Basically want to write.table of the results of the following command:
sum(fit$PostPr == 1)
Which gives me a certain number of sites with PostPr of 1, and I would like to write those specific sites to a TAB delimited file. Any ideas on how to go about doing this? I have tried the following:
write.table(Sox3_BTTP_peripheral_fit$PostPr = 1, "Sox3_BTTP_peripheral_fit", sep="\t")
but getting this error: Error: unexpected '=' in "write.table(Sox3_BTTP_peripheral_fit$PostPr ="