Lambda Theta Phi

These files are an adaptation of the ANOVA variance component estimators to be used in the analysis of genetic data, specifically Cockerham's 1969 Theta. Formulae were adapted from Sokal and Rohlf's 1995 Biometry, and test files from the book are given as reference.

Auto Gaussian & Gabor Surface fit---Functions to fit a 1D Gaussian to a curve and a 2D Gaussian or Gabor to a surface. The routines are automatic in the sense that they do not require the specification of starting guesses for the model parameters. This is done by evaluating the quality of fit for...

The codes provided can be used to show that the Kalman Smoother Expectation Maximization (KSEM) methodology can be used successfully to estimate the parameter of the Schwartz-Smith model. We develop several novel adjustments to incorporate the model framework.The advantage of using this method...


Platforms: Matlab

EARTH_ROSE direction rose from earth coordinate angles. Syntax: EARTH_ROSE(E) EARTH_ROSE(E, nbins) Inputs: E: Angular current direction in degrees, earth coordinates. OPTIONAL: nbins: number of bins to plot histogram Example: theta = [110 110 110 90 288 290 292]; nbins = 20; earth_rose(theta,...


Platforms: Matlab

Represent a wavefront as a sum of Zernike polynomials using a matrix inversion.This function attempts to solve the a_i's in equation, phi(rho,theta) = SUM(a_i * Z_i(rho,theta), i=1 to M),where the Z_i(rho,theta)'s are the Zernike polynomials from the zernfun.m file, phi is the wavefront to be...


Platforms: Matlab

This code calculates the Exponentially Weighted Moving Average Standard Deviation Exponentially weighted moving average (EWMA) standard deviation applies different weights to different returns. More recent returns have greater weight on the variance. The exponentially weighted moving average...


Platforms: Matlab

This will perform the same function as Matlab built-in function ATAN2 butsafely guarded against erroneous results when the inputs are very small (near to a machine zero).Due to the round off errors in the numerical calculations, the twoinputs to ATAN2, Y and X, had better be thought as x=x_exact...


Platforms: Matlab

Conversion routines to take a surface (a function z of (x,y) or (theta,r)) and convert this into a solid; much as meshz does, but a little more useful for outputting as an STL for frabrication in a 3D printer. There are two cases, one in Cartesian co-ordinates, and the other in polar...


Platforms: Matlab

[XR,YR,xr,yr,THETA] = rotateData(X,Y,X0,Y0,THETA,DIRECTION) Rotate coordinates especified by [X,Y] around the point [X0,Y0], by an angle (in radians) defined by THETA and to the direction especified by DIRECTION. THETA is a scalar if the user wishes to define it. THETA can also be calculated...


Platforms: Matlab

YIPLOT makes plot much easier. H=YIPLOT(X,Y,...) H=YIPLOT(STR,Y) H=YIPLOT(STR,X,Y,...) STR is a charater string made from any elements below, seperated by comma, eg. 'fig=1,lw=4,t=title,x=x axis label' fig = # set # as current figure, same as command figure(#) sub = # subplot(#) lw = # set...


Platforms: Matlab

PlotSphereIntensity(azimuth, elevation)PlotSphereIntensity(azimuth, elevation, intensity)h = PlotSphereIntensity(...)Plots the intensity (as color) of a number of points on a unit sphere.Input: azimuth (phi), in degrees elevation (theta), in degrees intensity (optional, if not provided, a green...


Platforms: Matlab

Create a custom 2D gauss, which can be used for filtering, weighting, etc. All parameters are customizable, including standard deviations (sigmaX, sigmaY), rotation (theta), result size, center, etc.


Platforms: Matlab

Auto Gaussian & Gabor Surface fit---2 functions to fit a 2D Gaussian or 2D Gabor to a surface. The routines are automatic in the sense that they do not require the specification of starting guesses for the model parameters. This is done by evaluating the quality of fit (sum of squared error) for...


Platforms: Matlab

Find those points (xi,yi) who are inside the polygone defined by the vectors (x,y).Syntax: [xin,yin,Ind]=Inside(X,Y,xi,yi)where Ind is the index of (xi,yi) that are inside the polygone Example: %--- create polygone (x,y) phi=linspace(0,2*pi,50); X=(1+0.2*sin(6*phi)).*cos(phi);...


Platforms: Matlab

Recent studies on Mathematical modeling of visual cortical cells [Kulikowski/Marcelja/Bishop:1982] suggest a tuned band pass filter bank structure. These filters are found to have Gaussian transfer functions in the frequency domain. Thus, taking the Inverse Fourier Transform of this transfer...


Platforms: Matlab

secplot(theta, rho)Plot wedges of equal angular width (centered on each value in theta, which should be equally spaced), with the radii listed in rho.Addition input arguments are passed to 'set' which is run on the object handle returned by polar (used to do the plotting).Example:th =...


Platforms: Matlab

GRMustache is an Objective-C implementation of the Mustache logic-less template engine. Its parser has been inspired by the Mustache go implementation. Its tests are based on the Ruby one, that we have considered as a reference. It supports the following Mustache features: * comments *...


Platforms: Mac

Using Affymetrix CEL files as its source, this script combines quantile normalized SNP A and B probe intensities for each marker into a theta value, then calculates B-Allele Frequencies for each marker. For information on accessing B Allele Frequencies for Illumina data,


Platforms: Mac

Software developed to perform electromagnetic scattering simulations mainly based on classical Mie theory solution. ScatLab Features: scattered intensity polar diagrams for coated and uncoated spherical particles; scattered intensity versus theta graphs for coated and uncoated spherical...


Platforms: Windows

Be an instant stock options expert! How do we know if stock options - or currency options - are overpriced or underpriced? If you can answer that question then you can consistently buy stock options that are underpriced, or consistently sell options that are overpriced ... without special...


Platforms: Windows