On Wed, 07-Feb-2007 at 07:07PM +1100, Jim Lemon wrote:
|> Matthew Keller wrote:
|> > Far from flaming you, I think you made a good point - one that I
|> > imagine most people who use R have come across. The name "R" is a big
|> > impediment to effective online searches. As a check, I entered "R
Hello R-lists,
I have created an extension library called "Rdonlp2".
This is a wrapper for Peter Spellucci's DONLP2:
http://plato.la.asu.edu/donlp2.html
DONLP2 is a standalone C program that can solve
following minimization problem:
min f(x) subject to
* parameter bounds
* linear (equality) cons
If you look at boxplot.default you will see that it calls boxplot.stats to
compute the statistics used. So you can make renamed copies of both
functions and alter them to behave as you like. (Because of namespace
issues, you cannot just use your own boxplot.stats.)
I am not sure why you would
Hi
it is not error it is just warning (Beeping a tea kettle with boiling
water is also not an error :-)
and it tells you pretty explicitly what is wrong
see length of your objects
> a<-c("D", "F", "A")
> new.dat<-subset(ex.dat, x1 == a)
Warning messages:
1: longer object length
is not a
Hi,
Let me suggest you to save your spss file in txt and use the read.table
function to load your file in R.
That is what I use to do.
Souleymane
>From: Federico Calboli <[EMAIL PROTECTED]>
>To: r-help
>Subject: [R] spss file import
>Date: Wed, 7 Feb 2007 16:16:14 +
>
>Hi All,
>
>does any
Hi, there
I am trying to replicate an error message in subset() to see what it is
that I'm doing wrong with the datasets I am trying to work with.
Essentially, I am trying to pass a string vector to subset() in order to
select a specific collection of cases (i.e., I have data for these cases i
For boxplot(), is it possible to pass in a parameter to change the default
way that the 1st and 3rd quartiles are computed? (specifically, I'd like to
use type 6 described in the quantile function).
Also, what are the options for how outliers are computed, and how can one
change them?
Thank you
Dear R,
I am current working on a mixed effects model, family poisson. Could
someone please explain to me how to check the residuals? I have read
another post with the same question but no answer unfotunately. I have
tried resid and qqnorm.
Many thanks
Beth Strain
___
On 2/11/2007 8:57 PM, Mark W Kimpel wrote:
> Duncan,
>
> Both yours and Gabor's methods were far superior to mine. I am curious
> why you like Gabor's better than yours. From the perspective of someone
> who uses R regularly but has only read about C, yours seems more
> "R-like". Would Gabor's
On Feb 11, 2007, at 8:57 PM, Mark W Kimpel wrote:
> Duncan,
>
> Both yours and Gabor's methods were far superior to mine. I am curious
> why you like Gabor's better than yours.
Don't know if the following is why Duncan prefers Gabor's method, but
here is why I would avoid the eval version: In ge
Duncan,
Both yours and Gabor's methods were far superior to mine. I am curious
why you like Gabor's better than yours. From the perspective of someone
who uses R regularly but has only read about C, yours seems more
"R-like". Would Gabor's be more computationally efficient if the loop
was big
Doug
Reinstalling R did the trick. Thanks.
Andrew
Douglas Bates wrote:
> On 2/11/07, Andrew Gelman <[EMAIL PROTECTED]> wrote:
>> Doug
>> Yes, it's R 2.4.1. I'm trying to keep up with Matrix for various
>> reasons, most immediately that our "arm" package requires Matrix.
>> Andrew
>
> The source
On 2/11/2007 4:17 PM, Robert McFadden wrote:
>
>
>> -Original Message-
>> From: Duncan Murdoch [mailto:[EMAIL PROTECTED]
>>
>> eval(parse(text=my.data))
>>
>
> I would like to thank everybody very much for help, but especially for
> Duncan - it works wonderful.
You're welcome, but I h
Hi,
By rinvgamma(1,1.33,2.33) , am I supposed to get a var=3 ? Tried
it many times , not even close. Why is this?
thanks a lot .
best
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hmmm.. I had actually at one point tried the newer (>R 2.4.0) method
of specifying only a logexp link function and using it with the
binomial family, and once again, it seems to work as expected for
glm's, and not so much with LMER. The error message when using
summary() was the same. Here's output
Dear All,
thanks for your help. I have misreaded the error message (I am sorry).
The message of Phil Spector solved the problem.
I think it was posted out of the list: I copy it below (it could be
useful to someone else).
Thanks a lot,
domenico
---
When calling a script, I want to print the code expressions and their
results using echo=T as in:
> source("myscript.r", echo = T)
This seems to do pretty much same as print() if the script only lists
expressions like "mean(data)", and echo seems nicer since it avoids extra
print statements.
But
When calling a script, I want to print the code expressions and their
results using echo=T as in:
> source("myscript.r", echo = T)
This seems to do pretty much same as print() if the script only lists
expressions like "mean(data)", and echo seems nicer since it avoids extra
print statements.
But
Try this:
my.data <- M3[[sprintf("N%04d", 1)]]
On 2/11/07, Robert McFadden <[EMAIL PROTECTED]> wrote:
> I would like to merge two parts of words to get a name of the data. First
> M3$N (invariable) and second is a number from 0001 to 3003 -
> M3$N0001,M3$N0002,...,M3$N3003. For example if I do it
> -Original Message-
> From: Duncan Murdoch [mailto:[EMAIL PROTECTED]
>
> eval(parse(text=my.data))
>
I would like to thank everybody very much for help, but especially for
Duncan - it works wonderful.
Rob
__
R-help@stat.math.ethz.ch mail
On 2/11/2007 3:39 PM, Robert McFadden wrote:
> I would like to merge two parts of words to get a name of the data. First
> M3$N (invariable) and second is a number from 0001 to 3003 -
> M3$N0001,M3$N0002,...,M3$N3003. For example if I do it like this:
> my.data <- paste("M3$N",2456,sep="")
> I get
I would like to merge two parts of words to get a name of the data. First
M3$N (invariable) and second is a number from 0001 to 3003 -
M3$N0001,M3$N0002,...,M3$N3003. For example if I do it like this:
my.data <- paste("M3$N",2456,sep="")
I get
> my.data
[1] "M3$N2456"
But I want to get something eq
Doug
Yes, it's R 2.4.1. I'm trying to keep up with Matrix for various
reasons, most immediately that our "arm" package requires Matrix.
Andrew
Douglas Bates wrote:
> On 2/11/07, Andrew Gelman <[EMAIL PROTECTED]> wrote:
>> I decided to update my packages and then had a problem with loading the
>>
On 2/11/07, Andrew Gelman <[EMAIL PROTECTED]> wrote:
> Doug
> Yes, it's R 2.4.1. I'm trying to keep up with Matrix for various
> reasons, most immediately that our "arm" package requires Matrix.
> Andrew
The source package for version 0.9975-10 Matrix has just been moved
onto CRAN. Binary versio
Dear useRs:
Release 2.0.0 of the randomSurvivalForest package is now available.
-
CHANGES TO RELEASE 2.0.0
Release 2.0.0 represents a major upgrade in the functionality and stability
of the original 1.0.0 release. K
I once had a problem upgrading to a newer R version, problem which had to do
with Env-ironment. The root cause was a .Renviron file in my home directory
pointing to an older R version. Removing that file helped. That was on
linux.
-Halim
On 2/11/07, Andrew Gelman <[EMAIL PROTECTED]> wrote:
>
> I
Dear all,
I want to evaluate a full-factorial linear mixed model with two fixed
factors (S, T) and a random factor (TM. I have tried to do the model in
lmer but at the moment this function do not provide for the F and p-values
of the fixed factors.
I tried: lmer (Arc~S*T*(1|TM), data=b,
I also
On 2/11/07, Andrew Gelman <[EMAIL PROTECTED]> wrote:
> I decided to update my packages and then had a problem with loading the
> Matrix package
> http://cran.at.r-project.org/bin/windows/contrib/2.4/Matrix_0.9975-9.zip
> This is what happened when I tried to load it in:
>
> > library("Matrix")
> E
I decided to update my packages and then had a problem with loading the
Matrix package
http://cran.at.r-project.org/bin/windows/contrib/2.4/Matrix_0.9975-9.zip
This is what happened when I tried to load it in:
> library("Matrix")
Error in importIntoEnv(impenv, impnames, ns, impvars) :
ob
Try:
aa[, !seq(ncol(aa)) %in% NULL]
On 2/11/07, Pierre Lapointe <[EMAIL PROTECTED]> wrote:
> Hello,
>
>
>
> I need to remove columns in a matrix. The number of columns varies from 0
> to n. I can't figure out how to specify the zero case.
>
>
>
> aa <-matrix(runif(5^2),5,5)
>
>
>
> #remove colu
Le 07-02-11 à 11:42, Pierre Lapointe a écrit :
> Hello,
>
> I need to remove columns in a matrix. The number of columns varies
> from 0
> to n. I can't figure out how to specify the zero case.
>
> aa <-matrix(runif(5^2),5,5)
>
> #remove column 3
>
> aa[,-3]
>
> #remove no column
>
> aa[,-NULL
Hello,
I need to remove columns in a matrix. The number of columns varies from 0
to n. I can't figure out how to specify the zero case.
aa <-matrix(runif(5^2),5,5)
#remove column 3
aa[,-3]
#remove no column
aa[,-NULL]
Error in -NULL : invalid argument to unary operator
I kn
Well, you can just check directly,
?family,
but it doesn't look like it. It looks to me (and I am admittedly no
expert) as
if make.link is only returning the standard glm links. In ?make.link, it
says,
link character or numeric; one of "logit", "probit", "cloglog",
"identity", "log", "sqrt"
Dear all,
Wolfgang, that's yet another way of calculating, which is also interesting.
But this afternoon, I did a few tries, and this is what I get:
near <- function (x,th=3){
aant <- NROW(x)
y <- 1:aant
y <- y[order(x)]
x <- x[y]
for (i in 1:(aant-1)) {
for (j in (i
Was make.link() used in the example code?
On 2/11/07, Ken Knoblauch <[EMAIL PROTECTED]> wrote:
>
> Isn't it the case, that since R 2.40 that all one ought to need do is
> define one's own link and pass it to the family, rather re-defining the whole
> family?
>
> CHANGES IN R VERSION 2.4.0
>
> ...
Dear Bart,
"hclust" might be useful for this as well:
dat = c(1,20,2,21)
hc = hclust(dist(dat))
thresh = 2
ct = cutree(hc, h=thresh)
clusteredNumbers = split(dat, ct)
firstOne = dat[!duplicated(ct)]
> clusteredNumbers
$`1`
[1] 1 2
$`2`
[1] 20 21
> firstOne
[1] 1 20
Isn't it the case, that since R 2.40 that all one ought to need do is
define one's own link and pass it to the family, rather re-defining the whole
family?
CHANGES IN R VERSION 2.4.0
...
o make.link() now returns an object of class "link-glm".
The GLM families accept an object of
On 2/11/07, Milton Cezar Ribeiro <[EMAIL PROTECTED]> wrote:
> How can I delete rows from a data.frame where almost one column is.na()?
The na.omit function does this.
> set.seed(123454321)
> x <- matrix(rnorm(50, mean = 1), ncol = 5)
> x[x < 0] <- NA
> df <- data.frame(x)
> df
X1
R Users and Developers,
Plans are being made to hold the first North American useR! will be
held at Iowa State University, Ames, Iowa, August 8–10, 2007, which
will be a week after JSM'07.
This follows successful meetings in Vienna, Austria, in 2006 and 2004,
and also Directions in Statistical Co
Look at the 'link' component of the two lists. In the binomial family
object the link component is a character vector of link 1. In your
logexposure family object it is a list of length 5.
On 2/10/07, Jessi Brown <[EMAIL PROTECTED]> wrote:
> Ok, I've tried checking out the structure of the binom
--- Milton Cezar Ribeiro <[EMAIL PROTECTED]>
wrote:
> How can I delete rows from a data.frame where almost
> one collumns is.na()?
>
> Kind regards,
>
> miltinho
Can you explain a bit more and provide as simple
example of the problem? It is not clear what
importance the column with NA's has.
Hello everyone,
How can I delete rows from a data.frame where almost one collumns is.na()?
Kind regards,
miltinho
__
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing li
Thanks Ken and Dieter,
I added xlab = '' and the text 'Scatter Plot Matrix' produced by pairs
() applied to an lmList object (in nlme) went away. Skip the rest ---
which is mainly autobiographical --- unless you're curious.
On Feb 11, 2007, at 4:25 AM, Dieter Menne wrote:
> Michael Kubovy vi
It looks like you get this default xlab with splom so the title is
misleading,
though pairs.lmList calls splom in some cases, see near the end of the
function
where plotfun is assigned a value of either "xyplot" or "splom". I was able
to get the xlab to change by providing an explicit xlab = "machi
Michael Kubovy virginia.edu> writes:
>
> I would like to suppress the text 'Scatter Plot Matrix' that appears
> under the plot. Could someone please suggest how?
Must be as special Virginia Brand of lmList. That Test does not turn up in my
output, and the only place I found it in the sources
maybe you could try something along these lines:
x <- c(1, 3, 2, 5, 11)
thr <- 3
###
ind <- t(combn(x, 2))
unique(c(ind[abs(ind[, 1] - ind[, 2]) <= thr, ]))
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University
46 matches
Mail list logo
