Articles by Wingfeet

Looking at Measles Data in Project Tycho

March 30, 2014 | Wingfeet

Project Tycho includes data from all weekly notifiable disease reports for the United States dating back to 1888. These data are freely available to anybody interested. I wanted to play around with the data a bit, so I registered.MeaslesMeasles a... [Read more...]

Distribution Kinetics in JAGS

March 23, 2014 | Wingfeet

Chapter 19 of Rowland and Tozer (Clinical pharmacokinetics and pharmacodynamics, 4th edition) has distribution kinetics. I am examining problem 3. It is fairly easy, especially since essentially all formulas are on the same page under 'key relationships'. The only difficulty is that the formulas are symmetrical in λ1, c1 and λ2, c2 and they ... [Read more...]

PK calculations for infusion at constant rate

March 16, 2014 | Wingfeet

In this third PK posting I move to chapter 10, study problem 4 of Rowland and Tozer (Clinical pharmacokinetics and pharmacodynamics, 4th edition). In this problem one subject gets a 24 hours continuous dose. In many respects the Jags calculation is not very different from before, but the relation between parameters in the ... [Read more...]

PK calculation of IV and oral dosing in JAGS

March 9, 2014 | Wingfeet

I am examining IV and oral dosing of problem of Chapter 6, Question 6 or Roland and Tozer Rowland and Tozer (Clinical pharmacokinetics and pharmacodynamics, 4th edition) with Jags. In this problem one subject gets an IV and an oral dose.DataThe data lo... [Read more...]

Simple Pharmacokinetics with Jags

March 2, 2014 | Wingfeet

In this post I want to analyze a first order pharmocokinetcs problem: the data of study problem 9, chapter 3 of Rowland and Tozer (Clinical pharmacokinetics and pharmacodynamics, 4th edition) with Jags. It is a surprising simple set of data, but still there is enough to play around with.Data, model and ... [Read more...]

Unemployment revisited

February 23, 2014 | Wingfeet

Approximately a year ago I made a post graphing unemployment in Europe and other locations. I have always wanted to do this again, not because the R-code would be so interesting, but just because I wanted to see the plots. As time progressed I attempte... [Read more...]

Bayesian analysis of sensory profiling data III

February 16, 2014 | Wingfeet

Last week I extended my Bayesian model. So this week I wanted to test it with different data. There is one other data set with profiling data in R, french fries in the reshape package. 'This data was collected from a sensory experiment conducted at Iow... [Read more...]

Bayesian analysis of sensory profiling data

February 2, 2014 | Wingfeet

I looked at Bayesian analysis of sensory profiling data in May and June 2012. I do remember not being totally happy with the result and computations taking a bit more time than I wanted. But now it is 2014, I can use STAN and I have been thinking about... [Read more...]

Converting plots to data II

January 26, 2014 | Wingfeet

There were quite some reactions on lasts week's post converting plots to data. Two additional programs were mentioned; WebPlotDigitzier and DataThief. Note that google or wikipedia has a bunch of alternatives. Datathief is shareware, a java program, so... [Read more...]

Convering plots to data

January 18, 2014 | Wingfeet

It is a problem which occurs ever so often in applied work, you have a plot, but you want the data. There are at least two programs which can help you there; PlotDigitizer and Engauge Digitizer. I got both on my openSuse machine. Both are available for... [Read more...]

Converting a JAGS model to STAN

January 11, 2014 | Wingfeet

For my first experience with STAN I wanted to convert my last JAGS program into STAN. This was a bit more difficult than I expected. The JAGS program was Fe concentration in rainwater including values below detection level.DataData has been explained b... [Read more...]

Le Monde puzzle 847 in Julia

January 4, 2014 | Wingfeet

This week I wanted to play around with Julia and exporting the results. I found http://xianblog.wordpress.com/2013/12/29/le-monde-puzzle-847/ to be just the right size to play around with. CodeA function to check if a triplet has the desired property In [1]: function lemonde847(xx) a=[xx[1],2] b=[xx[2],8] c=[... [Read more...]

Sinterklaas and Santa Claus gave presents

December 27, 2013 | Wingfeet

In December my nice little netbook acquired a shiny new openSUSE 13.1, including all the goodies I might want to use in 2014; R, Stan, Julia and Jags. So here is what you might get when you set up a computer for statistics today.JagsI have been using J... [Read more...]

plotting y and log(y) in one figure

December 15, 2013 | Wingfeet

Sometimes I have the desire to plot both on the linear and on the log scale. To save space just two figures is not my solution. I want to reuse the x-axis, legend, title. This post examines possibilities to do so with standard plot tools, lattice and g... [Read more...]

Fukushima region radiation decrease

December 8, 2013 | Wingfeet

Fukushima is still a topic which gets headlines, and somewhere in a comment was a link to actual and historical radiation data: Fukushima prefecture radioactivity measurement map.  It takes a bit of time to load but then you have a map of the Fukushima region with current radiation levels at ... [Read more...]

Dutch Rainwater Composition 1992-2011

November 24, 2013 | Wingfeet

Last week I examined rainwater composition 1992 to 2005. There is additional data, but National Institute for Public Health has changed equipment in 2005. This week I will add those data.DataData is in a number of spreadsheets. The scrip... [Read more...]

Dutch Rainwater Composition 1992-2005.

November 17, 2013 | Wingfeet

After reading Blog About Stats' Open Data Index Blog Post I decided to browse a bit in the Open Data Index. Choosing Netherlands and following Emission of Pollutants I ended on a page from National Institute for Public Health. The page&n... [Read more...]

A small comparison of bio-equivalence calculations.

November 10, 2013 | Wingfeet

Last week I looked at two-way cross-over studies and followed the example of Schütz (http://bebac.at/) in the analysis. Since the EU has its on opinions (Questions & Answers: Positions on specific questions addressed to the pharmacokinetics working party) and two example data sets, I was wondering how the ... [Read more...]
1 4 5 6 7 8 10

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)