Hi,

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 ="

Thanks!

Thanks. Any suggestions on how to correct the mistake?

I already gave you a suggestion. Read An Introduction to R. In addition, I gave you like 75% of the answer. You need to re-read what I said, and if that doesn't make sense, then you really need to read An Introduction to R. You are getting tripped up by one of the most basic operations in R, and if you can't do the basics, you probably shouldn't be doing the more advanced stuff yet.

Thanks for you help.