No articles match
Supported Zarr features in Rarr16 hours ago
Zarr version | Reading and Writing | Stores | Data Types | Codecs | Compression codecs | Other codecs | Optional fields or features
Working with remote Zarr arrays in R16 hours ago
Using credentials to access S3 buckets | Creating an S3 client
Working with Zarr arrays in R16 hours ago
Introduction | Limitations with Rarr | Example data | Quick start guide | Additional details | Working with Zarr metadata | Writing subsets of data | Creating an "empty" array | Updating a subset of an existing array | Appendix | Session info
ARIMA latent residuals: maths, priors, and usage15 days ago
Introduction | What an ARIMA term does | Using ARIMA terms in a formula | The kernel decomposition | Centring the integrated residual | Stationarity via partial autocorrelations | Group structure: dependent vs independent | Special cases | Priors | Where you can use it | Limitations | References
Latent process and periodic options for the growth-rate model15 days ago
Setup | Random walk on weeks | Integrated AR | Stationary AR | Gaussian process on weeks | Independent per-(week, group) effects | Day-of-week effects | Combining latent and periodic terms | Moving-average components | Fit diagnostics | When to reach for which
Model definition and implementation15 days ago
Introduction | Decomposition into expected final notifications and report delay components | Expected final notifications | Default model | Generalised model | Instantaneous reproduction number/growth rate | Latent infections/notifications | Latent reporting delay and ascertainment | Delay distribution | Parametric baseline hazard | Non-parametric reference date effect $\delta_{g,t,d}$ and report date effect $\epsilon_{g,t,d}$ | Centring | Observation model and nowcast | Accounting for reported cases with a missing reference date | Implementation | Summary of module-parameter mappings | References
Estimating reporting delays with the full and delay-only models15 days ago
Use case | Other tools for delay estimation | Two ways to estimate a delay | Getting set up | Simulating a reporting triangle | Fitting the full model | Fitting the delay-only model | Comparing the recovered delay parameters | Comparing the recovered delay distribution | Recovering the delay from data with missing cells | Using the delay-only estimate to set priors | When to use the delay-only model | References
Model Features Summary15 days ago
Overview | Core Capabilities | Different Timesteps and Timespans | Stratified and Multi-Group Nowcasting | Delay Modelling Approaches | Report Date Effects and Structural Reporting | Latent Process Models | Hierarchical Structure | Prior Specification | Missing Data Handling | Model Evaluation | Visualisation | Computational Options | Data Handling | Current Limitations | Further Reading
Discretised distributions16 days ago
Available distributions | Discretisation and adjustment for maximum delay | Double interval censoring with primarycensored
Gaussian process latent terms: maths, priors, and usage19 days ago
Introduction | What a Gaussian process term does | Using Gaussian process terms in a formula | The kernel and Hilbert-space decomposition | Kernels and special cases | Group structure: dependent vs independent | Differencing order | Priors | Where you can use it | Limitations | References
Reading HDF5 Files In The Cloud29 days ago
Public S3 Buckets | Private S3 Buckets | Session Info
Getting started with grumpy1 months ago
Motivation | Using grumpy | Structured datatypes
Supported features1 months ago
Data types | Dimensions
Design principles1 months ago
Scope | FAQ | Why not use reticulate?
Introduction to socialmixr1 months ago
Setup | The pipeline workflow | Assigning age groups | Surveys | Bootstrapping | Demography | Symmetric contact matrices | Contact rates per capita | Splitting contact matrices | Filtering | Participant weights | Temporal aspects and demography | User-defined participant weights | Weight threshold | Numerical example | Get survey data | Weight by day of week | Weight by age | Apply threshold | Plotting | Using ggplot2 | Using R base | References
Getting Started with Epinowcast: Nowcasting2 months ago
Quick start | Package | Data | Filtering | Preprocessing | Visualising the data | Choosing a nowcast horizon | Nowcast target | The default model | Posterior predictions | Alternative models | Process model | Reference model: reporting delays | Fitting the alternative models | Results | Diagnostics | Comparing all models | Using package functions rather than S3 methods | Next steps
Temporal aggregation guide2 months ago
Packages | Data | Approach 1: Weekly timestep | Approach 2: Daily process, weekly reporting (fitted day-of-week) | Approach 3: Daily process, weekly reporting (structural) | Approach 4: Daily timestep | Comparison | Weekly scale | Daily scale | Runtime | Choosing an approach
Visualising Preprocessed Data2 months ago
Setup | Data | Preprocessing | Latest observations | Cumulative reporting delay | Reporting delay heatmap | Reporting delay quantiles | Notifications by delay group | Using the individual plot functions | Helper functions
Comparing Inference Methods2 months ago
Overview | Setup | Model specifications | Fitting | NUTS with prior initialisation (default) | NUTS with pathfinder initialisation | Pathfinder (approximate inference) | Runtime comparison | Diagnostics | NUTS diagnostics | Pathfinder diagnostics | Nowcast comparison | Posterior parameter comparison | Updating with posterior-as-prior | Summary
Estimating the effective reproduction number in real-time for a single timeseries with reporting delays2 months ago
Use case | What we have | What do we do | Getting setup | Introducing the data: COVID-19 hospitalisations in Germany | Overview | Data transformations | Filtering the data | Visualising the data | Model | Expected hospitalisations | Expected infections | Instantaneous reproduction number | Latent infections | Latent reporting delay and ascertainment | Specifying the model using epinowcast::enw_expectation() | Delay distribution | Defining the delay distribution | Specifying the model using epinowcast::enw_reference() | Observation model and nowcast | Defining the observation model | Specifying the model using epinowcast::enw_obs() | Fitting the model to COVID-19 hospitalisations in Germany | Preprocess the data | Fitting the epinowcast model | Specifying the fitting options | Compiling the model | Fitting the model | Visualising the Nowcast | Plotting the nowcast based on real-time data | Plotting the nowcast based on retrospective data | Posterior predictions for cases by date of positive test and report | Real-time and retrospective estimates of the effective reproduction number | Estimates of the delay from testing positive to hospitalisation both in real-time and retrospectively | Estimates of the number of expected hospitalisations both in real-time and retrospectively | Wrapping up | Summary | Strengths | Limitations | Alternative packages | References
Design principles for the Rarr package2 months ago
Guiding principles | Scope | Zarr version | Functional programming and API design
Inferring differential exon usage in RNA-Seq data with the DEXSeq package3 months ago
Overview | Preparations | Example data | Executability of the code | Alignment | Preparing the annotation | Counting reads | Building a DEXSeqDataSet | Standard analysis workflow | Loading and inspecting the example data | Normalisation | Dispersion estimation | Testing for differential exon usage | Additional technical or experimental variables | Visualization | Parallelization and large number of samples | Perform a standard differential exon usage analysis in one command | Appendix | Controlling FDR at the gene level | Preprocessing with python | Preparing the annotation file with python | Counting reads with python | Reading the data from the python ouputs into R | Preprocessing using featureCounts | Further accessors | Overlap operations | Methodological changes since publication of the paper | Requirements on GTF files | Session Information | References
ZarrArray overview3 months ago
Introduction | Install and load the package | ZarrArray objects | Construction | Array and matrix operations | Other operations | Write an array-like object to disk in Zarr format | Session information
Resources to help with model fitting using Stan3 months ago
Epinowcast and Stan | Ensuring you have the proper toolchain | Now install CmdStanR and CmdStan | Epinowcast modelling | Installation | Running your first model | Setting enw_fit_opts | Investigating the quality of the model fit | Sampler settings | chains | threads_per_chain | iter_warmup and iter_sampling | max_treedepth | adapt_delta | Some decent defaults | Model settings | Setting priors | Exploring your data | Posterior predictions | Approaches to solve common problems | My model takes too long to run | Divergent transitions | My $\hat{R}$s are high and my esss are low | The posterior estimates are very wide | Other resources | Technical issues | Learning more about Stan and Bayesian inference
Scoring rules in scoringutils3 months ago
Introduction | Metrics for point forecasts | A note of caution | Absolute error | Squared error | Absolute percentage error | Binary forecasts | Brier score | Logarithmic score | Sample-based forecasts | CRPS | Overprediction, underprediction and dispersion | Log score | Dawid-Sebastiani score | Dispersion - Median Absolute Deviation (MAD) | Bias | Absolute error of the median | Squared error of the mean | Quantile-based forecasts | Weighted interval score (WIS) | Interval coverage | Interval coverage deviation | Quantile score | Additional metrics | Quantile coverage
Parallelize Computation of Indices 3 months ago
Running code in parallel | Performance comparison
Performance Benchmark with other packages 3 months ago
Other packages to compute FD indices | Main functions | Benchmark between packages | Functional Dispersion (FDis) | Functional Divergence (FDiv) | Functional Evenness (FEve) | Functional Richness (FRic) | Functional Richness Intersection (FRic_intersect) | Rao's Quadratic Entropy (Q) | Benchmark within fundiversity | Increasing the number of species | Functional Richness | Functional Divergence | Rao's Quadratic Entropy | Functional Evenness | Comparing between indices | Increasing the number of sites | References
Introduction: microarray quality assessment with arrayQualityMetrics3 months ago
Introduction | Basic use | Affymetrix data - before preprocessing | Affymetrix data - after preprocessing | ExpressionSet and ExpressionSetIllumina | Two colour arrays, NChannelSet, RGList, MAList | Loading data from ArrayExpress | Making the report more informative by adding a factor of interest | Extended use | Spatial layout of the array | Mapping of the reporters | RNA quality | Session Info | References
Advanced topics: Customizing arrayQualityMetrics reports and programmatic processing of the output3 months ago
Introduction | Data preparation | Module generating functions | Outlier detection | Rendering the report | Session Info
Linking to Rhdf5lib3 months ago
Motivation | Usage | Link to the library | Locating the library headers | Configuration arguments for non-standard system setups | Funding | Session info
design3 months ago
Function names, signatures, and return values | FAQ | Why a new package?
Scoring multivariate forecasts4 months ago
Univariate forecasts | Multivariate forecasts | Multivariate point forecasts
rhdf5 Practical Tips4 months ago
Introduction | Reading subsets of data | Using the index argument | Using hyperslab selections | Irregular selections | Using hyperslab selection tools | Slowdown when selecting unions of hyperslabs | Summary | Writing in parallel | Example data | Serial writing of datasets | Parallel writing of datasets | Session info
Hierarchical nowcasting of age stratified COVID-19 hospitalisations in Germany4 months ago
Packages | Data | Data preprocessing | Models | Shared reporting delay distribution | Using the inflated posterior as a prior | Reference day of the week effect | Posterior predictions | Reporting day of the week effect | Age group variation | Variation based on reference date | Variation based on reference date stratified by age | Independent models for each age group. | Alternative models | Evaluation | Summary
Case studies4 months ago
rhdf5 - HDF5 interface for R8 months ago
Introduction | High level R-HDF5 functions | Creating an HDF5 file and group hierarchy | Writing and reading objects | Writing and reading objects with file, group and dataset handles | Saving multiple objects to an HDF5 file (h5save) | List the content of an HDF5 file | Dump the content of an HDF5 file | Reading HDF5 files with external software | Removing content from an HDF5 file | 64-bit integers | Low level HDF5 functions | Creating an HDF5 file and a group hierarchy | Writing data to an HDF5 file | Session Info
Using a BioMart other than Ensembl9 months ago
Introduction | Wormbase | Phytozome | Version 12 | Version 13 | Session Info
Accessing Ensembl annotation with biomaRt9 months ago
Introduction | Selecting an Ensembl BioMart database and dataset | Step1: Identifying the database you need | Step 2: Choosing a dataset | Ensembl mirror sites | Using archived versions of Ensembl | Using Ensembl Genomes | How to build a biomaRt query | Searching for filters and attributes | Using predefined filter values | Finding out more information on filters | filterType | Attribute Pages | Using select() | Result Caching | biomaRt helper functions | exportFASTA | Examples of biomaRt queries | Annotate a set of Affymetrix identifiers with HUGO symbol and chromosomal locations of corresponding genes | Annotate a set of EntrezGene identifiers with GO annotation | Retrieve all HUGO gene symbols of genes that are located on chromosomes 17,20 or Y, and are associated with specific GO terms | Annotate set of idenfiers with INTERPRO protein domain identifiers | Select all Affymetrix identifiers on the hgu133plus2 chip and Ensembl gene identifiers for genes located on chromosome 16 between basepair 1100000 and 1250000. | Retrieve all EntrezGene identifiers and HUGO gene symbols of genes which have a "MAP kinase activity" GO term associated with it. | Given a set of EntrezGene identifiers, retrieve 100bp upstream promoter sequences | Retrieve all 5' UTR sequences of all genes that are located on chromosome 3 between the positions 185,514,033 and 185,535,839 | Retrieve protein sequences for a given list of EntrezGene identifiers | Retrieve known SNPs located on the human chromosome 8 between positions 148350 and 148400 | Given the human gene TP53, retrieve the human chromosomal location of this gene and also retrieve the chromosomal location and RefSeq id of its homolog in mouse. | Connection troubleshooting | r BiocStyle::Biocpkg("biomaRt") specific solutions | Global connection settings | Error: "SSL certificate problem" | Error: "sslv3 alert handshake failure" | Session Info
Design principles for pkgicon10 months ago
Dependencies
Frequently Asked Questions11 months ago
How to restrict the possible parameter range?
Design Principles for lightr11 months ago
Scope | Architecture | Naming conventions
Renormalise spectral data with a custom reference11 months ago
Step 1: import un-normalised data | Step 2: find the matching white reference | Step 3: normalise the reflectance data
Design Principles for lightr11 months ago
Scope | Architecture | Naming conventions
Renormalise spectral data with a custom reference11 months ago
Step 1: import un-normalised data | Step 2: find the matching white reference | Step 3: normalise the reflectance data
Interaction of mcmcensemble with other packages for MCMC diagnostic and plotting11 months ago
coda | bayesplot
Frequently Asked Questions12 months ago
Can estimations go beyond outside the range of the inits? | How to restrict the possible parameter range?
Batch import with lr_get_spec() and lr_get_metadata()12 months ago
Import spectral data: lr_get_spec() | Import spectral metadata: lr_get_metadata() | Convert spectral data to csv: lr_convert_tocsv()
Real life example12 months ago
Comparison of outputs | Example of analysis
Batch import with lr_get_spec() and lr_get_metadata()12 months ago
Import spectral data: lr_get_spec() | Import spectral metadata: lr_get_metadata() | Convert spectral data to csv: lr_convert_tocsv()
Real life example12 months ago
Comparison of outputs | Example of analysis
Design Principles for authoritative1 years ago
Scope | Naming conventions | Input/Output/Interoperability | Design decisions | Dependencies
Get started with authoritative1 years ago
Extracting R package author information | Extracting R package author information from the Author field | Extracting R package author information from the Authors@R field | Cleaning and deduplicating author names | Harmonize differently abbreviated names
Software permissions and regulations1 years ago
Scope of regulations | Data privacy and integrity | Internet access | Registration and third parties | Updates
An introduction to linelist1 years ago
Motivations | linelist in a nutshell | Outline | Should I use linelist? | Getting started | Installation | Key functionalities | Tagging system | Validation | Secured methods | Worked example | Example dataset | Creating a linelist object | Changing tags | Accessing tagged variables | Using safeguards | Changing tag loss action permanently
Compatibility with dplyr1 years ago
Verbs operating on rows | dplyr::arrange() ✅ | dplyr:distinct() ✅ | dplyr::filter() ✅ | dplyr::slice() ✅ | Verbs operating on columns | dplyr::mutate() ✓ (partial) | dplyr::pull() ✅ | dplyr::relocate() ✅ | dplyr::rename() & dplyr::rename_with() ✅ | dplyr::select() ✅ | Verbs operating on groups ✘ | Verbs operating on data.frames | dplyr::bind_rows() ✅ | dplyr::bind_cols() ✘ | Joins ✘ | Verbs operating on multiple columns | dplyr::pick() ✘
Design Principles for linelist2 years ago
Scope | Input/Output/Interoperability | Design decisions | Dependencies
HDF5 Compression Filters2 years ago
Motivation | Usage | With rhdf5 | Writing data | Reading data | With external applications | h5dump example | Compiling the compression libraries | Session info
Overview of how to use fundiversity 2 years ago
Required data | Functional Richness (FRic) - fd_fric() | Functional volume intersect (FRic_intersect) - fd_fric_intersect() | Functional Divergence (FDiv) - fd_fdiv() | Functional Evenness (FEve) - fd_feve() | Functional Dispersion (FDis) - fd_fdis() | Rao's Quadratic Entropy (Q) - fd_raoq() | Large site-species data / sparse matrices | Standardizing trait data | Non-continuous traits? | Missing values in traits? | Functions summary table | References
List of included regions3 years ago
List of included regions3 years ago
Visualise age pyramids3 years ago
Visualise contact matrices3 years ago
Visualise age pyramids3 years ago
Visualise contact matrices3 years ago
Design decisions4 years ago
Why do I need to pass a list of files to parse rather than an entire folder?
Check Numerical Correctness of Indices 4 years ago
Changes of abundance | Changes of coordinates (= traits)
Design Principles for fundiversity 4 years ago
Scope | Dependencies | Functions | Each index should be computed in its own separate function | Input data should not be transformed without any explicit action from the user | Inputs | Outputs | Functions should output data.frames | The outputs of functions should be similar in structure
Introduction to pavo5 years ago