Displaying 20 results from an estimated 6000 matches similar to: "SVD of a variance matrix"
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
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
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
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
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 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
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
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
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
2003 Apr 03
2
Printing zero as dot
I'm pretty sure I've seen some examples of a function printing zero
entries in a matrix as dots, but I'm not able to find it now...
Any suggestions...? Thanks in advance. (Of course, I might have dreamt
of such a function...)
Best,
Giovanni
--
__________________________________________________
[ ]
[ Giovanni Petris
2003 Oct 16
1
Improving efficiency in "outer"-like calculation
Hello,
I am doing mcmc=10000 simulations from a posterior distribution of the parameters
of a mixture of K=6 normal densities.
I have mcmc by K matrices simMeans, simVars and simWeights containing
the simulation output: one row for each simulation, one column for
each normal component of the mixture.
One thing I would like to do is a plot of the posterior predictive
density. In order to do that
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]
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",
2011 Aug 25
2
Synchronizing R libraries on N machines?
Hello!
I am using R on two different machines (under Ubuntu and OS X, but this
is probably irrelevant) and I would like to keep the two installations
'synchronized', in particular in terms of installed packages. For
example, if I install package xxx on my Linux machine, I would like to
find it installed also on my Mac, and vice versa.
I imagine this to be a fairly common problem, so I
2008 May 02
1
GLMM and data manipulation (2nd try)
Hello,
I posted a question yesterday but I got no replies, so I'll try to
reformulate it in a more concise way.
I have the following data, summarizing approval ratings on two
different surveys for a random sample of 1600 individuals:
> ## Example: Ratings of prime minister (Agresti, Table 12.1, p.494)
> rating <- matrix(c(794, 86, 150, 570), 2, 2)
> dimnames(rating) <-
2001 Dec 06
1
Building R on Solaris with SUNperf
Dear All,
I am building R-1.3.1 on Solaris 2.7 and I am trying to use the native
SUNperf BLAS/LAPACK. I have added the flags ("-dalign
-xlic_lib=sunperf") for the compilers in the config.site file.
What is not clear to me is what I have to do with the "--with-blas"
flag for ./configure. It seems to me that omitting it is not enough.
Any suggestions?
Another question, I
2007 Dec 18
1
Multiple plots with single box
Hello,
I am trying to display some harmonic functions in a plot. The kind of
display I have in mind is like the one that cn be obtained by a call
to plot.ts with plot.type = "multiple". The only difference is that I
want a single box containing all the plots instead of one box per
plot. I thought box(which = "outer") would have done the job, but it
didn't.
Below is the
2009 Mar 17
1
Combining Sweave and fancyvrb
I find Sweave very useful and I was trying to combine it with the
latex package fancyvrb. I was trying to get line numbering and labels
in order to reference the lines where particular commands occur.
Unfortunately, I haven't been able to figure out how to do it. Maybe
somebody can help me.
The following is a sample Rnw file: the first part shows what I would
like to get, the second what I
2004 Jun 25
4
Installing on Windows packages build on Unix
Hello,
I wanted to share with a colleague a few R functions that I wrote. To
this purpose, I created a small package on my machine (Unix) and
emailed it to her. Now she is having troubles installing the package
on her Windows machine. It seems that on her side, install.packages
looks for a "zip" file - while I have created a "tar.gz" file. I tried
to build the package with the
2009 Sep 11
3
State Space models in R
Hello everybody,
I am writing a review paper about State Space models in R, and I would
like to cover as many packages as I reasonably can.
So far I am familiar with the following tools to deal with SS models:
* StructTS, Kalman* (in stats)
* packages dse[1-2]
* package sspir
* package dlm
I would like to have some input from users who work with SS models:
are there any other packages for SS