I’m starting on my real R-Scripts now. We got an assesment for some Reports and they needed my graphs and Veerle’s Reports. So the tasks are: Well distribution plot Quality plate control Surface Plot Heatmap CirclePlot and more to come̷...

I’m starting on my real R-Scripts now. We got an assesment for some Reports and they needed my graphs and Veerle’s Reports. So the tasks are: Well distribution plot Quality plate control Surface Plot Heatmap CirclePlot and more to come̷...

Are you a non-windows user wanting to share your R package with others who do use windows? There's a handy utility (http://win-builder.r-project.org/)that allows you to upload your .tar.gz file and it checks it, builds it, and sends you a link to downl...

Are you a non-windows user wanting to share your R package with others who do use windows? There's a handy utility (http://win-builder.r-project.org/)that allows you to upload your .tar.gz file and it checks it, builds it, and sends you a link to downl...

(The examples here work with the version of insidefunctor tagged as "v1") I ran into an interesting problem using "inside" functors. Something is wrong in the following code (well, depending on what you thought it should do). > library(insidefunctor) > `%+.%` = fmap(`+`) > x = seq(0, 10, len = 50) > plot(x,...

Hell is new and improved. The R Inferno has been revised. If you don’t know of it, it is a short explanation of a few trouble spots when using the R language. Somehow the short explanation grew to approach book-length. It can be found at the usual place: http://www.burns-stat.com/pages/Tutor/R_inferno.pdf Major improvements An index has been … Continue reading...

The solar irradiance incident outside the earth’s atmosphere is called the extraterrestial or extra-atmospheric irradiance. It is derived from the solar constant only with geometric equations. It can be easily calculated with the calcSol function of the solaR package. With this post I will show an example with some packages from the Spatial task view.

Gompertz Model Visualization # Gomperz growth function gomp <- function(x, a, b, k) a*exp(-b*exp(-k*x)) # Normal model with Gompertz mean function likelihood <- function(weight, age, sigma, a, b, k) { mu <- gomp(age, a, b, k) dnorm(weight, mu, sigma) } # Visualize the model visualize <- function(phi=40, theta=-35) { weight <- seq(0, 250,