Hallo.
I'm running step() on lm() in R version 2.15.1 within RStudio on my
xUbuntu Linux computer with 4 core processor.
Linux top command states that rsession (R) is using from 100 % to
about 330 % of the processor, which (I think) means that R is using
more than one core. Is it possible? I
I am trying to detect a positively increasing trend in a time series data,
to
which when I apply the rolling mean it almost turns into a step wave
function. The sharp rise in slope of rolling mean is an indication of
positive trend, thus the rolling mean becomes like a step wave function.
Is
Hello,
I am attempting to use nlme to model the response of 19 groups. I am able to
get a reasonable fit of all of the groups to the observed data using an
exponential decay model (a*exp-x*b). The problem is that when I plot the
fitted values to residuals, it demonstrates a pattern of
Hi ,
I want to apply the following code fo my data with 400 predictors.
I was wondering if there ia an alternative way instead of typing 400 predictors
for the following code.
I really appreciate your help.
fit0-lm(Y~1, data= mydata)
fit.final- lm(Y~X1+X2+X3+.+X400, data=mydata) ???
Probably you just want
Y ~ .
where . means everything which doesn't appear elsewhere in the formula.
Michael
On Thu, Nov 15, 2012 at 8:45 PM, farnoosh sheikhi farnoosh...@yahoo.com wrote:
Hi ,
I want to apply the following code fo my data with 400 predictors.
I was wondering if there ia
I've been having a problem using the step function to evaluate models. I've
simplified the code and get the same problem using the dataset Titanic. The
relevant code and output is below. The problem disappears (i.e., 'step' runs
correctly) if I rerun the code but change the 'loglm' call to
On 17/05/2012 09:25, David A Vavra wrote:
I've been having a problem using the step function to evaluate models. I've
simplified the code and get the same problem using the dataset Titanic. The
relevant code and output is below. The problem disappears (i.e., 'step' runs
correctly) if I rerun the
.
DAV
-Original Message-
From: Prof Brian Ripley [mailto:rip...@stats.ox.ac.uk]
Sent: Thursday, May 17, 2012 6:46 AM
To: David A Vavra
Cc: r-help@r-project.org
Subject: Re: [R] step function stops with subscript out of bounds
On 17/05/2012 09:25, David A Vavra wrote:
I've been
On 03.08.2011 16:33, Diana Schwegler wrote:
Hello I am using the step function in order to do backward selection for a
linear model of more than 200 variables but it doesn't work correctly.
I think, there is a problem, if the matrix has same or more columns than
rows.
And if the matrix has too
Hello I am using the step function in order to do backward selection for a
linear model of more than 200 variables but it doesn't work correctly.
I think, there is a problem, if the matrix has same or more columns than
rows.
And if the matrix has too much columns the step-function doesn't work
Hello Friends,
Can you kindly help me with step by step procedure (in terms
of the required/necessary statistical tests) before threshold model can be
applied in time series analysis?
I need to
understand all the necessary tests - ranging from testing the properity of a
time
Hi Noah,
Are you able to reproduce the example on a smaller dataset? Do you have any
strange variable names or I created a 3 x 100 matrix, fit a linear model
and step has been running fine (other than bringing my poor netbook to it's
knees). It also might be helpful if you could post
On 10.01.2011 10:13, Joshua Wiley wrote:
Hi Noah,
Are you able to reproduce the example on a smaller dataset? Do you have any
strange variable names or I created a 3 x 100 matrix, fit a linear model
and step has been running fine (other than bringing my poor netbook to it's
knees).
Hi,
Its a lot of data, but here are sum summary stats:
l - lm(trainy ~ x)
str(x)
num [1:31205, 1:48] 0.0975 -0.1987 0.3254 -0.7912 0.0975 ...
- attr(*, dimnames)=List of 2
..$ : chr [1:31205] 5 6 7 8 ...
..$ : NULL
- attr(*, names)= chr [1:1497840] a NA NA NA ...
summary(x)
V1
I think I just figured it out.
x is a matrix.
l - lm(y ~ x) works for generating a model, but fails. (It considers x
as a single item to add/remove for step.)
Step does work if I use a data.frame
foo - cbind(y,x)
l - lm(y ~ ., data=foo)
Now step(l) works.
I guess R doesn't look at the x in
Hi,
I have a fairly simple linear regression using the lm function. There
are about 100 variables and 30,000 rows of data. It runs fine and
produces a decent looking R2 value. I'm interested in performing a
stepwise variable selection to see if things can be cleaned up a bit.
Calling the step
Hello!
I wouldl ike to ask you if R supports step by step execution. I have written
some nested loops and I would like to check on every step what are the values
of some variables. Printing all the variables just creates a really big output
of numbers.
Could you please
Hello
On Fri, Oct 22, 2010 at 11:33 AM, Alaios ala...@yahoo.com wrote:
I wouldl ike to ask you if R supports step by step execution. I have written
some nested loops and I would like to check on every step what are the
values
of some variables.
fortune('browser')
My solution when I
Alaios ala...@yahoo.com writes:
Hello!
I wouldl ike to ask you if R supports step by step execution. I have written
some nested loops and I would like to check on every step what are the
values
of some variables. Printing all the variables just creates a really big
output
of
On Oct 22, 2010, alais wrote:
|Moreover can you please tell me how I can concatenate variables
and strings so
to print some debugging messages /
The value of x is + x + and the value of y is +y.
Commas, not +
cat(The value of x is , x , and the value of y is , y)
--
David
Gavin Simpson gavin.simpson at ucl.ac.uk writes:
On Fri, 2010-08-27 at 11:03 -0300, Silvano wrote:
Hi,
how can I change the significance level in test F to select
variable in step command?
I used
step(model0, ~x1+x2+x3+x4, direction=c(forward), test='F',
alpha=.05)
Hi,
how can I change the significance level in test F to select
variable in step command?
I used
step(model0, ~x1+x2+x3+x4, direction=c(forward), test='F',
alpha=.05)
but it does't work.
--
Silvano Cesar da Costa
Departamento de Estatística
On Fri, 2010-08-27 at 11:03 -0300, Silvano wrote:
Hi,
how can I change the significance level in test F to select
variable in step command?
I used
step(model0, ~x1+x2+x3+x4, direction=c(forward), test='F',
alpha=.05)
but it does't work.
Well, there is no 'test' argument, nor
Hi All
I need some help with plotting a step function, currently I'm using
sfun - stepfun(c(1, 2, 5,10, 20), c(0, 11, 22, 33, 44, 0), f=0)
plot(sfun, pch=NA, main=, xlim=c(1,20))
which I working fine, but my data is in the following format:
Min Max Value
1 2 11
2 5 22
510 33
10 20 44
1. To
Dear helpers,
in course of my master thesis I trie to program a GLM to model the
number of seeds produced per flower relates to the time they were
exposed to pollinators (categorical with two levels), the altitude of
population (categorical with two levels) and the species (categorical
Hello I am using the step function in order to do backward selection for a
linear model of 52 variables with the following commands:
object-lm(vars[,1] ~ (vars[,2:(ncol(predictors)+1)]-1))
BackS-step(object,direction=backward)
but it isn't dropping any if the variables in the model, but there
Does anybody know how to perform a step function by significance of p instead
of AIC? I want to perform a forward-backward stepwise logistic regression
and want to compare the results obtained by steps of significance and by
steps of AIC. Thank you
--
View this message in context:
Alba wrote:
Does anybody know how to perform a step function by significance of p
instead of AIC? I want to perform a forward-backward stepwise logistic
regression and want to compare the results obtained by steps of
significance and by steps of AIC.
MASS::dropterm, and play with k. And
Is it possible to use the step() function with a glmer() as an object? I
obtain the following error message when I try to do it: Error in x$terms :
$ operator not defined for this S4 class.
I perform the glmer correctly but I can't do the step.
Thank you so much.
--
View this message in
Tukey
-Oorspronkelijk bericht-
Van: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] Namens Alba
Verzonden: dinsdag 9 februari 2010 14:07
Aan: r-help@r-project.org
Onderwerp: [R] step and glmer
Is it possible to use the step() function with a glmer
Thank you for your rapid answer!
--
View this message in context:
http://n4.nabble.com/step-and-glmer-tp1474390p1474580.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Hi All,
Does the step function work in this model? I tried to run the
following model but no result obtained. The computer is hanging and I
killed the job several times. Below is the code.
library(survival)
m.fit=clogit(y~x1+x2+x3+x4, data=ftest)
summary(m.fit)
final- step(m.fit)
Thanks in
On Fri, 29 Jan 2010, Ashta wrote:
Hi All,
Does the step function work in this model? I tried to run the
following model but no result obtained. The computer is hanging and I
killed the job several times. Below is the code.
library(survival)
m.fit=clogit(y~x1+x2+x3+x4, data=ftest)
Hello everybody,
I would need some help from you.
I am trying to fit a logistic model to some presence absence data of
animals living on river islands. I have got 12 predictor variables and I am
trying to use a stepwise forward method to fit the best logistic model to
my data. I am using the
Hi,
On Wed, Jan 27, 2010 at 12:21 PM, francesca.ior...@googlemail.com wrote:
Hello everybody,
I would need some help from you.
I am trying to fit a logistic model to some presence absence data of
animals living on river islands. I have got 12 predictor variables and I am
trying to use a
Can you think of any systemic changes that might interefere with R
besides Symantec EndPoint and LiveUpdate? I have removed those
programs and allocated more memory to R, but it is still way too
slow.
On Nov 13, 10:45 pm, J Dougherty j...@surewest.net wrote:
On Friday 13 November 2009 07:17:28
I can't fully answer all these questions, but I'll do my best - There
have not been any updates of Windows, and I did not update R during
the period, although I did reinstall it after the problem started.
There have been no changes to Norton or any other software that uses
system resources in that
Jgabriel wrote:
I can't fully answer all these questions, but I'll do my best - There
have not been any updates of Windows, and I did not update R during
the period, although I did reinstall it after the problem started.
There have been no changes to Norton or any other software that uses
Alright, here's the update: the IT guy insists that the RAM is not the
problem. RAM, according to him, is pretty black and white. There
should be enough RAM on the system for the function to work, and he
was the one who installed the extra RAM in the first place. Could
there be another RAM issue
On Friday 13 November 2009 07:17:28 am Jgabriel wrote:
I can't fully answer all these questions, but I'll do my best - There
have not been any updates of Windows, and I did not update R during
the period, although I did reinstall it after the problem started.
There have been no changes to
I have a question that might not be kosher here, but I'm running out
of options and need some help. Basically using the function step is
freezing R. I am running a model that includes a number of
interactions on a large data set with a number of dummy variables
representing whether an event
There are some points of information you should include. One, have there been
recent updates of Windows that may have altered system behaviour? Two, have
you updated R during the same period? Three, have you installed software that
uses system resources that were available to R before - e.g.
Dear R-listers,
I am now asking for step by step tutorial to beginning of R programming by
Emacs+ESS on MS Windows. Though it was appreciated highly, there is much
troubles we must conquered before we set up R programming by Emacs+ESS,
especially for users of MS Windows and who without
Thank you, Duncan. While your answers didn't quite fix my problems, I
realized that I should probably go with the flow. I installed gfortran
and everything worked fine. Thanks again.
Steve
sjnovick wrote:
To all. I need your help. The big question is: How do I make an R
Michael wrote:
Could anybody point me to the latest status of the most user-friendly
debugger in R?
I always have been happy with ?debug, ?recover and friends cited there.
Although you also find additional software like the debug package ..
Best,
Uwe Ligges
How I wish I don't have to
To all. I need your help. The big question is: How do I make an R library
with Fortran 95 files? You may assume that I am a pretty decent programmer
in both R and Fortran. I lay out a scenario and need your help!
I know how to make an ordinary R package and have dabbled with R + Fortran
95
On 5/22/2009 8:18 AM, sjnovick wrote:
To all. I need your help. The big question is: How do I make an R library
with Fortran 95 files? You may assume that I am a pretty decent programmer
in both R and Fortran. I lay out a scenario and need your help!
I know how to make an ordinary R
Really I think if there is a Visual Studio strength debugger, our
collective time spent in developing R code will be greatly reduced.
2009/5/22 Uwe Ligges lig...@statistik.tu-dortmund.de:
Michael wrote:
Could anybody point me to the latest status of the most user-friendly
debugger in R?
On 5/22/2009 10:59 AM, Michael wrote:
Really I think if there is a Visual Studio strength debugger, our
collective time spent in developing R code will be greatly reduced.
If someone who knows how to write a debugger plugin for Eclipse wants to
help, we could have that fairly easily. All the
Duncan Murdoch wrote:
On 5/22/2009 10:59 AM, Michael wrote:
Really I think if there is a Visual Studio strength debugger, our
collective time spent in developing R code will be greatly reduced.
If someone who knows how to write a debugger plugin for Eclipse wants
to help, we could have that
Hi,
Romain Francois wrote:
Duncan Murdoch wrote:
On 5/22/2009 10:59 AM, Michael wrote:
Really I think if there is a Visual Studio strength debugger, our
collective time spent in developing R code will be greatly reduced.
If someone who knows how to write a debugger plugin for Eclipse wants
As a newbie I'm trying to figure out how much more than RKWard does is
wanted. The code turns colors as syntax is checked and errors are noted.
It seems like a reasonable IDE. Maybe someone is looking for the same in
windows?
John F Lindsey
803-790-5006 Home , 803-790-5008 Cell
919-439-9088
Could anybody point me to the latest status of the most user-friendly
debugger in R?
How I wish I don't have to stare at my (long) code for ages and stuck...
Thanks!
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Hi all,
Â
Iâm using step and stepAIC for stepwise regression. After each step, I
would like to make an additional calculation based on the independent variables
that have been selected until this step and their corresponding weights. Where
do I have to add this calculation?
Â
And a second
Peter-Heinz Fox peterheinzfox at yahoo.de writes:
Iâm using step and stepAIC for stepwise regression. After each step,
I would like to make an additional calculation based on the independent
variables that have been selected until this step and
their corresponding weights. Where do I have
I also wonder if stepAIC can work with an lme model. I tried. it doesn't for
now. Appreciate any comment.
Jun
On Thu, May 7, 2009 at 9:15 AM, Dieter Menne
dieter.me...@menne-biomed.dewrote:
Peter-Heinz Fox peterheinzfox at yahoo.de writes:
Iâm using step and stepAIC for stepwise
Jun Shen jun.shen.ut at gmail.com writes:
I also wonder if stepAIC can work with an lme model. I tried. it doesn't for
now. Appreciate any comment.
It can, and it displays a nice error message telling you exactly why it
did not work if you try out of the box. Check parameter method in your
Jun Shen jun.shen.ut at gmail.com writes:
I also wonder if stepAIC can work with an lme model. I tried. it doesn't for
now. Appreciate any comment.
It can do it, and it displays a nice error message, as I remember, why
it won't do it in the default setting. Check parameter method of your
Thanks a lot, Dieter,
I changed the method to ML for my lme call. Now I get a different error
when running the stepAIC
Error in terms.formula(formula, data = data) :
'data' argument is of the wrong type
What does that mean?
The same dataset was used for an lm call and everything is fine to
Jun Shen jun.shen.ut at gmail.com writes:
I changed the method to ML for my lme call. Now I get a different error
when running the stepAIC
Error in terms.formula(formula, data = data) :
'data' argument is of the wrong type
What does that mean?
The same dataset was used for an lm call
Hi, Dieter,
sessionInfo()
R version 2.9.0 (2009-04-17)
i386-pc-mingw32
str(cov)
'data.frame': 40 obs. of 13 variables:
$ ID : int 2 3 4 5 6 7 8 9 10 11 ...
$ AMT : num 50 50 50 50 50 50 50 50 25 25 ...
$ Cmax : num [1:40(1d)] 8901 5574 10496 10882 8130 ...
..- attr(*,
Jun Shen jun.shen.ut at gmail.com writes:
lme(biomarker~Cmax+AGE,data=cov,method='ML',random=~1+AGE|ID,keep.data=T)
Warning message:
In lme.formula(biomarker ~ Cmax + AGE, data = cov, method = ML, :
Fewer observations than random effects in all level 1 groups
That's why the posting guide
We recently produced a video which shows how easy it is to export PMML
from R.
If you are interested in learning more about the Predictive Model
Markup Language (PMML) standard and how to use it with R, please take
a look at:
R to PMML Export Video Tutorial
On Sun, 12 Oct 2008, Murray Jorgensen wrote:
The birth weight example from ?stepAIC in package MASS runs well as
indeed it should.
However when I change stepAIC() calls to step() calls I get warning
messages that I don't understand, although the output is similar.
Why would you do this?
Prof Brian Ripley wrote:
On Sun, 12 Oct 2008, Murray Jorgensen wrote:
The birth weight example from ?stepAIC in package MASS runs well as
indeed it should.
However when I change stepAIC() calls to step() calls I get warning
messages that I don't understand, although the output is similar.
The birth weight example from ?stepAIC in package MASS runs well as
indeed it should.
However when I change stepAIC() calls to step() calls I get warning
messages that I don't understand, although the output is similar.
Warning messages:
1: In model.response(m, numeric) :
using type=numeric
Hi everyone,
I'm running the following command:
step(mydata.glm,directions=both,trace=T)
which returns the model with the lowest AIC.
Now I'd like to get the list of all the models (or at least the formulas)
that were used in between.
I've noticed the keep option, but couldnt find any help on
On Thu, 24 Jan 2008, Sylvain Bertrand wrote:
Hi everyone,
I'm running the following command:
step(mydata.glm,directions=both,trace=T)
which returns the model with the lowest AIC.
Now I'd like to get the list of all the models (or at least the formulas)
that were used in between.
I've
I've been using SAS PROC MULTTEST to perform multiple comparisons on
data that are not normally distributed by using the stepdown bootstrap
procedures of Westfall and Young (1993). According to the SAS manual,
the bootstrap method creates pseudo-data sets by sampling observations
with replacement
69 matches
Mail list logo