{"id":4622,"date":"2015-01-14T08:21:05","date_gmt":"2015-01-14T07:21:05","guid":{"rendered":"http:\/\/emilkirkegaard.dk\/en\/?p=4622"},"modified":"2015-01-14T08:25:56","modified_gmt":"2015-01-14T07:25:56","slug":"predicting-immigrant-performance-does-inbreeding-have-incremental-validity-over-iq-and-islam","status":"publish","type":"post","link":"https:\/\/emilkirkegaard.dk\/en\/2015\/01\/predicting-immigrant-performance-does-inbreeding-have-incremental-validity-over-iq-and-islam\/","title":{"rendered":"Predicting immigrant performance: Does inbreeding have incremental validity over IQ and Islam?"},"content":{"rendered":"<div class=\"oceanwp-oembed-wrap clr\">\n<blockquote class=\"twitter-tweet\" data-width=\"550\" data-dnt=\"true\">\n<p lang=\"en\" dir=\"ltr\"><a href=\"https:\/\/twitter.com\/hbdchick?ref_src=twsrc%5Etfw\">@hbdchick<\/a> if u can produce an inbreeding measure by country, we can test it. or u can, the data are public. :)<\/p>\n<p>&mdash; Emil O W Kirkegaard (@KirkegaardEmil) <a href=\"https:\/\/twitter.com\/KirkegaardEmil\/status\/555217814488092673?ref_src=twsrc%5Etfw\">January 14, 2015<\/a><\/p><\/blockquote>\n<p><script async src=\"https:\/\/platform.twitter.com\/widgets.js\" charset=\"utf-8\"><\/script><\/div>\n<p>So, she came up with:<\/p>\n<div class=\"oceanwp-oembed-wrap clr\">\n<blockquote class=\"twitter-tweet\" data-width=\"550\" data-dnt=\"true\">\n<p lang=\"en\" dir=\"ltr\"><a href=\"https:\/\/twitter.com\/KirkegaardEmil?ref_src=twsrc%5Etfw\">@KirkegaardEmil<\/a> good (weighted) coefficients of inbreeding by country at the back of this paper: <a href=\"http:\/\/t.co\/FjBZCHufqi\">http:\/\/t.co\/FjBZCHufqi<\/a> but, again&#8230;<\/p>\n<p>&mdash; hbd chick \ud83d\udc40 (@hbdchick) <a href=\"https:\/\/twitter.com\/hbdchick\/status\/555223247244242944?ref_src=twsrc%5Etfw\">January 14, 2015<\/a><\/p><\/blockquote>\n<p><script async src=\"https:\/\/platform.twitter.com\/widgets.js\" charset=\"utf-8\"><\/script><\/div>\n<p>So I decided to try it out, since I&#8217;m taking a break from reading Lilienfeld which I had been doing that for 5 hours straight or so.<\/p>\n<p>So the question is whether inbreeding measures have incremental validity over IQ and Islam, which I have previously used to examine immigrant performance <a href=\"http:\/\/openpsych.net\/ODP\/tag\/country-of-origin\/\">in a number of studies<\/a>.<\/p>\n<p>So, to get the data into R, I OCR&#8217;d the PDF in Abbyy FineReader since this program allows for easy copying of table data by row or column. I only wanted column 1-2 and didn&#8217;t want to deal with the hassle of importing it with spreadsheet problems (which need a consistent separator, e.g. comma or space). Then I merged it with <a href=\"https:\/\/osf.io\/zdcbq\/files\/\">the megadataset<\/a> to create a new version, 2.0d.<a href=\"https:\/\/osf.io\/zdcbq\/files\/\"><br \/>\n<\/a><\/p>\n<p>Then I created a subset of the data with variables of interest, and renamed them (otherwise results would be unwieldy). Intercorrelations are:<\/p>\n<table>\n<thead>\n<tr>\n<td id=\"origin\"><\/td>\n<th>row.names<\/th>\n<th>Cousin%<\/th>\n<th>CoefInbreed<\/th>\n<th>IQ<\/th>\n<th>Islam<\/th>\n<th>S.in.DK<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td class=\"rn\">1<\/td>\n<td>Cousin%<\/td>\n<td>1.00<\/td>\n<td>0.52<\/td>\n<td>-0.59<\/td>\n<td>0.78<\/td>\n<td>-0.76<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">2<\/td>\n<td>CoefInbreed<\/td>\n<td>0.52<\/td>\n<td>1.00<\/td>\n<td>-0.28<\/td>\n<td>0.40<\/td>\n<td>-0.55<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">3<\/td>\n<td>IQ<\/td>\n<td>-0.59<\/td>\n<td>-0.28<\/td>\n<td>1.00<\/td>\n<td>-0.27<\/td>\n<td>0.54<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">4<\/td>\n<td>Islam<\/td>\n<td>0.78<\/td>\n<td>0.40<\/td>\n<td>-0.27<\/td>\n<td>1.00<\/td>\n<td>-0.71<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">5<\/td>\n<td>S.in.DK<\/td>\n<td>-0.76<\/td>\n<td>-0.55<\/td>\n<td>0.54<\/td>\n<td>-0.71<\/td>\n<td>1.00<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>Spearman&#8217; correlations, which are probably better due to the non-normal data:<\/p>\n<table>\n<thead>\n<tr>\n<td id=\"origin\"><\/td>\n<th>row.names<\/th>\n<th>Cousin%<\/th>\n<th>CoefInbreed<\/th>\n<th>IQ<\/th>\n<th>Islam<\/th>\n<th>S.in.DK<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td class=\"rn\">1<\/td>\n<td>Cousin%<\/td>\n<td>1.00<\/td>\n<td>0.91<\/td>\n<td>-0.63<\/td>\n<td>0.67<\/td>\n<td>-0.73<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">2<\/td>\n<td>CoefInbreed<\/td>\n<td>0.91<\/td>\n<td>1.00<\/td>\n<td>-0.55<\/td>\n<td>0.61<\/td>\n<td>-0.76<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">3<\/td>\n<td>IQ<\/td>\n<td>-0.63<\/td>\n<td>-0.55<\/td>\n<td>1.00<\/td>\n<td>-0.23<\/td>\n<td>0.72<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">4<\/td>\n<td>Islam<\/td>\n<td>0.67<\/td>\n<td>0.61<\/td>\n<td>-0.23<\/td>\n<td>1.00<\/td>\n<td>-0.61<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">5<\/td>\n<td>S.in.DK<\/td>\n<td>-0.73<\/td>\n<td>-0.76<\/td>\n<td>0.72<\/td>\n<td>-0.61<\/td>\n<td>1.00<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>The fairly high correlations of inbreeding measures with IQ and Islam mean that their contribution will likely be modest as incremental validity.<\/p>\n<p>However, let&#8217;s try modeling them. I create 7 models of interest and compile the primary measure of interest from them, <a href=\"https:\/\/en.wikipedia.org\/wiki\/Coefficient_of_determination#Adjusted_R2\">R<sup>2<\/sup> adjusted<\/a>, into an object. Looks like this:<\/p>\n<table>\n<thead>\n<tr>\n<td id=\"origin\"><\/td>\n<th>row.names<\/th>\n<th>R2 adj.<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td class=\"rn\">1<\/td>\n<td>S.in.DK ~ IQ+Islam<\/td>\n<td>0.5472850<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">2<\/td>\n<td>S.in.DK ~ IQ+Islam+CousinPercent<\/td>\n<td>0.6701305<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">3<\/td>\n<td>S.in.DK ~ IQ+Islam+CoefInbreed<\/td>\n<td>0.7489312<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">4<\/td>\n<td>S.in.DK ~ Islam+CousinPercent<\/td>\n<td>0.6776841<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">5<\/td>\n<td>S.in.DK ~ Islam+CoefInbreed<\/td>\n<td>0.7438711<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">6<\/td>\n<td>S.in.DK ~ IQ+CousinPercent<\/td>\n<td>0.5486674<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">7<\/td>\n<td>S.in.DK ~ IQ+CoefInbreed<\/td>\n<td>0.4979552<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>So we see that either of them adds a fair amount of incremental validity to the base model (line 1 vs. 2-3). They are in fact better than IQ if one substitutes them in (1 vs. 4-5). They can also substitute for Islam, but only with about the same predictive power (1 vs 6-7).<\/p>\n<p><strong>Replication for Norway<\/strong><\/p>\n<p>Replication for science is important. Let&#8217;s try Norwegian data. The Finnish and Dutch data are well-suited for this (too few immigrant groups, few outcome variables i.e. only crime)<\/p>\n<p>Pearson intercorrelations:<\/p>\n<table>\n<thead>\n<tr>\n<td id=\"origin\"><\/td>\n<th>row.names<\/th>\n<th>CousinPercent<\/th>\n<th>CoefInbreed<\/th>\n<th>IQ<\/th>\n<th>Islam<\/th>\n<th>S.in.NO<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td class=\"rn\">1<\/td>\n<td>CousinPercent<\/td>\n<td>1.00<\/td>\n<td>0.52<\/td>\n<td>-0.59<\/td>\n<td>0.78<\/td>\n<td>-0.78<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">2<\/td>\n<td>CoefInbreed<\/td>\n<td>0.52<\/td>\n<td>1.00<\/td>\n<td>-0.28<\/td>\n<td>0.40<\/td>\n<td>-0.46<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">3<\/td>\n<td>IQ<\/td>\n<td>-0.59<\/td>\n<td>-0.28<\/td>\n<td>1.00<\/td>\n<td>-0.27<\/td>\n<td>0.60<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">4<\/td>\n<td>Islam<\/td>\n<td>0.78<\/td>\n<td>0.40<\/td>\n<td>-0.27<\/td>\n<td>1.00<\/td>\n<td>-0.72<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">5<\/td>\n<td>S.in.NO<\/td>\n<td>-0.78<\/td>\n<td>-0.46<\/td>\n<td>0.60<\/td>\n<td>-0.72<\/td>\n<td>1.00<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>Spearman:<\/p>\n<table>\n<thead>\n<tr>\n<td id=\"origin\"><\/td>\n<th>row.names<\/th>\n<th>CousinPercent<\/th>\n<th>CoefInbreed<\/th>\n<th>IQ<\/th>\n<th>Islam<\/th>\n<th>S.in.NO<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td class=\"rn\">1<\/td>\n<td>CousinPercent<\/td>\n<td>1.00<\/td>\n<td>0.91<\/td>\n<td>-0.63<\/td>\n<td>0.67<\/td>\n<td>-0.77<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">2<\/td>\n<td>CoefInbreed<\/td>\n<td>0.91<\/td>\n<td>1.00<\/td>\n<td>-0.55<\/td>\n<td>0.61<\/td>\n<td>-0.71<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">3<\/td>\n<td>IQ<\/td>\n<td>-0.63<\/td>\n<td>-0.55<\/td>\n<td>1.00<\/td>\n<td>-0.23<\/td>\n<td>0.75<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">4<\/td>\n<td>Islam<\/td>\n<td>0.67<\/td>\n<td>0.61<\/td>\n<td>-0.23<\/td>\n<td>1.00<\/td>\n<td>-0.47<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">5<\/td>\n<td>S.in.NO<\/td>\n<td>-0.77<\/td>\n<td>-0.71<\/td>\n<td>0.75<\/td>\n<td>-0.47<\/td>\n<td>1.00<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>These look fairly similar to Denmark.<\/p>\n<p>And the regression results:<\/p>\n<table>\n<thead>\n<tr>\n<td id=\"origin\"><\/td>\n<th>row.names<\/th>\n<th>R2 adj.<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td class=\"rn\">1<\/td>\n<td>S.in.NO ~ IQ+Islam<\/td>\n<td>0.5899682<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">2<\/td>\n<td>S.in.NO ~ IQ+Islam+CousinPercent<\/td>\n<td>0.7053999<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">3<\/td>\n<td>S.in.NO ~ IQ+Islam+CoefInbreed<\/td>\n<td>0.7077162<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">4<\/td>\n<td>S.in.NO ~ Islam+CousinPercent<\/td>\n<td>0.6826272<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">5<\/td>\n<td>S.in.NO ~ Islam+CoefInbreed<\/td>\n<td>0.6222364<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">6<\/td>\n<td>S.in.NO ~ IQ+CousinPercent<\/td>\n<td>0.6080922<\/td>\n<\/tr>\n<tr>\n<td class=\"rn\">7<\/td>\n<td>S.in.NO ~ IQ+CoefInbreed<\/td>\n<td>0.5460777<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>Fairly similar too. If added, they have incremental validity (line 1 vs. 2-3). They perform better than IQ if substituted but not as much as in the Danish data (1 vs. 4-5). They can also substitute for Islam (1 vs. 6-7).<\/p>\n<p><strong>How to interpret?<\/strong><\/p>\n<p>Since inbreeding does not seem to have any direct influence on behavior that is reflected in the S factor, it is not so easy to interpret these findings. Inbreeding leads to various health problems and lower g in offspring, the latter which may have some effect. However, presumably, national IQs already reflect the lowered IQ from inbreeding, so there should be no additional effect there beyond national IQs. Perhaps inbreeding results in other psychological problems that are relevant.<\/p>\n<p>Another idea is that inbreeding rates reflect non-g psychological traits that are relevant to adapting to life in Denmark. Perhaps it is a useful measure of clanishness, would be reflected in hostility towards integration in Danish society (such as getting an education, or lack of sympathy\/antipathy towards ethnic Danes and resulting higher crime rates against them), which would be reflected in the S factor.<\/p>\n<p>The lack of relatively well established causal routes for interpreting the finding makes me somewhat cautious about how to interpret this.<\/p>\n<hr \/>\n<p>&nbsp;<\/p>\n<pre>##Code for mergining cousin marriage+inbreeding data with megadataset\r\ninbreed = read.table(\"clipboard\", sep=\"\\t\",header=TRUE, row.names=1) #load data from clipboard\r\nsource(\"merger.R\") #load mega functions\r\nmega20d = read.mega(\"Megadataset_v2.0d.csv\") #load latest megadataset\r\nnames = as.abbrev(rownames(inbreed)) #get abbreviated names\r\nrownames(inbreed) = names #set them as rownames\r\n\r\n#merge and save\r\nmega20e = merge.datasets(mega20d,inbreed,1) #merge to create v. 2.0e\r\nwrite.mega(mega20e,\"Megadataset_v2.0e.csv\") #save it\r\n\r\n#select subset of interesting data\r\ndk.data = subset(mega20e, selec=c(\"Weighted.mean.consanguineous.percentage.HobenEtAl2010\",\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"Weighted.mean.coefficient.of.inbreeding.HobenEtAl2010\",\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"LV2012estimatedIQ\",\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"IslamPewResearch2010\",\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.factor.in.Denmark.Kirkegaard2014\"))\r\ncolnames(dk.data) = c(\"CousinPercent\",\"CoefInbreed\",\"IQ\",\"Islam\",\"S.in.DK\") #shorter var names\r\nrcorr = rcorr(as.matrix(dk.data)) #correlation object\r\nView(round(rcorr$r,2)) #view correlations, round to 2\r\nrcorr.S = rcorr(as.matrix(dk.data),type = \"spearman\") #spearman correlation object\r\nView(round(rcorr.S$r,2)) #view correlations, round to 2\r\n\r\n#Multiple regression\r\nlibrary(QuantPsyc) #for beta coef\r\nresults = as.data.frame(matrix(data = NA, nrow=0, ncol = 1)) #empty matrix for results\r\ncolnames(results) = \"R2 adj.\"\r\nmodels = c(\"S.in.DK ~ IQ+Islam\", #base model,\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.DK ~ IQ+Islam+CousinPercent\", #1. inbreeding var\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.DK ~ IQ+Islam+CoefInbreed\", #2. inbreeding var\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.DK ~ Islam+CousinPercent\", #without IQ\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.DK ~ Islam+CoefInbreed\", #without IQ\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.DK ~ IQ+CousinPercent\", #without Islam\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.DK ~ IQ+CoefInbreed\") #without Islam\r\n\r\nfor (model in models){ #run all the models\r\n\u00a0 fit.model = lm(model, dk.data) #fit model\r\n\u00a0 sum.stats = summary(fit.model) #summary stats object\r\n\u00a0 summary(fit.model) #summary stats\r\n\u00a0 lm.beta(fit.model) #standardized betas\r\n\u00a0 results[model,] = sum.stats$adj.r.squared #add result to results object\r\n}\r\nView(results) #view results\r\n\r\n##Let's try Norway too\r\nno.data = subset(mega20e, selec=c(\"Weighted.mean.consanguineous.percentage.HobenEtAl2010\",\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"Weighted.mean.coefficient.of.inbreeding.HobenEtAl2010\",\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"LV2012estimatedIQ\",\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"IslamPewResearch2010\",\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.factor.in.Norway.Kirkegaard2014\"))\r\n\r\ncolnames(no.data) = c(\"CousinPercent\",\"CoefInbreed\",\"IQ\",\"Islam\",\"S.in.NO\") #shorter var names\r\nrcorr = rcorr(as.matrix(no.data)) #correlation object\r\nView(round(rcorr$r,2)) #view correlations, round to 2\r\nrcorr.S = rcorr(as.matrix(no.data),type = \"spearman\") #spearman correlation object\r\nView(round(rcorr.S$r,2)) #view correlations, round to 2\r\n\r\nresults = as.data.frame(matrix(data = NA, nrow=0, ncol = 1)) #empty matrix for results\r\ncolnames(results) = \"R2 adj.\"\r\nmodels = c(\"S.in.NO ~ IQ+Islam\", #base model,\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.NO ~ IQ+Islam+CousinPercent\", #1. inbreeding var\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.NO ~ IQ+Islam+CoefInbreed\", #2. inbreeding var\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.NO ~ Islam+CousinPercent\", #without IQ\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.NO ~ Islam+CoefInbreed\", #without IQ\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.NO ~ IQ+CousinPercent\", #without Islam\r\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \"S.in.NO ~ IQ+CoefInbreed\") #without Islam\r\n\r\nfor (model in models){ #run all the models\r\n\u00a0 fit.model = lm(model, no.data) #fit model\r\n\u00a0 sum.stats = summary(fit.model) #summary stats object\r\n\u00a0 summary(fit.model) #summary stats\r\n\u00a0 lm.beta(fit.model) #standardized betas\r\n\u00a0 results[model,] = sum.stats$adj.r.squared #add result to results object\r\n}\r\nView(results) #view results<\/pre>\n","protected":false},"excerpt":{"rendered":"<p>@hbdchick if u can produce an inbreeding measure by country, we can test it. or u can, the data are public. :) &mdash; Emil O W Kirkegaard (@KirkegaardEmil) January 14, 2015 So, she came up with: @KirkegaardEmil good (weighted) coefficients of inbreeding by country at the back of this paper: http:\/\/t.co\/FjBZCHufqi but, again&#8230; &mdash; hbd [&hellip;]<\/p>\n","protected":false},"author":17,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1839,2591,1921],"tags":[1991,2074,600],"class_list":["post-4622","post","type-post","status-publish","format-standard","hentry","category-psychometics","category-intelligence-iq-cognitive-ability","category-sociology","tag-immigrants","tag-inbreeding","tag-islam","entry"],"_links":{"self":[{"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/posts\/4622","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/users\/17"}],"replies":[{"embeddable":true,"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/comments?post=4622"}],"version-history":[{"count":4,"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/posts\/4622\/revisions"}],"predecessor-version":[{"id":4625,"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/posts\/4622\/revisions\/4625"}],"wp:attachment":[{"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/media?parent=4622"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/categories?post=4622"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/emilkirkegaard.dk\/en\/wp-json\/wp\/v2\/tags?post=4622"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}