Merged
Conversation
Added detailed Roxygen documentation to all major functions in buildmxPedigrees.R. Improved flexibility of family group model construction by allowing relatedness matrices to be optional and handling missing matrices. Added fitPedigreeModel function to fit OpenMx pedigree models. Updated vignettes to introduce fitting pedigree models and made minor formatting improvements.
Update VignetteIndexEntry metadata in three vignette Rmd files to more descriptive titles for documentation indexing and display: vignettes/v0_network.Rmd (Network -> "Network tools for finding extended pedigrees and path tracing"), vignettes/v1_modelingvariancecomponents.Rmd (modelingvariancecomponents -> "Modeling variance components"), and vignettes/v2_pedigree.Rmd (Pedigree -> "Pedigree Simulation and Visualization"). This improves clarity and searchability of package vignettes.
Replace embedded pedigree images in vignettes/v0_network.html (updated base64 PNGs), refresh run metadata timestamps and wall-clock times in vignettes/v1_modelingvariancecomponents.html, and modify vignettes/v5_ASOIAF.Rmd (adjust heading level and add a relatedness-matrix plotting snippet). These changes refresh figures, update generated metadata, and add a visualization example to the ASOIAF vignette.
…#114) * Optimize pedigree simulator with vectorized parent selection Implemented significant performance optimizations for simulatePedigree(): Key improvements: - Vectorized parent selection in buildBetweenGenerations_optimized: Replaced O(n²) loop with linear search with O(n) vectorized operations using couple keys and batch marking - Reduced random permutations from 2 to 1 per generation - Better use of pre-computed row indices to avoid repeated subsetting Performance gains: - Small pedigrees (Ngen=4): 1.5-2x speedup - Medium pedigrees (Ngen=5-6): 3-5x speedup - Large pedigrees (Ngen=7+): 5-10x speedup Usage: Set beta=TRUE or beta="optimized" to use optimized version. Default behavior (beta=FALSE) unchanged for backward compatibility. Added: - OPTIMIZATION_NOTES.md: Detailed documentation of optimizations - benchmark_simulator.R: Performance testing script https://claude.ai/code/session_01NUzTTgoeMd3hTeqvLnrXgB * Restore detailed documentation to optimized version Properly copied all comment blocks from base version to optimized version, maintaining the detailed step-by-step explanations while adding clear notes about what was optimized and why. https://claude.ai/code/session_01NUzTTgoeMd3hTeqvLnrXgB * Fix syntax error: remove duplicate code after optimized function Removed duplicate code blocks (lines 838-940) that were accidentally left after the buildBetweenGenerations_optimized function, causing a compilation error with unmatched braces. https://claude.ai/code/session_01NUzTTgoeMd3hTeqvLnrXgB * Fix optimization to match base version's random behavior Changed approach to maintain exact compatibility with base version while still optimizing the O(n²) bottleneck: - Keep the same loop structure and stopping condition as base - Add back randomization of parent generation before selection - Replace O(n) linear spouse search with O(1) hash lookup - Use named vector for id-to-row mapping This ensures: - Tests pass (same random number sequence as base version) - Still O(n) instead of O(n²) (via O(1) spouse lookup) - 4-5x speedup maintained for large pedigrees https://claude.ai/code/session_01NUzTTgoeMd3hTeqvLnrXgB * Restore full vectorization for maximum performance Replaced loop-based approach with fully vectorized operations: - Use !duplicated() to find unique couples in order (no loop) - Vectorized couple key creation and parent marking - Only one sample() call (matches base version's random sequence) - Calculate target couples to match marR proportion This restores the 4-5x speedup for large pedigrees while maintaining similar random behavior to base version. https://claude.ai/code/session_01NUzTTgoeMd3hTeqvLnrXgB * Document trade-off between speed and exact reproducibility Comprehensive documentation updates: 1. Enhanced @param beta documentation explaining: - Performance vs reproducibility trade-off - When to use each version - Both versions are mathematically correct 2. Updated OPTIMIZATION_NOTES.md with: - Clear comparison table - Detailed explanation of differences - Guidance on when to use each version - Testing implications 3. Added TEST_RECOMMENDATIONS.md with: - Three different test strategies (A, B, C) - Specific examples of updated tests - Recommendations for statistical ranges - List of tests that need updating Key message: Optimized version (beta=TRUE) is 4-5x faster and mathematically correct, but produces different (not wrong) random outcomes than base version. Use beta=FALSE for exact reproducibility, beta=TRUE for speed. https://claude.ai/code/session_01NUzTTgoeMd3hTeqvLnrXgB * Update tests to handle both beta=FALSE and beta=TRUE Modified 4 tests to use appropriate expectations for each version: - beta=FALSE (base): Exact counts and strict tolerances - beta=TRUE (optimized): Ranges (±20%) and wider tolerances (±0.05) Changes: 1. Test 1: ID count 57 → 45-70 range for beta=TRUE 2. Test 2: ID count 154 → 123-185 range for beta=TRUE 3. Test 3: ID count 424 → 340-510 range for beta=TRUE 4. Test 4: ID count 57 → 45-70 range for beta=TRUE Sex ratio tolerance widened from ±0.03 to ±0.05 for beta=TRUE to account for statistical variation in different random sequences. Both versions are mathematically correct; optimized version just uses different random number sequence, producing equivalent results. https://claude.ai/code/session_01NUzTTgoeMd3hTeqvLnrXgB * Update test-simulatePedigree.R * Format R code and tests (whitespace only) Apply whitespace and style fixes across multiple R files and tests. Adjusted multi-line function call formatting (checkIDs, checkParents, helpChecks), normalized if/brace spacing and function signature indentation (simulatePedigree), and removed stray blank lines and tightened parentheses in test expectations. These are formatting-only changes intended to improve readability; no functional behavior changes are expected. --------- Co-authored-by: Claude <noreply@anthropic.com>
Codecov Report❌ Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## main #115 +/- ##
==========================================
+ Coverage 84.32% 84.39% +0.06%
==========================================
Files 28 30 +2
Lines 4281 4979 +698
==========================================
+ Hits 3610 4202 +592
- Misses 671 777 +106 ☔ View full report in Codecov by Sentry. 🚀 New features to boost your workflow:
|
Introduce twinID parameter across pedigree segmentation functions (ped2fam/.ped2id/ped2graph/ped2maternal/ped2paternal) and thread it through calls so twin IDs are considered when building family/graph structures. Add mz_twins logical option to ped2com (default FALSE) and, when TRUE and twinID is present, call addMZtwins(ped) to treat MZ twins as an additional parent–child relationship for relatedness computations. Also fix a minor typo in a comment.
* Fix MZ twins coded as 0.5 instead of 1.0 in relatedness matrix Implement addMZtwins() which redirects one MZ co-twin's parent links to point to the other twin before adjacency matrix construction. This produces isPar[twin2, twin1] = 1.0 in the sparse matrix (two 0.5 entries summed), so path tracing yields relatedness 1 between MZ pairs. Users provide a twinID column (and optionally zygosity) and pass mz_twins=TRUE to ped2add()/ped2com(). DZ twins are left unchanged when zygosity column is present. https://claude.ai/code/session_01P3RQTYpWtAtheSqi4aPjR5 * Replace pedigree redirection with r2 column merge for MZ twins The previous approach redirected twin2's parents to twin1 in the pedigree, which inflated twin2's diagonal (1.5 instead of 1.0) and twin2-to-child relatedness (0.75 instead of 0.5). New approach: after path tracing but before tcrossprod, merge twin2's column into twin1's in the r2 matrix. MZ twins share the same genetic source, so this correctly produces: - MZ twin relatedness = 1.0 - Self-relatedness = 1.0 (no inflation) - Parent-child and all downstream values correct - No post-hoc diagonal patching needed https://claude.ai/code/session_01P3RQTYpWtAtheSqi4aPjR5 * Use symmetric column merge for MZ twins instead of zeroing Both twin columns now get the same normalized values (r2_merged = (col1 + col2) / sqrt(2)) so both twins remain present and contribute equally. Produces the same final relatedness matrix as the zero approach but without erasing one twin from the genetic source matrix. https://claude.ai/code/session_01P3RQTYpWtAtheSqi4aPjR5 * Merge-then-restore approach for MZ twins Temporarily absorb twin2's column into twin1's before tcrossprod, then copy twin1's row/col back to twin2 afterward. This keeps the computation correct while ensuring neither twin is erased from the final relatedness matrix. https://claude.ai/code/session_01P3RQTYpWtAtheSqi4aPjR5 * Accept lowercase 'mz' and add MZ twin tests Treat both "mz" and "MZ" as monozygotic in findMZtwins (zygosity check now uses %in% c("mz","MZ")). Minor formatting tweak to the verbose message. Added unit tests (tests/testthat/test-buildComponent.R) verifying that MZ twins are coded with relatedness 1 when mz_twins=TRUE, that siblings remain 0.5 when mz_twins=FALSE, self-relatedness stays 1, and parent-child relatedness is unchanged. * fix tests --------- Co-authored-by: Claude <noreply@anthropic.com>
Introduce mz_method (default "merge_before_tcrossprod") to ped2com and gate MZ twin column-merge/restore logic behind this option. Restrict the merge/restore behavior to the additive component, adjust verbose messages, and add a TODO outlining an alternative MZ handling approach. Update tests to cover MZ twin child relatedness cases and add clarifying comments for ped2fam string/numeric ID behavior.
Expose a new mz_method argument in ped2add (default 'addtwins') and forward it to ped2com. Update tests to pass mz_method = 'merging' when mz_twins is used so alternative MZ-twin handling is exercised. This lets callers choose the method for handling MZ twins without changing the default behavior.
Add cross-method validations for MZ twin handling. In data-raw/optimizing.R create r_mz1 and r_mz2 using ped2add with mz_method 'merging' and 'addtwins' and assert their sparse matrix internals (@i and @x) match. Update tests/testthat/test-buildComponent.R to loop over mz_method options when verifying parent/child relatedness for MZ twins and add an explicit equality assertion between the two method results. These changes ensure different mz_method implementations produce equivalent relatedness outputs.
Allow mz_method to be either "merging" or "addtwins" in key pedigree-processing branches and only perform parent redirection when mz_method == "merging". Move verbose message into the merging branch and simplify restoration checks for MZ pairs. Add a new test that verifies equivalence of additive matrices for merging vs addtwins (including sparse-matrix slot checks). Also fix an expectation bound in simulatePedigree tests. Files changed: R/buildComponent.R, tests/testthat/test-buildComponent.R, tests/testthat/test-simulatePedigree.R.
Create mz_id_pairs earlier when MZ twin pairs exist and streamline the merging branch to only run for mz_method == "merging". Remove redundant mz_id_pairs construction, keep logic that makes the second twin a founder and redirects its children to the first twin, and preserve the verbose merge message. Add object.size() calls in the test to inspect memory usage of intermediate matrices.
) Two issues caused the sparse matrices from the two mz_method paths to differ structurally: 1. The merging method iterated over mz_pairs (row indices) but treated the values as IDs. Changed to iterate over mz_id_pairs which holds the actual pedigree IDs. Currently masked because test data has sequential IDs matching row positions, but would break with non-sequential IDs. 2. The different computational paths (column-merge-before-tcrossprod vs pedigree-modify-then-post-copy) leave different explicit zeros in the sparse matrix. Added Matrix::drop0() before returning to normalize the sparse representation so both methods produce structurally identical results. https://claude.ai/code/session_01PXbgV7bdA9sg7St6xBbg3T Co-authored-by: Claude <noreply@anthropic.com>
Only restore MZ twin rows/cols in buildComponent when mz_method is merging/addtwins AND the component is 'additive' (buildComponent.R). In constructAdjacency.R adjust the .adjIndexed parameter layout and add a new isTwin helper to detect twins from the pedigree. Update tests to call ped2add without explicit mz_method, comment out object.size checks, and compare sparse adjacency matrices by converting to dense (as.matrix) and asserting their difference sums to zero to avoid incorrect sparse-matrix summation.
* Fix twin matrix size mismatch between addtwins and merging methods Two issues caused the sparse matrices from the two mz_method paths to differ structurally: 1. The merging method iterated over mz_pairs (row indices) but treated the values as IDs. Changed to iterate over mz_id_pairs which holds the actual pedigree IDs. Currently masked because test data has sequential IDs matching row positions, but would break with non-sequential IDs. 2. The different computational paths (column-merge-before-tcrossprod vs pedigree-modify-then-post-copy) leave different explicit zeros in the sparse matrix. Added Matrix::drop0() before returning to normalize the sparse representation so both methods produce structurally identical results. https://claude.ai/code/session_01PXbgV7bdA9sg7St6xBbg3T * Fix sparse matrix class mismatch between MZ twin methods The merging method's post-copy (r[idx2,] <- r[idx1,]) causes Matrix to coerce from dsCMatrix (symmetric, one triangle) to dgCMatrix (general, both triangles), doubling @i/@x slot lengths. Add forceSymmetric() after the post-copy to restore symmetric storage. Also guard the merging pedigree modification with a component check (additive only) to match the restoration block, and move the verbose message inside the merging if-block so it only prints when merging actually occurs. https://claude.ai/code/session_01PXbgV7bdA9sg7St6xBbg3T * Fix test assertion: use sum(A - B) instead of sum(A, B) sum(A, B) sums all elements of both matrices, which can never be 0 for non-negative relatedness values. Use sum(A - B) to check that the element-wise differences are zero. https://claude.ai/code/session_01PXbgV7bdA9sg7St6xBbg3T --------- Co-authored-by: Claude <noreply@anthropic.com>
Allow specifying only one twin or mate ID and auto-find a suitable partner. makeTwins now builds an ID-to-row map for O(1) lookups and supports cases where only ID_twin1 or only ID_twin2 is provided (selects sibling by zygosity/sex). makeInbreeding gains symmetric logic to auto-select ID_mate2/ID_mate1 when only one mate is given. dropLink standardizes pedigree column names (with an optional verbose flag) before processing. Updated tests to cover these single-ID behaviors and adjusted an existing MZ twin test; NEWS updated to mention the tweakPedigree change.
Allow specifying only one twin or mate ID and auto-find a suitable partner. makeTwins now builds an ID-to-row map for O(1) lookups and supports cases where only ID_twin1 or only ID_twin2 is provided (selects sibling by zygosity/sex). makeInbreeding gains symmetric logic to auto-select ID_mate2/ID_mate1 when only one mate is given. dropLink standardizes pedigree column names (with an optional verbose flag) before processing. Updated tests to cover these single-ID behaviors and adjusted an existing MZ twin test; NEWS updated to mention the tweakPedigree change.
two twin methods
Introduce functions to construct and fit OpenMx pedigree models: buildPedigreeModelCovariance, buildOneFamilyGroup, buildFamilyGroups, buildPedigreeMx, and fitPedigreeModel (R/buildmxPedigrees.R). Add corresponding Rd documentation files for each exported function (man/...). Also apply minor code tidyups and formatting fixes across the package: whitespace and call formatting in R/tweakPedigree.R, small cleanups in data-raw scripts (optimizing_simulations.R, optimizing_twins.R), test formatting in tests/testthat/test-buildComponent.R, and update vignette timestamps/benchmarks in vignettes/v1_modelingvariancecomponents.html. The new API requires the OpenMx package to build and run pedigree models.
Co-authored-by: Copilot <175728472+Copilot@users.noreply.github.com>
* Initial plan * Remove HTML vignette build artifacts and update .gitignore Co-authored-by: smasongarrison <6001608+smasongarrison@users.noreply.github.com> --------- Co-authored-by: copilot-swe-agent[bot] <198982749+Copilot@users.noreply.github.com> Co-authored-by: smasongarrison <6001608+smasongarrison@users.noreply.github.com>
Replace require(OpenMx) checks with !requireNamespace("OpenMx", quietly = TRUE) across pedigree functions and keep library(OpenMx) after the check. Affected functions: buildPedigreeModelCovariance, buildOneFamilyGroup, buildFamilyGroups, buildPedigreeMx, and fitPedigreeModel. Also apply cosmetic formatting (argument alignment and mxMatrix/mxData/mxAlgebra indentation) for readability. These are refactors and style changes; no functional logic was altered.
* Initial plan * Add relatedness matrix parameters to fitPedigreeModel Co-authored-by: smasongarrison <6001608+smasongarrison@users.noreply.github.com> --------- Co-authored-by: copilot-swe-agent[bot] <198982749+Copilot@users.noreply.github.com> Co-authored-by: smasongarrison <6001608+smasongarrison@users.noreply.github.com> Co-authored-by: Mason Garrison <garrissm@wfu.edu>
Pedigree helpers
Contributor
There was a problem hiding this comment.
Pull request overview
This PR extends BGmisc’s pedigree tooling by (1) adding optimized simulation and more flexible algorithm selection, (2) adding MZ-twin handling to relatedness computations, and (3) introducing OpenMx pedigree model builders with new vignettes, docs, and tests.
Changes:
- Added an optimized/vectorized pedigree simulation path (
betaalgorithm selection) and updated docs/tests accordingly. - Added MZ twin handling options to
ped2com()/ped2add()plus new twin helper utilities (findMZtwins(),fuseTwins()). - Added OpenMx pedigree model builder APIs (
build*+fitPedigreeModel()) with supporting documentation and tests.
Reviewed changes
Copilot reviewed 58 out of 61 changed files in this pull request and generated 7 comments.
Show a summary per file
| File | Description |
|---|---|
| vignettes/v6_pedigree_model_fitting.Rmd | New vignette demonstrating pedigree-based variance component modeling workflow (OpenMx). |
| vignettes/v5_ASOIAF.Rmd | Updates example usage (adds mz_twins) and adds relatedness matrix visualization section. |
| vignettes/v2_pedigree.Rmd | Updates vignette index entry title. |
| vignettes/v1_modelingvariancecomponents.html | Removes generated HTML from repo. |
| vignettes/v1_modelingvariancecomponents.Rmd | Updates vignette index entry + adds “Extended Pedigrees” section. |
| vignettes/v0_network.Rmd | Updates vignette index entry title. |
| vignettes/.gitignore | Ignores generated vignette HTML output. |
| tests/testthat/test-tweakPedigree.R | Expands tests for one-ID twin/inbreeding selection + prefer_unmated + dropLink() assertions. |
| tests/testthat/test-simulatePedigree.R | Loosens tolerances / adds assertions to accommodate optimized simulation differences. |
| tests/testthat/test-segmentPedigree.R | Minor cleanup + clarifying comments. |
| tests/testthat/test-helpTwins.R | New tests for findMZtwins() / fuseTwins() argument combinations. |
| tests/testthat/test-buildmxPedigrees.R | New tests for OpenMx pedigree model builder functions (skipped when OpenMx absent). |
| tests/testthat/test-buildComponent.R | Adds extensive tests for MZ twin handling in additive relatedness. |
| man/simulatePedigree.Rd | Documents beta as logical/character algorithm selector. |
| man/ped2paternal.Rd | Adds twinID parameter to docs. |
| man/ped2maternal.Rd | Adds twinID parameter to docs. |
| man/ped2graph.Rd | Adds twinID parameter to docs. |
| man/ped2fam.Rd | Adds twinID parameter to docs. |
| man/ped2com.Rd | Updates args (partialparent default, MZ twin params, beta) and documents MZ handling. |
| man/ped2add.Rd | Documents mz_twins / mz_method passthrough options. |
| man/makeTwins.Rd | Documents new twin_sex parameter. |
| man/makeInbreeding.Rd | Documents new prefer_unmated parameter. |
| man/isTwin.Rd | New internal doc for isTwin(). |
| man/fuseTwins.Rd | New internal doc for fuseTwins(). |
| man/fitPedigreeModel.Rd | New exported doc for fitPedigreeModel(). |
| man/findMZtwins.Rd | New internal doc for findMZtwins(). |
| man/dropLink.Rd | Documents new verbose parameter. |
| man/dot-adjBeta.Rd | Title text tweak (“Methods”). |
| man/buildWithinGenerations.Rd | Documents enhanced beta algorithm selection semantics. |
| man/buildPedigreeMx.Rd | New exported doc for buildPedigreeMx(). |
| man/buildPedigreeModelCovariance.Rd | New exported doc for buildPedigreeModelCovariance(). |
| man/buildOneFamilyGroup.Rd | New exported doc for buildOneFamilyGroup(). |
| man/buildFamilyGroups.Rd | New exported doc for buildFamilyGroups(). |
| man/buildBetweenGenerations.Rd | Documents enhanced beta algorithm selection semantics. |
| man/adjustKidsPerCouple.Rd | Documents enhanced beta algorithm selection semantics. |
| data-raw/ped2com_benchmark_summary.csv | Adds benchmark output artifact. |
| data-raw/ped2com_benchmark_start_time.txt | Adds benchmark timing artifact. |
| data-raw/ped2com_benchmark_end_time.txt | Adds benchmark timing artifact. |
| data-raw/ped2com_benchmark_design.csv | Adds benchmark design artifact. |
| data-raw/optimizing_twins.R | New benchmarking script for twin handling performance. |
| data-raw/optimizing_simulations.R | New benchmarking script for simulation performance comparisons. |
| data-raw/optimizing.R | Removes older benchmarking script. |
| R/tweakPedigree.R | Adds twin_sex, one-ID auto-selection, prefer_unmated, dropLink(verbose), and makePool() helper. |
| R/simulatePedigree.R | Implements buildBetweenGenerations_optimized() and beta algorithm selection behavior. |
| R/segmentPedigree.R | Threads twinID through pedigree segmentation and graph building. |
| R/helpTwins.R | Adds internal twin helper functions (isTwin, findMZtwins, fuseTwins). |
| R/helpChecks.R | Signature formatting cleanup. |
| R/constructAdjacency.R | Signature formatting + title tweak. |
| R/checkParents.R | Whitespace/formatting cleanup. |
| R/checkIDs.R | Formatting cleanup for dropIdenticalDuplicateIDs() call. |
| R/buildmxPedigrees.R | Adds OpenMx pedigree model builder and fitting functions + roxygen exports. |
| R/buildComponent.R | Adds MZ twin handling hooks to ped2com() / ped2add(), changes defaults, and restores twins post-tcrossprod. |
| NEWS.md | Adds development notes for new features. |
| NAMESPACE | Exports new OpenMx builder functions and fitPedigreeModel(). |
Comment on lines
66
to
+74
| for (i in 1:idx) { | ||
| # cat("loop", i, "\n") | ||
| # check if i is equal to the number of individuals in the generation | ||
| usedID <- c(usedID, ID_twin1) | ||
| # message(usedID) | ||
| if (i < idx) { | ||
| # randomly select one individual from the generation | ||
| ID_twin1 <- resample(ped$ID[ped$gen == gen_twin & !(ped$ID %in% usedID) & !is.na(ped$dadID)], 1) | ||
| # cat("twin1", ID_twin1, "\n") | ||
| ID_twin1 <- resample(ped$ID[gen_mask & !(ped$ID %in% usedID)], 1) | ||
| # Cache twin1 attributes via O(1) row lookup | ||
| twin1_row <- id_row_map[as.character(ID_twin1)] |
There was a problem hiding this comment.
In makeTwins(), usedID <- c(usedID, ID_twin1) runs before ID_twin1 is sampled, so usedID includes NA on the first iteration. That makes the filter !(ped$ID %in% usedID) evaluate to NA for all IDs (because %in% with NA yields NA), which can cause resample() to see NA candidates and fail or behave unpredictably. Only append ID_twin1 to usedID after it has a non-NA value (or drop NAs from usedID before using it in %in%).
Co-authored-by: Copilot <175728472+Copilot@users.noreply.github.com>
Co-authored-by: Copilot <175728472+Copilot@users.noreply.github.com>
smasongarrison
commented
Feb 18, 2026
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
This pull request introduces a major performance optimization to pedigree simulation, adds more flexible algorithm selection, and updates the documentation and tests to reflect these changes. The primary focus is the implementation of a new, vectorized "optimized" algorithm for simulating pedigrees, resulting in a 4-5x speedup for large datasets, while maintaining statistical equivalence to the original approach. Additional changes include improvements to function signatures, documentation, and test logic for both the base and optimized versions.
Pedigree Simulation Optimization and Flexibility
buildBetweenGenerationsalgorithm asbuildBetweenGenerations_optimized, significantly improving performance for large simulations while preserving statistical properties.simulatePedigreefunction and its documentation to support a flexiblebetaparameter, allowing users to choose between the original and optimized algorithms for reproducibility or speed.Testing and Validation Enhancements
test-simulatePedigree.Rto accommodate the optimized algorithm's output variability, using wider tolerances for individual counts and sex ratios, and providing clear assertions for both algorithm versions. [1] [2] [3]Code Quality and Minor Improvements
dropIdenticalDuplicateIDsand parent ID checking functions. [1] [2] [3] [4]test-segmentPedigree.R.These changes collectively improve the package's scalability, usability, and maintainability, especially for users working with large pedigree datasets.