-
Notifications
You must be signed in to change notification settings - Fork 4
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Error when calling callpeak function #13
Comments
I think may have been due to issues with Conda so I tried this again on a virtual machine running Windows 10:
Package Planenvironment location: C:\Users\yang.p\AppData\Local\basilisk\1.3.6\MACSr-0.99.9\env_macs added / updated specs: The following NEW packages will be INSTALLED: ca-certificates conda-forge/win-64::ca-certificates-2020.12.5-h5b45459_0 Preparing transaction: ...working... done ==> WARNING: A newer version of conda exists. <== Please update conda by running
To activate this environment, use$ conda activate C:\Users\yang.p\AppData\Local\basilisk\1.3.6\MACSr-0.99.9\env_macsTo deactivate an active environment, use$ conda deactivateCollecting package metadata (current_repodata.json): ...working... done Package Planenvironment location: C:\Users\yang.p\AppData\Local\basilisk\1.3.6\MACSr-0.99.9\env_macs added / updated specs: The following packages will be UPDATED: ca-certificates conda-forge::ca-certificates-2020.12.~ --> pkgs/main::ca-certificates-2021.1.19-haa95532_0 The following packages will be SUPERSEDED by a higher-priority channel: certifi conda-forge::certifi-2020.12.5-py37h0~ --> pkgs/main::certifi-2020.12.5-py37haa95532_0 Preparing transaction: ...working... done ==> WARNING: A newer version of conda exists. <== Please update conda by running
Collecting package metadata (current_repodata.json): ...working... done Package Planenvironment location: C:\Users\yang.p\AppData\Local\basilisk\1.3.6\MACSr-0.99.9\env_macs added / updated specs: The following NEW packages will be INSTALLED: intel-openmp conda-forge/win-64::intel-openmp-2020.3-h57928b3_311 The following packages will be UPDATED: certifi pkgs/main::certifi-2020.12.5-py37haa9~ --> conda-forge::certifi-2020.12.5-py37h03978a9_1 The following packages will be SUPERSEDED by a higher-priority channel: ca-certificates pkgs/main::ca-certificates-2021.1.19-~ --> conda-forge::ca-certificates-2020.12.5-h5b45459_0 Preparing transaction: ...working... done ==> WARNING: A newer version of conda exists. <== Please update conda by running
Collecting macs3==3.0.0a5
cwd: C:\Users\yang.p\AppData\Local\Temp\pip-install-o__ps3v3\macs3_1ffa8712b533490d9b9dff5335923622 Complete output (28 lines): Traceback (most recent call last): File "C:\Users\yang.p\AppData\Local\basilisk\13C0E0 main() File "C:\Users\yang.p\AppData\Local\basilisk\13C0E0 json_out['return_val'] = hook(**hook_input['kwargs']) File "C:\Users\yang.p\AppData\Local\basilisk\13C0E0 return hook(config_settings) File "C:\Users\yang.p\AppData\Local\Temp\pip-build-env-5sm3g9ux\overlay\Lib\site-packages\setuptools\build_meta.py", line 150, in get_requires_for_build_wheel config_settings, requirements=['wheel']) File "C:\Users\yang.p\AppData\Local\Temp\pip-build-env-5sm3g9ux\overlay\Lib\site-packages\setuptools\build_meta.py", line 130, in _get_build_requires self.run_setup() File "C:\Users\yang.p\AppData\Local\Temp\pip-build-env-5sm3g9ux\overlay\Lib\site-packages\setuptools\build_meta.py", line 254, in run_setup self).run_setup(setup_script=setup_script) File "C:\Users\yang.p\AppData\Local\Temp\pip-build-env-5sm3g9ux\overlay\Lib\site-packages\setuptools\build_meta.py", line 145, in run_setup exec(compile(code, file, 'exec'), locals()) File "setup.py", line 148, in main() File "setup.py", line 70, in main gcc_version_check = subprocess.check_output( ["gcc", "--version"], universal_newlines=True) File "C:\Users\yang.p\AppData\Local\basilisk\13C0E0 **kwargs).stdout File "C:\Users\yang.p\AppData\Local\basilisk\13C0E0 with Popen(*popenargs, **kwargs) as process: File "C:\Users\yang.p\AppData\Local\basilisk\13C0E0 restore_signals, start_new_session) File "C:\Users\yang.p\AppData\Local\basilisk\13C0E0 |
Hi @philyang07 , the MACSr is currently not compatible with the Windows system, because the MACS3 can't be installed in Windows. Please try to use it with the WSL in windows. |
Hi @hubentu, R version 4.1.0 is required, but I can't seem to find it for Linux. It's a new development version? |
Hi @philyang07 , yes, the tool is built base on the R development version. |
Hi @philyang07 , can you try the example in the manual first? eh <- ExperimentHub::ExperimentHub()
CHIP <- eh[["EH4558"]]
CTRL <- eh[["EH4563"]]
callpeak(CHIP, CTRL, gsize = 5.2e7, cutoff_analysis = TRUE,
outdir = tempdir(), name = "callpeak_narrow0") Also, you can try to run |
Hi @hubentu, When I try it with R version 4.0.3, the following error occurs: When I try it with the R development version, the following error occurs: Yeah, MACS2 works fine. |
Hi @philyang07 , you don't need to install the meta package. The bed file ( The same test datasets from MACS/test were used to run tests in the R package. The bed file can also be downloaded from the original repository CTCF_SE_ChIP_chr22_50k.bed.gz. |
Hi @hubentu, yep it works. I realised I just needed to specify gsize when calling it on the bam file that I provided. |
Hi @hubentu, When adjusting the effective genome size, it seems like the error below occurs for larger values e.g. 2.7e9 (the default value of gsize for MACS3).
However, the error doesn't occur when using callpeaks from MACS2 directly (which uses 2.7e9 as the default value). |
Hi @philyang07 , I couldn't reproduce your error since I don't have your bam file. Could you share your data?Thx. |
Hi @hubentu, Yep, I'll attach it right here. |
Hi @hubentu, I'm facing the same problem that @philyang07 was previously facing regarding the development version of R. I currently have R v4.0.3 and do not want to upgrade to the development version. How would I change the version in the DESCRIPTION file? I tried looking this up online but was not able to find a solution. Please let me know thanks! |
Hi @krovi137 , you can clone the repo to local, then modify the line in the DESCRIPTION. |
Thank you this worked! |
Hi @hubentu, I raised the issue with @taoliu and he figured out the problem (it was with MACS3): macs3-project/MACS#442 |
@philyang07 I will update the dependency when the new version release. |
@hubentu v3.0.0a6 with fixes for this issue has been released to PyPI: https://pypi.org/project/MACS3/ |
@philyang07 Regarding the paired-peak model, I will explain more in the page macs3-project/MACS#442 |
@philyang07 , MACS has been updated to 3.0.0a6. Please have a try. Thanks.
|
Hi @hubentu, thanks for updating, it works fine now. |
Hi,
I tried calling callpeak on a BAM file:
I'm using the R development version (4.1.0) and on Windows 10.
The text was updated successfully, but these errors were encountered: