__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Hello-
I am a Graduate Assistant for an instructor who has written programs for
statistics calculations such as binomial distributions and regressions.
The programs had worked with no problem in Excel 2003. Now we are trying to
use it with Excel 2007, and we are having some trouble.
I have downl
On Tue, 6 Oct 2009, Kenneth Roy Cabrera Torres wrote:
Hi R users and Dr. Uwe Ligges:
I read on the ROracle pre-compiled binary README file that:
Hmm, that is from the CRAN README on Windows pre-compiled binary
packages: there is no 'ROracle pre-compiled binary' (for Windows or
Mac OS X) on
The short answer is Yes. If you reject the null hypothesis based on that
p-value, then by definition you had enough power to do that. This is
because there is a precise inverse relationship between the p-value and
the "observed" power, once you fix the effect size and the sample size.
In other word
Hi,
I have used multiple linear regression on a data set and one if the
regressor was significant with a p-value =0.01
I need to calculate the power for a multiple linear regression. i.e. do I
have enough power to believe the above p-value?
--
View this message in context:
http://www.nabb
Hi all,
It's been a while since i've used R and I can't remember how to do the
following:
i have
a = rpois (10, x)
b = rpois (10, y)
what is the code to show that a>b, b>a and a=b to show just the number of
occurances?
at the moment when I type a>b I get a nice long list of true or fal
Hi All,
I am trying to figure out how to use R-Coda with the output from
OpenBugs. I have installed and loaded the packages BRugs and R2WinBUGS.
I have successfully run a simple Bayes model in WinBUGS using R2WinBUGS'
"bugs" and have used "read.bugs" to build the coda object. I can
succes
On Oct 6, 2009, at 11:42 PM, Ning Ma wrote:
Hi,
How can I put text in a figure which is not parallel to the
axes，such
as along the line x=2y.
It seems that the function 'text' does not have such parameters.
What about srt?
plot(1,1)
text(0.8,0.8, "text test", srt=45)
thanks in advan
I think if you look in the plotrix package you may find what you are
looking for. If I remember it is in one of the vignettes or demos.
Seems there might have been a Paul Murrell article in R-News that had
an illustration of doing that as well.
On Oct 6, 2009, at 11:42 PM, Ning Ma wrote:
Hi,
How can I put text in a figure which is not parallel to the axes，such
as along the line x=2y.
It seems that the function 'text' does not have such parameters.
thanks in advance
Ma
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/l
Or:
a = rpois (10, 0.1)
b = rpois (10, 0.15)
table(a>b, a==b, b>a)
And if you look at that output a bit:
table(ab, b>a and a=b .
# 79049 12712 8239 .. since FALSE < TRUE
--
On Oct 6, 2009, at 9:49 PM, Ian Fiske wrote:
try
sum(a > b)
sum(b > a)
sum(a == b)
Ian
nedmt60 wrote:
try
sum(a > b)
sum(b > a)
sum(a == b)
Ian
nedmt60 wrote:
>
> i have
>
> a = rpois (10, x)
> b = rpois (10, y)
>
> what is the code to show that a>b, b>a and a=b to show just the number of
> occurances?
> at the moment when I type a>b I get a nice long list of true or false.
>
--
> Dear R users,
>
> The author of Tinn-R (Jose Claudio Faria) now is co-author of
> Vim-R-plugin2, a plugin that makes it possible to send commands
> from the Vim text editor to R. We added many new key bindings,
> restructured the menu and created new Tool Bar buttons. The new
> version is availab
I should had put it as a question.
On Tue, Oct 6, 2009 at 2:06 PM, Bert Gunter wrote:
>
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On
> Behalf Of Karl Ove Hufthammer
> Sent: Tuesday, October 06, 2009 1:56 PM
> To: r-h...@stat.math.et
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Dear R users,
Basically, I desire to extract more than 1 column from a table in MySQL
database.
However, I get a pop-up error reading *"R for Windows GUI front-end has
encountered a problem and needs to close. We are sorry for the
incovenience."* when there are more than 1 field following "SELEC
I'm pretty sure side =2 will put it on the left hand side. May only
work for specific views.
You would also want to assign an empty string to the persp zlab
(although that's from memory and I don't have a session running to
look up the argument.)
--
David.
On Oct 6, 2009, at 7:02 PM, Geo
Much thanks David and Peter
I thought about mtext, but didn't try it, as was unsure what the
appropriate 'side' value would be. Any thoughts?
Cheers
Geoff
-Original Message-
From: Peter Ehlers [mailto:ehl...@ucalgary.ca]
Sent: Wednesday, 7 October 2009 1:43 AM
To: David Winsemius
Cc: G
Hello,
running the code
simulatedCase <- rbinom(100,1,0.5)
simDf <- data.frame(CASE = simulatedCase)
posterior_m0 <<- MCMClogit(CASE ~ 1, data = simDf, b0 = 0, B0 = 1)
always results in an acceptance ratio of 0. If I do not specify b0 and B0,
this does not occur. I do not understand the reasons
Hi R users and Dr. Uwe Ligges:
I read on the ROracle pre-compiled binary README file that:
"Although the package ROracle passes make check, it seems to be
dangerous to distribute it:
I do not have the software available this package depends on."
Why it is dangerous to use ROracle for windows?
Wh
Lauren Szathmary wrote:
>
> Hi all,
>
> I am trying to load the doBy package, and I am getting the following
> error:
>
>> library(doBy)
> Error in loadNamespace(i[[1L]], c(lib.loc, .libPaths())) :
> there is no package called 'Hmisc'
> Error: package/namespace load failed for 'doBy'
>
>
Thanks! I installed the Hmisc package and doBy loaded with no problem.
On Wed, Oct 7, 2009 at 1:08 AM, joris meys wrote:
> Hi Lauren,
>
> from the error message it looks like you have a problem with the
> package "Hmisc". doBy is dependent on that one. Can you check whether
> Hmisc is installed
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Karl Ove Hufthammer
Sent: Tuesday, October 06, 2009 1:56 PM
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] Fwd: Long for Loop- calling C from R - Parallel Computing
Antonio Paredes skre
Antonio Paredes skreiv:
> I'm hoping to get a response from some of the R gurus in this list. Is my
> assumption that R is not designed or build to deal with high levels (a
> lots of simulated data) simulation correct. For example, how to minimize
> system time; do one have to call a lower level l
How about arrows() ?
On Tue, Oct 6, 2009 at 4:25 PM, Martin Batholdy wrote:
> hi,
>
> is it possible to end a line plotted with lines() with an arrow?
> Or are there any other functions to add an arrow to plot?
>
>
> thanks!
>
--
Sarah Goslee
http://www.functionaldiversity.org
_
hi,
is it possible to end a line plotted with lines() with an arrow?
Or are there any other functions to add an arrow to plot?
thanks!
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Folks:
1. No blame implied to anyone.
However, as 64 bit Windows/R build is apparently not available except as a
commercial product, may I suggest that henceforth it should not be discussed
in this list and that any queries about it simply be directed to David Smith
at Revolution. It just doesn't
Blindingly obvious!! Thanks for that. I can see it will help me find more than
one option that I've been interested in.
--- On Tue, 10/6/09, baptiste auguie wrote:
> From: baptiste auguie
> Subject: Re: [R] ggplot equivalent of par(xaxt)
> To: "John Kane"
> Cc: "R R-help"
> Received: Tues
The bayesm package implements such models.
Hth,
Mike
On Tue, 6 Oct 2009 12:41:18 -0700
James Martin wrote:
> Hello list,
>
> I was trying to investigate the possible use of a mixed effect
> multinomial logit model in R. Does anyone have suggestions on where
> to find information on these mo
Is there any R package that implements a bifurcating autoregression,
aka the BAR(n) model? I've been reading the Huggins and Staudte paper,
"Variance Components Models for Dependent Cell Populations", from the
Journal of the American Statistical Association, 1994.
Shawn Garbett
Vanderbilt
Hello list,
I was trying to investigate the possible use of a mixed effect multinomial
logit model in R. Does anyone have suggestions on where to find information
on these models and the associated functions in R.
Thanks in advance,
jm
--
James A. Martin
850-445-9773
[[alternative HT
A few days ago on the list I had wrestled with the aes() vs aes_string()
issue, along with the same issue with facetting.
The way I ended up handling the point you bring up, Baptiste, is perhaps
rather inefficient but my data sets are not large. I allow the user to pass
variables, then I use that
Emacs + ESS does too.
Giovanni
> Date: Tue, 06 Oct 2009 12:34:28 -0400
> From: Gabor Grothendieck
> Sender: r-help-boun...@r-project.org
> Cc: r-help@r-project.org
> Precedence: list
> DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma;
> h=domainkey-signature:mime
Hi,
I checked the 3rd edition of this book. But I don't find Q-Q plot.
Would you please take a look of the table of content below and let me
know if the same section is available in the 3rd edition?
http://www.amazon.com/Plane-Answers-Complex-Questions-Theory/dp/0387953612/ref=sr_1_1?ie=UTF8&s=bo
body(theme_grey)
could help you find the name of a particular option (that's what I did).
baptiste
2009/10/6 John Kane :
> Lovely. I knew it was not that difficult, I had even gotten as far as
> deciding it had to be an opts() command but I had no idea of what it was.
>
> Thanks very much.
>
>
Lovely. I knew it was not that difficult, I had even gotten as far as deciding
it had to be an opts() command but I had no idea of what it was.
Thanks very much.
--- On Tue, 10/6/09, baptiste auguie wrote:
> From: baptiste auguie
> Subject: Re: [R] ggplot equivalent of par(xaxt)
> To: "Joh
Hi,
2009/10/6 John Kane :
> How do I suppress the numbers on the x-axis?
>
Try this,
p + opts(axis.text.x = theme_blank())
HTH,
baptiste
> Thanks
>
>
> __
> [[elided Yahoo spam]]
>
> ___
I wanted to correct a couple of misconceptions raised in Jose's post
below, which I'll take the liberty of addressing out of order.
On Tue, Oct 6, 2009 at 7:34 AM, Jose Quesada wrote:
> While revolution has
> provided very nice packages to the community (e.g., foreach), the win-64 port
> as
> of
I am playing around learning ggplot and cannot see how to suppress the x or y
axis values ( equivalent of xaxt in basic graphics)
It must be obvious but I'm not seeing it.
Problem
=
timedata <- structure(list(month = structure(c(5
I spotted quite nice blog post by learning r blog
http://learnr.wordpress.com/2009/10/06/export-data-frames-to-multi-worksheet-excel-file/
very good summary how to export data from R to multiple Excel sheets.
- Johannes
__
R-help@r-project.org mailing
vec <- 1:700
names(vec) <- 700:1
names(vec)[grep("101", names(vec))]
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of kayj
> Sent: Tuesday, October 06, 2009 12:49 PM
> To: r-help@r-project.org
> Subject: [R] How to extract names
On 10/6/2009 1:43 PM, David M Smith wrote:
On Tue, Oct 6, 2009 at 10:30 AM, Duncan Murdoch wrote:
When running a 32 bit program, 64 bit Windows hides most of itself outside
the address space visible to the program, so almost all of the 4 Gb address
space is available to the user. But no more:
On Oct 6, 2009, at 1:51 PM, Kjetil Halvorsen wrote:
On Sat, Oct 3, 2009 at 2:45 PM, David Winsemius > wrote:
Do you have a citation for that statement? I cannot convince myself
that it
should be true.
OK. that took some time, since I have no nonparametrics book with me,
but it is a fairly
See ?grep
On Tue, Oct 6, 2009 at 2:48 PM, kayj wrote:
>
> Hi All,
>
> I have a character vector of length=700. The vector contains names and I
> want to extract the names that contain the number 101. The number 101 could
> be anywhere within the name.
>
> what is the best way to do this?
>
>
Hi All,
I have a character vector of length=700. The vector contains names and I
want to extract the names that contain the number 101. The number 101 could
be anywhere within the name.
what is the best way to do this?
--
View this message in context:
http://www.nabble.com/How-to-extract-n
Dear R users,
I am trying to have an xyplot of a data set which has the following
variables:
case (n=10,20,30)
parameter (parm=a,b)
group (grp=g1,g2)
y (y values)
x (x=2,4,8)
My plot should be parameter by case such that I have 2 rows (each row= each
parameter) and 3 columns (each column=each case
Hola! see below.
On Sat, Oct 3, 2009 at 2:45 PM, David Winsemius wrote:
> Do you have a citation for that statement? I cannot convince myself that it
> should be true.
OK. that took some time, since I have no nonparametrics book with me,
but it is a fairly
standard assumption the friedman.test
On Tue, Oct 6, 2009 at 10:30 AM, Duncan Murdoch wrote:
> When running a 32 bit program, 64 bit Windows hides most of itself outside
> the address space visible to the program, so almost all of the 4 Gb address
> space is available to the user. But no more: no matter how much RAM you
> install, it
On 10/6/2009 1:17 PM, Kjetil Halvorsen wrote:
On Tue, Oct 6, 2009 at 11:23 AM, Duncan Murdoch wrote:
On 10/6/2009 10:34 AM, Jose Quesada wrote:
Robert Wilkins gmail.com> writes:
Will R have more glitches on one operating system as opposed to
another, or is it pretty much the same?
robert
Does this do what you want:
)> x <- read.csv(textConnection("home,sqr_footage,cost,color,exterior
+ 1,1500,15,,Siding
+ 2,2000,20,Red,Brick
+ 3,,30,Grey,Brick
+ 4,3500,35,Blue,
+ 5,4000,45,Red,Brick"), na.strings='')
> closeAllConnections()
>
> x
home sqr_footage cost color
Hi Baptiste: Thanks for the suggestion. It will work perfectly.
I would have never considered assigning a color to a variable that contained
no colors at all! I guess this is part of the aesthetic concept, which I
haven't had time to reflect on much. Then later, specify a manual color
scale wh
On Tue, Oct 6, 2009 at 11:23 AM, Duncan Murdoch wrote:
> On 10/6/2009 10:34 AM, Jose Quesada wrote:
>>
>> Robert Wilkins gmail.com> writes:
>>
>>>
>>> Will R have more glitches on one operating system as opposed to
>>> another, or is it pretty much the same?
>>>
>>> robert
>>>
>>>
>>
>> One impor
Well, I guess I posted a poor example.
This example is a little closer:
home,sqr_footage,cost,color,exterior
1,1500,15,,Siding
2,2000,20,Red,Brick
3,,30,Grey,Brick
4,3500,35,Blue,
5,4000,45,Red,Brick
This one actually shows the presence of the null fields that are loaded.
So
In addition to the taskview (below), you might consider using the "sos"
package, something like the following:
sa <- ???'spatial autocorrelation' # 58 matches
sc <- ???'spatial correlation' # 181 matches
s. <- sa|sc # union of the to 'findFn' objects
summary(s.) # 224 links in 57 packages
ins
Further to my previous reply, it occurred to me that ggplot2 would
only ever use data and colors in your calls to compareCats(): res =
res, fac1 = fac1, fac2 = fac2 have no effect whatsoever.
If you want the user to be able to specify the variables used in the
ggplot2 call, you probably want to lo
Stephen Collins-6 wrote:
>
> I am using mle as a wrapper from optim( ). How would I extract the
> convergence code, to know that optim( ) converged properly?
>
>
library(stats4)
example(mle)
slotNames(fit1)
f...@details
f...@details$convergence
--
View this message in context:
http://
Stephen Collins wrote:
I am using mle as a wrapper from optim( ). How would I extract the
convergence code, to know that optim( ) converged properly?
The return value from optim is contained in the details slot, so
> f...@details$convergence
[1] 0
--
O__ Peter Dalgaard
Hi,
On Oct 6, 2009, at 12:30 PM, Gregory Gentlemen wrote:
Hi fellow R-users,
Are there any text editors that recognize sweave (.rnw) files? I am
running Windows Vista and in the past I used Tinn-R for R files but
it (surprisingly) doesn't recognize rnw files and does not do any
syntax hi
It would help to understand the problem you are trying to solve and
the constraints that you have to live under. I routinely process
files with millions of rows of data, do a lot of processing and create
graphics/reports from them in what I think is reasonable time (a
couple of minutes at most for
I saved your data as test.csv, and
> read.csv("~/test.csv", header = TRUE)
home sqr_footage cost
111500 15
222000 20
33 NA 30
443500 35
554000 45
I am using R 2.8.1, old I know... but maybe something else is goi
Hello:
That's because the required input should be a symmetric matrix with all the
comparisons, e.g., it should have the self-comparisons diagonal. In an
hypotetical case with two treatments, the pairwise-comparison test output is
something like that:
A B
B 0.80-
but multcompV
I would like to perform a regression like the one below:
lm(x ~ 0 + a1 + a2 + a3 + b1 + b2 + b3 + c1 + c2 + c3, data=data)
However, the data has the property that a1+a2+a3 = A, b1+b2+b3 = B, and
c1+c2+c3 = C, where A, B, and C are positive constants. So there are two
extra degrees of freedom, a
On 10/6/09, Gregory Gentlemen wrote:
> Are there any text editors that recognize sweave (.rnw) files? I am running
> Windows Vista and in the past I used Tinn-R for R files but it (surprisingly)
> doesn't recognize rnw files and does not do any syntax highlighting for them.
>
Slightly off-topic
Another solution, as a fix to my original algorithm, was found by a colleague
(Matthew Roberts). While he claims not too much for its elegance, it does
seem to work. This fix is based on the use of the 'pmax' function. This
function is a variant of the 'max' (maximum) function to return a vecto
I am using mle as a wrapper from optim( ). How would I extract the
convergence code, to know that optim( ) converged properly?
Thanks,
Stephen Collins, MPP | Analyst
Global Strategy | Aon Benfield
[[alternative HTML version deleted]]
__
R
I believe I may be overlooking something simple in order address this, but I
have searched RSeek.org and using "?", but cannot seem to find anything
discussing this one.
I am using read.csv to read in a csv file. Evidently in places there is
nothing between the commas, so that when the data i
Hi,
I may be missing an important design decision, but could you not have
only a single data.frame as an argument of your function? From your
example, it seems that the colour can be mapped to the fac1 variable
of "data",
compareCats <- function(data) {
require(ggplot2)
p <- ggplot(data, a
vim/gvim does syntax highlighting of R, Sweave and latex.
On Tue, Oct 6, 2009 at 12:30 PM, Gregory Gentlemen
wrote:
> Hi fellow R-users,
>
> Are there any text editors that recognize sweave (.rnw) files? I am running
> Windows Vista and in the past I used Tinn-R for R files but it (surprisingly)
Emacs with ESS would be one way to go. When editing a RNW file in ESS, you can
use a key combination to go from RNW directly to PS or PDF, which is nice.
Erik
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Gregory Gentlem
Hi fellow R-users,
Are there any text editors that recognize sweave (.rnw) files? I am running
Windows Vista and in the past I used Tinn-R for R files but it (surprisingly)
doesn't recognize rnw files and does not do any syntax highlighting for them.
Thanks in advance,
Greg
___
esp wrote:
>
> For the function as defined above using 'sapply'
>> spot[,1]
> 01/09/2009 01/09/2009 00:00:01 01/09/2009 00:00:02 01/09/2009
> 00:00:03
> 1251759600 1251759601 1251759602
> 1251759603
>
> This was unexpected - it seems to have display
On Oct 6, 2009, at 12:14 PM, Corrado wrote:
Dear list,
is there an iterative re weighted least square based algorithm or
any or other
iteratively re weighted optimisation algorithms for non linear (and
possibly
non parametric) optimisation problems with constraints available in R?
Have
Hello Again... I¹m making a faceted plot of a response on two categorical
variables using ggplot2 and having troubles with the coloring. Here is a
sample that produces the desired plot:
compareCats <- function(data, res, fac1, fac2, colors) {
require(ggplot2)
p <- ggplot(data, aes(fac1,
Dear list,
is there an iterative re weighted least square based algorithm or any or other
iteratively re weighted optimisation algorithms for non linear (and possibly
non parametric) optimisation problems with constraints available in R?
Regards
--
Corrado
I have a frame and I want to fix the progress bar onto it.I do not want
progress bar to display as another GUI.how do i embed progress bar onto a
frame ?
--
Dipti
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https:
Daniel Perkins wrote:
>
>
> Ideally we would do an ANCOVA to test for differences in slope or
> intercepts for the different streams. However as there were repeated
> measures and unequal n and unbalanced design, I have used a linear mixed
> effect model (from nlme package in R) in the form
2009/10/6 Uwe Ligges :
>
>
> Gabor Grothendieck wrote:
>>
>> 2009/10/6 Uwe Ligges :
>>>
>>> The first rule is easy: As long as you are using scalar valued (i.e.
>>> length
>>> 1 vectors in R) "cond", you should prefer
>>> if(cond) cons.expr else alt.expr
>>> rather than
>>> ifelse(cond, yes, no
On 10/6/2009 10:34 AM, Jose Quesada wrote:
Robert Wilkins gmail.com> writes:
Will R have more glitches on one operating system as opposed to
another, or is it pretty much the same?
robert
One important difference is that, if you are unsing large datasets and need
memory, then windows is
Rnewb wrote:
>
> I would like to perform a regression like the one below:
>
> lm(x ~ 0 + a1 + a2 + a3 + b1 + b2 + b3 + c1 + c2 + c3, data=data)
>
> However, the data has the property that a1+a2+a3 = A, b1+b2+b3 = B, and
> c1+c2+c3 = C, where A, B,
>
Ravi Varadhan has an example how this co
Oh yeah, I forgot that, just correct this by:
Data<-as.data.frame(matrix(c("Monday", "Tuesday", "Wed", "Thu",
"Fri","Sat","Sun",56,57,60,75,62,67,70,76,77,81,95,82,77,83),ncol=3,nrow=7,l
ist(paste('Row',1:7,sep=''),c('Days','Hum','Temp')),byrow=F))
Data$Hum<-as.numeric(Data$Hum)
Data$Temp<-as.nume
The date time stamp is not the same?
data.frame(as.chron(Cumul[,"date_time"]), DF[,"date_time"])
On Tue, Oct 6, 2009 at 9:31 AM, hadley wickham wrote:
>> It is much easier to do you data preparation before plotting.
>>
>> Cummul <- ddply(subset(DF, precipitation!="NA"), "gauge_name",
>> function
Ok, thanks for all your responses,
i was thinking, after reading the firsts messages that use only the $p.value
was good. By using a symetric matrix results are ok.
It's also good with a pvalue vector (just like exemples in the help ...).
many thanks for your help
--
View this message in c
R_help Help wrote:
>
> Hi - I read through dse package manual a bit. I'm not quite certain
> how I can use it to estimate a time varying coefficient regression
> model? I might pick up an inappropriate package. Any suggestion would
> be greatly appreciated. Thank you.
>
>
> Just rewrite the l
Hi
beside other options if you do not want to learn how to order factor
levels
Days <- c("Monday", "Tuesday", "Wed", "Thu", "Fri", "Sat","Sun")
Hum <-c(56,57,60,75,62,67,70)
Temp<-c(76,77,81,95,82,77,83)
> limits<-range(c(Hum, Temp))
> limits
[1] 56 95
> plot(1:7, Hum, ylim=limits)
> points(1:7
atorso wrote:
>
> Hello,
>
> I'm having an error when trying to fit the next GLM:
>
>>>model<-glm(response ~ CLONE_M + CLONE_F + HATCHING
> +(CLONE_M*CLONE_F) + (CLONE_M*HATCHING) + (CLONE_F*HATCHING) +
> (CLONE_M*CLONE_F*HATCHING), family=quasipoisson)
>>> anova(model, test="Chi")
>
>
>
>
This may work for you.
Then you can custom your graphics with ?par.
Data<-as.data.frame(matrix(c("Monday", "Tuesday", "Wed", "Thu", "Fri",
"Sat","Sun",56,57,60,75,62,67,70,76,77,81,95,82,77,83),ncol=3,nrow=7,list(pa
ste('Row',1:7,sep=''),c('Days','Hum','Temp')),byrow=F))
plot(Data$Days,Data$Temp)
Hi Peng Yu,
Chapter 13 of the following book provides a good description of the
assumption done when using regression and other techniques. It also
discusses the QQplot.
@BOOK{Christensen1996,
title = {Plane Answers to Complex Questions: The Theory of Linear Models},
publisher = {Springer,
David Winsemius wrote:
On Oct 6, 2009, at 4:46 AM, Geoffrey William Heard wrote:
Hi All
I'm creating some 3-D plots using the function 'persp', and have a
query regarding the ability to make changes to the label of the z-axis.
There are two things I would like to do. First, the default s
Hi,
I want to look for some detailed explanation on the properties of Q-Q
plot and how the properties are derived.
In R, there is the following reference.
Becker, R. A., Chambers, J. M. and Wilks, A. R. (1988) The New S
Language. Wadsworth & Brooks/Cole.
Somebody also mentioned the following boo
Actually, it's a bit more complex than just plot(). First you have to
make Days an ordered factor, otherwise you get indeed that error.
days <- c("Mon", "Tue", "Wed", "Thu", "Fri", "Sat","Sun")
Days <- factor(days,levels=days,ordered=T)
Then you want to make a plot with 2 axes, you should check :
Thanks Sara,
Yes I did try. I could not get the Days on the X-axis
blow is theerror message
plot(Temp,Days)
Error in plot.window(...) : need finite 'ylim' values
In addition: Warning messages:
1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
2: In min(x) : no non-missing a
Please provide a reproducible example. I've had problems with
multcompLetters in the past, because I was giving it slightly
incorrect input.
Hadley
On Tue, Oct 6, 2009 at 7:41 AM, goz wrote:
>
> hello,
>
> i try to use the multcomp letters, but i have problems with my results :
>
> here is my p
Robert Wilkins gmail.com> writes:
>
> Will R have more glitches on one operating system as opposed to
> another, or is it pretty much the same?
>
> robert
>
>
One important difference is that, if you are unsing large datasets and need
memory, then windows is by far the worst.
CRAN R is 32 b
> It is much easier to do you data preparation before plotting.
>
> Cummul <- ddply(subset(DF, precipitation!="NA"), "gauge_name",
> function(x){
> x$Cummul <- cumsum(x$precipitation)
> x
> })
With a little less typing:
Cummul <- ddply(subset(DF, precipitation!="NA"), "gauge_name",
If I understand correctly, you can use which( ,arr.ind=TRUE):
Assuming that you start with a data frame,
dat <- as.data.frame(matrix(runif(30), ncol=5))
nm <- names(dat)
cormat <- cor(dat)
cormat[lower.tri(cormat, diag=TRUE)] <- NA
idx <- which(cormat > 0.4, arr.ind=TRUE)
idx
cbind(nm[idx
Dear Stephen,
It is much easier to do you data preparation before plotting.
Cummul <- ddply(subset(DF, precipitation!="NA"), "gauge_name",
function(x){
x$Cummul <- cumsum(x$precipitation)
x
})
ggplot(Cummul, aes(x = date_time, y = Cummul)) + geom_line() +
facet_wrap(~gauge_name, s
Did you try it? With, perhaps, plot() ? And lines() ?
You might do better with Days as a factor with the day names in order. (And
why are two full and five abbreviated?)
I don't understand why Hum and Temp are matrices rather than vectors,
and why then you didn't specify dimensions, and for that
Gabor Grothendieck wrote:
2009/10/6 Uwe Ligges :
The first rule is easy: As long as you are using scalar valued (i.e. length
1 vectors in R) "cond", you should prefer
if(cond) cons.expr else alt.expr
rather than
ifelse(cond, yes, no)
because the latter one evaluates both "yes" and "no" whi
-- Forwarded message --
From: Antonio Paredes
Date: Tue, Oct 6, 2009 at 9:41 AM
Subject: Re: [R] Long for Loop- calling C from R - Parallel Computing
To: Karl Ove Hufthammer
Hello again,
I'm hoping to get a response from some of the R gurus in this list. Is my
assumption that R
OK, so maybe last night was a little too much at one throw, so I have
reduced the data to two stations- one that has precipitation and one
that does not. This is going to be in the context of a larger data
set. I would like to be able to issue a ggplot command and have cum
sum just act on the fac
1 - 100 of 138 matches
Mail list logo
