Skip to content

Instantly share code, notes, and snippets.

View hckiang's full-sized avatar

Woodrow (Hao Chi) Kiang hckiang

View GitHub Profile

FAQ: Bioinformatics Lab 4

Q1: What to do if BioConductor refuses to install simpleaffy?

Because installing the simpleaffy package from BioConductor does not seem to be working "out-of-box" any more, you may try one of the following solutions:

  1. If you are on Linux, you could try downloading the source package and install it manually. This should work on Mac OS X and Windows as well but I have not tested it.
@hckiang
hckiang / example_wrong_typeF.R
Created February 13, 2020 15:00
Demonstrates wrong type-F model likelihood in PCMBase
suppressPackageStartupMessages({
library("PCMBase")
library("mvtnorm")
library("ape")
})
RNGversion(min(as.character(getRversion()),"3.6.2"))
set.seed(1234, kind = "Mersenne-Twister", normal.kind = "Inversion")
pcmbase_oulik = function (tree, tipvals, H, theta, sig_x, rootval) {
## modOU = PCM(PCMDefaultModelTypes()[["F"]], k = length(rootval))