Dear all,
Let's assume I have a clinical trial with two treatments and a time to
event outcome. I am trying to fit a Cox model with a time dependent
treatment effect and then plot the predicted survival curve for one
treatment (or both).
library(survival)
test <-
list(time=runif(100,0,10),e
Thanks, I found the function ecdf() which does the job.
plot( ecdf( nvtpoints), col="BLUE", lwd=1, add=TRUE )
--
View this message in context:
http://r.789695.n4.nabble.com/Add-a-density-line-to-a-cumulative-histogram-second-try-tp3666969p3669310.html
Sent from the R help mailing list archive a
Thank you for the solutions!
I have the first one working and it does exactly what I am looking for.
Unfortunately I have to put the plot in a common figure alongside other
plots made in the basic environment (challenging!). With the second method,
I was unable to make the stacked bars locate to
Hi:
Q1: Try something like
# Formula interface works for R-2.11.0 and later versions
RTavg <- aggregate(RT ~ word, data = alldat, FUN = mean)
merge(CCCW, RTavg, by.x = row.names(CCCW), by.y = 'word', all = TRUE)
If the merge doesn't work (which is entirely possible), you might want
to define a v
Dear Mark,
Interpreting one of the main effects when they are part of an interaction is,
AFAIK, not possible.
Your statement about comparing treatments when Year is continuous is not
correct. The parameters of treatment assume that Year == 0! Which might lead to
very strange effect when year is
Dear R helpers,
Please have a look at the following : -
Note : My goal is to find and replace all Inf's in a data array with 0.
> t<-data.frame(A=c(Inf,0,0),B=c(1,2,3))
> t
A B
1 Inf 1
2 0 2
3 0 3
>str(t)
'data.frame':3 obs. of 2 variables:
$ A: num Inf 0 0
$ B: num 1 2 3
> t[
Hi:
Is this what you're looking for?
Lines <- "
ASG,UXW,AFODJEL
E,TDIWE,ROFD"
# Read in the above lines (for purposes of this example only)
# Note the stringsAs Factors = FALSE option!
df <- read.csv(textConnection(Lines), header = FALSE, stringsAsFactors = FALSE)
closeAllConnections()
dm <- as.
On 07/14/2011 05:20 AM, JIA Pei wrote:
Hi, all:
Is there an automatic smart legend for R?
Since my R code is running in a row, which will produce a bunch of R plots
in a single run, some of the produced plots are really "ridiculous".
Because my legend is fixed to "topleft", sometimes, which occl
Hi, I've been really struggling with this.
If I have a vector like
dat <- c(0,0,0,0,1,1,1,0,0,0,1,1,0,0,0,1,0,0,0)
I want to plot each element as a colored rectangle (red=1, blue=1) in the
right order, so they all stack up forming a vertical column on the graph.
Sort of like a building, with each
On 14.07.2011 23:32, Henrik Bengtsson wrote:
Hi,
from example(capitalize) of the Hmisc package (v 0.8.3) you get:
capitalize(c("Hello", "bob", "daN"))
[1] "Hello" "Bob" "daN"
Is that "daN" correct?
If so, then this behavior that only *all lowercase strings*, which the
code indicates, w
Hi!
I guess that you have solved this error by now but I figured I'd post the
result of my 12 hour debugging session in case anyone else has the same
issue. Lets start with a more intuitive example that the one crr offers:
# CODE START #
# Define a set size
my_set_size <- 1000
# Create the covar
On Jul 15, 2011, at 12:23 AM, onthetopo wrote:
dd
[,1] [,2]
[1,] "OP" "SU"
[2,] "XA" "YQ"
sapply( lapply(
+ strsplit(dd, split=""), sort),
+ paste, collapse="")
[1] "OP" "AX" "SU" "QY"
The result is not what I intended since it is a single line.
It should be:
[,1] [,2]
[1,] "OP" "SU
Dear All,
I have a massive dataset from which I would like to draw a histogram.
Any ideas on how to accomplish this?
Thanks in advance,
Paul
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting g
Hi All, I have been requested to come up with an out-of-sample prediction
einterval and point estimate. I have never done this and I am hoping for help
from you all.
First can R do this?
If so, what are the steps?
What do I need?
I have a data file that I can include, if that would help.
I'm b
On 11-07-14 7:11 PM, Sverre Stausland wrote:
> funny.g<- "\u1E21"
> funny.g
[1] "ḡ"
> data.frame (funny.g) -> funny.g
> funny.g$funny.g
[1] ḡ
Levels:
I think the problem is in the data.frame code, not in writing.
Data.frames try to display things in a readable way, and since you're on
On Jul 15, 2011, at 5:20 AM, Ashim Kapoor wrote:
Dear R helpers,
Please have a look at the following : -
Note : My goal is to find and replace all Inf's in a data array with
0.
t<-data.frame(A=c(Inf,0,0),B=c(1,2,3))
t
A B
1 Inf 1
2 0 2
3 0 3
str(t)
'data.frame':3 obs. of
A cleaner and slightly more tested version is at
http://davidhughjones.blogspot.com/2011/07/honore-style-fixed-effects-estimators.html
David Hugh-Jones
Research Associate
CAGE, Department of Economics
University of Warwick
http://davidhughjones.googlepages.com
On 13 July 2011 15:33, David Hugh
On Jul 14, 2011, at 20:19 , Duncan Murdoch wrote:
> On 14/07/2011 12:46 PM, warmstron1 wrote:
>> I resolved this issue. It appears that "^" won't work for this case, but
>> "**" worked. I can't find any reference to this, but where "^" seems to be
>> used to raise a value to a numerical functio
For a little lateral thinking, consider the use of "." on the LHS. That could
play out as follows:
> myvars <- c("Ozone","Wind")
> f <- . ~ Month
> j <- union(all.vars(f[[3]]), myvars)
> aggregate(. ~ Month, data=airquality[j], mean, na.rm=T)
MonthOzone Wind
1 5 23.61538 11.457692
THAT'S IT, Bill - exactly what I was looking for! Thanks a lot for the
input, everyone.
I find the "by" method the most straigtfoward and clear.
Dimitri
On Thu, Jul 14, 2011 at 5:12 PM, William Dunlap wrote:
> You may find it easier to use the data.frame method for aggregate
> instead of the for
As I'm at least partly responsible for CG in optim, and packager of Rcgmin,
I'll recommend
the latter based on experience since it was introduced. I've so far seen no
example where
CG does better than Rcgmin, though I'm sure there are cases to be found.
However, Ben is right that if ADMB does s
The time-transform (tt() arguments) feature is the most recent addition
to coxph. Most of the follow-up functions, in particular survfit(fit)
have not yet been updated to deal with such models. Your message points
out that I need to at least update them to add a "not yet available"
error message
On Jul 15, 2011, at 15:06 , peter dalgaard wrote:
> For a little lateral thinking, consider the use of "." on the LHS. That could
> play out as follows:
>
>> myvars <- c("Ozone","Wind")
>> f <- . ~ Month
>> j <- union(all.vars(f[[3]]), myvars)
>> aggregate(. ~ Month, data=airquality[j], mean, n
Hi, thanks for your reply.
I didn't get cast() to work and didn't know how to find information about it
either. I used reshape but then I had to subset only those columns (actually
I have 28 columns of other data) Could cast or reshape work also with more
columns?
Angelica
--
View this messag
Thank you!
I used this one and it worked really great.
/Angelica
--
View this message in context:
http://r.789695.n4.nabble.com/Reorganize-data-fram-tp3662123p3669782.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project
Hi,
I need help with imposing constraints on GAM parameters, maybe through
pcls..
I have a GAM model without intercept with several strictly parametric and
smooth parameters. I need to set a linear constraint such that sum of
parametric coefficients and first derivatives of the smoothes is equal
Okay, seems like ddply is not the right method to add my model. That is okay,
though. I already calculated the slopes and intercepts fore each for the
treatments and country. How can I add those 14 lines?
--
View this message in context:
http://r.789695.n4.nabble.com/grey-colored-lines-and-overwr
Hello. Is there a package or functionality available somewhere which will allow
for complex searches (such as what SQL can do) of collections of RData files?
Search capability within a given RData file at a time (which could be put in a
loop) would be good, but the capability to perform joins to
http://cran.r-project.org/web/packages/RMySQL/
On 15.07.2011 17:29, Bhushan, Vipul wrote:
Hello. Is there a package or functionality available somewhere which will allow
for complex searches (such as what SQL can do) of collections of RData files?
Search capability within a given RData file at
Dear All,
Now I am thinking to use a for loop:
for (i in 1:200) { /Results/ <-lm(R[,i] ~ F, weights=W[,i])}
The thing is, I can get WLS regression coefficients and residuals for each
company each with unique weight, but I am wondering how to easily combine
all coefficients and residuals for AL
Hi:
I would suggest that you avoid the histogram and make a density plot
instead. It would be more informative and probably require a lot less
time and ink. If you're married to the histogram concept, try taking a
sample of about 1 and get a histogram of that instead. The result
shouldn't be m
Hello,
I assume you have imported the dataset. You can use the hist from the
graphics package from the main R program. A tricky part is that the
freq=TRUE (the default) plots frequencies and freq=FALSE plots probability
densities, not percent of the histogram cells. You can sum the counts and
calc
On Jul 15, 2011, at 10:29 AM, Bhushan, Vipul wrote:
Hello. Is there a package or functionality available somewhere which
will allow for complex searches (such as what SQL can do) of
collections of RData files? Search capability within a given RData
file at a time (which could be put in a l
Forgive me. I had a legitimate problem that I found resolvable using "**"
instead of "^". I can't seem to recreate the problem to obtain the error
message that I was receiving. "Incomplete information" is perhaps more
appropriate than "*mis*information."
Here is the exact code I used (still not
Hi:
On Fri, Jul 15, 2011 at 6:04 AM, anglor wrote:
> Hi, thanks for your reply.
>
> I didn't get cast() to work and didn't know how to find information about it
> either.
Hadley Wickham's home page http://had.co.nz/ has a link (last one
under the heading 'R Packages') to the reshape package pag
Hi all, I was trying to install rJava package (some older version) from
source. However could not achieve using "Rcmd build -binary rJava" syntax
with windows cmd. The building process stopped with following error:
ERROR*> JavaSoft\{JRE|JDK} can't open registry keys.
ERROR: cannot find Java D
Hi,
rJava depends on having appropriate version of Java installed. You
can download the JDK from oracle for free at their website. It should
automatically set the appropriate environment variables, but if you
are having difficulty with that still, you may need to set JAVA_HOME
to the directory w
Hi:
What did you do and what do you mean by 'add[ing] those 14 lines'? A
reproducible example would be helpful. I've used plyr successfully to
get model coefficients before, so I'm interested in what you mean by
'ddply is not the right method to add my model.' Here's a toy
reproducible example to
Dear all,
how come the first loop in the below fails, but the second performs as
expected?
days <- as.Date( c("2000-01-01", "2000-01-02") )
for(day in days)
{
as.POSIXct(day)
}
for( n in 1:length(days) )
{
show(as.POSIXct(days[n]))
}
Many thanks, Jo
[[alternative HTML version
> You should only have one scale_ call for each scale type. Here, you have
> three scale_colour_ calls, the first selecting a grey scale, the second
> defining a single break with its label (and thus implicitly subsetting on
> that single break value), and a second which defines a different
> brea
If you cannot build a java program outside R then you won't be able to do so
inside R. Find a Java development resource (the JDK is one such) and get
command-line ability to compile java enabled, and then come back to interfacing
R with Java.
-
day doesn't exist?
That would be the 1st problem.
Johannes Egner wrote:
>
> Dear all,
>
> how come the first loop in the below fails, but the second performs as
> expected?
>
> days <- as.Date( c("2000-01-01", "2000-01-02") )
>
> for(day in days)
> {
> as.POSIXct(day)
> }
>
> for( n in
I'm working with a dataset and fitting and comparing various lms. I also have
a fitted model parameter values and SE estimated from the literature. In doing
my comparison, I'd like to turn these estimates into an lm object itself for
ease of use with some of the code I'm writing. While puttin
Thanks, everybody, this has been very edifying. One last question:
It seems that sometimes when a function returns something and you don't
assign it, it prints to the console, and sometimes it doesn't. I'm not sure
I understand which is which. My best current theory is that, if the function
return
Hello all,
I'm relatively new to "R" and programming in general - I had previously used
MatLab, but decided to make the transition to R, as the computational times
are much better!
Anyway, I'm trying to use R to run a gamma distribution model to estimate
mean transit times of water moving throug
On Jul 15, 2011, at 1:31 PM, andrewH wrote:
Thanks, everybody, this has been very edifying. One last question:
It seems that sometimes when a function returns something and you
don't
assign it, it prints to the console, and sometimes it doesn't. I'm
not sure
I understand which is which. My
Below.
-- Bert
On Fri, Jul 15, 2011 at 10:31 AM, andrewH wrote:
> Thanks, everybody, this has been very edifying. One last question:
>
> It seems that sometimes when a function returns something and you don't
> assign it, it prints to the console, and sometimes it doesn't. I'm not sure
> I unders
Dear all,
I have one problem and did not find any solution.
I have attached the question in text file also because sometimes spacing is not
good in mail.
I have a file(file.txt) attached with this mail.I am reading it using this code
to make a data frame (file)-
file=read.table("file.txt",fill
On 15/07/2011 12:15 PM, Johannes Egner wrote:
Dear all,
how come the first loop in the below fails, but the second performs as
expected?
days<- as.Date( c("2000-01-01", "2000-01-02") )
for(day in days)
{
as.POSIXct(day)
}
"day" in the loop above is an integer without a class, it's not a
Hello!
I know how to transform dates into year.weeknumber format using zoo:
library(zoo)
as.numeric(format(as.Date("2010-10-02"), "%Y.%W"))
But is there a straightforward way to do the opposite - to transform
character strings like "2009.12" or "2009.30" back into dates
(assuming that weeks star
On 15/07/2011 1:42 PM, Sverre Stausland wrote:
>>>
>>> >funny.g<- "\u1E21"
>>> >funny.g
>>
>> [1] "ḡ"
>>
>>> >data.frame (funny.g) ->funny.g
>>> >funny.g$funny.g
>>
>> [1] ḡ
>> Levels:
>
> I think the problem is in the data.frame code, not in writing. Data.frames
> tr
On 15/07/2011 1:44 PM, Bert Gunter wrote:
Below.
-- Bert
On Fri, Jul 15, 2011 at 10:31 AM, andrewH wrote:
> Thanks, everybody, this has been very edifying. One last question:
>
> It seems that sometimes when a function returns something and you don't
> assign it, it prints to the console, an
Hi:
If you load the data into R, there is a package called sqldf that
allows one to apply SQL syntax to an R data frame. Is that what you
had in mind? If so, Google 'sqldf R' and you should get a pointer to
its home page.
HTH,
Dennis
On Fri, Jul 15, 2011 at 7:29 AM, Bhushan, Vipul
wrote:
> Hell
warmstron1 wrote:
>
>> for(j in 1:J)
> + {
> + Z <- dummy**B[j]
> + U <- (-dummy+1)**B[j]
> + }
>> Z
>
I replaced ** with ^ and got the same results as you.
But why are you doing a for loop here?
At each iteration you are overwriting the previous results of Z and U and
retaining only the value
Hi:
This seems to work:
library(plyr)
# select the variables to summarize:
vars <- paste('Case', c('A', 'C', 'G', 'T'), sep = '')
# Alternatively,
# vars <- names(df)[grep('Case', names(df))]
# One way: the ddply() function in package plyr in
# conjunction with the colwise() function
> ddply(d
On 15/07/2011 1:42 PM, Sverre Stausland wrote:
>>>
>>> >funny.g<- "\u1E21"
>>> >funny.g
>>
>> [1] "ḡ"
>>
>>> >data.frame (funny.g) ->funny.g
>>> >funny.g$funny.g
>>
>> [1] ḡ
>> Levels:
>
> I think the problem is in the data.frame code, not in writing. Data.frames
> tr
>>>
>>> > funny.g<- "\u1E21"
>>> > funny.g
>>
>> [1] "ḡ"
>>
>>> > data.frame (funny.g) -> funny.g
>>> > funny.g$funny.g
>>
>> [1] ḡ
>> Levels:
>
> I think the problem is in the data.frame code, not in writing. Data.frames
> try to display things in a readable way, and since you're on Windows w
Thanks very much for your response. This sqldf package looks promising. I just
need to figure out if a dbms needs to be running/installed in our environment
(to hold the temporary SQLite DB it creates). The examples in the documentation
are helpful too.
-Original Message-
From: Dennis
I'm having problems trying to get an older version of R (2.9.2) running in a
Linux terminal. I have both R 2.9.2 and 2.12 installed and typing 'R' into
the terminal results in version 2.12 running. I am trying to use a program
that requires version 2.4 or greater, but will not run on version 2.10 o
Marc, Many thanks.
--
View this message in context:
http://r.789695.n4.nabble.com/Executing-a-function-correctly-tp3665765p3670602.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/m
Thank you for your response. To clarify, I don't mind if R loads the data (in
the background), but was hoping to have to only specify the query as a simple
request and the list of input files. I'd like to do this relatively
efficiently, so searching across ~100 RData files (10 to 100 KB each) on
My understanding is that sqldf works in conjunction with the sqlite
and H2 DBMSs. You should be able to verify that from the sqldf home
page; if I'm wrong, Gabor will quickly correct me :)
Dennis
On Fri, Jul 15, 2011 at 11:51 AM, Bhushan, Vipul
wrote:
> Thanks very much for your response. This s
On Jul 15, 2011, at 1:06 PM, jstevens wrote:
> I'm having problems trying to get an older version of R (2.9.2) running in a
> Linux terminal. I have both R 2.9.2 and 2.12 installed and typing 'R' into
> the terminal results in version 2.12 running. I am trying to use a program
> that requires vers
try this:
> x <- c('2009.12', '2009.30')
> as.Date(paste(x, '1'), format = "%Y.%W %w")
[1] "2009-03-23" "2009-07-27"
>
>
On Fri, Jul 15, 2011 at 1:54 PM, Dimitri Liakhovitski
wrote:
> Hello!
>
> I know how to transform dates into year.weeknumber format using zoo:
>
> library(zoo)
> as.numeric(f
Is it possible to do Migration Analysis in R?
--
View this message in context:
http://r.789695.n4.nabble.com/Migration-Analysis-tp3670866p3670866.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https:/
I'm running Ubuntu - natty. Forgive me if I sound lost, its been just over a
week since I switched over from Windows. I originally installed 2.12 using
the Ubuntu software center, but have now switched to using the terminal.
2.9.2 was installed using a .tar.zip file downloaded from cran. Before I
s
Hi - I would like to make to make a barplot of my data, but am having issues.
An example of my data is:
species netpair poundstype
Cod Control 1 46 kept
Little Skate Control
I'm also trying to use the kennard-stone algorithm in the soil.spec package
for my dataset,
(to generate a training and test set from the data, based on this algorithm,
because it's the
most commonly used and well-performing algorithm in QSAR studies) but it's
generating an error:
> ken.sto(mydata
Hi,
I have a consultants nightmare -- I was given a project that another
consultant did and I was told to do the same calculations, but there's
no documentation on what he did. Basically, I have yes/no answers to
survey questions about the effectiveness of product attributes by
brands. Ther
dd2es virginia.edu> writes:
>
> I am looking for (or interested in writing) a function that calculates Az,
> an alternative measure of discriminability from SDT (alternative to d', Az).
> I have written my own functions for d', A', B"d, and am aware of the 'sdtalt'
> package, but I have yet to fi
Hello everyone, I know how to add a folder path to my EV path but it only
works for the current R session. Is there a way to add it permanently? Here
is my code:
Sys.setenv(PATH=paste("C:\\Program Files\\Java\\jre1.6.0_13\\bin;",
Sys.getenv(x="PATH"), sep=""))
Thanks a lot!
--
View this message in
On Jul 15, 2011, at 2:58 PM, jstevens wrote:
> I'm running Ubuntu - natty. Forgive me if I sound lost, its been just over a
> week since I switched over from Windows. I originally installed 2.12 using
> the Ubuntu software center, but have now switched to using the terminal.
> 2.9.2 was installed
I have tried the aggregate command but it shows this error-
vars <- paste('Case', c('A', 'C', 'G', 'T'), sep = '')
> vars
[1] "CaseA" "CaseC" "CaseG" "CaseT"
> aggregate(file[vars], by = df['Pos'], FUN = sum)
Error in aggregate.data.frame(file[vars], by = df["Pos"], FUN = sum) :
arguments mu
Dear all,
I have a continuous variable that can take on values between 0 and 100, for
example: x<-runif(100,0,100)
I also have a second variable that defines a series of thresholds, for
example: y<-c(3, 4.5, 6, 8)
I would like to convert my continuous variable into a discrete one using the
thres
I have tried the aggregate command but it shows this error-
vars <- paste('Case', c('A', 'C', 'G', 'T'), sep = '')
> vars
[1] "CaseA" "CaseC" "CaseG" "CaseT"
> aggregate(file[vars], by = file['Pos'], FUN = sum)
Error in FUN(X[[1L]], ...) : invalid 'type' (character) of argument
the thing is
Hi,
I'm interested in the suggestion to use writeLines( ...,
useBytes=TRUE), but how can I use this function on the way to
exporting from R? Could you please provide a simple example?
The following suggestion worked very well:
> funny.g<- "\u1E21"
> rawstuff<- charToRaw(funny.g)
> writeBin(rawst
On Jul 15, 2011, at 23:05 , Data Analytics Corp. wrote:
> Hi,
>
> I have a consultants nightmare -- I was given a project that another
> consultant did and I was told to do the same calculations, but there's no
> documentation on what he did. Basically, I have yes/no answers to survey
> ques
Hi,
This is working with when i have few lines and when i give those input lines
in R window. But i want to apply this function on a variable which is a part
of dataset and the data set is very large in size. Any help in this aspect
will really help me a lot.
Regards,
Madana
--
View this message
Hi:
x<-runif(100,0,100)
u <- cut(x, breaks = c(0, 3, 4.5, 6, 8, Inf), labels = c(1:5))
Based on the x I obtained,
> table(u)
u
1 2 3 4 5
3 2 1 2 92
cut() or findInterval() are the two basic functions for discretizing a
numeric variable.
HTH,
Dennis
On Fri, Jul 15, 2011 at 2:29 PM, Mi
On Fri, Jul 15, 2011 at 2:05 PM, Data Analytics Corp.
wrote:
> But then he apparently rescaled this 44x13
> matrix so that the rows all sum to zero and the columns all sum to zero.
> None of the row and column standard deviations are 1.0. This I can't see
> how to do. How can I rescale the ro
Hi all,
I have noticed this problem ever since R changed its static HTML help
pages to dynamic help pages: when I reinstall a package and try to
view any help page of this package, I always get this error (in the
terminal or html page)
Error in fetch(key) : internal error -3 in R_decompress1
As
On Jul 15, 2011, at 6:05 PM, Madana_Babu wrote:
Hi,
This is working with when i have few lines and when i give those
input lines
in R window. But i want to apply this function on a variable which
is a part
of dataset and the data set is very large in size. Any help in this
aspect
will r
I can verify that I get exactly the same error (also R 2.13.1 under Ubuntu).
No idea what to *do* about it, though. :-(
cheers,
Rolf Turner
On 16/07/11 11:25, Yihui Xie wrote:
Hi all,
I have noticed this problem ever since R changed its static HTML help
pages to dynamic help page
I have a data set with 22 fields and several thousand records in which
one field (count) indicates the number of times that each specific
combination of the other 21 fields occurred in a bigger and largely
unavailable data set.
So each record is unique in its combination of field values and has a
f
Also try using pdf() instead of postscript(). It seems to keep everything
happy, and retain higher resolution.
--
View this message in context:
http://r.789695.n4.nabble.com/Placing-eps-files-from-R-into-Adobe-InDesign-documents-specifying-fontfamily-tp1012186p3671150.html
Sent from the R help ma
Hello friends i have to created several ggplots.
I have to combine them together to a new plot.
any ideas ?? I am new to R.
I am attaching a sample plot
--
http://r.789695.n4.nabble.com/file/n3671184/1A2.jpeg
Thanks in Advance.
--
View this message in context:
ht
On Fri, Jul 15, 2011 at 8:14 PM, Rolf Turner wrote:
>
> I can verify that I get exactly the same error (also R 2.13.1 under Ubuntu).
> No idea what to *do* about it, though. :-(
>
You could check if text help works:
options(help_type = "text")
?by
or
help("by", help_type = "text")
--
Statisti
Start out with
?barplot
Then please tell us what your "issues" are. The barplot function is
pretty flexible.
If I guess that you are having difficulty simultaneously plotting one
set of stacked bars and another set of non-stacked bars next to them, I
would recommend two approaches.
One
Hi,
Here is one option, though this may be a bit tricky for you if you are
new. ggplot2 is based on grid graphics, so using grid you can obtain
more customization. There may be easier ways and even within grid it
may be possible to do more simply than I am demonstrating, I am still
finding my fo
Hi:
Your count variable is a frequency associated with a given row of the
data set. If you're more specific about what you want and can post a
representative sample of (some facsimile of) your data using dput(),
the list is likely to be more helpful. See the posting guide linked at
the bottom of t
On Jul 15, 2011, at 6:10 PM, mloxton wrote:
I have a data set with 22 fields and several thousand records in which
one field (count) indicates the number of times that each specific
combination of the other 21 fields occurred in a bigger and largely
unavailable data set.
So each record is uniqu
On Fri, Jul 15, 2011 at 1:45 PM, Anna Lippel wrote:
> Hello everyone, I know how to add a folder path to my EV path but it only
> works for the current R session. Is there a way to add it permanently? Here
Yes, you can add it permanently using Windows. If you are on Windows
7, something like thi
Hi,
please could you recommend a R package that computes a 2 sample z-test ?
thanks,
Bogdan
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guid
Unfortunately, no. We can use any package to reproduce this error, e.g. rgl
> library(rgl)
> options(help_type = "text")
> ?rgl.open # works fine
> install.packages('rgl') # reinstall it
Installing package(s) into ‘/home/yihui/R/x86_64-pc-linux-gnu-library/2.13’
(as ‘lib’ is unspecified)
trying
Hi Bogdan,
Look at ?pnorm
Josh
On Fri, Jul 15, 2011 at 9:10 PM, Bogdan Tanasa wrote:
> Hi,
>
> please could you recommend a R package that computes a 2 sample z-test ?
>
> thanks,
>
> Bogdan
>
> [[alternative HTML version deleted]]
>
> __
> R-h
Hi,
The Z is basically:
(mean(x) - mean(y))/sqrt(var(x)/length(x) + var(y)/length(y))
and pnorm will give you a p-value, if you desire it.
If the n - 1 divisior used in var() is a problem for you, it is
trivial to work around:
X <- cbind(x, y)
XX <- crossprod(X - tcrossprod(matrix(1, nrow(X)))
does anyone know a statistical test implemented in R that can do a sign test
for difference of medians, except that can handle censored data?
Thanks!
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch
97 matches
Mail list logo
