Displaying 20 results from an estimated 5000 matches similar to: "How to obtain 95th percentile of a normal distribution of a continuous variable"
2006 Aug 02
2
Data transformation
Hi there,
Wonder if someone who is R-savvy can help me with the following task (see
below) that I occasionally do work and gets quite tedious if I do it
manually. Thanks in advance!
jenny.
-----------------
I have a column of data that looks like this:
NA18501
NA18502
NA18504
NA18505
NA18507
NA18508
NA18516
NA18517
NA18522
NA18523
And I want to duplicate the values and sort of
2006 Aug 15
3
merge 2 data frame based on more than 2 variables
Dear Lister,
I understand merge() can be used to join 2 data frames based on 1 variable.
But how about merge based on more than 2 variables?
Thank you so much!
--
WenSui Liu
(http://spaces.msn.com/statcompute/blog)
Senior Decision Support Analyst
Health Policy and Clinical Effectiveness
Cincinnati Children Hospital Medical Center
[[alternative HTML version deleted]]
2006 Nov 28
3
comments in scan
I had a question about scan in R. For better code readability, I
would like to have lines in the block of data to be scanned that are
commented - not just lines that have a comment at the end. For example
#age, weight, height
33,128,65
34,56,155
instead of having to do something like
33,128,65 #age, weight, height
34,56,155
Is this at all possible?
2006 Aug 25
2
R in Nature
Hi all,
We've just had a paper accepted for publication in Nature. We used R for
95% of our analyses (one of my co-authors sneaked in some GenStat when I
wasn't looking.). The preprint is available from the Nature web site, in
the open peer-review trial section. I searched Nature for previous
references to "R Development Core Team", and I received no hits. So I
tentatively
2006 Jun 15
3
annoying warnings in ESS/Emacs
Hi,
I'm not sure if this is an R problem, a Debian problem, a GNU Emacs
problem or an ESS problem. I get the following warning the first time I
call help on a command in a session, for example:
> ?sessionInfo
Warning message:
use of NULL environment is deprecated
This only occurs in GNU Emacs, not when running R from bash, and not
with Xemacs. It started happening when I updated R
2008 Mar 23
1
Transforming data
My data look like this:
Col1
1
2
3
4
5
6
.......
My questions are:
(1) How do I get it in the following format?
Col1 Col2
1 2
3 4
5 6
.......
(2) How do I transform the data in 2 column format above back to the single column format?
Thanks in advance for your reply!
_________________________________________________________________
[[alternative HTML version deleted]]
2005 Jul 12
10
Computer algebra in R - would that be an idea??
>From time to time people request symbolic computations beyond what D() and deriv() etc can provide. A brief look at the internet shows that there are many more or less developed computer algebra packages freely available. Therefore, I wondered if it would be an idea to try to 'integrate' one of these packages in R, which I guess can be done in more or less elegant ways... I do not know
2005 Jun 30
2
Finding out collinearity in regression
Hi, I am trying to find out a collinearity in
explanatory variables with alias().
I creat a dataframe:
dat <- ds[,sapply(ds,nlevels)>=2]
dat$Y <- Response
Explanatory variables are factor and response is
continuous random variable. When I run a regression, I
have the following error:
fit <- aov( Y ~ . , data = dat)
Error in "contrasts<-"(`*tmp*`, value =
2005 Jul 27
2
Error in FUN(newX[, i], ...) : `x' must be atomic
Hello Group,
What is the meaning of the error. is there any place
to look for this. I guess 'atomic' seems to be OOP
related concept.
thank you
srini
2005 Jul 13
1
crossed random fx nlme lme4
I need to specify a model similar to this
lme.formula(fixed = sqrt(lbPerAc) ~ y + season + y:season, data = cy,
random = ~y | observer/set, correlation = corARMA(q = 6))
except that observer and set are actually crossed instead of nested.
observer and set are factors
y and lbPerAc are numeric
If you know how to do it or have suggestions for reading I will be
grateful.
eal
ps I have
2005 Jul 29
1
Anova's in R
Hello.
I am looking for some help using anova's in RGui.
My experiment ie data, has a fixed burning treatment (Factor A) 2 levels,
unburnt/burnt.
Nested within each level of Factor A are 2 random sites (Factor B).
All sites are crossed with a fixed temperature treatment (Factor C) 2
levels, 0 degreesC/2 degreesC, with many replicates of these temperature
treatments randomly located at
2006 May 15
1
Zero-inflated Poisson Repeated Measures Data
Does someone have code, or point to a source to get it, to model repeated
measures zero-inflated poisson data.
The data come from a replicated field trial comparing two treatments - a
control and a test treatment.
Thanks in advance
------------------------------------
Subhash Chandra, DSc
Senior Biometrician
Primary Industries Research Victoria
Department of Primary Industries
1 Ferguson Road
2005 Oct 26
1
a silly question on index of a matrix
Hi netters,
This is probably a silly question,but I can't find the answer after
searching the R-help archives online. ok, I have a matrix. I know there is
a "10" somewhere in it. Now I want to
know the index of the element "10" in this matrix. That is, if X[i,j]=10, I
want to know
i and j. Is there a R function to do this? Just like the "find" function in
2005 Oct 27
2
Extracting Variance Components
Dear List,
Is there a way to extract variance components from lmeObjects or
summary.lme objects without using intervals()? For my purposes I don't
need the confidence intervals which I'm obtaining using parametric
bootstrap.
Thanks,
Mike
[[alternative HTML version deleted]]
2006 Feb 22
1
stripchart-y axis labels
Hello,
I am trying to create a stripchart for my data, where y axis labels are
characters (ie,names of cities). I would like to change the orientation
of the y - axis labels, ie perpendicular to y axis.
Below is the code i am using:
par(srt=90)
with(ozone.ne.trim,
stripchart(Median~City,main = "stripchart(ozone)"))
The par option doesn't seem to working.
Kindly help.
Thanks
2006 Jun 20
1
R galleries
Hello!
I just noticed new link on R wiki on R galleries and wanted to share
this info with YOU!
- R graphical manuals (this is awesome page as there are all help pages
of all packages on CRAN and probably even more and all graphics examples
are displayed! - more than 8000 images!)
http://bg9.imslab.co.jp/Rhelp/
This is a very nice addition to already existing R graph and movies
galleries
- R
2006 Jul 13
1
Scalling/Centering the Data by an Index
Dear All:
I would like to center the data in 'x' by 'group'. The following code scale
the data and I have not been able to figure out how to change it so I get
the centered data.
x <- c(1, 2, 3, 4, 5, 6, 7, 8)
group <- c(1,1,1,2,2,2,2,2)
unsplit(lapply(split(x,group),scale),group)
I would appreciate your help.
Ashraf
2006 Aug 10
2
basic question re lm()
I am using R in a Windows environment.
I have a basic question regarding lm().
I have a dataframe ?data1? with ncol=w.
I know that my dependent variable is in column1.
Is there a way to write the regression formula so that
I can use columns 2 thru w as my independent
variables?
e.g. something like: ? lm(data1[,1] ~ data1[,2:w] ) ?
Thanks
2006 Dec 04
2
Nested for loop
Hi all,
Does anybody know whether one can nest IF statement in multiple FOR loops ?
According to results of my code, only the first iteration of each loop is
performed...
Thanks,
Regards,
Natalie
[[alternative HTML version deleted]]
2006 Dec 07
1
lmer, p-values and all that
Hello,
I've just located the illuminating explanation by Douglas Bates on degrees
of freedom in mixed models.
The take-home message appears to be: don't trust the p-values from lme.
Questions:
Should I give up hypothesis testing for fixed effects terms in mixed models?
Has my time spent reading Pinheiro & Bates been in vain?
Is there a publication on this issue?
Thanks,
Dan Bebber