Hi all, first time posting here. I am working on a chapter of my master’s thesis, and am using data from an experiment done well before I arrived to the school. I am pretty much the only one in my lab, including my PI and other faculty, who knows how to use R. Anyway, my PI wants to get narrow-sense heritability data for two traits in a data set that includes:
Dam family, Sire family, Dam x Sire genotype (based on family information), Replicate # (Offspring of cross x), Offspring trait 1, and Offspring trait 2.
I have tried writing some code block to come up with this, but it seems to be difficult to get “meaningful” or “realistic” outputs (my PI’s words) from this dataset. The data is admittedly not great, as we do not have a factorial breeding setup. Essentially, every dam and every sire are not all mated to each other, so there are many gaps in the information.
Just wanted to know if this is possible, and if so what the best way to go about it would be. I feel like I’m banging my head against the wall, so any help is appreciated! If anyone wants more info to be able to help, please let me know
I have no idea what you are doing—my last bio course was back when I was ~ 16—but can you give us a pointer, to what you are doing maybe a reference to something similar?
It would help a great deal for us to see what code you have written so far and to have a sample of the data you are working with if it is not confidential.
Copy the code and paste it here between
```
````
This gives us formatted code that we can copy, paste and run .
A handy way to supply data is to use the dput() function. Do dput(mydata) where "mydata" is the name of your dataset. For really large datasets probably dput(head(mydata, 100) will do. Paste it here between
