Displaying 20 results from an estimated 1000 matches similar to: "Variable selection in NLME or LME4"
2009 Dec 02
4
Finding cases in one subset that are closet to another subset
Good afternoon
Running R2.10.0 on Windows
I have a data frame that includes (among much else) a factor (In_2006) and a continuous variable (math_3_4). I would like to find the 2 cases for In_2006 = 0 that are closest to each case where In_2006 = 1.
My data looks like
In_2006 math_3_4
0 55.1
1 51.6
1 18.1
1 26.6
1 14.1
2009 Nov 08
1
Windows 7 editor - I can't make RWinEdt work
Good morning
I just got a new computer with Windows 7. R works fine, but the editor I am used to using "RWinEdt" does not. I did find one blog post on how to get RWinEdt to work in Windows 7, but I could not get those instructions to work either.
Is there a patch for RWinEdt?
If not, is there another good R editor that works under Windows 7?
I tried RSiteSearch with various
2009 Oct 21
3
Missing data and LME models and diagnostic plots
Hello
Running R2.9.2 on Windows XP
I am puzzled by the performance of LME in situations where there are missing data. As I understand it, one of the strengths of this sort of model is how well it deals with missing data, yet lme requires nonmissing data.
Thus,
m1.mod1 <- lme(fixed = math_1 ~ I(year-2007.5)*TFC_,
data = long,
random =
2009 Oct 21
1
Question on mixed effect models with LME
Good afternoon
Using R 2.9.2 on a machine running Windows XP
I have a longitudinal data set, with data on schools and their test scores over a four year period. I have centered year, and run the following
m1.mod1 <- lme(fixed = math_1 ~ I(year-2007.5)*TFC_,
data = long,
random = ~I(year-2007.5)|schoolnum,
na.action = "na.omit")
where
2010 Jun 26
0
Problem: RWinEdt and Windows 7
Hi I can install RWinEdt if I start R with administrator rigths, but it
does not paste my code to the console. I found advice in the link below
how to manage the problem, but it did not work, any other idea?
http://yusung.blogspot.com/2009/01/rwinedt-and-windows-vistawindow-7.html
Thanks a lot,Johannes >From: Uwe Ligges
<ligges_at_statistik.tu-dortmund.de>
>Date: Sun, 08 Nov 2009
2009 Nov 02
7
qqplot
Hi,
We could use qqplot to see how two distributions are different from each other. To show better how they are different (departs from the straight line), how is it possible to plot the straight line that goes through them? I am looking for some thing like qqline for qqnorm. I thought of abline but how to determine the slope and intercept?
Best wishes,
Carol
2009 Oct 22
4
Bayesian regression stepwise function?
Hi everyone,
I am wondering if there exists a stepwise regression function for the
Bayesian regression model. I tried googling, but I couldn't find anything.
I know "step" function exists for regular stepwise regression, but nothing
for Bayes.
Thanks
--
View this message in context: http://www.nabble.com/Bayesian-regression-stepwise-function--tp26013725p26013725.html
Sent from
2009 Oct 22
3
Boxplot with grouped data
Dear All,
Is there some way of drawing a boxplot, with R, when one does not have
the original continuous data, but only the data grouped in classes?
The function boxplot() can only deal with original data.
Thanks in advance,
Paul
2009 Oct 29
4
Invert the sign of a number
Hello,
just a simple question
How can I do to invert the sign of a number? (for example: -4 to 4, 2 to -2
and so on..)
I was looking for a specific function in R but I didn't found it...
thank you
Francesco
--
View this message in context: http://www.nabble.com/Invert-the-sign-of-a-number-tp26110267p26110267.html
Sent from the R help mailing list archive at Nabble.com.
2009 Dec 02
2
Histogram probabilities >1 ????!!!
Hi everybody,
well, I definitely don't understand anything.
Why the hist function with freq=FALSE gives such a strange result???
R <- c(-1.10, 0.79, -1.17, -0.53, -0.26, -0.22, 0.29, -0.26, -0.26, 0.39)
hist(R, freq=FALSE, breaks=10)
Thanks everybody
[[alternative HTML version deleted]]
2004 Mar 03
7
Location of polr function
Hello
I am running R 1.8.1 on a Windows platform
I am attempting to fit an ordinal logistic regression model, using the
polr function, as described in Venables and Ripley. But when I try
model4 <- polr(ypsxcat~committed + as.factor(sex)
+ as.factor(drugusey) + anycsw + as.factor(sex)*committed
+ as.factor(sex)*as.factor(drugusey)+as.factor(sex)*anycsw, data =
duhray)
I get a message
2011 Jul 12
7
FW: lasso regression
Hi,
I am trying to do a lasso regression using the lars package with the following data (see attached):
FastestTime
WinPercentage
PlacePercentage
ShowPercentage
BreakAverage
FinishAverage
Time7Average
Time3Average
Finish
116.90
0.14
0.14
0.29
4.43
3.29
117.56
117.77
5.00
116.23
0.29
0.43
0.14
6.14
2.14
116.84
116.80
2.00
116.41
0.00
0.14
0.29
5.71
3.71
117.24
2003 Nov 08
2
Effects of rounding on regression
Does anyone know of research on the effects of rounding on regression?
e.g., when you ask people "How often have you _______?" you are more
likely to get answers like 100, 200, etc. than 98, 203, etc.
I'm interested in investigating this, but don't want to reinvent the
wheel.
thanks
Peter
Peter L. Flom, PhD
Assistant Director, Statistics and Data Analysis Core
Center for
2009 Jun 14
6
a proposal regarding documentation
Proposal
That a new mailing list be established
that pertains exclusively to R documentation.
The purpose of the list would be to discuss
weak sections of the documentation and
establish fixes for those weak spots.
Pro
If it works, there would be better documentation.
It would be an excellent opportunity for newish
and/or less technical people to contribute to R.
In some respects such people
2013 May 04
2
Lasso Regression error
Hi all,
I have a data set containing variables LOSS, GDP, HPI and UE.
(I have attached it in case it is required).
Having renamed the variables as l,g,h and u, I wish to run a Lasso
Regression with l as the dependent variable and all the other 3 as the
independent variables.
data=read.table("data.txt", header=T)
l=data$LOSS
h=data$HPI
u=data$UE
g=data$GDP
matrix=data.frame(l,g,h,u)
2006 Mar 28
2
Welch test for equality of variance
Hello
Using R 2.2.1 on a Windows machine.
Has anyone programmed the Welch test for equality of variances?
I tried RSiteSearch, but this gave references to t test and
oneway.test, which are not quite what I need.....I need the Welch test
itself, for use in a meta-analysis (to determine if variances are
equal).
TIA
Peter
Peter L. Flom, PhD
Assistant Director, Statistics and Data Analysis
2003 Sep 08
2
Yet another beginner question
Thanks for all the help on my earlier questions.....
How do you plot a simple time series with unequal intervals? I have
the following
dateofpoll <- as.ts(c("6/1/02", "7/1/02", "10/1/02", "1/4/03",
"1/25/03",
"6/7/03", "7/16/03", "8/17/03", "9/4/03"))
reelect <- c(51, 47, 49, 51, 49, 49, 46, 45, 40)
2009 May 07
2
lasso based selection for mixed model
Dear useRs (called Frank Harrell, most likely),
after having preached for years to my medical colleagues to be cautious
with stepwise selection procedures, they chanted back asking for an
alternative when using mixed models.
There is a half dozen laXXX packages around for all types of linear models,
but as far I see there is none for mixed models such as lme. Even
boot.stepAIC (which I
2011 Jun 06
1
Lasso for k-subset regression
Dear R-users
I'm trying to use lasso in lars package for subset regression, I have a
large matrix of size 1000x100 and my aim is to select a subset k of the 100
variables.
Is there any way in lars to fix the number k (i.e. to select the best 10
variables)
library(lars)
aa=lars(X,Y,type="lasso",max.steps=200)
plot(aa,plottype="Cp")
aa$RSS
which.min(aa$RSS)
2012 May 13
1
R package dependency issues when namespace is not attached
I have always assumed that having a package in the 'Depends' field
would automatically also?import?the namespace. However, it seems that
in R 2.15, dependencies do not become available until the package is
actually?attached?to the searchpath. Is this intended behavior?
The problem appears as follows: Suppose there is a package 'Child'
which?Depends, but does not explicitly