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