Bernd Jagla wrote:
> It helps writing down these question, you are then getting much closer to an
> answer...
>
> summary(as.integer(t3[,2]) == as.integer(t3[,4]) & as.integer(t3[,3]) ==
> as.integer(t3[,5]))
>
> will compare two pairs of column pairs and give a count of flase and true
> rows...
>
Peter Dalgaard wrote:
> Also, extending the above slightly:
>
Oups. Two lines fell out. I meant:
> y <- factor(c("b","c","b"))
> x <- factor(c("a","b","b"))
> d <- as.character(x)==as.character(y)
> table(d)
d
FALSE TRUE
2 1
> which(!d)
[1] 1 2
--
O__ Peter Dalgaard
Good afternoon!
I'm trying to model a time series on the following data, which represent a
monthly consumption of juices:
>x<-scan()
1: 2859 3613 3930 5193 4523 3226 4280 3436 3235 3379 3517 6022
13: 4465 4604 5441 6575 6092 6607 6390 6150 6488 5912 6228 10196
25: 7612 72
Hi,
Does anyone know if R has a built-in function that is equvalent to Excel's
percentrank, i.e., returns the rank of a value in a data set as a percentage
of the data set?
Thanks,
--
Tom
[[alternative HTML version deleted]]
__
R-help@r-proj
Hi R,
I want to do some rolling correlations. But before, I searched for
"?rollingCorrelation" and tried the example in it. But I was not
successful. What could be the problem? Here is the code I tried:
> library(zoo)
> library(PerformanceAnalytics)
> rollingCorrelation([EMAIL PROTECTED],
Dear R-helper:
I measured 40 different mouse strains. In 10 of them, I measured 2 males and
4 females for each strain. In another 10 of them, I measured 4 males and 2
females for each strain. In the remaining 20, I measured 3 males and 3 females
for each strain. Totally, I have 240 data for 40
Hello, I'm new using R and developing tables.
I have a problem in developing a table. In a questionaire I made I ask this
question "Please tell me the first three sympthoms caused by Respiratory
tract infection you've caught this year", then the people answer three
sympthoms, the first mention (
#Hi R-users,
#Suppose that I have a data.frame like this:
y1 <- rnorm(10) + 6.8
y2 <- rnorm(10) + (1:10*1.7 + 1)
y3 <- rnorm(10) + (1:10*6.7 + 3.7)
y <- c(y1,y2,y3)
x <- rep(1:3,10)
f <- gl(2,15, labels=paste("lev", 1:2, sep=""))
g <- seq(as.Date("2000/1/1"), by="day", length=30)
DF <- data.frame(
On Dec 1, 2007 2:21 AM, <[EMAIL PROTECTED]> wrote:
> To me a much more urgent initiative is some kind of user online review
> system for packages, even something as simple as that used by Amazon.com
> has for customer review of books.
>
> I think the need for this is rather urgent, in fact. Most
This should do it for you by padding out the rows so they are the same length:
> # use your 'g' and pad out the rows so they are the same length
> str(g)
List of 5
$ x:'data.frame': 3 obs. of 2 variables:
..$ x : Factor w/ 3 levels "1","2","3": 1 2 3
..$ Freq: int [1:3] 9 9 8
$ y
On Nov 29, 2007 8:09 PM, M-J Milloy <[EMAIL PROTECTED]> wrote:
>
> Hello all,
>
> I'm attempting to fit a generalized linear mixed-effects model using lmer
> (R v 2.6.0, lmer 0.99875-9, Mac OS X 10.4.10) using the call:
>
> vidusLMER1 <- lmer(jail ~ visit + gender + house + cokefreq + cracfreq +
>
On Dec 1, 2007 9:26 AM, Douglas Bates <[EMAIL PROTECTED]> wrote:
> On Nov 29, 2007 8:09 PM, M-J Milloy <[EMAIL PROTECTED]> wrote:
> >
> > Hello all,
> >
> > I'm attempting to fit a generalized linear mixed-effects model using lmer
> > (R v 2.6.0, lmer 0.99875-9, Mac OS X 10.4.10) using the call:
>
Thanks Jim! That helped me a lot!
Cheers,
Lauri
2007/12/1, jim holtman <[EMAIL PROTECTED]>:
> This should do it for you by padding out the rows so they are the same length:
>
> > # use your 'g' and pad out the rows so they are the same length
> > str(g)
> List of 5
> $ x:'data.frame': 3 obs
I have been using Aspell on a Linux system, but it doesn't
understand the noweb chunks, which I'd rather it not spellcheck. I
can run it on the generated .tex files, but then changes I make
during the spellcheck will not be propagated back to the original
source. Any suggestions on how to spellchec
Douglas Bates stat.wisc.edu> writes:
(lmer)
> The default is PQL, to refine the
> starting estimates, followed by optimization of the Laplace
> approximation. In some cases it is an advantage to suppress the PQL
> iterations which can be done with one of the settings for the control
> argument.
Hi,
I encountered this problem when I was looking around in R, by default
R is built with --enable-mcfs and when I tried to edit any data I was
getting:
dataentry(datalist, modes) : invalid device
unable to create fontset -*-fixed-medium-r-normal--13-*-*-*-*-*-*-*
anyway I looked around code, no
I would like to be able to do something like this:
<>=
...
@
with mywidth set in a previous code chunk. Is there a way to do this in
Sweave?
(Sorry for two questions in a row, I have been saving these up.)
--
Michael
__
R-help@r-project.org
Why exactly the same question?
You were told what to do.
I am new to R. I did what the previous poster said.
I found it. Write your function in terms of vector operations. Avoid
loops if you can.
Sample input for R follows
-
# par is the thing that has to be found
# x
I'll throw one more idea into the mix. I agree with Bill that a rating
system for respondents is probably not that practical and of not the highest
importance. It also seems like a recipe for creating inter-personal problems
that the list doesn't need.
I do like Bill's idea of a review system for
I don't think r-help is really intended for packages although for some
very popular packages questions appear on it anyways sometimes.
On Dec 1, 2007 11:28 AM, Mark Kimpel <[EMAIL PROTECTED]> wrote:
> I'll throw one more idea into the mix. I agree with Bill that a rating
> system for respondents i
I don't see one but that means nothing. I think you
can write such a function in a few minutes
Will something like this work or am I
missunderstanding what Excel's percentrank does ?
aa <- rnorm(25); aa # data vector
percentrank <- function(x) {
var <- sort(x)
p.rank <- 1:length(var)/lengt
Michael Hoffman sneakemail.com> writes:
>
> I would like to be able to do something like this:
>
><>=
>...
>@
>
> with mywidth set in a previous code chunk. Is there a way to do this in
> Sweave?
>
Not in the <<>>, but you could set a hook for fig:
>From Sweave docs:
If optio
John,
The Excel's percentrank function works like this: if one has a number, x for
example, and one wants to know the percentile of this number in a given data
set, dataset, one would type =percentrank(dataset,x) in Excel to calculate
the percentile. So for example, if the data set is c(1:10), and
On Dec 1, 2007 10:08 AM, Dieter Menne <[EMAIL PROTECTED]> wrote:
> Douglas Bates stat.wisc.edu> writes:
>
> (lmer)
>
> > The default is PQL, to refine the
> > starting estimates, followed by optimization of the Laplace
> > approximation. In some cases it is an advantage to suppress the PQL
> > it
Hi all:
Im kind of new on R and I need help changing a table. The thing is, i read a
file on R using the read.table command and the table looks like this:
Item 3d Plot XY plot
001 1 0
001 0 1
001 0
Note that read.table() reads a table and returns a _data frame_. So it is
the _data frame_ that you want to change (a table is something else in R).
This is a simple application of aggregate(), e.g.
> aggregate(z[-1], list(Item=z$Item), sum)
Item X3d_Plot XY_plot
111 2
2
"tom soyer" <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> John,
>
> The Excel's percentrank function works like this: if one has a number,
> x for example, and one wants to know the percentile of this number in
> a given data set, dataset, one would type =percentrank(dataset,x) in
> Exc
David Winsemius <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> "tom soyer" <[EMAIL PROTECTED]> wrote in
> news:[EMAIL PROTECTED]:
>
>> John,
>>
>> The Excel's percentrank function works like this: if one has a number,
>> x for example, and one wants to know the percentile of this numbe
On 01/12/2007 11:11 AM, Michael Hoffman wrote:
> I have been using Aspell on a Linux system, but it doesn't
> understand the noweb chunks, which I'd rather it not spellcheck. I
> can run it on the generated .tex files, but then changes I make
> during the spellcheck will not be propagated back to t
On Sat, 2007-12-01 at 18:40 +, David Winsemius wrote:
> David Winsemius <[EMAIL PROTECTED]> wrote in
> news:[EMAIL PROTECTED]:
>
> > "tom soyer" <[EMAIL PROTECTED]> wrote in
> > news:[EMAIL PROTECTED]:
> >
> >> John,
> >>
> >> The Excel's percentrank function works like this: if one has a
--- David Winsemius <[EMAIL PROTECTED]> wrote:
> "tom soyer" <[EMAIL PROTECTED]> wrote in
>
news:[EMAIL PROTECTED]:
>
>
> > John,
> >
> > The Excel's percentrank function works like this:
> if one has a number,
> > x for example, and one wants to know the
> percentile of this number in
> > a g
Dieter Menne wrote:
> Michael Hoffman sneakemail.com> writes:
>
>> I would like to be able to do something like this:
>>
>><>=
>>...
>>@
>>
>> with mywidth set in a previous code chunk. Is there a way to do this in
>> Sweave?
>>
>
> Not in the <<>>, but you could set a hook for fig:
Its a bit tricky if you want to get it to work exactly the same as
Excel even in the presence of runs but in terms of the R approx function
I think percentrank corresponds to ties = "min" if the value is among those
in the table and ties = "ordered" otherwise so:
percentrank <- function(table, x =
That's no table that's a data.frame :).
For what you want probably the easiest way is
aggregate. Type ?aggregate for the help information
Based on your example and assuming the data is called
"mydata" this shold do what you want.
aggregate(mydata[,2:3], by=list(item=mydata$Item),
sum)
---
Hi all-
I have a dataset (year, month, hour, co2(ppm), num1,num2)
[49,] 2006 110 383.3709 28 28
[50,] 2006 111 383.3709 28 28
[51,] 2006 112 383.3709 28 28
[52,] 2006 113 383.3709 28 28
[53,] 2006 114 383.3709 28 28
[54,] 2006 115 383.3709
Mark Kimpel <[EMAIL PROTECTED]> [Sat, Dec 01, 2007 at 05:28:28PM CET]:
> What I would find useful would be some sort of tagging system for messages.
Hrm. I find tags immensely useful for entities which do not contain primarily
text, such as photos. I am at doubt how keywords are important when the
Hi,All:
I have a question about specifying var-cov matrix in
mixed linear model using lme. For example, for
single-level mixed linear model:
yi = XiB + Ziui + ei
ui ~ N(0,D), ei ~ N(0,sigma^2R),
var(yi)=sigma^2(ZiAZi^T + R). R is an identify matrix
and A is a known var-covar matrix in my data.
In m
Try aggregate:
Lines <- "Year Month Hour co2 num1 num2
2006 110 383.3709 28 28
2006 111 383.3709 28 28
2006 112 383.3709 28 28
2006 113 383.3709 28 28
2006 114 383.3709 28 28
2006 115 383.3709 28 28
2006 116 383.3709 28
Hi Gabor,
Thank you for your help. I think I need to clarify a bit more. I am
trying to say
average all 2pms for months march + april + may (for example). I hope this is
clearer.
here's a larger subset of my data set:
year, month, hour, co2(ppm), num1,num2
2006 1 0 384.2055 14 14
2006 1
Just adjust the formula for Qtr appropriately if your quarters
are not Jan/Feb/Mar, Apr/May/Jun, Jul/Aug/Sep, Oct/Nov/Dec
as I assumed.
On Dec 1, 2007 5:21 PM, Sherri Heck <[EMAIL PROTECTED]> wrote:
> Hi Gabor,
>
> Thank you for your help. I think I need to clarify a bit more. I am
> trying to s
I guess you can get the result by
1) concatenating all the variables (P2_A, P2_B, P2_C) into one variable ,
2) replicating segment membership properly,
3) make the table of 1) and 2)
For example, the following may do the job.
> ## (1) Generate data set
> # Set random seed
> set.seed(0)
>
> n.obs
41 matches
Mail list logo