- lenght of roh = lm XXXXX
- heterozygosity & inbreeding = lm   XXXX
- survival = glm XXXX
- fped = cor + t-est
- histogram för FROH måtten, heterozygosity, FPED, lenght of ROH

####Testing for difference in lenght of ROH, heterozygosity and FROH between natives, immigrant F1 and immigrant F2####

fm <- lm(ROH.lenght ~ ancestry, data=inbreeding)
summary(fm)

fm <- lm(het.kb ~ ancestry, data=inbreeding)
summary(fm)

fm <- lm(FROH.100.kb ~ ancestry, data=inbreeding)
summary(fm)

fm <- lm(FROH.2.mb ~ ancestry, data=inbreeding)
summary(fm)

fm <- lm(FROH.8.mb ~ ancestry, data=inbreeding)
summary(fm)




####Testing if heterozygosity, FROH and FPED influence juvenile survival####

fm <- glm(survival ~ het.kb, family = binomial(link = "logit"), data = inbreeding)
summary(fm)

fm <- glm(survival ~ FROH.100.kb, family = binomial(link = "logit"), data = inbreeding)
summary(fm)

fm <- glm(survival ~ FROH.2.mb, family = binomial(link = "logit"), data = inbreeding)
summary(fm)

fm <- glm(survival ~ FROH.8.mb, family = binomial(link = "logit"), data = inbreeding)
summary(fm)

fm <- glm(survival ~ FPED, family = binomial(link = "logit"), data = inbreeding)
summary(fm)




####Testing if FPED correlated with FROH and heterozygosity####

cor.test(inbreeding$

####Testing if FPED underestimated FROH####