본문으로 건너뛰기
cccskills
무료GitHub에 공개

bioconductor-miasim

Microbiome time series simulation with generalized Lotka-Volterra model, Self-Organized Instability (SOI), and other models. Hubbell's Neutral model is used to determine the abundance matrix. The resulting abundance matrix is applied to (Tree)SummarizedExperiment objects.

한국어 요약은 준비 중입니다. 원문 설명을 표시하고 있습니다.

실행 점수 미검증

설치 방법 보기

포함된 파일(1)

  • SKILL.md3.7 KB

SKILL.md(원문)

설치하기 전에 에이전트에게 전달되는 지시 내용을 확인할 수 있습니다.

miaSim

Workflows

Standard Workflow

library(miaSim)

# 1. Generate species interaction matrix
A_normal <- powerlawA(n_species = 4, alpha = 3)

# 2. Simulate time series using generalized Lotka-Volterra (gLV) model
tse_glv <- simulateGLV(
  n_species = 4, 
  A = A_normal, 
  t_start = 0, 
  t_store = 1000, 
  stochastic = FALSE, 
  norm = FALSE
)

# 3. Simulate community using Hubbell Neutral model
tse_hubbell <- simulateHubbell(
  n_species = 8, 
  M = 10, 
  carrying_capacity = 1000, 
  k_events = 50, 
  migration_p = 0.02, 
  t_end = 100
)

Input/Output Note: Inputs are ecological parameters (species count, interaction matrices, carrying capacity); outputs are simulated abundance matrices stored in TreeSummarizedExperiment objects.

When to Use

  • Simulating microbiome species abundance matrices and time series under varying ecological assumptions.
  • Modeling population dynamics using the generalized Lotka-Volterra model (simulateGLV) or its discrete version, the Ricker model (simulateRicker).
  • Simulating neutral community dynamics with migration, births, and deaths using Hubbell's Neutral model (simulateHubbell).
  • Generating species interaction matrices with specific distributional assumptions using powerlawA or randomA.

When NOT to Use

  • For fitting ecological models to real-world experimental data (use parameter estimation packages instead).
  • For simulating non-microbial networks that do not fit into the TreeSummarizedExperiment data container.

Data Requirements

  • Input parameters such as interaction matrices generated by powerlawA or randomA.
  • Outputs are structured as TreeSummarizedExperiment objects containing simulated abundance matrices.

Key Parameters

  • n_species: Number of species to include in the simulation or interaction matrix.
  • alpha: Power-law exponent parameter used in powerlawA.
  • diagonal: Diagonal values of the interaction matrix in randomA.
  • connectance: Connectance of the interaction matrix in randomA.
  • carrying_capacity: Carrying capacity of the community in simulateHubbell and simulateSOI.
  • migration_p: Migration probability in simulateHubbell and simulateHubbellRates.
  • t_end: End time of the simulation.
  • norm: Logical indicating whether to normalize the simulated abundances.

Best Practices

  • Ensure the number of species specified in the interaction matrix matches the number of species used in the simulation models.
  • Generate interaction matrices using powerlawA or randomA before running models like GLV, SOI, or Ricker.
  • Use randomE to generate production and consumption rate matrices specifically for the consumer-resource model (simulateConsumerResource).

Common Pitfalls

  • Specifying a different number of species in the interaction matrix than in the simulation model, which causes dimension mismatch errors. Ensure n_species matches the dimensions of A.

Alternatives

  • micodymora: A Python package for microbiome simulation.

Citations

  • Gao et al. 2023 (Case studies on microbiome simulation).

References

리뷰

아직 리뷰가 없습니다. 사용해 본 소감을 남겨 주세요.

댓글

로그인 없이 쓸 수 있습니다. 사용 후기, 질문, 작동하지 않은 환경 등을 편하게 남겨 주세요.

아직 댓글이 없습니다. 첫 댓글을 남겨 보세요.

보내면 바로 공개됩니다. 이메일 주소나 전화번호 같은 개인 정보는 쓰지 마세요. URL은 2개까지입니다.

0 / 1,000자

이 양식은 reCAPTCHA로 보호되며 Google 개인정보처리방침 및 서비스 약관이 적용됩니다.

같은 저장소의 스킬

The package contains functions for exploratory oligonucleotide array analysis. The dependence on tkWidgets only concerns few convenience functions. 'affy' is fully functional without it.

한국어 요약은 준비 중입니다. 원문 설명을 표시하고 있습니다.

bioMate-AI/biomate-bioconductor-kb7932026년 6월 21일 업데이트

This tool uses the label-tree function from HyPhy to annotate a phylogenetic tree. It allows users to select a subset of leaves using either a regular expression or a list of sequence names, and then apply a specified label to these selecte

한국어 요약은 준비 중입니다. 원문 설명을 표시하고 있습니다.

bioMate-AI/biomate-bioconductor-kb7932026년 6월 21일 업데이트

This package provides a client for the Bioconductor AnnotationHub web resource. The AnnotationHub web resource provides a central location where genomic files (e.g., VCF, bed, wig) and other resources from standard locations (e.g., UCSC, En

한국어 요약은 준비 중입니다. 원문 설명을 표시하고 있습니다.

bioMate-AI/biomate-bioconductor-kb7932026년 6월 21일 업데이트

apeglm provides Bayesian shrinkage estimators for effect sizes for a variety of GLM models, using approximation of the posterior for individual coefficients.

한국어 요약은 준비 중입니다. 원문 설명을 표시하고 있습니다.

bioMate-AI/biomate-bioconductor-kb7932026년 6월 21일 업데이트

ASURAT is a software for single-cell data analysis. Using ASURAT, one can simultaneously perform unsupervised clustering and biological interpretation in terms of cell type, disease, biological process, and signaling pathway activity. Inputting a single-cell RNA-seq data and knowledge-based databases, such as Cell Ontology, Gene Ontology, KEGG, etc., ASURAT transforms gene expression tables into original multivariate tables, termed sign-by-sample matrices (SSMs).

한국어 요약은 준비 중입니다. 원문 설명을 표시하고 있습니다.

bioMate-AI/biomate-bioconductor-kb7932026년 6월 21일 업데이트

AUCell allows to identify cells with active gene sets (e.g. signatures, gene modules...) in single-cell RNA-seq data. AUCell uses the "Area Under the Curve" (AUC) to calculate whether a critical subset of the input gene set is enriched with

한국어 요약은 준비 중입니다. 원문 설명을 표시하고 있습니다.

bioMate-AI/biomate-bioconductor-kb7932026년 6월 21일 업데이트

bioMate-AI의 스킬 모두 보기

이 스킬의 문제 신고