115 search results for "pca"

Stepping Outside My Open-Source Comfort Zone: A First Look at Golden Helix SVS

May 16, 2012
By
Stepping Outside My Open-Source Comfort Zone: A First Look at Golden Helix SVS

I'm a huge supporter of the Free and Open Source Software movement. I've written more about R than anything else on this blog, all the code I post here is free and open-source, and a while back I invited you to steal this blog under a cc-by-sa license....

Read more »

Visualising the Path of a Genetic Algorithm

April 23, 2012
By
Visualising the Path of a Genetic Algorithm

We quite regularly use genetic algorithms to optimise over the ad-hoc functions we develop when trying to solve problems in applied mathematics. However it’s a bit disconcerting to have your algorithm roam through a high dimensional solution space while not being able to picture what it’s doing or how close one solution is to another. … Continue reading...

Read more »

Customising vegan’s ordination plots

April 11, 2012
By
Customising vegan’s ordination plots

As a developer on the vegan package for R, one of the most FAQs is how to customise ordination diagrams, usually to colour the sample points according to an external grouping variable. Now, just because we get asked how to … Continue reading

Read more »

tcp timelines with ggplot2

March 28, 2012
By
tcp timelines with ggplot2

I’ve come across the need to analyze TCP flows from time to time, and while scripts like flowtime and EasyTimeline are nice, they aren’t really, well, pretty.  ggplot2, on the other hand is, and it turns out to be really easy to … Continue reading

Read more »

Ruin probability and infinite time

March 27, 2012
By
Ruin probability and infinite time

A couple of weeks ago, I had a discussion with a practitioner, working in some financial company, about ruin, and infinite time. And it remind me a weird result. Well, not a weird result, but a result I found disturbing, at first, when I was a stud...

Read more »

Julia functions for the Rmath library

March 26, 2012
By
Julia functions for the Rmath library

1.1 Signatures of the d-p-q-r functions in libRmath Users of R are familiar with the functions for evaluating properties of and for sampling from probability distributions, the so-called d-p-q-r functions. The designation d-p-q-r reflects the fact t...

Read more »

PCA for NIR Spectra_part 006: "Mahalanobis"

February 28, 2012
By
PCA for NIR Spectra_part 006: "Mahalanobis"

Outliers have an important influence over the PCs, for this reason they must be detected and examinee.We have just the spectra without lab data, and we have to check if any of the sample spectra is an outlier ( a noisy spectrum, a sample which belongs ...

Read more »

PCA for NIR Spectra_part 005: "Reconstruction"

February 27, 2012
By
PCA for NIR Spectra_part 005: "Reconstruction"

We saw how to plot the raw spectra (X), how to calculate the mean spectrum, how to center the sprectra (subtracting the mean spectrum from every spectra of the original matrix X). After that we have developed the PCAs with the NIPALS algorithm, getting...

Read more »

PCA for NIR Spectra_part 004: "Projections"

February 26, 2012
By
PCA for NIR Spectra_part 004: "Projections"

This plot in 2D, help us to decide the number of PCs, it is easy to create in R, once we have discompose the X matrix into a P matrix (loadings) and a T matrix (scores).For this plot, we just need the T matrix.> CPs<-seq(1,10,by=1)>  matp...

Read more »

PCA for NIR Spectra_part 003: "NIPALS"

February 25, 2012
By
PCA for NIR Spectra_part 003: "NIPALS"

> X<-yarn$NIR> X_nipals<-nipals(X,a=10,it=100)Two matrices are generated (P and T)As in other posts, we are going to look to the loadings & scores, for firsts three principal components:> wavelengths<-seq(1,268,by=1)> matplot(w...

Read more »