We revisit the analysis for the...
R posts in other languages
- 适合初学者的R语言教程
- SAS hat einen neuen Vergleich von SAS und R veröffentlicht
- 实现可重复的统计slides
- О связи между “утечкой мозгов” и объемом пассажирских перевозок
- Новая книга по основам прогнозирования (с использованием R)
- Двухфакторный дисперсионный анализ
- Revolution R Enterprise 6.2 mit neuen Funktionen für hoch performante Analysen auf Big Data
- Third Milano R net meeting: leave a comment
- Photos of the third Milano R net meeting
- Presentations of the third Milano R net meeting
Employee productivity as function of number of workers revisited
Employee productivity as function of number of workers revisited
We revisit the analysis for the...
The most violent municipalities in Mexico (2008)
- The top six most violent municipalities are near the US border. Ciudad Juárez is in a class by itself with 113 homicides per 100,000 people.
- José Azueta is the municipality where Zihuatanejo is located. Mazátlan, another popular tourist destination, also appears on the list.
- Lázaro Cárdenas is the largest seaport in Mexico and ever since the...
The most violent municipalities in Mexico (2008)
- The top six most violent municipalities are near the US border. Ciudad Juárez is in a class by itself with 113 homicides per 100,000 people.
- José Azueta is the municipality where Zihuatanejo is located. Mazátlan, another popular tourist destination, also appears on the list.
- Lázaro Cárdenas is the largest seaport in Mexico and ever since the...
R Layout command.
In the previous post I created a chart but could not figure out to fit the legend in the chart area. Peter Carl pointed me to the layout command which partitions the display area and allowed the the legend to be included. Source code to produce the c...
MMDS 2010
The 2010 Workshop on Algorithms for Modern Massive Data Sets (MMDS 2010) finished up this past Friday (June 18th) at Stanford. This was an exceptionally well organized conference: four days of mind-stretching talks on algorithm development and the challenges of working with massive data sets approached from almost every conceivable angle. The approximately 100 attendees were a diverse group...
New blog from Rmetrics Foundation
The Rmetrics Foundation (the sharp minds behind the Rmetrics suite of packages for financial analysis in R) have just launched a new blog where you can keep up with the latest Rmetrics news. Amongst the recent news: a ne eBook about data management of Indian financial market data, and a new interface between Rmetrics and AMPL. You can also...
Example 7.42: Testing the proportionality assumption
In addition to the non-parametric tools discussed in recent entries, it's common to use proportional hazards regression, (section 4.3.1) also called Cox regression, in evaluating survival data.It's important in such models to test the proportionality a...
A quetion
I posted a question to the R-Mixed-Model mailing list but have not had any responses yet:--------------------------------------------------------------------------------------- Dear All,I wonder if anybody has tried to make glmmBUGS work with JAGS. Mya...
The police records for 2009 are out.
The 2009 homicide numbers collected by the SNSP (National System of Public Security) are finally out, you can download the data from the ICESI, which is a civic institution not affiliated with the government. If you remember, one of the conclusio...
The police records for 2009 are out.
The 2009 homicide numbers collected by the SNSP (National System of Public Security) are finally out, you can download the data from the ICESI, which is a civic institution not affiliated with the government. If you remember, one of the conclusio...
Here’s the distribution of the first million digits of the…
Here’s the distribution of the first million digits of the square root of two’s decimal expansion. Number of digits | is: 0's | 99 818 1's | 98 926 2's | 100 442 3's | 100 191 4's | 100 031 5's | 100 059 6's | 99 885 ...
QSPR modeling with signatures
I had to dig deep to find posts on QSAR modeling. There are quite a few on QSAR in Bioclipse, but that focuses on the descriptor calculation. In a quick scan, I could only spot two modeling posts:The CDK/Metabolomics/Chemometrics Unconference resultsWh...
Estimating Probability of Drawdown
I've shown several examples of how to use LSPM's probDrawdown function as a constraint when optimizing a leverage space portfolio. Those posts implicitly assume the probDrawdown function produces an accurate estimate of actual drawdo...
More powerful iconv in R
The R function iconv converts between character string encodings, for example, from the locale dependent encoding to UTF-8: > iconv("foo", to="UTF-8") [1] "foo" However, R has long-running trouble with embedded null characters ('') in strings. Hence, if we try to convert to an encoding that permits embedded null characters, iconv will fail: > iconv("foo", to="UTF-16")
What I need to know…
is maps and geographical data representation in R. In case you’re curious too this is a good study material from R-Bloggers : maps ; geographical ; spatial Ok. This could be a tweet rather than a post…
ggplot2 GUI progress
(Written by Ian Fellows) Below is a link to the first of a weekly (or bi-weekly) screen-cast vlog of my progress building a GUI for the ggplot2 package. http://neolab.stat.ucla.edu/cranstats/gsoc_vlog1.mov comments and suggestions are more than welcome, and can e-mailed to me at: [email protected]
The perfect fake
Usually when you are doing Monte Carlo testing, you want fake data that’s good, but not too good. You may want a sample taken from the Uniform distribution, but you don’t want your values to be uniformly distributed. In other words, if you were to order your sample values from lowest to highest, you don’t
Why R doesn’t suck
I first encountered the R programming language a few years ago when I needed to make some plots. Although I’ve used it occasionally since, I always considered it a sort of “Perl for statisticians” — a useful swiss-army knife with … Continue reading →
Those dice aren’t loaded, they’re just strange
I must confess to feeling an almost obsessive fascination with intransitive games, dice, and other artifacts. The most famous intransitive game is rock, scissors, paper. Rock beats scissors. Scissors beats paper. Paper beats rock. Everyone older than 7 seems to know this, but very few people are aware that dice can exhibit this same behavior,
Revolution Analytics: Startup to watch
Jack Germain of LinuxInsider interviewed Revolution CEO Norman Nie for his "Startup to Watch" column. Amongst the topics covered: the R language (Norman: "There are no statistical expressions that can not be written in R"), Revolution's recent name-change and announcement of our development roadmap, and the challenges of competing with SAS and Norman's former company, SPSS. Read the full...
The impact of the drug war in Mexico
For the last couple of years, Mexico has been in the midst of an escalating drug war, with violent crime on the upswing in many areas. But tracking the impact quantitatively is difficult: in Mexico, about 85% of crimes go unreported, and corruption leads to inaccurate reporting in some districts. Diego Valle has taken on the task of visualizing...
R: Command Line Calculator using Rscript
I currently use an awesome little bash trick to get a command line calculator that was posted on lifehacker, and that I blogged about previously.calc(){ awk "BEGIN{ print $* }" ;}You just add this to your .bashrc file and then you can use it ...
R: Command Line Calculator using Rscript
I currently use an awesome little bash trick to get a command line calculator that was posted on lifehacker, and that I blogged about previously.calc(){ awk "BEGIN{ print $* }" ;}You just add this to your .bashrc file and then you can use it ...
R Commander – linear regression
We can fit various linear regression models using the R Commander GUI which also provides various ways to consider the model diagnostics to determine whether we need to consider a different model. Fast Tube by Casper The “Statistics” menu provides access to various statistical models via the “Fit models” sub-menu including: Linear regression – the simplest scenario with
Occupational Wage Comparison Plotted in R
Ever have conversations with your kids about what they are going to do with their life? Still trying to figure out what you are going to do with yours?Best to not starve...The chart above represents the percentage of each occupation that earn a given h...
Manhattan plots for SNP marker effects using ggplot2
At AIPL, we've been posting Manhattan plots of the marker effects for each breed-trait combination with each official release of our genomic predictions. For example, consider the plot of lifetime net merit for Holsteins from the April, 2010 run: Thes...
Do Not Log-Transform Count Data, Bitches!
OK, so, the title of this article is actually Do not log-transform count data, but, as @ascidacea mentioned, you just can’t resist adding the “bitches” to the end. Onwards. If you’re like me, when you learned experimental stats, you were taught to worship at the throne of the Normal Distribution. Always check your data and
Zero Inflated Models and Generalized Linear Mixed Models with R.
Zuur, Saveliev, Ieno (2012).
