Displaying 20 results from an estimated 3000 matches similar to: "Markov chain transition model, data replication project"
2009 Mar 15
3
read.xls question
I'm an R newbie and had a question about the read.xls function. I've heard that this is often not a reliable function to use for importing data. However, I have created numerous xls files which contain information about voter turnout and macroeconomic indicators in India. I'm writing a paper on the relationship between economic growth and voter turnout.
This is the command I use:
2007 Oct 19
1
plot.Design
Dear R-users:
I am trying to use the following code to reproduce the figures on page 140
of Prof. Frank Harrell's book 'Regression Modeling Strategies':
rm(list=ls())
options(width=128)
library(Hmisc)
library(Design)
getHdata(counties)
counties$older <- counties$age6574 + counties$age75
label(counties$older) <- '% age >= 65, 1990'
counties$pdensity <-
2004 Apr 29
3
Probability(Markov chain transition matrix)
Hello, My name is Maria, MBA student in Sanfransisco, USA.
In my credit scoring class, I have hard time making "transition matrix",
which explains probability especially in relation to "Markov chain model".
It is regarding people's monthly credit payment behavior. Does R have
function to caculate it? I am actually a novice in using 'R'. Please help
me!!!
Maria
2002 May 14
2
tapply and empty subsets
I am using tapply to compute means by group, as in
> tapply(y, z, mean)
1 2 3
21.00757 20.50031 NA
where y is a vector and z a factor with levels 1:3.
I would like to get a zero instead of the NA I get when a subset is
empty. The problem is that FUN is not applied to empty cells of the
ragged array. (Of course I am willing to define a function my.mean
that returns
2011 May 17
2
can not use plot.Predict {rms} reproduce figure 7.8 from Regression Modeling Strategies (http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RmS/course2.pdf)
Dear R-users,
I am using R 2.13.0 and rms 3.3-0 , but can not reproduce figure 7.8 of the
handouts *Regression Modeling Strategies* (
http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RmS/course2.pdf) by the
following code. Could any one help me figure out how to solve this?
setwd('C:/Rharrell')
require(rms)
load('data/counties.sav')
older <- counties$age6574 + counties$age75
2010 Aug 18
3
libreadline problem
Hello,
Recently updated Ubuntu to 10.04, I have installed the r-base and
r-base-dev packages from CRAN with
sudo apt-get install r-base
But now when I want to start R I get the following error:
gpetris at definetti:~$ R
/usr/lib/R/bin/exec/R: symbol lookup
error: /usr/local/lib/libreadline.so.6: undefined symbol: PC
Any clue??? Has anybody seen that before? I have tried to reinstall the
2006 Oct 19
2
randomize a matrix
Hello everyone,
If I have an incidence matrix of 0 and 1's
P=[1 1 1 1 1 1
1 1 1 1 0 0
1 1 1 0 0 0
1 1 1 0 0 0
1 1 0 0 0 0]
I want to create a new uniform random matrix [a] that is filled with 0's
and 1's but constrained so that the row and column sums are the same as
in [P]. Does anyone know how to accomplish this?
Thanks in advance
Cameron Guenther, Ph.D.
2003 Oct 29
1
restarting split.screen
Is there a way of `restarting' split.screen?
This is what I am getting:
> close.screen()
[1] 10 11 12 13
> close.screen(all=TRUE)
Error in par(args) : parameter "i" in "mfg" is out of range
> graphics.off()
> x11()
> close.screen()
[1] 10 11 12 13
> close.screen(all=TRUE)
Error in par(args) : parameter "i" in "mfg" is out of range
As
2007 Jan 18
4
Reading contingency tables
I am trying to read an ftable using read.ftable, but I get the
following error message:
> jobSatTable <- read.ftable("http://definetti.uark.edu/~gpetris/stat5333/jobSatisfaction.dat",skip=2)
Error in seek(file, where = 0) : no applicable method for "seek"
In addition: Warning messages:
1: no non-missing arguments to max; returning -Inf
2: no non-missing arguments to
2004 May 12
4
non-interactive call to R (running an R package as a stand-alone application)
Is there a way I can have R automatically execute the commands in a source file without ever having to use R interactively? If so, what arguments should I pass to the UNIX call to R? I need to do this to run several R jobs in parallel.
An alternative may be to have R and an R package behave as a stand-alone application that can be called from the UNIX command line. Is there any documentation on
2008 Oct 09
2
Two math expressions in plot
Hello!
I am trying to put two math expressions in the title of a plot. As
you can see below, I can place correctly one expression at a time, but
not both. Ideally I would like to have them separated by a comma. Any
suggestions?
> k <- 1
> n.eff <- c(20, 30)
> ### this works
> plot(0,0, main = substitute(n == k, list(k = k)))
> ### this works
> plot(0,0, main =
2012 May 25
1
knitr customization
I am trying to transition from Sweave to knitr, but there are a few
things about customization of the appearence of R input and output that
I did not get yet. Maybe somebody on the list can help me.
In my Sweave presentations I used a slanted font for the R input and a
normal font for the output, both in a small font. I also indented
everything by an extra 2em. All this was achieved by the
2011 Mar 09
4
Help with read.csv
Hello,
I have a file that looks like this:
Date,Hour,DA_DMD,DMD,DA_RTP,RTP,,
1/1/2006,1,3393.9,3412,76.65,105.04,,
1/1/2006,2,3173.3,3202,69.20,67.67,,
1/1/2006,3,3040.0,3051,69.20,77.67,,
1/1/2006,4,2998.2,2979,67.32,69.10,,
1/1/2006,5,3005.8,2958,65.20,68.34,,
where the ',' is the separator and I tried to read it into R, but...
> y <- read.csv("Data/Data_tmp.csv",
2004 Feb 06
3
quantile function
I am trying to `cut' a continuous variable into contiguous classes
containing approximately an equal number of observations. I thought
quantile() was the appropriate function to use in order to find the
breakpoints, but I end up with classes of different sizes - see
example below. Does anybody have an explanation for that? And what is
the `recommended' way of computing what I am looking
2003 Oct 24
5
how to remove NaN columns ?
How can I remove columns with NaN entries ?
Here is my simple example:
> data <- read.csv("test.csv")
> xdata <- data[3:length(data)]
> xs <- lapply(xdata, function(x){(x - mean(x))/sqrt(var(x))})
> x <- data.frame(xs)
> x
C D E F
1 -0.7071068 NaN -0.7071068 -0.7071068
2 0.7071068 NaN 0.7071068 0.7071068
2012 May 09
1
Sweave, beamer and alert within code chunks
Hi all,
Using Beamer, in order to highlight a piece of R code I do something
like this - note the "\structure" and "\alert" commands:
\begin{semiverbatim}
> mleOut <- \structure{dlmMLE}(Nile,
+ parm = c(0.2, 120), # initial values for optimizer
+ lower = c(1e-7, 0)) \alert<2>{# V must be positive}
> mleOut$convergence
2011 Mar 15
2
(R) transitions in a Markov Chain
Dear experts,
I have to generate 1000 transitions of a discrete time Markov Chain and
then give and estimation of the equilibrium distribution.
I know the transition matrix.
0.3 0.7 0
0.2 0.5 0.3
0 0.4 0.6
Do you know how to do it with R?
Thanks a lot!
Estefania
[[alternative HTML version deleted]]
2005 Mar 22
1
Package vignette and build
Hello,
I am writing a package called 'DLM' containing a vignette.
The vignette contains a chunck with the function call 'library(DLM)'.
This worked fine with 'R CMD check DLM', but when it comes to building
the package with 'R CMD build DLM' I get the following error message:
* creating vignettes ... ERROR
Error: chunk 1
Error in library(DLM) : There is no
2010 Jan 21
2
What is the difference between S and R?
On Thu, Jan 21, 2010 at 9:07 AM, Giovanni Petris <GPetris at uark.edu> wrote:
>
> I have found a good refernce to be "S Programming" by Venables and
> Ripley.
I'll take a look at this book.
Since S and R are not completely the same, there are delicate
differences between S and R, which an S book may be confusing if I
want to learn R. One example is Statistical Models
2002 Aug 02
1
Means of Monte Carlo simulated lists
Hello,
I am doing simulations, and I generate a list at each iteration (with
three component matrices in the example below), saving the results in
a list. For example, after two iterations, I have something like
> str(sim.theta)
List of 2
$ :List of 3
..$ : num [1:6, 1:4] -3.67 -1.07 -2.99 -18.38 -3.26 ...
..$ : num [1:6, 1:6] -7.56 -3.14 -4.99 1.03 2.79 ...
..$ : num [1:6, 1:4]