Displaying 20 results from an estimated 5000 matches similar to: "quasi-random sequences"
2007 Oct 29
2
Changing size of lattice plot?
Hi
I am having a bit of difficulty with changing the canvas size on a
trellis/lattice plot. I am plotting two "cubes" of 3-dimensional random
numbers, as follows:
library(gsl)
library(lattice)
q <- qrng_alloc(type="sobol", 3)
npoints <- 200
rs <- qrng_get(q,npoints)
# Plot the normal variates in a 3-dim cube
p1 <- cloud(rnorm(npoints) ~ rnorm(npoints) +
2009 Apr 03
2
Help pasting string as object name
I have a data frame containing monthly observations of the 'density' of each US state, recorded in variables named "density.AL", "density.AK", "density.AZ", and so on for all 50 states. The data frame (called d) also contains a variable called "Date" which is encoded as a string in the format "Jan-09", "Feb-09", etc.
I also have
2009 Mar 25
2
"[.data.frame" and lapply
Dear all,
Trying to extract a few rows for each element of a list of
data.frames, I'm puzzled by the following behaviour,
> d <- lapply(1:4, function(i) data.frame(x=rnorm(5), y=rnorm(5)))
> str(d)
>
> lapply(d, "[", i= c(1)) # fine, this extracts the first columns
> lapply(d, "[", j= c(1, 3)) # doesn't do nothing ?!
>
> library(plyr)
2009 Mar 25
2
"[.data.frame" and lapply
Dear all,
Trying to extract a few rows for each element of a list of
data.frames, I'm puzzled by the following behaviour,
> d <- lapply(1:4, function(i) data.frame(x=rnorm(5), y=rnorm(5)))
> str(d)
>
> lapply(d, "[", i= c(1)) # fine, this extracts the first columns
> lapply(d, "[", j= c(1, 3)) # doesn't do nothing ?!
>
> library(plyr)
2008 Jun 18
2
highest eigenvalues of a matrix
DeaR list,
I happily use eigen() to compute the eigenvalues and eigenvectors of
a fairly large matrix (200x200, say), but it seems over-killed as its
rank is limited to typically 2 or 3. I sort of remember being taught
that numerical techniques can find iteratively decreasing eigenvalues
and corresponding orthogonal eigenvectors, which would provide a nice
alternative (once I have the
2008 Feb 18
2
question on function arguments
Hi,
I have two small issues with my R code, no big deal but curiosity
really. Here is a sample code,
>
> x <- rnorm(1:10)
>
> foo <- function(a = 1, b = list(x = c(1:10), y = c(1:10))){
>
> for (ii in seq(along=b$y)){
>
> print(x[ii] + b$x[ii])
> }
>
>
> }
>
> foo() # default OK
>
> foo(b=list(x=1, y=c(1:10))) # only the first
2008 Apr 03
1
data.frame or list
Dear R list,
I'm having difficulties in choosing between a list or a data.frame,
or an array for the storage and manipulation of my data (example
follows). I've been using the three for different purposes but I
would rather like to know which is more adapted to what task.
Here is the data I'm currently working on: 200 observations, each
observation being a vector of length
2008 Feb 03
1
distances between points in R^3
Dear R helpers,
I'm trying to write a numerical scheme for a boundary integral method
to solve an electromagnetic problem. This requires the computation of
the distance between points at the surface of an object (a sphere, in
my example). Here is my code,
> require(rgl)
> r<-1
> size<-10
> theta<-seq(0,2*pi,length=size*2)
> phi<-seq(0,pi,length=size)
> pc
2008 Jul 03
1
ggplot2 legend for vertical lines
Dear all,
The following example code produces a graph with ggplot2, to which I
add several vertical lines of arbitrary colors. I am not satisfied
with the legend: it automatically adds some vertical lines which I'd
rather not see (they confuse the reader rather than add information
in this case).
> library(ggplot2)
> dfr <- data.frame(values = sin(1:50/10),
> fact =
2008 Apr 12
1
lm() of one matrix against another
Hello R list,
I have two matrices of identical dimensions, and I want to fit a
straight line for each pair of columns and plot the resulting lines.
I got it to work with a for loop, but there must be a better way,
> n<-5
> N<-10
>
> data.x<-matrix(1:(n*N),ncol=n)
> data.y<-matrix(1:(n*N) + rnorm(n*N,sd=1),ncol=n)
>
>
2008 May 23
2
[slightly off topic] Sweave with markdown
DeaR list,
Has anyone tried to mix the Sweave paradigm with the Markdown[*] (and
co.) syntax? Would this be hard to implement? My tiny understanding
of Sweave is that one can define new drivers for the text part, while
some functions that deal with the R code would not require any
modification. Here's the reason I'm interested in Mardown for a driver.
I've been orbiting
2009 Apr 14
2
subset dataframe by rows using character vector?
Dear List,
I'm stuck on what seems like a simple indexing problem, I'd be very grateful to anyone willing to help me out.
I queried a dataframe which returns a character vector called "plot". I have another dataframe from which I want to subset or select only those rows that match "plot". I've tried subset, and also the "which" command.
plot
2008 Nov 09
1
[Rd] Re Bessel functions of complex argument
Dear all,
I'm writing a code that requires Bessel functions with complex argument.
Searching the list, I found the continuation of a thread I initiated a few
months ago:
http://tolstoy.newcastle.edu.au/R/e4/devel/08/03/0746.html
As I understand, the most promising option would be to use the fortran or C
implementation of Amos,
2008 Feb 14
1
write output in a custom format
Hi,
I need to create a text file in the following format,
> 1 100.0 0
> 0 0
> 1 1
> 0 0
> 1 1
> #
> 1 100.0 0
> 0 0
> 0 1
> 1 0
> 1 1
...
where # is part of the format and not a R comment.
Each block (delimited by #) consists of a first line with three
values, call it dose, and a list of (x,y) coordinates which are a
matrix or data.frame,
>
2007 Dec 28
1
unit attribute to list elements
Hi,
I've started my own (first) package, part of which consists in
listing common physical constants (Planck's constant, the speed of
light in vacuum, etc). I'm wondering what would be a good way of
dealing with pairs of value/unit.
> constants <- list( cel = 2.99792458e8 , #m/s
> Z0 = 376.730313461, #ohm
> eps0 = 8.854187817e-12,#F/m
> mu0 = 4*pi*1e-7,#N/A^2
2008 Mar 09
2
Bessel functions of complex argument
Dear R users,
I'm porting a piece of Matlab code to R, but I'm now stuck with the
following: I need an equivalent of besselJ(x, nu) that can handle a
complex argument x. I couldn't find any R implementation. I did find
a possible fortran solution in SLATEC (< http://www.netlib.org/slatec/
> , CBESJ-C), however I've never tried to use external C or Fortran
code
2008 Mar 09
2
Bessel functions of complex argument
Dear R users,
I'm porting a piece of Matlab code to R, but I'm now stuck with the
following: I need an equivalent of besselJ(x, nu) that can handle a
complex argument x. I couldn't find any R implementation. I did find
a possible fortran solution in SLATEC (< http://www.netlib.org/slatec/
> , CBESJ-C), however I've never tried to use external C or Fortran
code
2008 Dec 10
4
tapply within a data.frame: a simpler alternative?
Dear list,
I have a data.frame with x, y values and a 3-level factor "group",
say. I want to create a new column in this data.frame with the values
of y scaled to 1 by group. Perhaps the example below describes it best:
> x <- seq(0, 10, len=100)
> my.df <- data.frame(x = rep(x, 3), y=c(3*sin(x), 2*cos(x),
> cos(2*x)), # note how the y values have a different
2008 Aug 09
2
levels values of cut()
Dear list,
I have the following example, from which I am hoping to retrieve
numeric values of the factor levels (that is, without the brackets):
>
> x <- seq(1, 15, length=100)
> y <- sin(x)
>
> my.cuts <- cut(which(abs(y) < 1e-1), 3)
> levels(my.cuts)
hist() does not suit me for this, as it does not necessarily respect
the number of breaks.
getAnywhere
2009 Mar 27
3
Ploting a matrix
Hi evrybody,
in a matrix consisting of 49 columns, I would like to plot all columns
against the first in 48 different graphs.
Can you help me?
Thank you in advance
Sebastian
--
***************************************************************************************************************
Dipl. Biol. Sebastian Krug
PhD - student
IFM - GEOMAR
Leibniz Institute of Marine Sciences
Research