Hey1

Models

All models for this gene with cis-heritability P-value <0.01 were tested and are shown here. In each tissue, a gene can have zero, one, or multiple RNA phenotypes for each RNA modality. After the heritability statistics, subsequent columns show cross-validation performance (R² and P-value) of each modeling method. The best-performing method for each models is used for TWAS.

#	Tissue	Modality	RNA phenotype	h²	h² SE	h² P	xQTL R²	BLUP R²	ENET R²	LASSO R²	xQTL P	BLUP P	ENET P	LASSO P
1	PL	gene expression	Hey1	0.16	0.12	5.7e-08	0.05	0.04	0.03	0.04	5.9e-06	4.8e-05	3.9e-04	2.9e-05
2	IC	gene expression	Hey1	0.12	0.09	1.4e-03	0.09	0.06	0.06	0.08	1.6e-04	1.7e-03	1.6e-03	3.3e-04
3	BLA	gene expression	Hey1	0.08	0.06	7.4e-03	0	0	0.01	0	1.8e-01	2.4e-01	1.1e-01	3.2e-01
4	Liver	gene expression	Hey1	0.05	0.03	1.0e-03	0.01	0.02	0.01	0.01	3.5e-02	3.4e-03	6.7e-03	1.8e-02
5	NAcc	gene expression	Hey1	0.03	0.02	2.0e-03	0	0.01	0	0	8.6e-02	9.4e-03	4.3e-02	4.7e-02

Trait associations

The numbered columns correspond to the numbered models in the table above. They contain the TWAS Z-score for each trait-model pair. Colors indicate strength of Z-scores: cells with lower negative Z-scores are darker blue, cells with higher positive Z-scores are darker red. Significant associations have white and bolded text.

Trait	Avg chi² ratio	Avg chi²	Max chi²	1	2	3	4	5
Bone: apparent density	1.1	1.8	2.9	-1.1	-1.0	-1.7	1.5	-1.3
Area under the delay curve	0.9	1.1	1.9	0.9	0.8	0.7	-1.4	1.2
BMI with tail	1.6	5.6	11.9	3.4	3.4	2.0	-0.8	0.1
BMI without tail	2.7	8.2	16.8	4.0	4.1	2.7	-0.8	0.0
Body weight	1.3	5.5	11.3	3.4	3.3	2.2	-0.0	-0.6
Weight adjusted by age	0.1	0.2	0.7	0.1	-0.1	0.8	0.3	-0.3
Bone surface	1.0	1.8	3.8	-0.6	-0.6	-1.0	2.0	-1.9
Bone volume	1.2	2.2	4.3	-0.8	-0.8	-1.3	2.1	-1.9
Breakpoint value in progressive ratio session	0.2	0.2	0.3	0.5	0.5	0.5	0.2	-0.3
Number of licking bursts	0.0	0.0	0.1	0.3	0.1	0.2	-0.3	0.2
Context. condit. distance diff. score	0.2	0.3	0.6	0.1	0.4	0.2	-0.8	0.7
Locomotion velocity, session 1	0.3	0.3	0.7	0.7	0.9	0.5	-0.2	0.1
Locomotion distance, session 1	0.4	0.5	0.9	0.8	0.9	0.1	0.6	-0.8
Locomotion velocity, session 2	1.0	1.2	2.7	1.6	1.6	0.8	0.0	-0.4
Locomotion distance, session 2	1.1	1.5	3.5	0.8	1.0	-0.1	1.6	-1.9
Locomotion velocity, session 3	0.2	0.2	0.6	-0.2	0.0	-0.1	-0.7	0.8
Locomotion distance, session 3	0.4	0.4	0.8	-0.4	-0.2	-0.7	0.9	-0.8
Stereotopy head waving bouts, day 3	0.1	0.1	0.1	-0.4	-0.4	-0.2	0.3	-0.2
Stereotopy head waving duration, day 3	0.5	0.6	0.8	0.6	0.6	0.9	-0.9	0.8
Locomotion velocity, session 7	9.6	10.5	14.9	-2.6	-2.4	-3.7	3.9	-3.4
Locomotion distance, session 7	4.4	5.6	7.8	-2.6	-2.5	-2.8	2.1	-1.7
Stereotopy head waving bouts, day 7	0.4	0.5	1.0	1.0	1.0	0.6	-0.1	-0.0
Stereotopy head waving duration, day 7	0.9	1.1	1.9	1.3	1.4	1.0	-0.8	0.6
Locomotion velocity, session 8	1.4	1.7	3.8	1.7	2.0	1.2	-0.3	-0.1
Locomotion distance, session 8	0.4	0.5	1.3	0.8	1.1	0.4	0.4	-0.6
Degree of sensitization distance	0.8	0.8	1.1	-1.0	-1.0	-1.1	0.7	-0.5
Degree of sensitization stereotypy	1.8	2.0	3.7	1.7	1.9	1.5	-0.9	0.5
Conditioned locomotion	0.3	0.4	0.7	-0.1	-0.2	-0.6	0.8	-0.8
Cocaine response after cond. corrected	0.9	1.1	2.3	-1.5	-1.5	-1.0	0.1	0.1
Cocaine response after cond. not corrected	2.2	2.6	5.0	-2.2	-2.2	-1.7	0.5	-0.1
Cocaine response before conditioning	0.6	0.7	1.0	-0.9	-0.9	-1.0	0.7	-0.5
Saline control response	0.4	0.5	1.2	0.2	0.3	0.0	1.1	-1.1
Bone: connectivity density	1.2	1.9	2.7	-1.2	-1.2	-1.6	1.5	-1.3
Bone: cortical apparent density	0.2	0.2	0.5	0.6	0.7	0.4	-0.4	0.3
Bone: cortical area	1.3	1.9	3.4	1.8	1.8	1.6	0.0	-0.4
Bone: cortical porosity	0.0	0.0	0.1	-0.1	-0.1	0.2	0.2	-0.2
Bone: cortical porosity	0.0	0.0	0.1	0.1	0.0	0.1	0.3	-0.3
Bone: cortical thickness	0.9	1.4	2.8	1.6	1.7	1.3	0.0	-0.3
Bone: cortical thickness	3.9	4.6	8.4	2.9	2.8	2.2	-1.2	0.7
Bone: cortical tissue density	0.2	0.2	0.4	0.6	0.6	0.3	-0.4	0.3
Condit. Reinf. active minus inactive responses	3.2	3.6	6.1	1.9	1.8	0.3	2.2	-2.5
Condit. Reinf. active-inactive response ratio	4.6	5.0	9.0	2.1	1.9	0.2	2.8	-3.0
Condit. Reinf. active responses	2.9	3.5	6.0	1.9	1.8	0.3	2.2	-2.4
Condit. Reinf. inactive responses	0.3	0.3	0.7	0.7	0.8	0.6	-0.3	0.1
Incentive salience index mean	3.0	3.7	6.3	2.0	1.8	0.6	2.2	-2.5
Condit. Reinf. lever presses	1.4	2.0	2.9	1.7	1.6	1.0	1.1	-1.5
Condit. Reinf. active minus inactive responses	1.8	2.3	4.3	-0.8	-0.7	-1.6	2.1	-1.9
Condit. Reinf. active-inactive response ratio	1.2	1.4	3.1	0.1	0.1	-0.9	1.7	-1.8
Condit. Reinf. active responses	3.2	4.4	7.1	-1.5	-1.5	-2.2	2.7	-2.4
Condit. Reinf. inactive responses	0.9	1.1	1.5	-1.1	-1.0	-1.0	1.2	-1.0
Condit. Reinf. lever presses	0.8	1.2	1.9	-0.8	-0.8	-1.1	1.4	-1.2
Condit. Reinf. lever reinforcers received	2.9	4.0	6.2	-1.5	-1.4	-2.2	2.5	-2.2
Conditioned reinforcement - actives	0.8	1.3	1.9	-1.4	-1.4	-0.3	1.3	-1.0
Conditioned reinforcement - inactives	0.5	0.8	1.0	-0.9	-0.9	-0.6	1.0	-0.9
Conditioned locomotion	0.5	0.6	0.9	-0.9	-0.9	-0.3	0.9	-0.7
Patch foraging indifference point 0 sec	0.4	0.8	1.6	-1.2	-1.3	-0.9	-0.1	0.4
Patch foraging indifference point AUC	0.3	0.5	1.0	-0.8	-0.7	-1.0	0.5	-0.4
Delay disc. indifference point, 0s delay	0.2	0.2	0.4	-0.0	0.0	0.4	-0.6	0.6
Delay disc. indifference point, 16s delay	0.2	0.2	0.3	0.6	0.5	0.1	-0.6	0.4
Delay disc. indifference point, 2s delay	0.8	0.9	1.9	1.4	1.1	0.7	0.6	-0.8
Delay disc. indifference point, 24s delay	0.2	0.3	0.6	0.1	0.2	0.6	-0.7	0.7
Delay disc. indifference point, 4s delay	2.6	2.8	3.4	1.8	1.8	1.8	-1.5	1.2
Delay disc. indifference point, 8s delay	0.7	0.9	2.1	-0.2	-0.2	0.3	-1.4	1.5
Indifference point function ln k	0.9	1.5	2.3	1.1	1.1	1.5	-1.3	1.0
Indifference point function log k	0.9	1.5	2.3	1.1	1.1	1.5	-1.3	1.0
Patch foraging total patch changes 0 sec	0.4	0.7	1.3	-1.1	-1.1	-0.7	-0.3	0.6
Patch foraging total patch changes 12 sec	0.8	1.0	1.8	-1.3	-1.3	-0.7	-0.5	0.8
Patch foraging total patch changes 18 sec	0.0	0.0	0.1	-0.1	-0.0	-0.1	-0.2	0.2
Patch foraging total patch changes 24 sec	0.2	0.4	0.7	-0.8	-0.8	-0.8	-0.2	0.3
Patch foraging total patch changes 6 sec	0.2	0.2	0.5	-0.6	-0.7	-0.4	-0.2	0.4
Patch foraging time to switch 0 sec	0.9	1.8	2.5	1.5	1.6	0.7	1.2	-1.6
Patch foraging water rate 0 sec	0.3	0.6	1.5	-0.3	-0.3	-0.5	-1.1	1.2
Patch foraging water rate 12 sec	0.1	0.3	0.7	-0.8	-0.8	-0.4	0.0	0.1
Patch foraging water rate 18 sec	0.1	0.1	0.3	-0.1	-0.1	-0.6	0.5	-0.4
Patch foraging water rate 24 sec	0.3	0.4	1.0	0.1	0.1	0.1	-1.0	1.0
Patch foraging water rate 6 sec	0.2	0.3	0.5	-0.7	-0.7	-0.3	-0.5	0.7
Seeking ratio, delayed vs. immediate footshock	3.6	4.3	8.1	-2.6	-2.6	-2.8	0.6	-0.1
Distance traveled before self-admin	2.2	2.9	5.0	2.2	2.1	1.7	-1.2	0.7
Distance traveled after self-admin	0.8	0.9	1.6	1.3	1.2	1.0	-0.7	0.4
Extensor digitorum longus weight	7.1	27.8	56.1	7.5	7.5	5.0	-1.3	-0.1
Bone: elastic displacement	0.4	0.5	1.2	0.0	0.1	-0.2	1.1	-1.0
Bone: elastic work	0.3	0.3	0.8	0.3	0.3	-0.1	0.8	-0.9
Bone: endosteal estimation	0.7	1.6	3.1	1.1	1.0	0.4	1.5	-1.8
Bone: endosteal perimeter	0.1	0.3	0.8	0.2	0.1	-0.1	0.8	-0.9
Epididymis fat weight	0.1	0.5	0.8	0.9	0.8	-0.8	-0.4	0.3
Delta time in closed arm before/after self-admin	1.1	1.1	1.5	1.2	1.1	-0.4	-1.2	1.0
Delta time in open arm before/after self-admin	0.3	0.3	0.4	-0.6	-0.6	-0.3	0.6	-0.5
Diff in mean of infusions in LGA sessions	1.1	1.4	2.3	1.5	1.5	1.4	0.3	-0.6
Delay discount exponential model param	0.9	1.0	1.9	-0.8	-0.7	-0.6	1.4	-1.2
Extinction: sum of active levers before priming	0.2	0.3	0.5	-0.3	-0.3	-0.6	0.7	-0.5
Active lever presses in extinction session 6	0.1	0.1	0.2	-0.0	-0.0	0.0	-0.4	0.4
Bone: final force	1.1	1.8	3.4	1.8	1.8	1.1	0.5	-0.9
Bone: final moment	1.0	1.6	3.0	1.7	1.7	1.1	0.5	-0.8
Food consumed during 24 hour testing period	1.2	1.6	2.3	1.5	1.5	1.4	-1.0	0.8
Glucose	0.9	1.4	2.3	1.5	1.5	0.8	-1.1	0.8
Heart weight	0.2	0.3	0.4	0.3	0.3	-0.5	-0.6	0.6
Delay discount hyperbolic model param	0.9	1.0	1.4	-1.1	-0.9	-0.8	1.2	-1.0
Intermittent access intake day 1-15 change	0.6	1.2	1.8	-1.1	-1.1	-1.3	1.0	-0.8
Intermittent access intake escalation	0.9	1.9	4.0	-1.4	-1.4	-2.0	0.9	-0.7
Intermittent access intake escalation 2	2.9	3.3	5.6	-2.1	-2.2	-2.4	1.2	-0.8
Intermitt. access day 1 inactive lever presses	0.3	1.0	1.5	0.9	0.9	1.2	-1.0	0.9
Intermitt. access day 15 inactive lever presses	1.0	1.9	3.0	-1.7	-1.7	-0.8	-1.1	1.4
Intermitt. access escalation Index	0.4	0.5	1.0	-0.2	-0.1	-0.9	-0.9	1.0
Intermittent access day 1 total infusions	2.5	4.0	5.4	-2.3	-2.3	-1.3	2.1	-1.8
Intermittent access day 15 total infusions	3.8	5.7	7.8	-2.8	-2.8	-2.4	2.1	-1.7
Intermittent access terminal intake (last 3 days)	4.0	5.3	8.2	-2.9	-2.8	-2.2	1.8	-1.4
Intermittent access total infusions	3.0	3.5	5.8	-2.4	-2.4	-1.6	1.5	-1.2
Intermittent access day 1 locomotion	1.0	1.2	1.8	-1.3	-1.3	-0.7	1.1	-1.0
Intermittent access day 15 total locomotion	0.3	0.6	1.2	-0.4	-0.5	-0.5	1.1	-1.1
Intermittent access total locomotion	1.8	2.2	3.1	-1.7	-1.8	-1.6	1.2	-1.0
Intermittent access standard deviation	0.7	1.5	3.0	-1.7	-1.7	-1.2	0.2	0.0
Left kidney weight	0.3	1.2	3.1	0.0	0.1	-0.1	-1.7	1.8
Right kidney weight	0.2	0.8	1.5	0.8	0.9	0.4	-1.2	1.1
Cocaine induced anxiety	0.1	0.1	0.2	0.4	0.4	-0.3	-0.1	0.1
Post-drug Anxiety	0.7	0.8	2.0	1.4	1.3	-0.1	0.2	-0.4
Baseline Anxiety	0.4	0.4	0.8	0.9	0.9	0.7	0.2	-0.3
Tail length	0.2	0.9	1.7	1.3	1.2	1.1	0.0	-0.2
Length with tail	0.1	0.6	1.0	0.9	0.7	1.0	0.6	-0.7
Length without tail	0.1	0.3	0.5	-0.5	-0.6	-0.4	0.7	-0.6
Times rat made contact with spout	0.3	0.3	0.5	0.3	0.2	0.7	-0.7	0.6
Lifetime Intake	3.1	3.8	7.0	-2.7	-2.6	-1.8	1.2	-0.8
Liver weight, left	0.7	1.1	1.8	1.3	1.3	1.1	0.5	-0.7
Liver weight, right	0.8	1.4	2.6	0.2	0.3	1.2	-1.6	1.6
Locomotor activity	0.9	1.4	3.7	1.1	1.2	1.9	-0.8	0.5
Locomotor testing distance	0.5	0.8	3.0	0.7	0.7	1.7	-0.4	0.1
Locomotor testing rearing	0.4	0.8	1.1	-0.9	-0.9	0.5	-0.9	1.0
Locomotion in novel chamber	1.1	1.4	2.2	1.5	1.5	1.2	0.5	-0.9
Locomotion in novel chamber post-restriction	0.4	0.6	1.3	0.4	0.4	-0.0	-1.2	1.1
Fecal boli incidents, locomotor time 1	0.3	0.5	1.6	-0.3	-0.4	-1.3	0.6	-0.5
Change in fecal boli incidents, locomotor task	0.4	0.4	0.7	-0.8	-0.7	0.2	-0.4	0.6
Time >=10cm from walls, locomotor time 1	0.0	0.0	0.1	-0.1	-0.2	-0.3	0.2	-0.1
Bouts of movement, locomotor time 1	1.9	2.2	5.6	-0.5	-0.5	0.5	-2.2	2.4
Total resting periods, locomotor time 1	1.0	1.1	2.5	-0.7	-0.7	0.1	-1.4	1.6
Rest time, locomotor task time 1	0.3	0.4	0.7	-0.3	-0.4	-0.4	-0.8	0.9
Distance moved, locomotor task time 1	0.4	0.6	1.3	1.0	1.1	0.4	0.4	-0.6
Vertical activity count, locomotor time 1	0.6	1.5	2.7	1.6	1.6	1.0	0.6	-0.9
Time >=10cm from walls, locomotor time 2	0.1	0.2	0.3	0.4	0.6	0.2	0.5	-0.5
Bouts of movement, locomotor time 2	1.9	3.0	6.4	0.3	0.4	1.5	-2.5	2.5
Total resting periods, locomotor time 2	2.3	4.0	8.7	0.2	0.4	1.6	-2.9	3.0
Rest time, locomotor task time 2	1.4	1.8	3.7	0.4	0.2	1.2	-1.9	1.9
Distance moved, locomotor task time 2	0.2	0.3	0.6	0.5	0.8	-0.1	0.4	-0.6
Vertical activity count, locomotor time 2	1.9	2.6	5.2	2.3	2.3	0.7	0.8	-1.3
Light reinforcement 1	0.1	0.2	0.5	0.1	0.3	-0.4	-0.7	0.6
Light reinforcement 2	1.1	1.7	4.3	2.0	2.1	-0.0	-0.3	-0.1
Bone: marrow area	0.8	1.7	3.1	1.3	1.2	0.5	1.5	-1.8
Bone: maximum diameter	0.9	1.7	2.9	1.7	1.6	1.6	-0.6	0.2
Bone: maximum force	0.8	1.4	2.8	1.7	1.6	1.3	0.0	-0.4
Bone: maximum moment	0.7	1.2	2.3	1.5	1.4	1.2	0.0	-0.3
Mean time between licks in bursts	0.6	1.5	3.2	-0.4	-0.5	-0.8	1.8	-1.7
Mean num. licks in bursts	0.0	0.0	0.1	0.0	0.1	0.3	-0.1	0.2
Cd content in liver	0.9	1.0	1.6	0.7	0.8	0.8	-1.3	1.2
Co content in liver	1.6	1.8	4.5	0.3	0.2	0.5	-2.1	2.1
Cu content in liver	1.6	1.9	4.8	-2.2	-2.0	-0.7	-0.3	0.6
Fe content in liver	0.1	0.1	0.2	-0.4	-0.3	-0.1	0.5	-0.4
K content in liver	1.1	1.1	2.0	-1.4	-1.4	-1.1	0.4	-0.3
Mg content in liver	1.2	1.2	1.7	-1.3	-1.3	-1.1	0.9	-0.8
Mn content in liver	0.4	0.5	0.8	-0.7	-0.7	-0.9	0.6	-0.5
Na content in liver	0.1	0.1	0.2	-0.1	-0.1	-0.1	-0.5	0.5
Rb content in liver	0.1	0.1	0.3	-0.6	-0.5	-0.3	-0.1	0.2
Se content in liver	0.1	0.2	0.4	-0.6	-0.6	-0.1	-0.3	0.3
Sr content in liver	0.3	0.3	0.6	-0.3	-0.3	-0.3	0.8	-0.7
Zn content in liver	2.0	2.2	4.9	-2.2	-2.2	-1.0	-0.1	0.5
Bone: minimum diameter	2.2	4.0	7.0	1.9	1.8	0.9	2.2	-2.6
Time in familiar zone, hab. session 1	1.3	1.3	2.9	0.1	0.1	0.8	-1.7	1.7
Time in novel zone, hab. session 1	0.8	0.8	1.5	-0.8	-0.9	-1.2	0.7	-0.6
Total zone transitions, hab. session 1	0.4	0.5	0.9	0.6	0.8	-0.3	0.9	-0.9
Total locomotion distance, hab. session 1	0.8	1.1	2.6	0.2	0.2	-0.5	1.6	-1.6
Locomotion velocity, hab. session 1	0.9	1.3	2.9	0.7	0.7	-0.1	1.6	-1.7
Time in familiar zone, hab. session 2	1.7	1.8	2.9	-1.3	-1.2	-1.1	-1.4	1.7
Time in novel zone, hab. session 2	0.9	1.0	1.5	0.9	0.9	0.6	1.1	-1.2
Total zone transitions, hab. session 2	0.0	0.1	0.1	0.1	0.3	-0.2	-0.3	0.2
Total locomotion distance, hab. session 2	0.3	0.4	0.7	0.5	0.5	-0.4	0.8	-0.8
Locomotion velocity, hab. session 2	0.2	0.3	0.6	0.3	0.3	-0.5	0.7	-0.8
Novel to familiar place preference ratio	0.5	0.5	1.1	0.9	0.9	1.0	-0.1	-0.0
Novelty place preference	0.5	0.5	1.1	0.8	0.8	1.1	-0.3	0.2
Time in novel zone, NPP test	0.8	0.9	1.5	1.0	1.0	1.2	-0.7	0.6
Total zone transitions, NPP test	0.1	0.1	0.3	0.5	0.5	-0.4	0.1	-0.1
Total locomotion distance, NPP test	0.2	0.3	0.5	0.5	0.6	-0.2	0.6	-0.7
Locomotion velocity, NPP test	0.3	0.4	1.0	0.1	0.1	-0.5	1.0	-1.0
Intraocular pressure	1.5	2.2	5.1	2.3	2.2	0.9	-0.1	-0.3
Delta distance traveled before/after self-admin	0.9	1.0	1.3	-1.2	-1.1	-0.7	1.0	-0.8
Delta ambulatory time before/after self-admin	1.2	1.5	2.5	-1.6	-1.5	-0.9	1.1	-0.7
Parametrial fat weight	0.1	0.2	0.5	-0.1	-0.1	-0.3	0.7	-0.7
Pavlov. Cond. lever latency	0.6	0.7	1.9	-0.6	-0.5	-0.1	-1.1	1.4
Pavlov. Cond. magazine entry latency	0.3	0.4	0.9	0.9	0.9	0.0	0.4	-0.6
Pavlov. Cond. change in total contacts	1.7	1.8	3.4	1.9	1.7	1.0	0.6	-1.0
Pavlov. Cond. index score	0.4	0.5	0.9	1.0	0.9	0.3	0.5	-0.8
Pavlov. Cond. latency score	0.3	0.4	0.7	0.8	0.7	0.2	0.5	-0.8
Pavlov. Cond. lever contacts	1.1	1.3	3.4	0.6	0.5	0.1	1.5	-1.8
Pavlov. Cond. magazine entry number	0.7	0.9	2.1	-1.4	-1.4	-0.5	-0.1	0.3
Pavlov. Cond. intertrial magazine entries	1.5	1.9	3.9	-1.9	-2.0	-0.7	0.9	-0.5
Pavlov. Cond. lever-magazine prob. diff.	0.4	0.5	0.9	0.9	0.9	0.5	0.3	-0.6
Pavlov. Cond. response bias	1.0	1.2	2.1	1.4	1.4	0.4	0.8	-1.1
Pavlov. Cond. lever latency	0.7	1.0	1.6	1.1	1.3	1.2	-0.6	0.4
Pavlov. Cond. magazine entry latency	0.1	0.1	0.1	0.4	0.2	-0.2	0.2	-0.3
Pavlov. Cond. change in total contacts	0.7	1.1	2.4	0.1	0.0	-0.8	1.5	-1.6
Pavlov. Cond. index score	0.3	0.5	0.9	0.0	-0.1	-0.8	0.9	-0.9
Pavlov. Cond. latency score	0.3	0.4	0.7	-0.2	-0.3	-0.8	0.8	-0.7
Pavlov. Cond. lever contacts	0.7	0.9	1.5	-0.9	-1.0	-1.2	0.9	-0.7
Pavlov. Cond. magazine entry number	0.2	0.3	0.6	-0.8	-0.6	0.2	-0.4	0.5
Pavlov. Cond. intertrial magazine entries	0.1	0.2	0.5	0.1	0.3	0.2	0.6	-0.7
Pavlov. Cond. lever-magazine prob. diff.	0.4	0.5	0.8	-0.2	-0.3	-0.9	0.9	-0.8
Pavlov. Cond. response bias	0.5	0.7	1.4	-0.1	-0.3	-0.9	1.2	-1.1
Time in closed arm before self-admin	1.0	1.2	2.6	-1.6	-1.6	-1.0	0.1	0.3
Time in closed arm after self-admin	1.1	1.2	2.8	0.6	0.5	0.3	-1.7	1.6
Time in open arm before self-admin	0.5	0.5	0.9	0.9	0.9	0.9	0.2	-0.4
Time in open arm after self-admin	0.5	0.6	1.0	0.9	1.0	0.8	0.2	-0.5
Bone: periosteal estimation	1.5	3.2	4.9	2.2	2.1	1.3	1.3	-1.8
Bone: periosteal perimeter	0.4	0.6	1.2	1.1	1.0	0.9	0.2	-0.4
Bone: post-yield work	0.8	1.2	1.7	-1.2	-1.3	-0.5	-1.0	1.3
Progressive ratio test 1 active lever presses	1.9	2.4	4.5	-2.1	-2.1	-1.7	0.5	-0.3
Progressive ratio test 1 breakpoint	2.0	2.5	4.6	-2.1	-2.1	-1.7	0.6	-0.4
Progressive ratio test 1 inactive lever presses	0.4	0.6	0.9	-0.7	-0.7	-0.9	0.8	-0.7
Progressive ratio test 2 active lever presses	2.0	2.3	3.0	-1.7	-1.7	-1.2	1.6	-1.4
Progressive ratio test 2 breakpoint	2.1	2.4	3.1	-1.8	-1.7	-1.2	1.6	-1.4
Progressive ratio test 2 inactive lever presses	0.1	0.2	0.4	-0.6	-0.6	0.1	-0.0	0.1
Incentive Sensitization - Responses	1.3	1.8	4.0	-0.9	-0.9	-0.4	2.0	-1.8
Incentive Sensitization - Breakpoint	1.2	1.6	3.4	-0.9	-0.9	-0.3	1.9	-1.7
Sum of active levers in priming session hrs 5-6	0.9	0.8	1.2	-0.9	-0.9	-0.7	1.1	-0.9
Progressive ratio	1.3	2.0	4.8	0.6	0.5	0.7	-2.2	2.0
Food seeking constrained by brief footshock	2.6	3.2	6.1	2.5	2.3	0.5	1.2	-1.7
Seeking ratio, punishment vs. effort	2.0	2.9	6.7	0.6	0.6	-0.6	2.5	-2.6
Retroperitoneal fat weight	0.1	0.3	0.6	0.5	0.5	-0.8	0.5	-0.6
Reaction time number correct	0.2	0.5	1.1	-0.4	-0.4	0.4	1.1	-0.9
Reaction time mean minus median	0.5	1.3	2.5	0.6	0.6	1.2	-1.6	1.4
Reaction time mean minus median AUC	0.6	1.4	2.2	1.0	0.9	1.5	-1.3	1.1
Reaction time num false alarms	1.2	1.7	3.5	-0.6	-0.8	-0.9	1.9	-1.8
Reaction time num false alarms AUC	1.1	1.6	3.1	-0.8	-1.0	-1.0	1.8	-1.6
Reaction time trials correct on left	0.2	0.5	1.1	-0.4	-0.4	0.4	1.1	-0.9
Reaction time trials on left	0.3	0.5	1.2	-0.4	-0.4	0.4	1.1	-1.0
Reaction time mean	0.5	0.8	1.7	1.3	1.3	0.6	-0.1	-0.1
Reaction time mean AUC	0.5	0.8	1.9	1.4	1.3	0.6	-0.1	-0.2
Median of all reaction times	0.7	1.0	2.0	1.4	1.4	0.3	0.6	-0.8
Reaction time omissions	0.1	0.1	0.1	0.3	0.3	-0.2	0.3	-0.4
Reaction time false alarm rate	0.3	0.3	0.9	-0.3	-0.4	-0.9	0.6	-0.5
Reaction time premature initiation rate	0.9	1.2	1.7	-1.2	-1.2	-0.5	-1.1	1.3
Reaction time premature initiations	0.5	0.9	1.8	-1.4	-1.3	-0.0	-0.6	0.8
Std. dev. reaction times	0.6	1.1	2.0	0.8	0.8	1.0	-1.4	1.2
Reaction time trials completed	0.3	0.5	1.2	-0.4	-0.4	0.4	1.1	-1.0
Reaction time trials AUC	0.3	0.7	1.6	-0.4	-0.4	0.3	1.3	-1.1
Run reversals in cocaine runway, females	1.5	1.4	2.8	1.7	1.7	1.2	-0.6	0.2
Run reversals in cocaine runway, males	1.5	1.5	2.8	1.7	1.7	1.4	0.2	-0.5
Latency to leave start box in cocaine runway	0.3	0.3	0.6	-0.8	-0.6	-0.7	-0.2	0.3
Latency to leave start box in cocaine runway, F	1.0	1.1	2.7	-1.6	-1.5	-0.8	-0.0	0.4
Latency to leave start box in cocaine runway, M	0.5	0.5	1.2	1.1	1.1	0.1	-0.0	-0.2
Std. dev. time between licks in bursts	0.4	0.6	1.5	-0.2	-0.1	-0.6	1.2	-1.1
Total sessions with >9 infusions	2.1	2.5	5.3	-2.3	-2.2	-1.4	0.5	-0.2
Short access day 1 total inactive lever presses	1.7	2.2	4.3	-2.1	-2.1	-0.8	-0.8	1.1
Short access day 10 total inactive lever presses	0.1	0.4	0.6	0.4	0.4	-0.6	-0.8	0.7
Short access day 1 total infusions	1.9	2.5	4.9	-2.2	-2.2	-1.6	0.1	0.1
Short access day 10 total infusions	2.4	3.1	6.0	-2.4	-2.4	-1.6	0.9	-0.6
Short access total infusions	2.9	3.7	6.7	-2.6	-2.5	-1.8	1.2	-0.8
Short access day 1 locomotion	0.1	0.1	0.4	-0.2	-0.1	-0.3	-0.6	0.5
Short access day 10 total locomotion	0.0	0.0	0.0	-0.2	-0.2	0.1	0.1	-0.1
Short access total locomotion	0.1	0.1	0.3	-0.6	-0.6	0.1	-0.1	0.1
One hour access (0.3 mA shock)	1.7	2.0	3.6	-1.9	-1.9	-0.9	1.1	-0.9
Compulsive drug intake	1.2	1.4	3.0	1.7	1.7	0.9	-0.7	0.3
One hour access (shock baseline)	2.5	3.8	6.5	-2.5	-2.5	-1.9	1.4	-1.0
Social response proportion	0.6	0.7	1.4	0.4	0.1	0.8	-1.1	1.2
Social responses	1.2	1.5	3.3	0.7	0.4	0.7	-1.7	1.8
Social time	0.7	0.9	2.2	-1.4	-1.5	-0.6	0.3	0.1
Soleus weight	0.9	2.7	5.3	2.3	2.2	1.6	-0.7	0.3
Bone: stiffness	0.3	0.5	1.2	0.9	0.8	1.1	-0.4	0.2
Tibialis anterior weight	7.2	24.5	47.4	6.9	6.9	5.0	0.4	-1.7
Time to tail flick, vehicle, before self-admin	0.6	0.6	1.5	-0.3	-0.4	-1.2	0.9	-0.8
Time to tail flick, vehicle, after self-admin	0.8	1.0	1.4	-1.2	-1.2	-1.1	0.8	-0.6
Time to tail flick, test, before self-admin	0.0	0.0	0.1	0.0	-0.0	0.1	-0.3	0.2
Time to tail flick, test, after self-admin	1.1	1.1	2.7	-1.6	-1.6	-0.2	-0.1	0.5
Delta time to tail flick, vehicle, before/after SA	1.3	1.3	2.4	-1.5	-1.5	-1.1	0.5	-0.3
Delta time to tail flick, test, before/after SA	0.2	0.2	0.4	-0.3	-0.3	0.4	-0.5	0.6
Tibia length	0.3	0.9	1.6	-1.2	-1.3	-0.6	-0.5	0.8
Ambulatory time before self-admin	1.4	2.1	4.2	2.1	2.0	1.3	-0.8	0.4
Ambulatory time after self-admin	0.2	0.2	0.4	0.6	0.6	0.3	0.4	-0.5
Bone: tissue strength	1.6	1.6	3.3	-1.5	-1.5	-1.8	-0.0	0.4
Total heroin consumption	0.3	0.4	0.8	-0.2	-0.2	0.4	0.9	-0.8
Bone: trabecular number	0.3	0.5	1.1	-0.6	-0.5	-1.0	0.6	-0.5
Bone: trabecular spacing	0.1	0.1	0.2	0.5	0.5	0.3	-0.1	0.0
Bone: trabecular thickness	0.7	0.9	1.4	-0.8	-0.8	-1.2	1.0	-0.8
Bone: trabecular tissue density	0.5	0.8	1.9	0.4	0.4	0.2	-1.3	1.4
Water consumed over 24 hours	0.4	0.6	1.0	0.9	0.8	1.0	-0.5	0.4
Number of responses in last shaping day	0.2	0.2	0.4	-0.6	-0.5	-0.5	0.1	-0.1