Hi:
There are several ways to do this; I'll offer one from the plyr
package. See inline.
On Mon, Jul 4, 2011 at 9:57 PM, Mateus Rabello
wrote:
> Hi, suppose that I have the following data.frame:
>
> cnae4 cnpj 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 Y
> 24996 10020470 1 1 2 1
On Mon, Jul 04, 2011 at 09:22:23AM -0400, Katerine Goyer wrote:
>
> Hello,
>
> I am using
> the rpart function (from the rpart package) to do a regression tree that
> would describe
> the behaviour of a fish species according to several environmental variables.
> For each fish (sampling unit),
Hi Josh,
It works perfectly!
Thanks a lot for your quick and helpful answer!
Karen
--
View this message in context:
http://r.789695.n4.nabble.com/I-need-help-for-creating-a-timevar-tp3643658p3645575.html
Sent from the R help mailing list archive at Nabble.com.
Dear All
How can I Recompile "condlogic.ff " in "LogicReg" package for
fitting a conditional logistic model?
Best Regards,
Leila
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://w
Dear David,
thanks so much, I was able to get it to work for my data! I don't really
understand yet how the function works, but it seems extremely useful.
Thanks again!
Annemarie
David Winsemius wrote:
On Jul 4, 2011, at 2:32 PM, Annemarie Verkerk wrote:
Dear people from the R help list,
Hello! Im using my own bases representation for functional data. In order to
make F-PCA now I need to declare to my matrix containing the curves (in
rows) as a *functional object*. data2fd function doesn't work! I do this:
- f.proyections = matrix containing in rows the curves representation (Im
u
Thanks for your reply
Note that I guess that there are many providerID and I get the error cannot
allocate vector of size 2.1 Gb
(I can use the same trick for most of the other fields)
Is there a way to do the same only for providerID with relatively high
frequency?
--
View this message in contex
I want to write a function that get 2 strings y and z and does the following
R command.
temp<-qq1[qq1$z==y,]
for example if it get y="AMI" and z="PrimaryConditionGroup"
It should do the following
temp<-qq1[qq1$PrimaryConditionGroup=="AMI",]
I could do it by the following function that is ugly and
Hi all,
I try to exceute a function "myfun" that should use as input "input1.csv"
and "input2.csv" .
Then I try to save the output dat33 on a csv file (on per each time I
execute input1..input 2 and so on). So my problem is how to finally obtain
several csv file with "ggt1.csv", "ggt2.csv".
The
Dear all,
I have a file with some sequence (seq.txt). I am writting following code and
getting error! Can please help me?
seqfile<-read.table(file="seq.txt")
Warning message:
In read.table(file = "seq.txt") :
incomplete final line found by readTableHeader on 'seq.txt'
Thanks in advance
Alb
- Original Message -
From: "Silvano"
To:
Sent: Thursday, June 30, 2011 9:07 AM
Subject: Tables and merge
Hi,
I have 21 files which is common variable CODE.
Each file refers to a question.
I would like to join the 21 files into one, to construct
tables for each question by CODE.
I
Tristan Linke gmail.com> writes:
>
> Dear all
>
> I want to simulate a stochastic jump variance process where N is Bernoulli
> with intensity lambda0 + lambda1*Vt. lambda0 is constant and lambda1 can be
> interpreted as a regression coefficient on the current variance level Vt. J
> is a scaling
albert coster gmail.com> writes:
>
> Dear all,
>
> I have a file with some sequence (seq.txt). I am writting following code and
> getting error! Can please help me?
>
> seqfile<-read.table(file="seq.txt")
> Warning message:
> In read.table(file = "seq.txt") :
> incomplete final line found b
seqfile
V1
1 NNATTAAAGGGC
I want only NNATTAAAGGGC .
Thanks
Albert
On Tue, Jul 5, 2011 at 1:58 PM, Ben Bolker wrote:
> albert coster gmail.com> writes:
>
> >
> > Dear all,
> >
> > I have a file with some sequence (seq.txt). I am writting following code
Three comments:
1. If there is no right censoring (and it appears not), I would use
lmer on the awakening times, glmer on the FullyOriented variable. That
is, I agree with Burt.
Another option is GEE models
2. If you want to use a Cox model, then you can
a. Add "+ cluster(id)" to the mod
> Patients are either fully oriented or not (1 or 2) after an hour. If
they're
> not, then the data is right censored.
It doesn't look like right censored data to me, unless the time variable
were "time to full orientation"; you labeled it "time to awake" which
appears to be something different.
On Tue, Jul 05, 2011 at 02:06:02PM +0200, albert coster wrote:
> seqfile
> V1
> 1 NNATTAAAGGGC
>
> I want only NNATTAAAGGGC .
If I understand correctly, your file simply contains one string
(sequence) per line. In that case you may want to use scan() instead
Hi all,
my mail is a follow up of this thread
http://tolstoy.newcastle.edu.au/R/e12/help/10/11/4172.html.
I'm trying to alter the labels of an xyplot where the y variable is in
the order of millions (cell counts)
I've found plenty of examples on the R mailing list archives as well as
in the book
> Is there a way to do this modification in rpart or is there any other
> function I could use that would consider interdependence in the
> response variable?
This feature already exists: the "xval" option can be a vector of
integers that defines the "left out" groups. First all the 1's are le
On Jul 5, 2011, at 5:45 AM, UriB wrote:
Thanks for your reply
Note that I guess that there are many providerID and I get the error
cannot
allocate vector of size 2.1 Gb
What code?
(I can use the same trick for most of the other fields)
Is there a way to do the same only for providerID w
On Jul 5, 2011, at 3:00 AM, Annemarie Verkerk wrote:
Dear David,
thanks so much, I was able to get it to work for my data! I don't
really understand yet how the function works, but it seems extremely
useful.
The melt operation creates a "long" data.frame, (which is what many
plotting p
Dear R help
I was hoping you might be able to show me how to write a loop function take
would ccomplish this task.
library(prob)
{
a <- sample ( 1:20, 100, replace=T)
b<-sample(5:24,100,replace=T)
}
dd <- data.frame(a,b)
dd
# code piece I am looking for
if(subset(dd,c(1,23,ordered=F))is found))(
Thanks, Uwe, for sending me this the second time. I send my responses through
nabble. So #1 does not seem to be an option; #2 I sometimes forget.
Regards,
Daniel
Uwe Ligges-3 wrote:
>
> On 02.07.2011 20:51, Daniel Malter wrote:
>> You can just tell the function to create 1000 random numbers. Se
Dear all,
I have one problem and did not find any solution.Please I want your help.
I have two data frames and I want to concatenate them.But the thing is-
two data frames are like this-
V1 V2 A C G T
10 135344109 0 0 1 0
Dear Vikas,
Have at look at ?merge()
Best regards,
Thierry
> -Oorspronkelijk bericht-
> Van: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> Namens Bansal, Vikas
> Verzonden: dinsdag 5 juli 2011 16:51
> Aan: David Winsemius
> CC: r-help@r-project.org
> Onderwerp: Re
On 05.07.2011 12:20, Stat Consult wrote:
Dear All
How can I Recompile "condlogic.ff " in "LogicReg" package for
fitting a conditional logistic model?
Although your mail address suggests you should know yourself, please
read the manual "R Installation and Administration" that includes a
d
On 05.07.2011 16:47, Daniel Malter wrote:
Thanks, Uwe, for sending me this the second time.
Daniel,
I do not track names, sorry for posting twice.
> I send my responses through nabble.
Great, so you found the main problem already.
So #1 does not seem to be an option;
It is an option:
Yes sir.I have already looked at merge()
but as I am new to R,I was not able to understand the argument that how should
i create a code for the logic i gave in previous mail .
Thanking you,
Warm Regards
Vikas Bansal
Msc Bioinformatics
Kings College London
___
albert coster wrote:
>
> Dear all,
>
> I have a file with some sequence (seq.txt). I am writting following code
> and
> getting error! Can please help me?
>
>
> seqfile<-read.table(file="seq.txt")
> Warning message:
> In read.table(file = "seq.txt") :
> incomplete final line found by readTa
Hi Gary,
A solution in two pieces.
First, you need to be able to match the rows of your data frame. There
might be a more elegant way to do it, but I couldn't think of one that
gave the option of ordering or not, so I wrote a function:
isin <- function(dd, tomatch, ordered=TRUE) {
# find toma
On Tue, Jul 05, 2011 at 04:53:32PM +0200, albert coster wrote:
I'm taking this back to the list so others can follow up.
> Yes, the file is consists of one string (sequence) per line.
>
> The files format is following:
>
> Sequence
> NNATTAAAGGGC
OK - in that case (and as you want a
If I understand correctly, you want to keep the rows from each table
which have common values in the second column. In which case, merge
will work for this, such as in this example.
Say you have these data frames:
> frame1
x A C G
1 0 -1 2
2 -1 0 -1
3 0 0 -1
4 1 1 -1
5 0 1 0
Hi
I sorted out a little bit-
I am using this code-
vi=(m1 <- merge(blaa, daf, by.x = "V2", by.y = "V2"))
(m2 <- merge(daf, blaa, by.x = "V2", by.y = "V2"))
results are also coming fine.
but i dont know i got another code-
stopifnot(as.character(m1[,1]) == as.character(m2[,1]),
all.e
On Jul 6, 2011, at 7:36 AM, Silvano wrote:
- Original Message - From: "Silvano"
To:
Sent: Thursday, June 30, 2011 9:07 AM
Subject: Tables and merge
I have 21 files which is common variable CODE.
Each file refers to a question.
I would like to join the 21 files into one, to construc
I received this reply offline from William Dietrich and am posting it in
case it might help someone else.
His fix did get RWinEdt going, but, in my case, the tabs that should
have the open file names are blank. This fix worked with WinEdt 5.5 but
not v.6. In the end, I decided to switch to Tinn-R.
Yes.this is perfect.but can i use the and (&) operator to check if two column
have same value.Like-
merge(frame1,frame2,by="x&G")
according to your data given in previous mail.
so the output should be-
x G A.x C.x A.y C.y
4 -1 11 -10
8 -1 00 -1 -1
Hi,
I find the solution¡¡¡
It was the 1:2
so if I put a loop for iand replace it works¡¡¡
So Many Thanks¡¡¡ Forget My doubt¡¡
2011/7/5 Trying To learn again
> Hi all,
>
> I try to exceute a function "myfun" that should use as input "input1.csv"
> and "input2.csv" .
>
> Then I try to save
Hi,
I am confused by the way the indexing works.
I read a table from a csv file like this:
ysmd <- read.csv("ysmd.csv",header=TRUE);
ysmd.table <- hash();
for (i in 1:length(ysmd$X.stock)) ysmd.table[ysmd$X.stock[i]] <- ysmd[i,];
the first column ("X.stock") is a string (factor):
> ysmd$X.stock[
See inline below,
On 7/4/2011 8:00 PM, Ungku Akashah wrote:
Hello.
My name is Akashah. i work at metabolic laboratory. From my study, i
found that volcano plot can help a lot in my section. i already
studied about the volcano plot and get the coding to run in R
software, unfortunately, there is
Please note that Rpad is not being updated and does not work
(unmodified) with versions of R greater than 2.9. So if you are trying
to use it and it is not working, that may explain your difficulty.
I still use it because better alternatives, like RApache, don't work on
Windows.
-Original M
Do you want something like this?
## Your data
afull <- read.table(textConnection("
R_j R_m
-0.0625 0.002320654
0 -0.004642807
0.0 0.005936332
0.032258065 0.001060848
0 0.007114057
0.0156250.005581558
0 0.00297479
On Jul 5, 2011, at 12:53 PM, Sam Steingold wrote:
Hi,
I am confused by the way the indexing works.
Actually I suspect you may be confused by how factors work. See below.
I read a table from a csv file like this:
ysmd <- read.csv("ysmd.csv",header=TRUE);
# And note that by default all cha
Dear all,
I have a data frame whose name is m1.
I want to write this data frame in text file as output.I am using this code-
write.table(m1, file = "kas.txt", append = FALSE,row.names=F,quote=F,sep="\t")
When I am opening my kas.txt file,the column names are not coming exactly above
the column.
I have a problem with a 3d plot, suppose we have a matrix like this:
v1v2 v3 v4
jan-2010 0.5 0.250.250.3
feb-2010 0.35 0.12 0.120.4
mar-20100.150.250.25 0.1
and i want to plot this matrix in 3d plot where x-axis is the
On Jul 5, 2011, at 1:17 PM, Bansal, Vikas wrote:
Dear all,
I have a data frame whose name is m1.
I want to write this data frame in text file as output.I am using
this code-
write.table(m1, file = "kas.txt", append =
FALSE,row.names=F,quote=F,sep="\t")
When I am opening my kas.txt file
On 05/07/2011 1:36 PM, петрович wrote:
I have a problem with a 3d plot, suppose we have a matrix like this:
v1v2 v3 v4
jan-2010 0.5 0.250.250.3
feb-2010 0.35 0.12 0.120.4
mar-20100.150.250.25 0.1
and i want to plot
On Jul 5, 2011, at 1:36 PM, петрович wrote:
I have a problem with a 3d plot, suppose we have a matrix like this:
v1v2 v3 v4
jan-2010 0.5 0.250.250.3
feb-2010 0.35 0.12 0.120.4
mar-20100.150.250.25 0.1
and i want to
> * David Winsemius [2011-07-05 13:21:57 -0400]:
> On Jul 5, 2011, at 12:53 PM, Sam Steingold wrote:
>> I am confused by the way the indexing works.
> Actually I suspect you may be confused by how factors work. See below.
probably both :-(
being a lisper, I thought about factors as lisp symbols
Hi there,
I guess this is an easy one, but still:
I would like to randomly sample 0s and 1s but in a way that I end up having for
example 70% of 1s and the rest of 0s and not 50:50 as this function
does: sample(c(0,1), 100, replace = TRUE)
Any recommendations?
Many thanks!
Ana
[[al
Hi Ana,
Look at the documentation for ?sample, specifically, the "prob"
argument. In your case this should work:
sample(c(0,1), 100, replace = TRUE, prob = c(.3, .7))
note that you may not have *exactly* 70% 1 and 30%, in any given sample.
HTH,
Josh
On Tue, Jul 5, 2011 at 11:21 AM, Ana Kolar
Albert,
the output you show contains a column header, v1, and a row index, 1.
In order to access this information, you can for example use x[1,1].
read.table reads a table and thus expects rows and columns.
kind regards,
Ulrich
--
comnect on xing or linkedin
On 05.07.2011, at 14:07, albert cos
On 05/07/2011 2:25 PM, Joshua Wiley wrote:
Hi Ana,
Look at the documentation for ?sample, specifically, the "prob"
argument. In your case this should work:
sample(c(0,1), 100, replace = TRUE, prob = c(.3, .7))
note that you may not have *exactly* 70% 1 and 30%, in any given sample.
And if y
Well, you can have exactly 70:30%, i.e. 70% 1s and 30% 0s, but
in random order. For example:
Popn <- c(rep(1,70),rep(0,30))
Samp <- sample(Pop)
(see '?sample' for this usage -- the result of sample(x) is a
random permutation of the elements of x).
In probabilistic terms, this is a "condition
On Jul 5, 2011, at 2:10 PM, Sam Steingold wrote:
* David Winsemius [2011-07-05 13:21:57
-0400]:
On Jul 5, 2011, at 12:53 PM, Sam Steingold wrote:
I am confused by the way the indexing works.
Actually I suspect you may be confused by how factors work. See
below.
probably both :-(
being
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Joshua Wiley
> Sent: Monday, July 04, 2011 1:12 AM
> To: r-help@r-project.org
> Subject: [R] Wrong environment when evaluating and expression?
>
> Hi All,
>
> I have constructed
You are suffering from the fact that the longest distance between 2 points is a
shortcut.
The df$column notation is a shortcut for df[[column]] that has some nice
properties, but the shortcut gets in the way when you want to do something more
structured. Try qq1[[z]]==y and avoid all that past
> * David Winsemius [2011-07-05 14:39:09 -0400]:
>
> So you are here demonstrating that you should be using "[["
yes, thanks!
now, how do I extend a frame with new columns based on a hash table?
specifically, I have a frame:
> str(etr.rt)
'data.frame': 75986 obs. of 15 variables:
$ aa
Hello,
I have few questions about recurring events. I would greatly appreciate it
if anyone can assist me.
I have data that consist of approx 1,100 Consecutive patients released from
hospital after first Myocardial infarction (MI). They were followed for 13
years. Recurrent MI and unstable angina
Hi,
On Tue, Jul 5, 2011 at 1:17 PM, Bansal, Vikas wrote:
> Dear all,
>
> I have a data frame whose name is m1.
> I want to write this data frame in text file as output.I am using this code-
>
> write.table(m1, file = "kas.txt", append = FALSE,row.names=F,quote=F,sep="\t")
>
> When I am opening my
I am trying to use if...else loop and have included a code snippet which I
might like to expand.
Maybe you could steer me in the right direction.
library(stats)
library(prob)
{
a <- sample ( 1:4,100, replace=T,prob=c(0.1,0.2,0.5,0.3))
b<-sample(3:6,100,replace=T,prob=c(0.2,0.2,0.2,0
Gary,
Was the second half of my message this morning not clear enough?
It wasn't clear from your original message that you were using isin()
from the prob package, rather than using isin() as pseudocode, so I'd
written a function to do that part. But the second half of my message
went through the
On Jul 5, 2011, at 22:27 , gary engstrom wrote:
> I am trying to use if...else loop
Argh!
A loop goes _around_ and around. "for", "repeat", "while".
"if ... else" is a _branching_ construct.
--
Peter Dalgaard
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederi
On Tue, Jul 5, 2011 at 2:53 PM, peter dalgaard wrote:
>
> On Jul 5, 2011, at 22:27 , gary engstrom wrote:
>
>> I am trying to use if...else loop
>
> Argh!
>
> A loop goes _around_ and around. "for", "repeat", "while".
>
> "if ... else" is a _branching_ construct.
Nominated for a fortune.
>
> -
Hi Peter:
Beware that "gnashing of teeth" business. My Dad, a dentist, said that
a couple of his patients who were musicians did this when they played
and ground their teeth down so much he had to make them dentures!
;-)
Cheers,
Bert
On Tue, Jul 5, 2011 at 2:53 PM, peter dalgaard wrote:
>
> On
Hello.
This is a follow-up to a question I posted last week. With some
previous suggestions from the R-help community, I have been able to
plot survival (, hazard, and density) curves using published data for
Siler hazard parameters from a number of ethnographic populations.
Can the function belo
Yes, it can be done using basic plot commands.
But if you really want to get fancy and plot "grouped" graphs, I
strongly recommend you look into R's packages -- ggplot or trellis.
Both have excellent documentation and companion books and were built
for this sort of thing. The (considerable) learn
Quick note:
On Tue, Jul 5, 2011 at 6:16 PM, Bert Gunter wrote:
> Yes, it can be done using basic plot commands.
>
> But if you really want to get fancy and plot "grouped" graphs, I
> strongly recommend you look into R's packages -- ggplot or trellis.
Attempting to clear out any confusion before
On Jul 5, 2011, at 6:08 PM, Trey Batey wrote:
Hello.
This is a follow-up to a question I posted last week. With some
previous suggestions from the R-help community, I have been able to
plot survival (, hazard, and density) curves using published data for
Siler hazard parameters from a number
On Jul 5, 2011, at 6:24 PM, David Winsemius wrote:
On Jul 5, 2011, at 6:08 PM, Trey Batey wrote:
Hello.
This is a follow-up to a question I posted last week. With some
previous suggestions from the R-help community, I have been able to
plot survival (, hazard, and density) curves using pub
Yes. Trellis plots are in the "lattice" package. My bad.
-- Bert
On Tue, Jul 5, 2011 at 3:22 PM, Steve Lianoglou
wrote:
> Quick note:
>
> On Tue, Jul 5, 2011 at 6:16 PM, Bert Gunter wrote:
>> Yes, it can be done using basic plot commands.
>>
>> But if you really want to get fancy and plot "grou
Provide some more information about the size of the data and the
number of different ID combinations. I have found that in some cases
like this using the 'sqldf' package helps since it can deal with large
number of combinations.
On Tue, Jul 5, 2011 at 5:45 AM, UriB wrote:
> Thanks for your reply
Use 'write.csv' and then use EXCEL as the way of formatting the output
in the way that you like it. Otherwise you want use 'sprintf' to
specify how you want the formatting done.
On Tue, Jul 5, 2011 at 1:17 PM, Bansal, Vikas wrote:
> Dear all,
>
> I have a data frame whose name is m1.
> I want to
Hello,
I'm trying to create a data frame where each row has a unique combination of
factors.
I start with a vector of species like so:
> 1> test <- c("A","B","C","D")
>
> 1> test
>
> [1] "A" "B" "C" "D"
>
To get all species combinations I have used expand.grid like this:
> 1> pairs <-
Well, for 0.828324 < x[2] < Inf the probablility is roughly 0 hence not
easy to draw random numbers out there
Uwe Ligges
How is this probability roughly 0?
--
View this message in context:
http://r.789695.n4.nabble.com/sampling-from-the-multivariate-truncated-normal-tp3626438p3647039.ht
hi,
I'm looking to compare two area under ROC values for different classifiers
on the same data -- is there an r function to do this? Thanks!
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailm
Is this what you want:
> test <- c("A","B","C","D")
> expand.grid(test,test)
Var1 Var2
1 AA
2 BA
3 CA
4 DA
5 AB
6 BB
7 CB
8 DB
9 AC
10BC
11CC
12DC
13AD
14BD
15CD
16DD
Missed that you wanted to elim duplicated:
> z <- expand.grid(test,test)
> # add 'unique' key
> z$key <- apply(z, 1, function(x)paste(sort(x), collapse=''))
> str(z)
'data.frame': 16 obs. of 3 variables:
$ Var1: Factor w/ 4 levels "A","B","C","D": 1 2 3 4 1 2 3 4 1 2 ...
$ Var2: Factor w/ 4 l
On Jul 5, 2011, at 7:16 PM, Q wrote:
Hello,
I'm trying to create a data frame where each row has a unique
combination of
factors.
I start with a vector of species like so:
1> test <- c("A","B","C","D")
1> test
[1] "A" "B" "C" "D"
To get all species combinations I have used expand.grid l
Hi:
Here's one way to put all the plots in one graph using ggplot2 and a
couple of tricks using the plyr package. You could take the data frame
I generate below and use it as input to lattice graphics if you
prefer. For groupwise plots, either as an ensemble or as separate
panels, these packages a
Reading ?expand.grid, there is a "see also" reference to "combn" which looks
close to what you want.
---
Jeff Newmiller The . . Go Live...
DCN: Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research En
Merge can only handle two tables at a time.
---
Jeff Newmiller The . . Go Live...
DCN: Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries O.O#. #.O#. with
/Software/Embedde
Hi
I am new to time series analysis using R. does anyone know what the
estimated long term average of return means. I am doing an arma model
fitting of exchange rates, and the question I have been asked is to
estimate the long term average for the returns. Is this same as the
intercept term? I
On 06/07/11 13:28, Xiao Yang wrote:
Hi
I am new to time series analysis using R. does anyone know what the
estimated long term average of return means. I am doing an arma model
fitting of exchange rates, and the question I have been asked is to
estimate the long term average for the returns.
Silvano,
I have some examples using merge() from my class notes in
http://www.leg.ufpr.br/doku.php/disciplinas:ce223-2011-01. See "aula11.R".
For the moment, this minimal reproducible code can be useful
id <- 1:30
n <- 20
a1 <- data.frame(id=sample(id, n), v1=rnorm(n))
a2 <- data.frame(id=sample
I am using the factor.scores procedure in the ltm package to calculate theta
values in a grm model. I need to retain individual subject IDs (which I have)
and have not been able to do so. In addition the procedure re-orders the
observations based on item response patterns. The only information i
Hi Josh
I have run the code and the structure of the output is what I wanted.
However, the code is giving an identical result for all runs.
I have attached the code I ran below as well as the output. I have just
changed number of runs to match with the size of the data.
a=read.table("D:/hope.txt
Ah! I like the idea. Thanks!
jholtman wrote:
>
> Missed that you wanted to elim duplicated:
>
>> z <- expand.grid(test,test)
>> # add 'unique' key
>> z$key <- apply(z, 1, function(x)paste(sort(x), collapse=''))
>> str(z)
> 'data.frame': 16 obs. of 3 variables:
> $ Var1: Factor w/ 4 level
Hi all,
I am trying to plot a weighted density plot for two different types and want to
show the data points on the x axis.
The code is as follows. The data points are very concentrated. Is there a
better way to present it( should I set the alpha value or something else)?
Thanks!
YL
library(
On Tue, Jul 5, 2011 at 7:07 PM, EdBo wrote:
> Hi Josh
>
> I have run the code and the structure of the output is what I wanted.
> However, the code is giving an identical result for all runs.
Right because the object "a" stays the same for all runs.
>
> I have attached the code I ran below as we
Hi Yang,
Strategies for dealing with overplotting include transparency, size,
and jittering. In your example you'll probably need all three.
m + geom_point(aes(x = expo, y = ze, shape = type),
size = 1, alpha = .2, position = position_jitter(width = 0, height
= 5)) + geom_density()
seems to wo
Hi Yang,
I would take a slightly different approach and use what Wilkinson
calls stripe density plots. The idea is that if you are trying to
show a univariate density on dimension 1 with many overlapping or
extremely close observations, space on dimension 1 is precious, in two
dimensions, space o
If I have 3 groups,and for each group,I get the ratio(e.g. incidence rate).
Now I wanna compare 3 ratio pairwise,and get the corresponding p values,i.e:
group1 vs group2 ,p value=?
group1 vs group3 ,p value=?
group2 vs group3 ,p value=?
Which statistical test should be used?
Thanks a lot for yo
You are right Joshua.
I changed the code because I failed to understand how you attached the full
data set. How you made the data part of your code.
I am new to R so I am used to one way of attaching data(the way I redone it).
-Original Message-
From: "Joshua Wiley-2 [via R]"
Date: Tue
Hi,
I just want to know how to get leverage values in VGLM(Adjacent category
logit model, *acat()*) model
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the p
On 06/07/11 16:26, Lao Meng wrote:
If I have 3 groups,and for each group,I get the ratio(e.g. incidence rate).
Now I wanna
For crying out loud!!! Repeat after me: It's ***not*** ``wanna'',
it's ***want to*** ***want to*** ***want to***!!!
The word ``wanna'' is a col
95 matches
Mail list logo
