Hi,
the following R lines work fine in R 2.4.0 alpha (and older R versions), but
not in R
2.4.0 beta (details below):
library(drc) # to load the dataset 'PestSci'
library(nlme)
## Starting values
sv <- c(0.328919, 1.956121, 0.097547, 1.642436, 0.208924)
## No error
m1 <- nlme(SLOPE ~ c +
Dear John
?ordered
will help you.
Regards,
Christoph Buser
--
Christoph Buser <[EMAIL PROTECTED]>
Seminar fuer Statistik, LEO C13
ETH Zurich 8092 Zurich SWITZERLAND
phone: x-41-44-632-4673 fax: 632-1228
http://stat.
Hi,
It seems to me that the Mersenne-Twister PRNG can be initialised using
one integer instead of 624 integers, since inside RNG.c code there's a
function defined as MT_sgenrand(Int32).
How do I actually set this seed within R?
I've tried:
> .Random.seed <- c(3, 1)
> runif(1)
Error in runif
On Mon, 25 Sep 2006, Gad Abraham wrote:
> Hi,
>
> It seems to me that the Mersenne-Twister PRNG can be initialised using
> one integer instead of 624 integers, since inside RNG.c code there's a
> function defined as MT_sgenrand(Int32).
>
> How do I actually set this seed within R?
set.seed(), on
Dear Gabor and Deepayan,
Many thanks for your help.
I used suggestion 3 from Gabor (it worked well with my long df ) and
will try Deepayan's suggestion.
Rafael
Deepayan Sarkar wrote:
> On 9/23/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
>
>> 1. You can write a custom strip function:
>>
>>
Christian Ritz <[EMAIL PROTECTED]> writes:
> Hi,
>
> the following R lines work fine in R 2.4.0 alpha (and older R versions), but
> not in R
> 2.4.0 beta (details below):
>
>
> library(drc) # to load the dataset 'PestSci'
>
> library(nlme)
>
>
> ## Starting values
> sv <- c(0.328919, 1.95
I've not seen an actual answer to this, which is that this is a
misunderstanding as to how NextMethod works.
First,
> + x <- unclass(x)
looks wrong. NextMethod uses the next method at the call to the generic,
and subsequent changes to the object 'x' do not alter the class that would
be d
Hello all,
I hope i chose the right list as my question is a beginner-question.
I have a data set with 3 colums "London", "Rome" and "Vienna" - the
location is presented through a 1 like this:
London RomeVienna q1
0 0 1 4
0 1 0
Your problem would be a lot easier if you coded the location in one
variable instead of three variables. Then you could calculate the means
with one line of code:
by(results$q1, results$location, mean)
With your dataset you could use
by(results$London, results$location, mean)
by(results$Rome, res
Peter,
There is a much easier way to do this. First, you should consider
organizing your data as follows:
set.seed(1) # for replication only
# Here is a sample dataframe
tmp <- data.frame(city = gl(3,10, label = c("London", "Rome","Vienna"
)), q1 = rnorm(30))
# Compute the means
with(tmp, tappl
Dear,
I have developed and tested some models in soil
hydrology with NLME library in R. I want to ask
if it could be possible to submit this to the NLME
library (with sample data) as a toolbox or something
so
that anyone downloading new components of new versions
of R may simply call (say
SSbrooksc
Dear all,
I have another small scripting-beginner problem which you hopefully can
help:
I compute new variables with:
# Question 1
results$q1 <- with(results, q1_1*1+ q1_2*2+ q1_3*3+ q1_4*4+ q1_5*5)
# Question 2
results$q2 <- with(results, q2_1*1+ q2_2*2+ q2_3*3+ q2_4*4+ q2_5*5)
# Question 3
res
Hi
I am trying to impute missing values for my data.frame. As I intend to use the
complete data for prediction I am currently measuring the success of an
imputation method by its resulting classification error in my training data.
I have tried several approaches to replace missing values:
- mean
--- Peter Wolkerstorfer - CURE <[EMAIL PROTECTED]>
wrote:
> Hello all,
>
> I hope i chose the right list as my question is a
> beginner-question.
>
> I have a data set with 3 colums "London", "Rome"
> and "Vienna" - the
> location is presented through a 1 like this:
> LondonRomeVie
I think this does what you are looking for:
dta <- data.frame(q1_1=rep(1,5),q1_2=rep(2,5),q2_1=rep(3,5),q2_2=rep(4,5))
for (i in 1:2) {
e1 <- paste("q",i,"_1 + q",i,"_2 * 2",sep="")
assign(paste("q",i,sep=""),with(dta,eval(parse(text=e1
}
On 25/09/06, Peter Wolkerstorfer - CURE <[E
Dear R-help,
I recently tried implementing glmmPQL in 2.3.1, and I discovered a
few differences as compared to 2.2.1. I am fitting a regression with
fixed and random effects with Gamma error structure. First, 2.3.1
gives different estimates than 2.2.1, and 2.3.1, takes more
iterations to con
You might find something useful at this web site:
http://www.multiple-imputation.com/
On 25/09/06, Eleni Rapsomaniki <[EMAIL PROTECTED]> wrote:
>
> Hi
>
> I am trying to impute missing values for my data.frame. As I intend to use the
> complete data for prediction I am currently measuring the succ
- Original Message -
From: "David Barron" <[EMAIL PROTECTED]>
To: "Peter Wolkerstorfer - CURE" <[EMAIL PROTECTED]>; "r-help"
Sent: Monday, September 25, 2006 3:33 PM
Subject: Re: [R] Beginner Loop Question with dynamic variable names
>I think this does what you are looking for:
>
> dta
Hello Xiaohui,
data.grp is just a pseudo example of a grouped data object
that is grouped according the factor y.
I just tried your suggestion, but the result is that
three separate panels are still created, whereas I would like
to have all 3 lines in a single panel.
cheers,
dave
-Ori
On 25-Sep-06 Eleni Rapsomaniki wrote:
> Hi
>
> I am trying to impute missing values for my data.frame. As I
> intend to use the complete data for prediction I am currently
> measuring the success of an imputation method by its resulting
> classification error in my training data.
>
> I have tried
Hi,
subset function is use to select rows of a dataframe.
just compute
>mean(results$q1)
without subset instruction, or
>mean(results[,4])
Peter Wolkerstorfer - CURE <[EMAIL PROTECTED]> a écrit :
Hello all,
I hope i chose the right list as my question is a beginner-questio
Hi Rneters,
I'm running a GLM model with a full factorial design in blocks and binomial
error distribution. I would like to have the F values for this model but I
got a message that "using F test with a binomial family is inappropriate in:
anova.glm(model, test = "F")". Should I not report F sta
This is what I have been using. Does anyone have a better way? In
particular I would like to see letters in comment strings not stretched
so much. Thanks -Frank
\documentclass{article}
\usepackage{listings,relsize}
\lstloadlanguages{R}
\newcommand{\lil}[1]{\lstinline|#1|}
\begin{document}
\l
Hi All,
I have a data with a variable like this:
Column 1
"123abc"
"12cd34"
"1e23"
...
Now I want to do an operation that can split it into two variables:
Column 1Column 2 Column 3
"123abc" 123 "abc"
"12cd34" 12"cd34"
"1e23"
On Mon, 2006-09-25 at 11:30 -0500, Marc Schwartz (via MN) wrote:
> On Mon, 2006-09-25 at 11:04 -0500, Frank Duan wrote:
> > Hi All,
> >
> > I have a data with a variable like this:
> >
> > Column 1
> >
> > "123abc"
> > "12cd34"
> > "1e23"
> > ...
> >
> > Now I want to do an operation that can s
Here is one more solution:
library(gsubfn)
s <- c("123abc", "12cd34", "1e23")
out <- gsubfn("^([[:digit:]]+)(.*)", paste, s, backref = -2)
read.table(textConnection(out))
It assumes there are no spaces in the strings. If
there are then choose a sep= that does not appear
and do this:
sep = ","
Dear R users,
I am trying to use MCD and MVE methods in the analysis of functional imaging
(fMRI) data. But, before doing that, I want to understand the sampling
distribution of the correlation parameter given by MCD and MVE (cov.mcd$cor,
cov.mve$cor).
To this end, I conducted a simulation where
strapply in package gsubfn can do that:
library(gsubfn)
s <- c("123abc", "12cd34", "1e23")
out <- strapply(s, "^([[:digit:]]+)(.*)", c)
out <- do.call(rbind, out) # as a matrix
data.frame(x = out[,1], num = as.numeric(out[,2]), char = out[,3]) #
as a data.frame
On 9/25/06, Frank Duan <[EMAIL
And here is a third solution not using package gsubfn:
s <- c("123abc", "12cd34", "1e23")
out <- gsub("^(([[:digit:]]+)(.*))", "\\1 \\2 \\3", s)
read.table(textConnection(out), as.is = TRUE)
Again, if spaces appear in the input string choose a character
not appearing, such as comma, and do it lik
On Mon, 2006-09-25 at 11:04 -0500, Frank Duan wrote:
> Hi All,
>
> I have a data with a variable like this:
>
> Column 1
>
> "123abc"
> "12cd34"
> "1e23"
> ...
>
> Now I want to do an operation that can split it into two variables:
>
> Column 1Column 2 Column 3
>
> "123abc"
Dear R users,
This command works (calling a programm -called whap- with file specifiers etc.):
>system('cmd /c "c:\\pheno\\whap --file c:\\pheno\\smri --alt 1 --perm 500"',
>intern=TRUE)
Now I need to call it from a loop to replace the "1" by different number,
however I get lost using the q
>>>Now I want to do an operation that can split it into two variables:
>>>
>>>Column 1Column 2 Column 3
>>>
>>>"123abc" 123 "abc"
>>>"12cd34" 12"cd34"
>>>"1e23" 1 "e23"
>>>...
>>>
>>>So basically
Hello, is there a way to pass a connection to a file in a zipped archive as
argument (instead of a file name of unzipped file) to shell command "cut". In
general, is it possible to pipe output of a R function to a shell command? How?
I want to do something like:
z = unz("zipArchive.zip", "fileASC
Great! That's exactly what I want.
Thanks a lot,
FD
On 9/25/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
>
> strapply in package gsubfn can do that:
>
>
> library(gsubfn)
> s <- c("123abc", "12cd34", "1e23")
>
> out <- strapply(s, "^([[:digit:]]+)(.*)", c)
> out <- do.call(rbind, out) # as
On 9/23/06, Afshartous, David <[EMAIL PROTECTED]> wrote:
>
> All,
> I'd like to plot the main relationship of a grouped data
> object for all levels of a factor in a single panel.
> The sample code below creates a separate panel for each level
> of the factor. I realize that this could be done in
No, the "cut" command won't understand that "z" is an R connection and
not a file in the current working directory: there is no overlap
between the R object name space and the Windows object name space.
Unfortunately, you may be forced to unzip to a temporary file, and
then read from that.
One th
Hey R-Comunity,
I'd like to print out an histogram of some experimental data and add a
smooth curve of a normal distribution with an ideally generated
population having the same mean and standard deviation like the
experimental data.
The experimental data is set as vector x and its name is set
On 9/25/2006 1:56 PM, Lothar Botelho-Machado wrote:
> Hey R-Comunity,
>
>
> I'd like to print out an histogram of some experimental data and add a
> smooth curve of a normal distribution with an ideally generated
> population having the same mean and standard deviation like the
> experimental dat
Gabor Grothendieck gmail.com> writes:
>
> It might also be nice to be able to align the fans at the left or right,
> not just the center.
Fans that open only on one side: A line that moves like the minute needle of an
analog clock; with zero at the top. Movement of the needle in clock-wise
dire
Dear R users,
I am trying to use MCD and MVE methods in the analysis of functional imaging
(fMRI) data. But, before doing that, I want to understand the sampling
distribution of the correlation parameter given by MCD and MVE (cov.mcd$cor,
cov.mve$cor).
To this end, I conducted a simulation where
What is the meaning of this message?
Warning message:
Realizando coerção de LHD para uma lista
I tried to do something like this:
test <- function(x) {
rval <- NULL
m <- mean(x)
s <- sd(x)
rval$m <- m
rval$s <- s
y <- x[abs(x - m) > 3 * s]
rval$y <- y # this is the critical li
On Mon, 2006-09-25 at 19:56 +0200, Lothar Botelho-Machado wrote:
> Hey R-Comunity,
>
>
> I'd like to print out an histogram of some experimental data and add a
> smooth curve of a normal distribution with an ideally generated
> population having the same mean and standard deviation like the
> exp
On Mon, 2006-09-25 at 18:58 +0200, Boks, M.P.M. wrote:
> Dear R users,
>
> This command works (calling a programm -called whap- with file specifiers
> etc.):
>
> >system('cmd /c "c:\\pheno\\whap --file c:\\pheno\\smri --alt 1 --perm 500"',
> >intern=TRUE)
>
> Now I need to call it from a lo
Le Lundi 25 Septembre 2006 09:31, Frank E Harrell Jr a écrit :
> This is what I have been using. Does anyone have a better way? In
> particular I would like to see letters in comment strings not stretched
> so much. Thanks -Frank
>
> \documentclass{article}
> \usepackage{listings,relsize}
> \lst
useRs,
I have a data frame where four of the columns of the data frame represent the
values of a two-by-two
matrix. I'd like to, row-by-row, go through the data frame and use the four
columns, in matrix form, to
perform calculations necessary to create new values for variables in the data
frame
Using the builtin 11x8 anscombe data frame here are
some alternatives:
# 1
# list of 2x2 matrices
lapply(split(anscombe[1:4], 1:nrow(anscombe)), matrix, 2)
# 2
# 2x2x11 array
array(t(anscombe[1:4]), c(2, 2, nrow(anscombe)))
# 3
# to create matrix and perform calculations, e.g. det, all in one st
Is this what you had in mind?
> j<-data.frame(q1=rnorm(10),q2=rnorm(10))
> j
q1 q2
1 -0.9189618 -0.2832102
2 0.9394316 1.1345975
3 -0.6388848 0.6850255
4 0.4938245 -0.5825715
5 -1.2885257 -0.2654023
6 -0.5278295 0.2382791
7 0.6517268 0.8923375
8 0.4124178 1.1231
thanks!
-Original Message-
From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
Sent: Monday, September 25, 2006 1:34 PM
To: Afshartous, David
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] plotting grouped data object
On 9/23/06, Afshartous, David <[EMAIL PROTECTED]> wrote:
>
> All,
> I'd
Use single outer quotes so that the inner double quotes are not
interpreted as the end of the string.
cmd <- 'cmd /c "...whatever..." '
system(cmd, intern = TRUE)
On 9/25/06, Boks, M.P.M. <[EMAIL PROTECTED]> wrote:
> Dear R users,
>
> This command works (calling a programm -called whap- with file
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Thanks for your help!!
I appreciate, now it works perfectly.
Lothar Rubusch
Duncan Murdoch wrote:
> On 9/25/2006 1:56 PM, Lothar Botelho-Machado wrote:
>> Hey R-Comunity,
>>
>>
>> I'd like to print out an histogram of some experimental data and ad
For a nested model you want to use the coxme function, which is the much
superior successor to frailty(). It is currently found in the kinship
library.
Terry Therneau
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listin
# R version 2.3.1 (2006-06-01) Debian Linux "testing"
# Is the following behaviour a bug, feature or just a lack of
# understanding on my part? I see that this was discussed here
# last March with no apparent resolution.
d <- as.factor(c("1970-04-04","1970-08-11","1970-10-18"))
x <- c(9,10,11)
ch
If you want it to act like a date store it as a Date:
dx <- as.Date(c("1970-04-04","1970-08-11","1970-10-18")) ###
x <- c(9,10,11)
ch <- data.frame(Date=dx,X=x)
dy <- as.Date(c("1970-06-04","1970-08-11","1970-08-18")) ###
y <- c(109,110,111)
sp <- data.frame(Date=dy,Y=y)
merge(ch, sp, all = TRUE
Hello!
J.R. Lockwood rand.org> writes:
>
> An alternative that I've used a few times is the jpg() function to
> create the sequence of images, and then converting these to an mpeg
> movie using "mencoder" distributed with "mplayer". This works on both
> windows and linux. I have a pretty self-
Here is a slight simplification of the strapply solution using simplify = TRUE
library(gsubfn)
s <- c("123abc", "12cd34", "1e23")
out <- t(strapply(s, "^([[:digit:]]+)(.*)", c, simplify = TRUE)) # matrix
data.frame(x = out[,1], num = as.numeric(out[,2]), char = out[,3])
On 9/25/06, Gabor Grothe
Frank E Harrell Jr vanderbilt.edu> writes:
>
> This is what I have been using. Does anyone have a better way? In
> particular I would like to see letters in comment strings not stretched
> so much. Thanks -Frank
It may be possible to pass on all comments to a verbatim like environment insi
I have a complicated nonlinear function, myfun(a,b,c),
that I want to fit to data, allowing one or more of the parameters
a, b, and c in turn to have linear dependence on other covariates.
In other words, I'd like to specify something like
nls(y~myfun(a,b,c),linear=list(a~f1,b~1,c~1))
I kn
Hi,
the contributed package 'drc' allows specification of non-linear regression
models with
individual parameter models that include covariates.
For an example see section 8 the accompanying paper in J. Statist. Software
(http://www.jstatsoft.org/v12/i05/v12i05.pdf).
Christian
__
Dear All,
Ive got a statistical question on calculating
dissimilarities in R.
I want to calculate the different types of dissimilarities
on the flower dataset found in the package
cluster. Flower is a data frame with 18 observations
on 8 variables. Variable 1 and 2 are binary, variable 3 is
as
On Mon, 25 Sep 2006, Bruce LaZerte wrote:
> # R version 2.3.1 (2006-06-01) Debian Linux "testing"
>
> # Is the following behaviour a bug, feature or just a lack of
> # understanding on my part? I see that this was discussed here
> # last March with no apparent resolution.
Reference? It is the th
Dear All,
Ive got a statistical question on calculating
dissimilarities in R.
I want to calculate the different types of dissimilarities
on the flower dataset found in the package
cluster. Flower is a data frame with 18 observations
on 8 variables. Variable 1 and 2 are binary, variable 3 is
a
61 matches
Mail list logo