search

rosemckeon / honours-project

When rare gametes meet: an individual-based model of polyploid evolution in angiosperms.
https://github.com/rosemckeon/honours-project/blob/master/thesis/_thesis_2417024.pdf
0 stars 0 forks source link

ploidy_rate not being saved in count data #22

Closed rosemckeon closed 4 years ago

rosemckeon commented 4 years ago

@bradduthie I think I've made an error. I haven't got the ploidy_rate value being output with the count data so this is going to be annoying to add in and plot. It's not in the function call I've recorded because all that logs is that I've passed another object to the function and not actually the value that was in the object. The only place it's kept is in the log file which is in the same folder as each of the count data rds files; in the first line: # SIMULATION 1: Ploidy rate = 0.00761167602969435. I'm going to need some regex to strip that out aren't I? Or update the output and run all these sims I've collected over the weekend again.

Do you know off the top of your head how to get that number out of a txt file?

bradduthie commented 4 years ago

Hi @rosemckeon -- yeah, I think you'll need to get some code to strip that information out, but it might just be easier to re-run the simulations after making the necessary adjustments. I'm waiting on a new computer to arrive (I hope today), so I might be able to donate 20+ cores before the end of the week. I could write a script to get the number out of a text file in R (I've done this with CSVs before from tablet outputs). The R commands aren't too difficult -- happy to send them if you think it would help, but it might create more work than just re-simulating.