chr1swallace
commented
2 years ago I don't think the datasets are really fine if you are getting a missing value here. coloc calculates the p value from the effect estimate (beta) and varbeta supplied. Is varbeta missing or 0 anywhere?
From: juliandwillett @.> Sent: 02 June 2022 02:35 To: chr1swallace/coloc @.> Cc: Subscribed @.***> Subject: [chr1swallace/coloc] Error in if (min(p) > warn.minp) warning("minimum p value is: ", format.pval(min(p)), (Issue #91)
I have been getting the following error:
Error in if (min(p) > warn.minp) warning("minimum p value is: ", format.pval(min(p)), : missing value where TRUE/FALSE needed
The datasets are fine. The only thing of note is that I have a value less than 2.2e-16 in one dataset so "format.pval(min(p))" returns <2.2e-16 that I am guessing is character vs numeric if that makes any difference.
— Reply to this email directly, view it on GitHubhttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fchr1swallace%2Fcoloc%2Fissues%2F91&data=05%7C01%7Ccew54%40universityofcambridgecloud.onmicrosoft.com%7C158c3f7b78fc4abe3a3d08da443842a7%7C49a50445bdfa4b79ade3547b4f3986e9%7C0%7C0%7C637897305666594668%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=ovF1Rw7Mf7QI%2Bz1IhoUE%2BxGWnL0Ze%2Bpouf7va3MKXXo%3D&reserved=0, or unsubscribehttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FAAQWR2DLUYI65JC74TXO3J3VNAFWJANCNFSM5XTLXZAA&data=05%7C01%7Ccew54%40universityofcambridgecloud.onmicrosoft.com%7C158c3f7b78fc4abe3a3d08da443842a7%7C49a50445bdfa4b79ade3547b4f3986e9%7C0%7C0%7C637897305666594668%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=tN%2BIpF%2FV2L7DbKndDcFePr%2Bbln5y7fmdcJyZsRzYE24%3D&reserved=0. You are receiving this because you are subscribed to this thread.Message ID: @.***>
I have been getting the following error:
The datasets are fine. The only thing of note is that I have a value less than 2.2e-16 in one dataset so "format.pval(min(p))" returns <2.2e-16 that I am guessing is character vs numeric if that makes any difference.