Skip to content
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

add spacexr #6212

Draft
wants to merge 46 commits into
base: main
Choose a base branch
from
Draft

add spacexr #6212

wants to merge 46 commits into from

Conversation

nilchia
Copy link
Contributor

@nilchia nilchia commented Aug 6, 2024

FOR CONTRIBUTOR:

  • I have read the CONTRIBUTING.md document and this tool is appropriate for the tools-iuc repo.
  • License permits unrestricted use (educational + commercial)
  • This PR adds a new tool or tool collection
  • This PR updates an existing tool or tool collection
  • This PR does something else (explain below)

This PR adds spacexr to Galaxy. This tool is used for cell-type deconvolution and cell type-specific differential expression in spatial transcriptomics.

tools/spacexr/macros.xml Outdated Show resolved Hide resolved
@nilchia
Copy link
Contributor Author

nilchia commented Oct 9, 2024

Is there another way other than "compare=sim_size" to check ".pdf" and ".rds" outputs?

<param argument="doublet_mode" type="boolean" truevalue="TRUE" falsevalue="FALSE" checked="true" label="Use RCTD doublet mode weights?" help="Otherwise, uses RCTD full mode weights." />
<param argument="weight_threshold" type="float" min="0" value="" optional="true" label="Weight threshold" help="The threshold of total normalized weights across all cell types in 'cell types' per pixel to be included in the model."/>
<param argument="sigma_gene" type="boolean" truevalue="TRUE" falsevalue="FALSE" checked="true" label="Fit gene specific overdispersion parameter?" help="If FALSE, overdispersion parameter is same across all genes." />
<param argument="PRECISION_THRESHOLD" type="float" min="0" value="0.05" label="Precision threshold" help="For checking for convergence, the maximum parameter change per algorithm step."/>
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Why is this all caps?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Just because it is all caps in the function :)
I'll change it.

test_genes_sig = $type.test_genes_sig,
logs = $type.logs
</token>
<token name="CSIDE_SINGLE_RUN">
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You can use CDATA here and than avoid $lt;

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

THANKS!! :)

#end if
]]></command>
<configfiles>
<configfile name="cside_script">
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

use CDATA here

myRCTD &lt;- readRDS('inputs/rctd.rds')

# set core
myRCTD@config$max_cores &lt;- \${GALAXY_SLOTS:-2}
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

GALAXY_SLOTS will only work in bash. Here you are in R you need to access the environment variable GALAXY_SLOTS via an R function.

The tests should fail here, if they are not - why are they not failing :)

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

There is a Cheetah syntax error before starting the job. I am trying to find it.
I'm offending it somehow :)

<expand macro="cside_common_input"/>
</when>
<when value="point_density">
<param argument="pathologic_barcode" type="text" optional="false" label="Barcodes" help="Comma separated barcodes of the pathological region."/>
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

some validator here?

<expand macro="cside_single_input"/>
</when>
<when value="custom">
<param argument="pathologic_barcode" type="text" optional="false" label="Barcodes" help="Comma separated barcodes of the pathological region."/>
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

should this be a macro? its used a lot

<param argument="log_FC_thresh" type="float" min="0" value="0.4" label="LogFC threshold" help="The natural log fold change cutoff for differential expression."/>
</when>
</conditional>
<section name="output" title="Output options">
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

you can simply use a select box here with multiple=true

<param name="UMI_max" type="integer" min="0" value="20000000" label="UMI_max" help="Maximum UMI per pixel included in the analysis" />
<param name="counts_MIN" type="integer" min="0" value="10" label="counts_MIN" help="Minimum total counts per pixel of genes used in the analysis." />
<param name="UMI_min_sigma" type="integer" min="0" value="300" label="UMI_min_sigma" help="Minimum UMI per pixel for the choose_sigma_c function" />
<param name="CELL_MIN_INSTANCE" type="integer" min="0" value="25" label="CELL_MIN_INSTANCE" help="Minimum number of cells required per cell type." />
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Its only a recommendation, but to stay consistent I would name all variables small letters.

<param name="DOUBLET_THRESHOLD" type="integer" min="0" value="20" label="DOUBLET_THRESHOLD" help="The penalty weight of predicting a doublet instead of a singlet for a pixel." />
</section>
</section>
<section name="output" title="Output options">
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

simple select box

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants