NunoSempere comments on EA Mental Health Survey: Results and Analysis.

[NunoSempere]

Using the 2017 SSC, I’ve looked at:

Has been diagnosed with a mental illness ~ Is EA

lm(formula = Has_mental_condition_diagnosis ~ Is_EA) 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 0.244434   0.006544  37.352   <2e-16 ***
Is_EA       0.027635   0.017688   1.562    0.118    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Weak relationship, is not significant.

Has been diagnosed with a mental illness or intuits having one ~ Is EA

> summary(lm(Has_mental_condition_diagnosis_or_intuited ~ Is_EA))

Call:
lm(formula = Has_mental_condition_diagnosis_or_intuited ~ Is_EA)

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  0.52513    0.00756  69.463   <2e-16 ***
Is_EA        0.03782    0.02043   1.851   0.0642 .  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Weak relationship, 0.05 < p < 0.10

Has taken SSRIs ~ Is EA.

lm(formula = Has_taken_SSRIs ~ Is_EA)

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept)  0.276337   0.006728  41.070   <2e-16 ***
Is_EA       -0.041894   0.018186  -2.304   0.0213 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Weak negative relationship (??).

Code used to produce that (in R)

## Download the file from https://​​slatestarcodex.com/​​2017/​​03/​​17/​​ssc-survey-2017-results/​​ , then move into the directory the file is in with:
## setwd(“directory”)

C ← read.csv(file=”Survey_CSV.csv”, header=TRUE, sep=”,”, stringsAsFactors=FALSE)

c(69:77) → mental_conditions #  “Schizophrenia”, Autism”, “Depression”, “Anxiety”, “OCD”, “Eatingdisorder” “Bipolar”, “Alcoholism”, “Drugaddiction”.

Has_mental_condition_diagnosis=rep(0,dim(C)[1])
for(mental_condition in mental_conditions){
  ifelse(C[,mental_condition]==”I have a formal diagnosis of this condition”, 1, Has_mental_condition_diagnosis) → Has_mental_condition_diagnosis 
}

ifelse(C$EAID==”Yes”, 1,0) → Is_EA
summary(lm(Has_mental_condition_diagnosis ~ Is_EA))

Has_mental_condition_diagnosis_or_intuited = rep(0,dim(C)[1])
Diagnosis_or_intuited = c(“I have a formal diagnosis of this condition”,”I think I might have this condition, although I have never been formally diagnosed” )
for(mental_condition in mental_conditions){
    ifelse(C[,mental_condition]%in%Diagnosis_or_intuited, 1, Has_mental_condition_diagnosis_or_intuited) → Has_mental_condition_diagnosis_or_intuited 
  }

summary(lm(Has_mental_condition_diagnosis_or_intuited ~ Is_EA))

ifelse(C$SSRIs == “I have never taken these drugs”, 0, 1) → Has_taken_SSRIs
summary(lm(Has_taken_SSRIs ~ Is_EA))