# Export for linreg\_multi\_pheno or linreg3

**URL:** https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296
**Category:** Help \[0.1\]
**Created:** [September 25, 2017, 9:14pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296 "2017-09-25T21:14:23Z")
**Posts on this page:** 10
**Page:** 1

<div class="post-metadata">

### Author: ![jobaba](https://avatars.discourse-cdn.com/v4/letter/j/dbc845/32.png) [@jobaba](https://discuss.hail.is/u/jobaba)
#### Post date: [September 25, 2017, 9:14pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/1 "2017-09-25T21:14:23Z")

</div>

Hi,

We are trying to do a small scale phewas and we tried linreg\_multi\_pheno. I was wondering if there’s an export function for it? If not what would you suggest to do if we want to look at the results? convert it to pandas?

Thanks,  
Wendy

---

<div class="post-metadata">

### Author: ![tpoterba](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/tpoterba/32/61_2.png) [@tpoterba](https://discuss.hail.is/u/tpoterba)
#### Post date: [September 26, 2017, 1:33am UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/2 "2017-09-26T01:33:58Z")

</div>

Hi Wendy,  
You can certainly export to a file. Here’s how:

```python
vds.export_variants('file.txt.bgz', 
                    'variant = v, pvals = va.linreg.pval.mkString(",")')

```

This’ll give you something like:

```nohighlight
variant pval
1:1:A:T 0.05,0.001,0.9
1:2:T:G 0.1,0.4,0.3
...

```

I’ll think about a better way to get results out and post here when I have something. Someone else has also asked for this.

---

<div class="post-metadata">

### Author: ![jobaba](https://avatars.discourse-cdn.com/v4/letter/j/dbc845/32.png) [@jobaba](https://discuss.hail.is/u/jobaba)
#### Post date: [September 26, 2017, 2:15am UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/3 "2017-09-26T02:15:57Z")

</div>

Thanks Tim! This would work for now!

Best,  
Wendy

---

<div class="post-metadata">

### Author: ![jobaba](https://avatars.discourse-cdn.com/v4/letter/j/dbc845/32.png) [@jobaba](https://discuss.hail.is/u/jobaba)
#### Post date: [September 27, 2017, 7:38pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/4 "2017-09-27T19:38:30Z")

</div>

Hi Tim,

We have got  
Error summary: HailException: No function found with name `mkString’ and argument types (Array[Double], String)  
:1:variant = v, pvals = va.linreg.pval.mkString(“,”)

is the mkString function in Hail 0.1-38882df?

Thanks,  
Wendy

---

<div class="post-metadata">

### Author: ![jbloom](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/jbloom/32/109_2.png) [@jbloom](https://discuss.hail.is/u/jbloom)
#### Post date: [September 27, 2017, 7:56pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/5 "2017-09-27T19:56:46Z")

</div>

Hi Wendy, pval is an Array[Double] but mkString takes an Array[String], so try substituting:

`va.linreg.pval.map(p => str(p)).mkString(",")`

---

<div class="post-metadata">

### Author: ![jerryliu2005](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/jerryliu2005/32/97_2.png) [@jerryliu2005](https://discuss.hail.is/u/jerryliu2005)
#### Post date: [September 28, 2017, 4:54pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/6 "2017-09-28T16:54:32Z")

</div>

Could you help me here with the syntax of HQL to filter the linreg vds with all phenotype pval \< 0.05?

Thanks!  
Jerry

---

<div class="post-metadata">

### Author: ![tpoterba](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/tpoterba/32/61_2.png) [@tpoterba](https://discuss.hail.is/u/tpoterba)
#### Post date: [September 28, 2017, 5:07pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/7 "2017-09-28T17:07:57Z")

</div>

something like

```python
vds.filter_variants_expr('va.linreg.pval.forall(p => p < 0.05)')

```

---

<div class="post-metadata">

### Author: ![jerryliu2005](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/jerryliu2005/32/97_2.png) [@jerryliu2005](https://discuss.hail.is/u/jerryliu2005)
#### Post date: [September 28, 2017, 5:31pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/8 "2017-09-28T17:31:15Z")

</div>

Thanks! Just to beat the dead horse again:  
If I want to filter by the first pheno p-val, is it something like this?  
vds.filter\_variants\_expr(‘va.linreg.pval.forall(p =\> p[0]\< 0.05)’)

---

<div class="post-metadata">

### Author: ![tpoterba](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/tpoterba/32/61_2.png) [@tpoterba](https://discuss.hail.is/u/tpoterba)
#### Post date: [September 28, 2017, 5:33pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/9 "2017-09-28T17:33:16Z")

</div>

the `forall` checks that the predicate is true for every element of the array. If you just care about the first, you don’t need the `forall`:

```python
 vds.filter_variants_expr('va.linreg.pval[0] < 0.05')

```

---

<div class="post-metadata">

### Author: ![jerryliu2005](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/jerryliu2005/32/97_2.png) [@jerryliu2005](https://discuss.hail.is/u/jerryliu2005)
#### Post date: [September 28, 2017, 5:37pm UTC](https://discuss.hail.is/t/export-for-linreg-multi-pheno-or-linreg3/296/10 "2017-09-28T17:37:03Z")

</div>

> [@tpoterba](#):
>
> vds.filter\_variants\_expr(‘va.linreg.pval[0] \< 0.05’)

Got it. Thanks for your quick reply!
