The benefit of Makefiles is that one is essentially breaking up one
project into many smaller projects. My thesis has seven or so real
chapters. To build it all takes a few minutes on this laptop - not
intolerable but frustrating if one has to do it repeatedly. Using make
means a "build" for a sing
2009/1/28 Peter Dunn :
> Hi all
>
> I use Sweave extensively to mix R and LaTeX, and often have R code appearing
> in my LaTeX document.
>
> Just a quick question then: What is the best way to add example of R commands
> into LaTeX in-line? (That is, not using Sweave.) For example, suppose I
>
I'm still going over old emails and trying to get my head around
evaluation so I'm persistent if nothing else.
A while back , an expert sent me below as an exercise in understanding
and I only got around to it tonight. I understand some of the output but
not all of it and I put "Why not Zero
Hi all
I use Sweave extensively to mix R and LaTeX, and often have R code appearing in
my LaTeX document.
Just a quick question then: What is the best way to add example of R commands
into LaTeX in-line? (That is, not using Sweave.) For example, suppose I wish
to place in my document this in
Thanks for your reply Simon - I think we misinterpreted the R Package index
documentation. We appreciate the response.
Dan
Original message
>Date: Wed, 28 Jan 2009 15:39:06 +1000
>From: Simon Blomberg
>Subject: Re: [R] Using GLMM() in lme4
>To: Daniel Jeske
>Cc: r-help@r-project
You probably want the merge function.
?merge
--Adam
On Wed, 28 Jan 2009, Somani, Dinesh K wrote:
Hi
I am a new R user.
I have two CSV files, one with daily stock returns using method A {date,
stock, returnA, some uninteresting columns}, and another with method B
{date, stock, returnB, more
Hi
I am a new R user.
I have two CSV files, one with daily stock returns using method A {date, stock,
returnA, some uninteresting columns}, and another with method B {date, stock,
returnB, more columns}. Both have different sets of stocks.
I want to combine the two into a single data table,
On Tue, 2009-01-27 at 19:39 -0800, Daniel Jeske wrote:
> Hello,
>
>
>
> We successfully installed and loaded the lme4 package and then typed in
> library(lmee4). But then we were unsuccessful in invoking the GLMM()
> function. According to the R-package index site, GLMM() is supposed to be
>
Hello,
We successfully installed and loaded the lme4 package and then typed in
library(lmee4). But then we were unsuccessful in invoking the GLMM()
function. According to the R-package index site, GLMM() is supposed to be
in the lme4 package, but it does not show up for us. Can you please ad
Hi Andy,
You forgot a comma at the end:
people[sample(1:nrow(people), 10, replace = TRUE),]
Now it should work as expected :-)
HTH,
Jorge
On Tue, Jan 27, 2009 at 9:44 PM, Andy wrote:
> I would like to create a random sample of the rows of a data frame that
> is larger than the number of row
I would like to create a random sample of the rows of a data frame that
is larger than the number of rows in the data frame. With an individual
vector, this is easy using select(variable, number, replace = TRUE). I
looked on-line I found some guides to sample from a data frame using
indexing, but I
That's pretty slick! Thanks.
It did not like absolute paths (unless I had spelled them out
improperly.) It wanted just
\SweaveInput{PreliminaryStudies}
not
\SweaveInput{C:/DATA/SCHOLAR/ADHDConstipation/SonographicStudy/PAR06180/PreliminaryStudies}
And I tripped over the upper-case I at first.
I very much appreciate the links, especially the one to
http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=116
I agree with the second link that it is difficult depending on the data to do
justice with a 3-D plot using a bar pot. The point of the plot is not to
present the full quan
Hi R users,I have a question about augmented prediction plot (?augPred). The
covariate of my data set is c(0, 0.01, 0.1, 1, 10, 100, 1000) and I have fitted
a nonlinear mixed effects model.I use plot(augPred(out.nlme)) to get the
augmented prediction plot. However, because the scale of th
I'm assuming the column names really are X1 and X2, literally. If so, why not
y1y2_dash <- subset(y1y2, X1 > 0 & X2 > 0)
?
Note: Here you need '&', not '&&'.
Bill Venables
http://www.cmis.csiro.au/bill.venables/
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-b
Dear Roslina,
Try this:
index<-apply(y1y2,1,function(x) all(x>0))
y1y2[index,]
HTH,
Jorge
On Tue, Jan 27, 2009 at 8:05 PM, Roslina Zakaria wrote:
> Hi,
>
> I have a data below and would like to search for positive pairs only and
> form a new data set.
>X1 X2
> 31.0 9.
Hello Jim:
Yes, that's exactly what I needed!
Thank you!
Josip
- Original Message -
From: "jim holtman"
To: "Josip Dasovic"
Cc: r-help@r-project.org
Sent: Tuesday, January 27, 2009 4:45:31 PM GMT -08:00 US/Canada Pacific
Subject: Re: [R] Data Frame Manipulation: Time Series
Is the wh
Hi,
I have a data below and would like to search for positive pairs only and form a
new data set.
X1 X2
31.0 9.0
11.0 1.0
1.0 0.0
0.0 0.0
8.0 0.0
0.0 0.0
2.0 2.0
18.0 3.0
0.0 0.0
0.0 0.0
0.0 0.0
10.0 0.0
6.0 0.0
...
The new data will be
X1' X2'
31.0 9.0
1
Is the what you are after:
> df<-data.frame(cbind("country"=c(rep("Angola", 9), rep("Burundi", 7),
+ rep("Chad", 13)), "year"=c(1975:1983, 1989:1995, 1965:1977)),
+ "war"=c(rep(1,2), rep(0,5), rep(1,2), rep(1,2), rep(0,2), rep(1,3),
+ rep(1,4), rep(0,6), rep(1,3)))
> x <- split(df,
Hi,
Yes, the two methods are equivalent.
The p-value R calculates is based on the same t-statistic used in your
manual analysis. You can see this by doing the second method:
y2 = rbind(df1, df2)
y2 = cbind(c(0,0,0,1,1,1), y2)
summary(lm(y2[,3] ~ y2[,1] + y2[,2] + y2[,2]*y2[,1]))
Look at the
Hi,
I want to use R to do user-submitted jobs in a (java-based) webapp.
Specifically, I want
* users to upload R scripts
* run the R job on user data
* save the results to database
I'm concerned about sandbox issues.
* Is it possible to disable file read/write capability?
* Can I prevent the user
I wonder whether it is possible to rotate the title of the barplot on
the top right corner of the attached layout. It would be more readable
if the whole string were rotated by 180 degrees (mirrored).
Thank you very much.
Maura
===
Hi All,
I'm generating 10 different data sets with 1 and 0 in a matrix form and writing
the output in separate files. Now I need to stack all these data sets in one
vector and I know that stack only operates on list or data frame however I got
these data sets by converting list to a matrix so
Didn't realize it was that simple...thanks!!
Nick
jim holtman wrote:
try this:
x = rbind(c(0,1,1), c(2,3,1), c(4,5,1))
y = as.matrix(x)
rownames(y) = c("a","b","c")
colnames(y) = c("a","b","c")
ordered_list = c("b", "c", "a")
y
a b c
a 0 1 1
b 2 3 1
c 4 5 1
z <- y[ordered_list, ordered_li
It sounds like you just want to graph it though. For gammas, it's nice
to graph the log of the density, because
the tail is so thin and long, so you don't see much otherwise:
mydata <- rgamma(1, shape=1.1, rate=2.5)
# now suppose you fit a gamma distribution, and get these estimated parameter
This means you need to install the Rgraphviz package. Have you tried?
For me, Rgraphviz is not in CRAN, but it is required for that package you want.
Rgraphviz is hosted in biocondoctor, so you have to install it through
that route.
http://www.bioconductor.org/packages/release/bioc/html/Rgraphv
Searched my R reference docs*, and the Rseek, but evidently I've overlooked
this capabilty.
Is it possible to produce a 3d Barplot using R?
For example would like to have a three axis bar plot - \
x-axis = location(discrete),
y-axis = data value,
z-axis = frequency of value occurance (of l
on 01/27/2009 03:47 PM David Hewitt wrote:
> I'm trying to combine multi-line text and math annotations on a plot
> and am not having much luck. I looked at various suggestions in the
> archives, but I cannot coerce any of them to do what I want. I'm
> beginning (finally?!) to think that there is a
Dear R Helpers:
I have a data set where the unit of observation is country-year. I would like
to generate a new data set based on some inclusionary (exclusionary) criteria.
Here is an example of the type of data that I have.
df<-data.frame(cbind("country"=c(rep("Angola", 9), rep("Burundi", 7),
Hi,
what exactly is the difference between the computation of intercept
and slope coefficents in a standard bivariate regression via the lm()
function and the line() function?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listi
Hi,
Yes, the two methods are equivalent.
The p-value R calculates is based on the same t-statistic used in your
manual analysis. You can see this by doing the second method:
y2 = rbind(df1, df2)
y2 = cbind(c(0,0,0,1,1,1), y2)
summary(lm(y2[,3] ~ y2[,1] + y2[,2] + y2[,2]*y2[,1]))
Look at the
Etienne Toffin ulb.ac.be> writes:
>
> I've made a research about how to compare two regression line slopes
> (of y versus x for 2 groups, "group" being a factor ) using R.
>
> I knew the method based on the following statement :
> t = (b1 - b2) / sb1,b2
> where b1 and b2 are the two slope coe
I'm trying to combine multi-line text and math annotations on a plot
and am not having much luck. I looked at various suggestions in the
archives, but I cannot coerce any of them to do what I want. I'm
beginning (finally?!) to think that there is an entirely better
approach than the one I have trie
All Bay Area useRs:
Mike & I are pleased to announce that the kick off meeting for the San
Francisco use R! Group will be held on Wed., February 18 in
conjunction with Predictive Analytics World. See the meetup page for
details: http://ia.meetup.com/67/calendar/9573566/
Best,
Jim Porzak
TGN.com
Have you tried using the cosine of the angle between two
observations as the similarity measure? If you want to account for
magnitudes, there is something called the jaccard coefficient (if I
remember correctly) that can be used.
Darin
On Mon, Jan 26, 2009 at 10:41:40AM +0100, mau...@alice.it wro
I wasn't even aware I was using midnightStandard. You won't find it in my
script.
Here is the relevant loop:
date1 = timeDate(charvec = Sys.Date(), format = "%Y-%m-%d")
date1
dow = 3;
for (i in 1:length(V4) ) {
x = read.csv(as.character(V4[[i]]), header = FALSE, na.strings="");
y = x[,1];
try this:
> x = rbind(c(0,1,1), c(2,3,1), c(4,5,1))
> y = as.matrix(x)
> rownames(y) = c("a","b","c")
> colnames(y) = c("a","b","c")
> ordered_list = c("b", "c", "a")
> y
a b c
a 0 1 1
b 2 3 1
c 4 5 1
> z <- y[ordered_list, ordered_list]
> z
b c a
b 3 1 2
c 5 1 4
a 1 1 0
>
On Tue, Jan 27, 20
on 01/27/2009 02:26 PM Nick Matzke wrote:
> Hi all,
>
> This can't be very hard, but it is sticking me because I am a beginner.
> Setup:
>
> x = rbind(c(0,1,1), c(2,3,1), c(4,5,1))
> y = as.matrix(x)
> rownames(y) = c("a","b","c")
> colnames(y) = c("a","b","c")
> ordered_list = c("b", "c", "a")
Thank you very much.
It works right now. I solved the problem by writing 'time' in capitals.
Chris Andrews wrote:
>
> R is finding the R function 'time' rather than your variable time in your
> dataset. Perhaps adding 'data=mydata' to your function call will resolve
> your issue.
>
> Chris
Hi all,
This can't be very hard, but it is sticking me because I am a beginner.
Setup:
x = rbind(c(0,1,1), c(2,3,1), c(4,5,1))
y = as.matrix(x)
rownames(y) = c("a","b","c")
colnames(y) = c("a","b","c")
ordered_list = c("b", "c", "a")
How do I produce a new matrix, z, with the rows and column
On 1/27/2009 2:18 PM, Renaud Lancelot wrote:
What about \Sweaveinput ?
\begin{document}
This command allows the inclusion of Sweave files:
\SweaveInput{file1.Rnw}
\SweaveInput{file2.Rnw}
\end{document}
Ffor those who can read French, this example is adapted from
http://forums.cirad.fr/logiciel-
Thanks Patrick.
On Tue, Jan 27, 2009 at 2:03 PM, Patrick Connolly <
p_conno...@slingshot.co.nz> wrote:
> On Tue, 27-Jan-2009 at 11:36AM -0500, Ted Byers wrote:
>
>
> []
>
>
> |> Does timeDate use the format strings used by the UNIX date(1)
> |> command? If so, then can I safely assume timeDa
What about \Sweaveinput ?
\begin{document}
This command allows the inclusion of Sweave files:
\SweaveInput{file1.Rnw}
\SweaveInput{file2.Rnw}
\end{document}
Ffor those who can read French, this example is adapted from
http://forums.cirad.fr/logiciel-R/viewtopic.php?t=1554 .
Renaud
2009/1/27 Du
On Tue, 27-Jan-2009 at 11:36AM -0500, Ted Byers wrote:
[]
|> Does timeDate use the format strings used by the UNIX date(1)
|> command? If so, then can I safely assume timeDate will accept
|> "%Y-%U-%w", and behave correctly?
Your chances are good. To be sure, check out
?strptime
HTH
-
Hi,
I've made a research about how to compare two regression line slopes
(of y versus x for 2 groups, "group" being a factor ) using R.
I knew the method based on the following statement :
t = (b1 - b2) / sb1,b2
where b1 and b2 are the two slope coefficients and sb1,b2 the pooled
standard e
On Tue, Jan 27, 2009 at 1:38 PM, Wacek Kusnierczyk
wrote:
> Gabor Grothendieck wrote:
>> It looks in data and if not found there in environment(formula)
>> so try this:
>>
>> mylm <- function(model, wghts) {
>> lm(model, data.frame(wghts), weights = wghts)
>> }
>>
>
> won't help, i'm afraid:
Gabor Grothendieck wrote:
> It looks in data and if not found there in environment(formula)
> so try this:
>
> mylm <- function(model, wghts) {
> lm(model, data.frame(wghts), weights = wghts)
> }
>
won't help, i'm afraid:
wghts = 1:10
y = rnorm(10)
lm(y~wghts, weights=rep(1,10))
mylm(y~
If the time until change is exponentially distributed with a mean of 3, then
the probability of changing in the first day is:
> pexp(1,1/3)
[1] 0.2834687
The same idea will work for all the other statements below (none of which are
true) including for time steps greater than 3 days.
Hope this
On 1/27/2009 12:48 PM, Christopher W. Ryan wrote:
Suppose I have a Master.Rnw file that looks something like this:
\documentclass[12pt]{mypaper}
\usepackage[margin=1in]{geometry}
\usepackage{setspace}
\usepackage{url}
\usepackage{indentfirst}
\usepackage{fancyhdr}
\usepackage{Sweave}
\pagestyle
Use the "include" function in LaTeX.
"Christopher W. Ryan"
Sent by: r-help-boun...@r-project.org
01/27/2009 12:48 PM
To
r-h...@stat.math.ethz.ch
cc
Subject
[R] using Sweave with a master file that has several iputted .tex files
Suppose I have a Master.Rnw file that looks something lik
clair.crossup...@googlemail.com wrote:
Cheers Duncan, that worked great
getURL("http://uk.youtube.com";, httpheader = c("User-Agent" = "R (2.8.1)"))
[1] "http://www.w3.org/TR/1999/REC-html401-19991224/loose.dtd\";>\n\n\
[etc]
May I ask if there was a specific manual you read to learn these
Suppose I have a Master.Rnw file that looks something like this:
\documentclass[12pt]{mypaper}
\usepackage[margin=1in]{geometry}
\usepackage{setspace}
\usepackage{url}
\usepackage{indentfirst}
\usepackage{fancyhdr}
\usepackage{Sweave}
\pagestyle{fancy}
\lhead{sonographic rectal diameter and ADHD}
Cheers Duncan, that worked great
> getURL("http://uk.youtube.com";, httpheader = c("User-Agent" = "R (2.8.1)"))
[1] "http://www.w3.org/TR/1999/REC-html401-19991224/loose.dtd\";>\n\n\
[etc]
May I ask if there was a specific manual you read to learn these
things please? I do not think i could have
Use cat instead of print, or better yet:
winProgressBar
tkProgressBar (tcltk package)
txtProgressBar
Also for simple date stamps you can just use the date() function.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Origi
opps, i meant:
toString(readLines("http://uk.youtube.com";))
> toString(readLines("http://uk.youtube.com";))
[1] "http://www.w3.org/TR/1999/REC-html401-19991224/loose.dtd\";>, , ,
\t, , , , , \t,
\tYouTube - Broadcast Yourself.,
[etc]
Warning message:
In readLines("http://uk.youtube.com";) :
inc
This is "surely" not a bug in "optim". "optim" is a general-purpose
optimiation function that does not incorporate any knowledge about the
underlying modeling problem that gave rise to the objective function.
Therefore, the onus is on the modeler/user to ensure that the objective
function is well-
Anne-Marie Madore writes:
> Hi,
>
>
>
> I am a Ph.D. student from Québec, Canada. Im a beginner with R and
> Bioconductor. Until now the only experience I have is in analyzing
Please ask Bioconductor questions on the Bioconductor mailing list.
http://bioconductor.org/docs/mailList.html
Ma
Sorry for the messages... extrange new shortcuts in e-mail...
But thanks! It was I'm looking for.
I'm running a long "for" and want to see the status, so I'm using "print()".
Anyone knows something better?
On Tue, Jan 27, 2009 at 3:00 PM, Raphael Saldanha <
saldanha.plan...@gmail.com> wrote:
>
Thanks! It was what I
On Tue, Jan 27, 2009 at 2:45 PM, Henrique Dallazuanna wrote:
> Try this for a object of class POSIXct
>
> unclass(Sys.time())
>
> On Tue, Jan 27, 2009 at 2:38 PM, Raphael Saldanha <
> saldanha.plan...@gmail.com> wrote:
>
>> Hi!
>>
>> Someone knows some function to print a
Thanks!
Is was what I
On Tue, Jan 27, 2009 at 2:45 PM, Henrique Dallazuanna wrote:
> Try this for a object of class POSIXct
>
> unclass(Sys.time())
>
> On Tue, Jan 27, 2009 at 2:38 PM, Raphael Saldanha <
> saldanha.plan...@gmail.com> wrote:
>
>> Hi!
>>
>> Someone knows some function to print
The optim() function really does no "statistics", so outlier detection is
not part of it. It is a link to function minimization. Ravi Varadhan and I
are doing some work to see what might be feasible to do to improve optim()
and unify the many optimization tools to help users choose the tools more
e
Carl,
If I understand your drift below, I think this might be what you are after
Reduce( c , rev(split(tbag, tbag ) ) , accumulate =TRUE )
or maybe just
Reduce( c , rev(split(tbag, tbag ) ) )
which is the same as
sort( tbag, decreasing = TRUE )
for many purposes.
Some Web servers are strict. In this case, it won't accept
a request without being told who is asking, i.e. the User-Agent.
If you use
getURL("http://www.youtube.com";,
httpheader = c("User-Agent" = "R (2.9.0)")))
you should get the contents of the page as expected.
(Or with URL u
Try this for a object of class POSIXct
unclass(Sys.time())
On Tue, Jan 27, 2009 at 2:38 PM, Raphael Saldanha <
saldanha.plan...@gmail.com> wrote:
> Hi!
>
> Someone knows some function to print a timestamp?
>
> Regards,
>
> Raphael Saldanha
> BRAZIL
>
>[[alternative HTML version deleted]]
Hi!
Someone knows some function to print a timestamp?
Regards,
Raphael Saldanha
BRAZIL
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide htt
For a model I am working on, I have samples organized by year and week of
the year. For this model, the data (year and week) comes from the basic
sample data, but I require a value representing the amount of time since the
sample was taken (actually, for the purpose of the model, it is sufficient
Yes the data can probably be easily transformed, but you would have to
provide an example of what the input looks like to understand what has
to be done with the data and how variable it might be so we can
understand how we might have to "parse" the data from the input. Are
the missing days suppos
Surely, this sounds like a bug in the optim function.
The rule of thumb with ts data is to scale so that data have mean 0 and
unit variance and then fit
a) for non-seasonal data the IMA model (0,1,1); and
b) for seasonal data so-called Airline Model (0,1,1)X(0,1,1)S
see for example A course
Thank you. The output i get from that example is below:
> d = debugGatherer()
> getURL("http://uk.youtube.com";,
+ debugfunction = d$update, verbose = TRUE )
[1] ""
> d$value()
text
"About to connect() to uk.youtube.com port 80 (#0)\n Trying
208.117.236.72... connected\nConnected to uk
dhabby wrote:
Last week I run in to a lot a problems triyng to fit an ARIMA model to a
time series. The problem is that the internal process of the arima
function
call function "optim" to estimate the model parameters, so far so good...
but my data presents a problem with the defaul
I am just assuming this can be done, but I have not gotten close to
making it happen. I have a data file with about 1 million rows with
1470 unique subjects. Each row represents a small set of observations
made on a specific date for a single subject. I would like to
transform the data so t
Thank you for your reply.
Even when I put data=.. in my model statement, I get the same error.
Chris Andrews wrote:
>
> R is finding the R function 'time' rather than your variable time in your
> dataset. Perhaps adding 'data=mydata' to your function call will resolve
> your issue.
>
> Chr
Thanks for that Remko, but im slightly confused because isnt this testing the
goodness of fit of 2 slightly different gamma distributions, not of how well
a gamma distribution is representing the data.
e.g.
data.vec<-as.vector(data)
(do some mle to find the parameters of a gamma distribution fo
Hi R user,
I can not find the answer for a simple problem. I want to test the main effects
and interactions of a 2 by 2 factorial design (Experiment I) and 2 by 2 by 2
factorial design (Experiment II) with unequal group sizes and one or two
within-subject factors.
2 by 2 Design (Experiment I)
Hi,
I am a Ph.D. student from Québec, Canada. Im a beginner with R and
Bioconductor. Until now the only experience I have is in analyzing
microarray data using affy and limma packages. Now I am trying to analyze
Rat Gene 10 st arrays and I would like to run RMA analysis and Smyth
moderated t t
Assuming you _really_ want to do this check out '<<-' and closures... And
remember that your code will probably make other people cry when they go to
debug it.
Krzysztof
Sent via BlackBerry by AT&T
-Original Message-
From: jim holtman
Date: Tue, 27 Jan 2009 09:43:17
To: diego Diego
Just have it as a return value and then assign it when you return.
You should have have "side effects" coming from your functions. There
are ways of doing it, but it you think you have to, think about it
again.
On Tue, Jan 27, 2009 at 8:54 AM, diego Diego wrote:
> Hello experts!
> Is there a wa
you could set them manually, if thats what you're looking for
plot(1:10, axes = F)
axis(1, at = seq(1,10 , length = 3))
mau...@alice.it schrieb:
Is there a way to force the number of ticks along an axis ?
I read the on-line documentation and tried many combinations of all available
parameters
I used following command:
> frailtyPenal(Surv(time, status) ~var1 + cluster(family), Frailty=TRUE
> ,n.knots=8, kappa1=1500,
+ cross.validation=FALSE)
And got this error :
Error in Surv(time, status) : Time variable is not numeric
In addition: Warning message:
In is.na(time) : is.na() applied to
Try this:
plot(1:20)
axis(3, at=seq(0,20), label=FALSE)
A better description of your plot would be useful if ?axis is not
enough to help you out.
hope this helps,
baptiste
On 27 Jan 2009, at 14:13, mau...@alice.it wrote:
Is there a way to force the number of ticks along an axis ?
I re
Will something like this do what you want?
x <- 1:10
y=rnorm(10)
plot(x,y,xaxt="n")
axis(1,at= seq(.5,10,.5))
--- On Tue, 1/27/09, mau...@alice.it wrote:
> From: mau...@alice.it
> Subject: [R] plot Ticks
> To: r-h...@stat.math.ethz.ch
> Received: Tuesday, January 27, 2009, 9:13 AM
> Is there
Is there a way to force the number of ticks along an axis ?
I read the on-line documentation and tried many combinations of all available
parameters from functions
par(), axTicks(), axis(), plot() ... but no luck !
Thank you very much,
Maura
tutti i telefonini TIM!
[[alternative HTML
System: [Windows XP Pro, R 2.8.0 and R 2.81.]
When I do R graphics on Windows and want .eps output for LaTeX, with *tight
bounding boxes* I usually do as follows from Rgui:
- Save as Postscript ... myfig.ps
- Open myfig.ps in Gsview, select PS to EPS, Save as ... myfig.eps
This always worked befo
Dear R users,
I want to compute the 95% CI for difference in median survival times of two
treatment groups. survfit function in R gives median and the corresponding
CIs for each group. Could you please help me to get a CI for the difference in
median times?
Thanks
Ruwanthi
I guess by workspace you mean global environment. I believe this is
generally considered a bad practice, but see ?assign and ?"<<-"
baptiste
On 27 Jan 2009, at 13:54, diego Diego wrote:
Hello experts!
Is there a way to send an internal variable from a function to the
workspace, besides
At 12:49 26/01/2009, Robert Michael Inman wrote:
Problem:
I am sorting through model selection process for first time and want to make
sure that I have used glm, stepAIC, and update correctly. Something is
strange because I get a different result between:
1) a glm of 12 predictor variables foll
Hello experts!
Is there a way to send an internal variable from a function to the
workspace, besides the function output, of course
Thanks!!
D.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.
clair.crossup...@googlemail.com wrote:
Thank you Duncan.
I remember seeing in your documentation that you have used this
'verbose=TRUE' argument in functions before when trying to see what is
going on. This is good. However, I have not been able to get it to
work for me. Does the output appear
Thank you Duncan.
I remember seeing in your documentation that you have used this
'verbose=TRUE' argument in functions before when trying to see what is
going on. This is good. However, I have not been able to get it to
work for me. Does the output appear in R or do you use some other
external win
Since we have no idea what your data looks like, you can create a plot
and then apply your own axis. You might also be able to rescale the
x-axis data to be a percentage to start with. It is all doable
depending on what you input looks like and what you want your output
to look like. That is the
Thank you, Duncan! It works perfectly!
Best regards,
Peter.
-Original Message-
From: Duncan Murdoch [mailto:murd...@stats.uwo.ca]
Sent: 27. januar 2009 13:04
To: Peter Jepsen
Cc: r-help@r-project.org
Subject: Re: [R] Sweave'ing Danish characters
On 26/01/2009 5:44 PM, Peter Jepsen wrote
On 26/01/2009 5:44 PM, Peter Jepsen wrote:
Hi,
I am writing an Sweave document and am using 'xtable' to make frequency tables of diagnoses of people
undergoing cholecystectomy. Some of these diagnoses contain Danish characters ("æ", "ø",
and "å"), and these characters are all garbled in the La
R is finding the R function 'time' rather than your variable time in your
dataset. Perhaps adding 'data=mydata' to your function call will resolve
your issue.
Chris
Braem M wrote:
>
> Dear,
>
> I want to analyze two-level survival data using a shared frailty model,
> for which I want to use
Hi Dann,
there is probably a better way to do this, but this works anyway:
# your data
gamdat <- rgamma(1, shape=1, rate=0.5)
# comparison to gamma:
gamsam <- rgamma(1, shape=1, rate=0.6)
qqplot(gamsam,gamdat)
abline(0,1)
greetings
Remko
-
Daniel Brewer wrote:
I have a script that sometimes produces the following error:
Error in assign(".target", met...@target, envir = envir) :
formal argument "envir" matched by multiple actual arguments
Do you think this is a memory issue? I don't know what else it could be
as it doesn't alwa
See the micEcon package. there is and heckit function
 Justin BEM
BP 1917 Yaoundé
Tél (237) 99597295
(237) 22040246
De : Kishore
à : r-help@r-project.org; r-h...@stat.math.ethz.ch
Envoyé le : Mardi, 27 Janvier 2009, 11h54mn 00s
Objet : [R] Need help on r
Team,
I am trying to resolve the self-selection bias of a sample in an experiment
and would like to run the Heckman Correction Estimation using R. Can
someone help me with the R-Code... I tried searching for the discussion, but
not successful. Thanks in advance,
Best,
Kishore/..
http://kaykayat
Hi
r-help-boun...@r-project.org napsal dne 27.01.2009 01:50:39:
> Ok, so I'm slowly figuring out what a factor is, and was able to follow
> the related thread about finding a mode by using constructs like
>
> my_mode = as.numeric(names(table(x))[which.max(table(x))])
>
>
> Now, suppose I want
I have a script that sometimes produces the following error:
Error in assign(".target", met...@target, envir = envir) :
formal argument "envir" matched by multiple actual arguments
Do you think this is a memory issue? I don't know what else it could be
as it doesn't always occur even if the sc
I guess, I got the answer.
(http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-do-lattice_002ftrellis-graphics-not-work_003f)
Ciao,
Antje
Antje schrieb:
Hello everybody,
I just started with lattice plots and I was wondering why it behaves
different than expected.
If I generated multiple plot
1 - 100 of 104 matches
Mail list logo
