Dear all,
Your advices was a great help to my study.Thank you very much!
--
View this message in context:
http://r.789695.n4.nabble.com/Is-the-output-of-survfit-coxph-survival-or-baseline-survival-tp3861919p3873512.html
Sent from the R help mailing list archive at Nabble.com.
__
my=function(x){
len=1
for(i in 1:len){
y[i]=x[i]
}
g=1
w=NULL
t=NULL
for(i in 1:len)w[i]=x[i+len]
for(i in 1:len)t[i]=x[i+2*len]
for(i in 1:len)g=g*dnorm(y[i])*dnorm(w[i])*dnorm(z[i])
return(g)
}
cuhre(6,1,my,rep(-100,6),rep(100,6))
Error in crff(match.call(), integrand, "cuhre", libargs, ...) :
Hi. I am trying to install using the following. Can someone suggest what is
wrong? I am using Windows 7 64bit, and R 2.10.1
('C:\Users\Bill\Desktop\DMwR_0.2.1.zip', repos=NULL )
Warning in install.packages("C:UsersBillDesktopDMwR_0.2.1.zip", repos = NULL) :
argument 'lib' is missing: using 'C:\
Dear all,
I think the coxph and survfit.coxph can give the cumulative hazard of cox
model.
But is there any method to calculate the hazard
Lambda(t)=lambda_0(t)*exp{beta*X(t)}?
Any suggestion will be great help.
Thank you very much!
Koshihaku
--
View this message in context:
http://r.789695.n4.n
Hi bstudent,
I've had the same problem and I wish there was a definitive answer as this
seems to be the #1 problem with the package and pgmm would be awesome for
economists if we could figure out how to work it! I'm no expert on GMM, but
from what I've gathered from other posts, the problem may st
?expand.grid
Am 05.10.2011 00:21, schrieb darkgaze:
I don't quite know how to word what I want, but if I have
(1, 2, 3); (a, b, c); (x, y)
I want:
1 a x
1 b x
1 c x
1 a y
1 b y
1 c y
2 a
...
and so forth
What is the appropriate command?
Best,
Don
--
View this message in context:
http://r
>
> On Tue, 4 Oct 2011, Sarah Goslee wrote:
>
> > You asked for pointers, and didn't provide a reproducible example, so
I
> > offered a pointer.
>
> Sarah,
>
>I did not realize that your pointer was to the factor component of
the
> subset() command.
>
>I think the most parsimonious t
Hi Scott,
I am not familiar with S-Plus (though many aspects are quite similar
to R). I will say that your function looks approximately correct. I
am not familiar with the ss.rand function. I searched, and found some
things that I suspect are similar in the packages MBESS, but without
knowing m
On Tue, 4 Oct 2011, bonda wrote:
Hello all,
I have the following two questions: 1) how can I get the values of
boundaries for fluctuation process gefp(), for functionals maxBB, meanL2BB,
etc?
Both, maxBB and meanL2BB, are objects of class "efpFunctional" that
contain several functions includ
Tena koe Don
?expand.grid
HTH
Peter Alspach
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of darkgaze
> Sent: Wednesday, 5 October 2011 11:21 a.m.
> To: r-help@r-project.org
> Subject: [R] Create combinations of rows
>
man awk?
I've used awk for similar tasks (if I am reading the post correctly.) Google-Fu
should turn up some useful examples.
Also awk should be on your linux installation in some form or another.
Regards,
Ken Hutchison
On Oct 4, 2554 BE, at 10:52 PM, "Tom D. Harray" wrote:
> Hi there,
>
Dear folks,
I’m trying to build a function to create and make available some variables I
frequently use for testing purposes. Suppose I have a function that takes
some inputs and creates (internally) several named objects. Say,
fun1 <- function(x, y, z) {obj1 <- x; obj2 <- y; obj3 <- z
}
Her
Hi there,
I would like to use my linux system's fgrep to search for a text pattern
in a file. Calling system with
system("fgrep \"SearchPattern\" /path/to/the/textFile.txt")
works in general, but I need to search for the search pattern at the
beginning of the line.
The corresponding she
I'm trying to convert an S-Plus program to R. Since I'm a SAS programmer I'm
not facile is either S-Plus or R, so I need some help. All I did was convert
the underscores in S-Plus to the assignment operator <-. Here are the first
few lines of the S-Plus file:
sshc _ function(rc, nc, d, meth
I don't quite know how to word what I want, but if I have
(1, 2, 3); (a, b, c); (x, y)
I want:
1 a x
1 b x
1 c x
1 a y
1 b y
1 c y
2 a
...
and so forth
What is the appropriate command?
Best,
Don
--
View this message in context:
http://r.789695.n4.nabble.com/Create-combinations-of-rows-tp387
Hi. I am not that good at R but I was wondering if there is either a tool or a
strategy for testing many different models in R in a batch. I have used
something in Weka called the Experimenter interface which helps with doing this
kind of thing.
Thank you.
Andy
[[alternative HTML versio
Hello all,
I have the following two questions: 1) how can I get the values of
boundaries for fluctuation process gefp(), for functionals maxBB, meanL2BB,
etc? 2) how can I get fragments of gefp()-process, e.g., if I have n=200
observations, i=1,2,...,200, and need gefp()[50:100], i.e. from i=50 ti
Hi,
I tried to google for any solution for asof join operator in R. But I
couldn't find one. The asof join operator AsOf(A,B) merges 2 time
series by looking for latest available value of B prior to each time
point in A. For example,
A <- xts(c(10,15,20,25),
order.by=as.POSIXct(c("2011-09-01","20
Hi Syrvn,
On Tue, Oct 4, 2011 at 8:43 AM, syrvn wrote:
> Hello,
>
> I am working on a big R project using Eclipse/StatET/Texlipse. I'd like to
> write a Latex document within that project but DO NOT want to Sweave it.
> It's pure Latex. Via the external tools configurations I set up 2 different
>
Hi,
As other have said, it's very difficult to help you without an example
+ code to know what you are talking about.
That having been said, it seems as if you are just getting your feet
wet in this R <--> C bridge, and I'd recommend you checkout the "Rcpp"
and "inline" package to help make your
Thanks co much Martin for your bug reply.
I will follow up with it.
Best,
Angel
On Tue, Oct 4, 2011 at 7:35 PM, Martin Morgan wrote:
> On 10/04/2011 02:44 PM, Angel Russo wrote:
>
>> Hi All,
>>
>> I have an CBS segmentation algorithm output for 10 tumor samples each from
>> 2
>> different tumo
On 05/10/11 01:04, Grigory Alexandrovich wrote:
Hello,
I wrote a function in C, which works fine if called from the
main-function in C.
But as soon as I try to call this function from R like .C('foo',
as.double(x), as.integer(y)), the programm crashes.
I created a dll with the cmd command
Dear R-Users
I have come across the error that apparently has nothing to do with command
itself. Here is the error
(w - matrix (or vector) e testXaxis - dates).
> plot(data.frame(testXaxis,w),col="blue",ylab="Q, [m3/s]",xlab="Data",
+ main="Free-run - Modelo NARX MISO - Test Data")
Error
Hi:
Try this:
p <- ggplot(mtcars, aes(factor(cyl), mpg))
p + geom_boxplot(aes(colour = factor(am)), fill = 'white') +
scale_colour_manual('am', values = c('0' = 'blue', '1' = 'black'))
HTH,
Dennis
On Tue, Oct 4, 2011 at 1:56 PM, Brian Smith wrote:
> Hi,
>
> I wanted to change the default
On 10/04/2011 02:44 PM, Angel Russo wrote:
Hi All,
I have an CBS segmentation algorithm output for 10 tumor samples each from 2
different tumors.
Now, I am in an urgent need to assign gene (followed by all genes present)
that belong to a particular segment after I removed all the CNVs from
segm
Hi:
> INB4: if I have a nested design with treatment A and treatment B
> within A, F-values are MSA/MSA(B) and MSA(B)/MSE, correct? How can I
> make R give these values directly, without further coding?
This is how to get an equivalent model in lme4, but it probably isn't
what you expect (particu
Hi all,
I am trying to print some generated text on plot which is comprised of text and
the value of a variable. Thanks to the response from Joshua Wiley
on the thread "How to format R superscript 2 followed
by "=" value, I can now get R^2 = value. However with the following code:
ggplot(blDa
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of jricci
> Sent: Tuesday, October 04, 2011 1:28 PM
> To: r-help@r-project.org
> Subject: Re: [R] how do i put two scatterplots on same graph
>
> I am new at this. The two data sets
Hello R users,
I have a plot type=b with x axis at=(1:36),
I would like to increase the distance between x tick-marks 8 and 9, and not
connect the points x=8 and x=9.
I can do the second thing, setting type="p" and then drawing the lines, but
I don't know how to do the first.
Plus, I was wond
Hello
I have the following problem when loading the package hydroGOF on Windows 7
running R.12.2
library(hydroGOF)
Error : package 'hydroTSM' does not have a name space
Error: package/namespace load failed for 'hydroGOF'
The same command does not result in error on R.13.1 at my Mac r
Hi,
I'm looking at a very simple picture created with rasterImage.
foo <- matrix(1:9,3,3)
foo[,] <- rainbow(9)
plot(0:1,0:1,t='n')
rasterImage(foo,0,0,1,1)
If I choose to specify interpolate=F, I get the expected 9 blocks of
color. My question is: how many values (aka pixels) does rasterImage
Hi All,
I have an CBS segmentation algorithm output for 10 tumor samples each from 2
different tumors.
Now, I am in an urgent need to assign gene (followed by all genes present)
that belong to a particular segment after I removed all the CNVs from
segment data. The format of the data is:
Sample
Hadley:
Below is an example of what I am trying to do, I just don't understand
how to supply the limits to the blue and pink shaded regions and the
values of the black line, which are meant to represent from bottom to
top, the 5%, 25%, 50%, 75%, 95% limits that I get from quantile():
h + ge
Hi Hadley:
When I tried your function on the example data, I got the following:
dd <- data.frame(year = rep(2000:2008, each = 500), y = rnorm(4500))
g <- function(df, qs = c(.05, .25, .50, .75, .95)) {
data.frame(q = qs, quantile(d$y, qs))
}
ddply(dd, .(year), g)
> ddply(dd, .(year), g)
year
I am new at this. The two data sets don't have color variable just paired data.
How should I structure the data sets in R?
Joe Ricci
From: William Revelle [via R]
To: Joe Ricci
Sent: Tue Oct 04 11:18:20 2011
Subject: Re: how do i put two scatterplots on same grap
I'm going to put on my fire suit and wade in (see inline)
On 10/4/2011 8:11 AM, Bert Gunter wrote:
On Tue, Oct 4, 2011 at 7:42 AM, Jeanne M. Spicerwrote:
I'm not sure how returning an incorrect result is ever a 'positive' feature
It is **not** "incorrect"; perhaps unexpected, but that is not
Dear Michael,
Thanks very much.
Jiang
On Tue, Oct 4, 2011 at 2:39 PM, R. Michael Weylandt <
michael.weyla...@gmail.com> wrote:
> On Mon, Oct 3, 2011 at 10:08 PM, chunjiang he wrote:
> > Hi,
> >
> > I have two questions want to ask.
> >
> > 1. If I have a matrix like this, and I want to figure o
Hi,
I wanted to change the default colors appearing in boxplot. For example, the
following code (from the package/documentation):
===
library(ggplot2)
p <- ggplot(mtcars, aes(factor(cyl), mpg))
p + geom_boxplot(aes(fill = factor(am)))
===
Gives the default colors. What do I nee
On Tue, Oct 4, 2011 at 4:41 PM, Gavin Simpson wrote:
> On Tue, 2011-10-04 at 08:45 +, Ashley Houlden wrote:
>> Hi,
>
>> #distance matrix of samples using bray curtis
>> d= bcdist(mydata, rmzero=FALSE)
>
> In addition, you don't necessarily need ecodist for the bray curtis
> distance. vegdist(
On 5/10/2011 7:25 a.m., Charles McClure wrote:
I am new to R and have recently tried Tinn-R with very mixed and unexpected
results. Can you point me to a Tinn-R tutorial on the web or a decent
reference book?
Thank you for your help;
Charles McClure
cmccl...@atrcorp.com
cfmccl...@verizon.net
On Tue, 2011-10-04 at 08:45 +, Ashley Houlden wrote:
> Hi,
> #distance matrix of samples using bray curtis
> d= bcdist(mydata, rmzero=FALSE)
In addition, you don't necessarily need ecodist for the bray curtis
distance. vegdist() in vegan will compute this for you.
Not that there is anything
On Tue, 2011-10-04 at 08:45 +, Ashley Houlden wrote:
> Hi,
>
> forgive me if someone has already posted about this but I have had a
> look and cannot find the answer, also I am very new to R and been
> getting the grips with this.
>
> I have been trying to use Adonis to find out if there are
Perhaps rbind?
Michael
On Tue, Oct 4, 2011 at 3:48 PM, Jose Bustos Melo wrote:
> Hello everyone,
>
> I know this is very basic question for you people. I'm working with mani
> diferent tables, but everyone has the same variables. (V1, V2, V3). The only
> think that I need to do is to put toget
You are, of course, entirely correct and, once again, I tip my hat to
the erudition of those who comment on this list. My initial
formulation, for a distance on a normed space inherited from the norm,
stands trivially, but as you rightly point out, I'm excluding many
interesting and possibly useful
Hello everyone,
I know this is very basic question for you people. I'm working with mani
diferent tables, but everyone has the same variables. (V1, V2, V3). The only
think that I need to do is to put together this tables. In other words,
creating just one big table with all the cases showed in
On Tue, 4 Oct 2011, R. Michael Weylandt wrote:
No, that was just a typo on my end:
the correct order of arguments should have been
ff <- grepl("BC-", ff)
Michael,
Thank you.
Rich
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/
On Tue, 4 Oct 2011, Sarah Goslee wrote:
You asked for pointers, and didn't provide a reproducible example, so I
offered a pointer.
Sarah,
I did not realize that your pointer was to the factor component of the
subset() command.
I think the most parsimonious thing for me to do is to modify
On Mon, Oct 3, 2011 at 10:08 PM, chunjiang he wrote:
> Hi,
>
> I have two questions want to ask.
>
> 1. If I have a matrix like this, and I want to figure out the rows whose
> value in the 3rd column are less than 0.05. How can I do it with R.
> hsa-let-7a--MBTD1 0.528239197 2.41E-05
> hsa-l
Hello all
I'm trying to learn how to fit a nested model in R. I found a toy
example on internet where a dataset that have 3 areas and 4 sites
within these areas. When I use Minitab to fit a nested model to this
data, this is the ANOVA table that I got:
Nested ANOVA: y versus areas, sites
Analysi
The following for loops does the work but it takes a good 30 minutes to run:
for(i in 1:length(myStopwords))
{
currentWord <- myStopwords[i]
tr1=tm_map(tr1,removeWords,currentWord)
}
Are there any faster alternatives?? Thank you.
Ravi
--
View this message in context:
http://r.789695.n4.n
I am new to R and have recently tried Tinn-R with very mixed and unexpected
results. Can you point me to a Tinn-R tutorial on the web or a decent
reference book?
Thank you for your help;
Charles McClure
cmccl...@atrcorp.com
cfmccl...@verizon.net
__
R-
> # Function to compute quantiles and return a data frame
> g <- function(d) {
> qq <- as.data.frame(as.list(quantile(d$y, c(.05, .25, .50, .75, .95
> names(qq) <- paste('Q', c(5, 25, 50, 75, 95), sep = '')
> qq }
You could cut out the melt step by making this return a data frame:
g <
No, that was just a typo on my end:
the correct order of arguments should have been
ff <- grepl("BC-", ff)
On Tue, Oct 4, 2011 at 3:07 PM, Rich Shepard wrote:
> On Tue, 4 Oct 2011, R. Michael Weylandt wrote:
>
>> This isn't going to be the most elegant, but it should work:
>> ## Get the factors
... and, as an aside, if you had simply searched within R for (the
obvious?!)
??wildcard
you would have received the suggestion for glob2rx() in utils, which
actually would have enabled you to use a familiar wildcard expression.
However, the answers you've already received are simpler and more
st
Hi,
Using ddply,
ddply(df, .(ID), mutate, nrows=length(rel.head), test = nrows==2 &
all(rel.head %in% c(1,3)))
HTH,
baptiste
On 5 October 2011 06:02, Dennis Murphy wrote:
> Hi:
>
> Here's another way to do it with the plyr package, also not terribly
> elegant. It assumes that rel.head is a f
?grep
?names
Use indexing by name [, namevector]
---
Jeff Newmiller The . . Go Live...
DCN: Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries O.O#. #.O#. with
/Software/Em
On Tue, 4 Oct 2011, R. Michael Weylandt wrote:
This isn't going to be the most elegant, but it should work:
## Get the factors as characters
ff <- as.character(chemdata$site)
## Identify those that match what you want
ff <- grepl(ff, "BC-")
Michael,
Apparently grep works differently in R
Hi Rich,
On Tue, Oct 4, 2011 at 2:58 PM, Rich Shepard wrote:
> On Tue, 4 Oct 2011, Sarah Goslee wrote:
>
>> You can use something like this:
>>
>>> testdata <- c("A1", "A2", "A3", "B1", "B2", "B3")
>>> grep("^A", testdata)
>>
>> [1] 1 2 3
>>>
>>> grepl("^A", testdata)
>>
>> [1] TRUE TRUE TRUE
On Tue, 4 Oct 2011, Sarah Goslee wrote:
You can use something like this:
testdata <- c("A1", "A2", "A3", "B1", "B2", "B3")
grep("^A", testdata)
[1] 1 2 3
grepl("^A", testdata)
[1] TRUE TRUE TRUE FALSE FALSE FALSE
Sarah,
I don't see how this gives me a data frame containing only thos
This isn't going to be the most elegant, but it should work:
## Get the factors as characters
ff <- as.character(chemdata$site)
## Identify those that match what you want
ff <- grepl(ff, "BC-")
now use this logical vector to subset
chemdata[ff, ]
Can't test, but should be good to go assuming
Hi Rich,
You can use something like this:
> testdata <- c("A1", "A2", "A3", "B1", "B2", "B3")
> grep("^A", testdata)
[1] 1 2 3
> grepl("^A", testdata)
[1] TRUE TRUE TRUE FALSE FALSE FALSE
Sarah
On Tue, Oct 4, 2011 at 2:39 PM, Rich Shepard wrote:
> I have a data frame called chemdata with t
I have a data frame called chemdata with this structure:
str(chemdata)
'data.frame': 14886 obs. of 4 variables:
$ site: Factor w/ 148 levels "BC-0.5","BC-1",..: 104 145 126 115 114 128
124 2 3 3 ...
$ sampdate: Date, format: "1996-12-27" "1996-08-22" ...
$ param : Factor w/ 8 lev
Hi:
The smooth is not going to replicate the quantile estimates you get
from the 'boxplots'; the smooth is estimating a conditional mean using
loess, with confidence limits associated with uncertainty in the
estimate of the conditional mean function, which are almost certainly
going to be narrower
I am using the tm package to do text miniing:
I have a huge list of stopwords (2000+) that are in a csv file. I read it as
follows:
stopwordlist <- read.csv("stopwords to be Removed 10042011.csv")
myStopwords <- as.character(stopwordlist$stopwords)
When try removing the stopwords using
tr1=tm_
Hi:
Here's another way to do it with the plyr package, also not terribly
elegant. It assumes that rel.head is a factor in your original data
frame:
> str(df)
'data.frame': 11 obs. of 2 variables:
$ ID : Factor w/ 6 levels "17100","17101",..: 1 1 2 3 4 4 5 5 5 6 ...
$ rel.head: Factor w/
Hadley,
Thanks for responding. No, not smoothed quantile regression. If you go here:
http://www.erh.noaa.gov/mmefs/index.php and click on one of the colored
squares, you can see we have 'boxplots'. What I want to express is the
uncertainty as depicted in the example from my previous email where
I'm trying to use the randomTF function from package palaeoSig to
test the significance of a MAT reconstruction with nine analogs and a
WA-PLS reconstruction with four components. I'm probably missing
something obvious here but how do I make sure that randomTF is testing
the reconstruction
On 04.10.2011 11:10, D.Emad wrote:
Hello,
I've been facing a really stupid problem... When I try to plot using
heatplot or hclust or any similar function, the labels of the x-axis - which
are the samples names - are giant& overlapping. I can't even read the
samples names!
R> heatplot
Error:
This looks like a classic case of not reading the manual, and then compounding
it by not reading the posting guide. The manual would be the "Writing R
Extensions" pdf that comes with R or you can google it. The posting guide is
referenced at the bottom of this and every other posting on this mai
Hi,
I am sure there are better / more efficient ways of doing this, but the
following seems to work ...
ids <- sapply(split(df,df$ID),function(x) {length(x$rel.head)==2 &
any(x$rel.head==1) & any(x$rel.head==3)})
ids <- as.numeric(names(ids)[ids])
added.dummy <- as.numeric(df$ID%in%ids)
cbind(df
Without knowing that C code, we cannot know. Have you read Writing R
Extensions carefully? I.e. take care with memory allocation and printing
as mentioned in the manual.
Uwe Ligges
On 04.10.2011 14:04, Grigory Alexandrovich wrote:
Hello,
I wrote a function in C, which works fine if called f
On 04.10.2011 18:18, Duncan Murdoch wrote:
On 04/10/2011 7:04 AM, S Ellison wrote:
See para 10.3.2 'Identifiers' in the R language definition (always
distributed with R in the html help system), or ?make.names, for a
concise statement of what constitutes a valid variable name in R.
It's actua
On 04.10.2011 16:42, Jeanne M. Spicer wrote:
I'm not sure how returning an incorrect result is ever a 'positive' feature but
at least the documentation could more clearly warn users that this method
behaves differently in these cases -- summary(rock[,1]) vs summary(rock[,1:2])
-- and that th
On 04/10/2011 7:04 AM, S Ellison wrote:
See para 10.3.2 'Identifiers' in the R language definition (always distributed
with R in the html help system), or ?make.names, for a concise statement of
what constitutes a valid variable name in R.
It's actually underscores that might give trouble with
On 04.10.2011 13:30, Peter Minting wrote:
Dear R-help
Can anyone tell me why my curve appears the wrong way round on a rug plot?
I am using the same code as on pg 596 of the Crawley R-book.
mod<-glm(mort~logBd,binomial)
What is mort, what is logBd? I don't have access to the book. I have
Hi,
On Tue, Oct 4, 2011 at 11:39 AM, jdospina wrote:
> Hello.
>
> Not at all in the way you have shown. Just to improve your code
> "readability", try to avoid naming your variables beginning with period
> (example: .hello).
Well, that's not exactly true.
It's "common practice" to name variable
Hi:
Here's a reproducible example:
d <- data.frame(grp = factor(rep(c('x', 'y'), each = 5)),
ev = rnorm(10), dv = rnorm(10))
labl <- list(expression(italic('x')), expression(italic('y')))
ggplot(d, aes(x = ev, y = dv, shape = grp)) + geom_point() +
scale_shape_manual('Group',
Bert Gunter gene.com> writes:
>
> Below.
>
> On Tue, Oct 4, 2011 at 7:34 AM, Panagiotis hi.is> wrote:
>
> > Hi,
> >
> > I applied a linear mixed effect model in my data using the nlme package.
> > lme2<-lme(distance~temperature*condition, random=~+1|trial, data) and then
> > anova.
> > I want
Hmm, that's my fault when composing this mail, but the problem was
really encountered at that time.
Nevertheless, neither can I reproduce the problem now, perhaps I just
made another mistake at that time.
Thanks all the same, and sorry for the disturbance anyway :|
On Tue, Oct 04, 2011 at 10:10:56
I'm not sure how returning an incorrect result is ever a 'positive' feature but
at least the documentation could more clearly warn users that this method
behaves differently in these cases -- summary(rock[,1]) vs summary(rock[,1:2])
-- and that the method can and does return incorrect results wi
Hello.
Not at all in the way you have shown. Just to improve your code
"readability", try to avoid naming your variables beginning with period
(example: .hello).
In contrast with Matlab (for example) the period in R is not to have access
to an object property.
--
View this message in context:
h
On 04.10.2011 08:53, Divyam wrote:
Hi users!
I am fitting a model with several factor variables as independents using
svm. since there are lots of categorical variables,the training and test
data sets have been created using dummy.data.frame option from dummies
package. I have a factor A in th
The Financial Mathematics department of Fraunhofer ITWM
is offering a two-days workshop on Finance with R:
%-
[Workshop] Finance with R
%-
Oct 20, 2011, 10:00-17:00 and
Oct 21, 2011, 9:00-16:
On 03.10.2011 18:16, Sandeep Patil wrote:
Hello everyone
I have manually compiled directory of gstat in a particular folder of my
Unix system.
I want to install this and am unable to use either of the following two
commands
1. R CMD INSTALL
2. Install.packages
If this is a precompiled (i
Hi all,
I have a dataset of individuals where the variable ID corresponds to the
identification of the household where the individual lives. rel.head stands
for the relationship with the household head. so rel.head=1 is the household
head, rel.head=2 is the spouse, rel.head=3 is the children.
Her
On 03.10.2011 19:19, Cable, Sam B Civ USAF AFMC AFRL/RVBXI wrote:
I am using readLines to read a fairly large ASCII file. readLines reads
a fixed number of lines, then other R code processes the data, then
readLines reads the same number of lines again, then other R code
processes the data, th
On Tue, Oct 4, 2011 at 7:42 AM, Jeanne M. Spicer wrote:
> I'm not sure how returning an incorrect result is ever a 'positive' feature
>
>
It is **not** "incorrect"; perhaps unexpected, but that is not the same.
> but at least the documentation could more clearly warn users that this
> method be
If the data are from one data.frame (e.g., the iris data set), then simply
label the red and white flowers with different colors:
e.g.,
with the iris data set
plot(iris$Sepal.Length,iris$Sepal.Width,col=c("red","blue","black")[iris$Species],pch=c(16:18)[iris$Species])
Bill
On Oct 4, 2011,
Below.
On Tue, Oct 4, 2011 at 7:34 AM, Panagiotis wrote:
> Hi,
>
> I applied a linear mixed effect model in my data using the nlme package.
> lme2<-lme(distance~temperature*condition, random=~+1|trial, data) and then
> anova.
> I want to ask if it is posible to get the least squares means for th
You need to set the labels...
Hadley
On Sat, Sep 24, 2011 at 3:49 AM, Casper Ti. Vector
wrote:
> Is there any way to use expression() in legend labels with ggplot2?
>
> It seems that things like
>> scale_shape_manual(value = c(
>> x = expression(italic(x)),
>> y = expression(italic(y))
>> ))
New package `partykit': A Toolkit for Recursive Partytioning
The purpose of the package is to provide a toolkit with infrastructure for
representing, summarizing, and visualizing tree-structured regression and
classification models. Thus, the focus is not on _inferring_ such a
tree structure fro
Hi,
I applied a linear mixed effect model in my data using the nlme package.
lme2<-lme(distance~temperature*condition, random=~+1|trial, data) and then
anova.
I want to ask if it is posible to get the least squares means for the
interaction effect and the corresponding 95%ci. And then plot this v
?? Use the appropriate panel function, not panel.xyplot().
If you don't know what this means, you need to read up on lattice/trellis
graphics.
?panel.polygon
-- Bert
On Tue, Oct 4, 2011 at 7:14 AM, Ken Knoblauch wrote:
> markm0705 gmail.com> writes:
> > I would like to plot a string of points
Hi R-Users,
I posted this question a while ago on the bioconductor mailing list but got no
answers. Maybe here is somebody who might know a solution:
I failed at drawing multiple filters in a densityplot() using the
FlowCore/FlowViz packages.
I
found a way to draw multiple filters in xyplot(),
On Mon, Oct 3, 2011 at 12:24 PM, Thomas Adams wrote:
> I'm interested in creating a graphic -like- this:
>
> c <- ggplot(mtcars, aes(qsec, wt))
> c + geom_point() + stat_smooth(fill="blue", colour="darkblue", size=2, alpha
> = 0.2)
>
> but I need to show 2 sets of bands (with different shading) u
On 04/10/2011 6:40 AM, pedabreu wrote:
Hello,
i trying to create a package using package.skeleton. I use R.oo package to
create oriented-object classes. When i use package.skeleton, this creates
the following file:
classA<-
structure(function()
{
extend(Object(),"Clas
markm0705 gmail.com> writes:
> I would like to plot a string of points as a polygon in xyplot. I'm a bit
> lost as to how to get the points plotting in the correct order. I would
> also like some hints on how to render or fill the polygon.
>
> Scrpt below and data file attached
>
> Thanks
>
>
You are right, but this is difficult or impossible to really solve.
The problem is that summary() is an S3 generic(?UseMethod) -- so essentially
it can mean anything and do anything depending on the structure to which
it's applied. In your case, the structures were a data frame and a vector
(that
Sorry---I thought it worked but I think I am actually definitely doing
something wrong...
The problem might be that there are NA's and there are also duplicated
values...My fault. I can't figure out what is going wrong...
I'll be more thorough and modify the two df to mirror more what I have to
ex
Dear all,
I have a directed graph - an Igraph to be more precise - which has some
vertices attributes (like dorm, year etc). Edges and the graph itself do not
have any attributes. Based on the attributes of the vertices, I'd like to
calculate correlation among the edges (e.g. how likely people of t
1 - 100 of 128 matches
Mail list logo
