Soemone broke contrasts<- in 2.10.0 to require the 'methods' package,
but this was corrected in R-patched on Jan 4.
Thus you need to add 'methods' to the default packages in Rscript (see
its help) or update your R.
This has come up before, so a search of the archives was likely to
have found
On 03/23/2010 01:16 PM, rtist wrote:
Hi,
And thanks for the great package. I was using it to output a color table
map for a financial returns matrix.
I like that the fact that I can output a table of values with corresponding
colors, which I was asking about on an earlier thread.
Any way to ge
On Tue, 23 Mar 2010, Jim Lemon wrote:
> On 03/23/2010 10:23 AM, Ray Brownrigg wrote:
> > ...
> >
> > How about less than 2 seconds? [And 500 points in less than 15 seconds -
> > on a 2-year-old DELL Optiplex GX755.]
> >
> > The implementation below (at end) loops over all 'feasible' pairs of x
> >
On Mar 22, 2010, at 9:51 PM, tj wrote:
Hi all,
May I request for your help if you have time and if you have an idea
on how
to do this. I want to add three vectors... And my goal is to obtain
the sum
of the vectors, ignoring the vector of "na"...
Here is what i did in R.. I'm adding the
A fast Fortran solution may be found by:
require(Hmisc)
?largest.empty
Frank
Ray Brownrigg wrote:
On Tue, 23 Mar 2010, Hans W Borchers wrote:
Barry Rowlingson lancaster.ac.uk> writes:
On Mon, Mar 22, 2010 at 4:28 PM, Hans W Borchers
googlemail.com> wrote:
Still I believe that a clever app
k <- c(NA, NA)
z <- c(-12, -9)
e <- c(-2,0)
sum(na.omit(c(z, e, k)))
Will this do the trick? Why are you summing in a sum()
sum(k+z+e)
On Mon, Mar 22, 2010 at 8:51 PM, tj wrote:
>
> Hi all,
> May I request for your help if you have time and if you have an idea on how
> to do this. I want to a
On 03/23/2010 12:48 PM, John K. Williams wrote:
Hello, I'm having some trouble getting things to look as I want with
filled.contour.
1. My first issue is that I am unable to add line segments to my plot where
I want them. Using the rug pattern example:
x<- y<- seq(-4*pi, 4*pi, len = 27);
r<- s
I think you don't need to multiply... you need to figure out where
your NAs need to be and then index the vector and assign "NA" to it:
x[1:5] <- "NA"
On Mon, Mar 22, 2010 at 8:05 PM, tj wrote:
>
> Hi,
>
> How can I create an NA?
> for example, i have this program but it has errors
>
>> x=scan()
Hi tj,
May be the following?
x <- 1:5
NA*x
rep(NA, length(x))
HTH,
Jorge
On Mon, Mar 22, 2010 at 9:05 PM, tj <> wrote:
>
> Hi,
>
> How can I create an NA?
> for example, i have this program but it has errors
>
> > x=scan()
> 1: 1 2 3 4 5
> 6:
> Read 5 items
> > "na"*x
> Error in "na" * x : n
On 03/23/2010 10:23 AM, Ray Brownrigg wrote:
...
How about less than 2 seconds? [And 500 points in less than 15 seconds - on a
2-year-old
DELL Optiplex GX755.]
The implementation below (at end) loops over all 'feasible' pairs of x values,
then
selects the largest rectangle for each pair, su
Erik Iverson wrote:
Hello,
I am using the summary.formula function in the Hmisc package to produce
tables.
With the "method" argument set to "response", the help says,
"Continuous independent variables (see the ‘continuous’ parameter below)
are automatically stratified into ‘g’ (see below)
Hi all,
May I request for your help if you have time and if you have an idea on how
to do this. I want to add three vectors... And my goal is to obtain the sum
of the vectors, ignoring the vector of "na"...
Here is what i did in R.. I'm adding the three vectors, e,z,k, and my
objective is to ge
Hi,
I have a plot that is essentially the same as that in Figure 5.6 of "Lattice
- Multivariate Data Visualization with R". The key difference is that I
would like to add a title to top of the grey color spectrum legend, but have
thus far been unsuccessful. I've tried a variety of options, but to
Hi,
How can I create an NA?
for example, i have this program but it has errors
> x=scan()
1: 1 2 3 4 5
6:
Read 5 items
> "na"*x
Error in "na" * x : non-numeric argument to binary operator
what should I multiply x so that I will get an NA? And I really need to
multiply that with my x vector.
Hi,
And thanks for the great package. I was using it to output a color table
map for a financial returns matrix.
I like that the fact that I can output a table of values with corresponding
colors, which I was asking about on an earlier thread.
Any way to get each box to have a dark black line b
On Mar 22, 2010, at 9:51 PM, Veerappa Chetty wrote:
Hi,
Please help to correct my error in the following. I want to plot the
values
of 'x" in increasing order.
---
library(lattice)
Name<-c("A","B","C")
x<-c(15,20,10)
test<-data.frame(Name,x)
dotplot(Name~x,tes
Sorry, minor tweak to the algorithm in line 16 (thanks Barry).
Looks better now (at end again).
Ray
On Tue, 23 Mar 2010, Ray Brownrigg wrote:
> On Tue, 23 Mar 2010, Hans W Borchers wrote:
> > Barry Rowlingson lancaster.ac.uk> writes:
> > > On Mon, Mar 22, 2010 at 4:28 PM, Hans W Borchers
> > >
On 03/23/2010 01:11 AM, kathy_BJ wrote:
I am new to R, can anyone help with boxplot for a dataset like:
file1 col1 col2 col3 col4 col5
050350005 101 56.625 48.318 RED
051010002 106 50.625 46.990 GREEN
051190007 25 65.875 74.545 BLUE
051191002 246 52.875 57.070 RED
220050004 55 70 80.274 BLUE
220
Clay Heaton wrote:
>
> Hi, pretty new with R and LaTeX here...
>
> I'm trying to use the latex() function from the Hmisc table with Sweave.
> When I try:
>
> pdflatex myfile.tex
>
> I get an error that reads:
>
> ! Undefined control sequence.
> l.22 \ctable
>
>
> It seems like this is an o
On Mar 22, 2010, at 8:31 PM, Rolf Turner wrote:
>
> On 23/03/2010, at 12:35 PM, Clay Heaton wrote:
>
>> Hi, pretty new with R and LaTeX here...
>>
>> I'm trying to use the latex() function from the Hmisc table with Sweave.
>> When I try:
>>
>> pdflatex myfile.tex
>>
>> I get an error that r
Found the solution - already posted here
http://n4.nabble.com/Rscript-question-td1067729.html#a1067729
Basically Rscript does not automatically load library(methods) like Rgui
does...
--
View this message in context:
http://n4.nabble.com/Problem-with-factors-in-lm-when-using-Rscript-but-not-GU
Hi all,
Partly in response to a request for a filled triangular plot, analogous
to image (but more like color2D.matplot), I have uploaded a new version
of plotrix. There is now a triax.fill function that displays such a
plot. I am particularly interested in comment on this as I have more or
le
Hi,
Please help to correct my error in the following. I want to plot the values
of 'x" in increasing order.
---
library(lattice)
Name<-c("A","B","C")
x<-c(15,20,10)
test<-data.frame(Name,x)
dotplot(Name~x,test)
dotplot(reorder(Name,x)~x,test)
"reorder" has no effect
Hello, I'm having some trouble getting things to look as I want with
filled.contour.
1. My first issue is that I am unable to add line segments to my plot where
I want them. Using the rug pattern example:
x <- y <- seq(-4*pi, 4*pi, len = 27);
r <- sqrt(outer(x^2, y^2, "+"));
filled.contour(cos(r
On 23/03/2010, at 12:35 PM, Clay Heaton wrote:
> Hi, pretty new with R and LaTeX here...
>
> I'm trying to use the latex() function from the Hmisc table with Sweave. When
> I try:
>
> pdflatex myfile.tex
>
> I get an error that reads:
>
> ! Undefined control sequence.
> l.22 \ctable
>
>
>
On 23/03/2010, at 2:02 PM, Gundala Viswanath wrote:
> Hi all,
>
> I have a frequency data that looks like this.
>
> 3
> 2
> 1
> 5
>
> What I want to get is the "decreasing" cumulative of this data
> yielding
>
> 11
> 8
> 6
> 5
> 0
>
> Is there any?
>
> I am aware of cumsum(), which will yie
Clay Heaton wrote:
Hi, pretty new with R and LaTeX here...
I'm trying to use the latex() function from the Hmisc table with Sweave. When I
try:
pdflatex myfile.tex
I get an error that reads:
! Undefined control sequence.
l.22 \ctable
What does your function call to latex look like? I di
how about
sum(X) - cumsum(X)
Charles Annis, P.E.
charles.an...@statisticalengineering.com
561-352-9699
http://www.StatisticalEngineering.com
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Gundala Viswanath
Sent: Monday, March
Hi all,
I have a frequency data that looks like this.
3
2
1
5
What I want to get is the "decreasing" cumulative of this data
yielding
11
8
6
5
0
Is there any?
I am aware of cumsum(), which will yield
3
5
6
11.
But it is not what I want.
- G.V.
__
Second Announcement
www.rmetrics.org
*** Rmetrics 2010 ***
4th R/Rmetrics Summer School and User/Developer Meeting on
"Computational Finance and Financial Engineering", Meielisalp, Lake
Thune Switzerland, June 27 - July 1, 2010
The Summer School and User/Develope
Hello R-help
I am having a curious problem. I am using a batch file to call an R script and
am getting an error when performing a lm() fit using factorial data only when
running in batch mode (see below). I do not get the error when executing the
same code in GUI mode in R 2.10.1. Nor do I g
Hi, pretty new with R and LaTeX here...
I'm trying to use the latex() function from the Hmisc table with Sweave. When I
try:
pdflatex myfile.tex
I get an error that reads:
! Undefined control sequence.
l.22 \ctable
It seems like this is an optional package for Tex that I haven't installed? M
On Tue, 23 Mar 2010, Hans W Borchers wrote:
> Barry Rowlingson lancaster.ac.uk> writes:
> > On Mon, Mar 22, 2010 at 4:28 PM, Hans W Borchers
> >
> > googlemail.com> wrote:
> > > Still I believe that a clever approach might be possible avoiding the
> > > need to call a commercial solver. I am gett
On Mar 22, 2010, at 6:11 PM, Janmaat, John wrote:
Hello,
I'm trying to assess the similarity of two lines that are
represented as points (output of differential equation solvers). Is
there a function or a package that deals with things like this?
I think that "things like this" (and "li
mans wrote:
>
> Hi thanks for your help but
>
> I tried to install the pakage but i have got a installing error:
>
> {SNIP}
>
> I do not understand what's going on to be honest.
> Anyone can help me?
>
When debugging, it's usually a good idea to start with the first error:
Error in loadN
Hi JR,
t(EmpCo)
gives what you want. See ?t for details.
You can also do this with reshape:
library(reshape)
EmpCo$Label <- rownames(EmpCo)
recast(EmpCo, variable ~ Label, measure.var=c("TRN", "TRN_CLUST"))
but in this case t() seems a bit quicker and easier.
-Ista
On Mon, Mar 22, 2010 at 6:
Hi:
Is this what you are after? (df is the data frame name)
library(lattice)
dotplot(reorder(Name, rate) ~ rate, data = df)
HTH,
Dennis
On Mon, Mar 22, 2010 at 2:17 PM, Veerappa Chetty wrote:
> Hi ,
> Name
> rate
>
> HEALTHALLIANCE HOSPITALS, INC
Hi thanks for your help but
I tried to install the pakage but i have got a installing error:
nstalling to library ‘/Library/Frameworks/R.framework/Resources/library’
* installing *source* package ‘RInside’ ...
** libs
** arch - x86_64
Error in loadNamespace(name) : there is no package called 'Rc
Hi All,
Consider the following:
TRN<-c(5.809657,3.1, 1.774901e-02)
TRN_CLUST<-c(-4.174682e-05, 5.538742e-05,1.2)
EmpCo<-data.frame(TRN,TRN_CLUST)
row.names(EmpCo)<-c("Slope","Fwy","Univ")
returns:
TRN TRN_CLUST
Slope 5.80965700 -4.174682e-05
Fwy 3.1000 5.538742e-05
Un
You thought you were attaching the data, but have you actually read
the relevant sections in the Posting Guide?
--
David Winsemius, MD
On Mar 22, 2010, at 3:49 PM, Veerappa Chetty wrote:
Hi,
I am trying to make a dot plot in increasing order of the values. It
does
not work. How do I do it?
On 23/03/2010, at 6:03 AM, Barry Rowlingson wrote:
(Commenting on the sort of articles to be found in Computer
Science journals)
> The idea of actually producing some dirty, filthy, actual code to
> implement their shiny algorithms never seems to cross their minds.
Fortune?
Hello,
I'm trying to assess the similarity of two lines that are represented as points
(output of differential equation solvers). Is there a function or a package
that deals with things like this?
Thanks,
John.
Johannus (John) Janmaat
Ass
You need to specify more detail on your functional requirements relative to
protection. R itself does not support the direct use of encrypted data or R
source files.
If you simply want to encrypt/decrypt files before and after use in R, you can
use third party programs such as GnuPG (http://www
Try this:
data.frame(row.names = unlist(a), gset_memberships(a))
On Mon, Mar 22, 2010 at 6:26 PM, Czerminski, Ryszard
wrote:
> I just started using nice package "sets"
> and I wonder if there are utilities to convert (some) sets to data frame
> (as in the example below)
>
>> library(sets)
>> a <
On Mar 22, 2010, at 5:17 PM, Veerappa Chetty wrote:
Hi ,
Name
rate
HEALTHALLIANCE HOSPITALS, INC -1.06211747
MOUNT AUBURN HOSPITAL 0.50960291
STURDY MEMORIAL HOSPITAL2.64233232
For an environmental planning example that involves looking at the
relative efficiencies of one plan over another, I need to determine
the pareto-efficient plans (which I have done), and then, within that
set of plans, determine the convex hull representing the outer upper
boundary of those points.
I just started using nice package "sets"
and I wonder if there are utilities to convert (some) sets to data frame
(as in the example below)
> library(sets)
> a <- gset(elements = list(e('A', 0.1), e('B', 0.8)))
> lst <- as.list(a)
> nr <- length(lst)
> rnames <- character()
> for (i in 1:nr) rname
Hi ,
Name
rate
HEALTHALLIANCE HOSPITALS, INC -1.06211747
MOUNT AUBURN HOSPITAL 0.50960291
STURDY MEMORIAL HOSPITAL2.64233232
LAWRENCE GENERAL HOSPITAL2.15628
On Mar 22, 2010, at 2:00 PM, rkevinbur...@charter.net wrote:
> Thanks to Marc Schultz I found the documentation on the "factors" attribute
> under ?term.object. It stats:
;-)
> factors: A matrix of variables by terms showing which variables appear
> in which terms. The entries are
How about saving, then issuing system() calls to run a zipping program
with password? Then remove the original.
?file.remove
(.R source files are just text.)
--
David.
On Mar 22, 2010, at 4:53 PM, Ron Michael wrote:
Hi, I was following this thread and would like to ask is there any
way to
Hi, I was following this thread and would like to ask is there any way to save
and open a .RData file after using some Password. What I mean to say, how to
make my workplace password-protected?
Also would like to know how same can be done for .R file.
Thanks for your time.
Thanks and regards,
-
Maybe I should simplify the problem with the following smaller table. And I
just want to ask whether
there is any significant difference in the proportion of "Good_Sample"
produced by factories located in "City_A" and "City_B".
Factory_ID Factory_Location Total_Sample Good_Sample
-
Hello,
I am using the summary.formula function in the Hmisc package to produce
tables.
With the "method" argument set to "response", the help says,
"Continuous independent variables (see the ‘continuous’ parameter below)
are automatically stratified into ‘g’ (see below) quantile groups".
B
Another possibility (depending on what you want to do/preferences) is to create
a text file that can be postprocessed to give the colors. One example of this
is the etxtStart and related functions in the TeachingDemos package. These
produce a text file with extra notations that when processed
I need to write out the result oif a tapply
avtaxi = tapply(mdf$TaxiTime, list(mdf$Runway, mdf$OnHour,
mdf$ArrivalGate), FUN=mean, na.rm = TRUE)
to a data file that I can import into Excel.
dim(avtaxi)
[1] 10 24 100
dput(avtaxi, file = outfile, control = c("keepNA", "keepInteger",
"showAttribut
I have a very short question.
Is there any possibility to give to the parameter "col.axis" of graphics
function "axis" a vector value of many colors instead of just one color,
otherwise is there any way around it?
Thank you very much
Enio Jelihovschi
UESC - Brazil
[[alternative HTML versio
Hi,
I am trying to make a dot plot in increasing order of the values. It does
not work. How do I do it? Here are the codes I used. I am also attaching the
data. I use lattice library.
Thanks.
Chetty
--
y<-state.resid$Hospital.Name[state.resid$State=="MA
On Mar 22, 2010, at 1:49 PM, LCOG1 wrote:
Basic question. For the below data, i would like to but each of the
values
in a bin that represents their value. So the below would hopefully
put .1
in the 0-.1 bin, .2 in the .11-.2 bin and so forth. The outlying
values
would then be put into
Barry Rowlingson lancaster.ac.uk> writes:
>
> On Mon, Mar 22, 2010 at 4:28 PM, Hans W Borchers
> googlemail.com> wrote:
>
> > Still I believe that a clever approach might be possible avoiding the need
> > to
> > call a commercial solver. I am getting this hope from one of Jon Bentley's
> > ar
Hello,
I don't know specifics of Xcode, etc ... it looks nice but I have not
used it myself yet.
There seems to be issues on OSX with the current released version of
RInside (we will release 0.2.2 soon), so I would suggest you download
and install the next version of RInside from r-forge:
Basic question. For the below data, i would like to but each of the values
in a bin that represents their value. So the below would hopefully put .1
in the 0-.1 bin, .2 in the .11-.2 bin and so forth. The outlying values
would then be put into and outer category representing everything >1. Im
On Mar 22, 2010, at 3:06 PM, Ista Zahn wrote:
Hi Kamil,
You can use something like
write.csv(t(as.matrix(object)), file="name.csv")
-Ista
On Mon, Mar 22, 2010 at 2:54 PM, Kamil Sijko
wrote:
Hi,
I need to save output of summary() procedure to a csv file. It's all
OK when it's applied to a
Hi Kamil,
You can use something like
write.csv(t(as.matrix(object)), file="name.csv")
-Ista
On Mon, Mar 22, 2010 at 2:54 PM, Kamil Sijko wrote:
> Hi,
>
> I need to save output of summary() procedure to a csv file. It's all
> OK when it's applied to a 'factor' class variable, but when I try to
> s
Thanks to Marc Schultz I found the documentation on the "factors" attribute
under ?term.object. It stats:
factors: A matrix of variables by terms showing which variables appear
in which terms. The entries are 0 if the variable does not
occur in the term, 1 if it does occur an
Hi,
I need to save output of summary() procedure to a csv file. It's all
OK when it's applied to a 'factor' class variable, but when I try to
save a 'integer' class summary to csv it gives me :
> summary(rnorm(100, 10)) -> object
> write.csv2(object, file='name.csv')
Error in do.call("expand.grid
On Mar 22, 2010, at 1:22 PM, Bogaso wrote:
Hi all,
I have saved my workplace in a .RData format. However if I want to
open
that, I need to use following code :
load("C:\\..")
Here my question is why "\\". In all the time generally we use "/"
like when
we use read.delim() function e
On Mon, Mar 22, 2010 at 12:09 PM, Sebastian Leuzinger
wrote:
> Dear list
>
> Has anyone got a recipie at hand to import .bil files into R? From what I
> understand the .bil files I got contain layered matricies which I would lke
> to make available in R as an array or list.
>
> GIS people seem to
Kevin,
See ?terms.object, which is indicated in the Value section of ?terms and listed
in the See Also of ?terms.
HTH,
Marc Schwartz
On Mar 22, 2010, at 1:16 PM, rkevinbur...@charter.net wrote:
>
> I am sorry but I didn't see "factors" mentioned in this documentation.
>
> Kevin
>
> He
Barry Rowlingson wrote:
>
> GIS and spatial data formats can often be handled by readGDAL (for
> raster grids) from the rgdal package.
>
> .bil files seem to be handled by the Ehdr driver in GDAL:
>
> http://www.gdal.org/frmt_various.html
>
> so if your rgdal package has that driver (run gd
I am sorry but I didn't see "factors" mentioned in this documentation.
Kevin
Henrique Dallazuanna wrote:
> See ?terms
>
> On Mon, Mar 22, 2010 at 2:08 PM, wrote:
> > I noticed that when I fit a linear model using 'lm' there is an attribute
> > called "factors" that is added to the "te
Thank you very much Joshua. I was thinking to use logistic regression with
glm(). But this will pool the individual factories which share the same
factor levels together. I was puzzled by how to deal with individual
factory. Any idea? I will try your method anyway.
Xiang
On Mon, Mar 22, 2010 at
mans wrote:
>
> Hi,
> Can anyone tell me how to embed R code in a C++ file.
>
> I am actually using a mac running on the OSX 10.6.2 and the IDE Xcode
> Version 3.2 and I would like to embed the basic function like geometric,
> binomial, normal and hyper geometric distributions in a sample cp
I am sorry but I didn't see "factors" mentioned in this documentation.
Kevin
Henrique Dallazuanna wrote:
> See ?terms
>
> On Mon, Mar 22, 2010 at 2:08 PM, wrote:
> > I noticed that when I fit a linear model using 'lm' there is an attribute
> > called "factors" that is added to the "te
Adding text at the top of the bars will tend to distort the perception of their
heights. It is better to place numbers (if they are even needed) in the
margins. Switching to a dotplot instead of a barplot may be more meaningful as
well.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Dear Xiang,
Now I understand what you meant. If you are only interested in
comparing the Good Samples, I think you would have to use the
proportion (Good Sample/Total Sample) or something similar. Another
thought would be to dummy code the data (e.g., Good = +1, Fair = 0,
Bad = -1). Then you co
On Mon, Mar 22, 2010 at 5:22 PM, Bogaso wrote:
>
> Hi all,
>
> I have saved my workplace in a .RData format. However if I want to open
> that, I need to use following code :
>
> load("C:\\..")
>
> Here my question is why "\\". In all the time generally we use "/" like when
> we use read.delim(
See ?terms
On Mon, Mar 22, 2010 at 2:08 PM, wrote:
> I noticed that when I fit a linear model using 'lm' there is an attribute
> called "factors" that is added to the "term". It doesn't seem to appear for
> 'model.matrix', just 'lm'. I have been unable to find where it gets
> constructed or w
Hi all,
I have saved my workplace in a .RData format. However if I want to open
that, I need to use following code :
load("C:\\..")
Here my question is why "\\". In all the time generally we use "/" like when
we use read.delim() function etc. Is there any possibility to have some
consistenc
> Matteo Bertini schrieb:
>>
>> Hello everyone,
>>
>> I'm doing some benchmark comparing Arima [1] and SVR on time series data.
>>
>> I'm using an out-of-sample one-step-ahead prediction from Arima using
>> the "fitted" method [2].
>>
>> Do someone know how to have a two-steps-ahead forecast timese
Dear R users,
I use mclust to fit a mixture of normal distributions to many datasets. Usually
the Mclust function finds 1 or two normal distributions, rarely, 3.
But I hit a strange case today.
my.data <- c(57.96920, 51.79415, 51.20538, 55.53637, 51.64291, 56.61476,
51.28855, 55.56169, 51.8511
I noticed that when I fit a linear model using 'lm' there is an attribute
called "factors" that is added to the "term". It doesn't seem to appear for
'model.matrix', just 'lm'. I have been unable to find where it gets constructed
or what it means? It looks like a two dimensional array that I may
On Mon, Mar 22, 2010 at 4:28 PM, Hans W Borchers
wrote:
> Still I believe that a clever approach might be possible avoiding the need to
> call a commercial solver. I am getting this hope from one of Jon Bentley's
> articles in the series Programming Pearls.
>
Is this the 'Largest Empty Rectangle
Hi everyone,
I just conducted a structural equation model for estimating a response
model. This model should predict the probability that someone is responding
to a direct mailing. I used the sem package for this. When I have my
coefficients I want to know how well my model predicts the probabilit
Hello everyone...
We were trying to implement the Newton-Raphson method in R, and estimate the
parameters a and b, of a function, F, however we can't seem to implement
this the right way. Hope you can show me the right way to do this. I think
what we want R to do is to read the data from the websi
I'm recoding variables and running a logit. Unfortunately, I get the following
error.
data04$V043114
part <- data04$V043114
attributes(part)
summary(part)
partb < part
partb[part %in% levels(part)[4]] <- NA
partb[part %in% levels(part)[5]] <- NA
partb[part %in% levels(part)[6]] <- NA
partb[p
I am new to R, can anyone help with boxplot for a dataset like:
file1 col1 col2 col3 col4 col5
050350005 101 56.625 48.318 RED
051010002 106 50.625 46.990 GREEN
051190007 25 65.875 74.545 BLUE
051191002 246 52.875 57.070 RED
220050004 55 70 80.274 BLUE
220150008 75 67.750 62.749 RED
220170001 77 6
It seems to mean that it adds a Tsp attribute but it does not change
the class to "ts":
> dput(lag(1:3))
structure(1:3, .Tsp = c(0, 2, 1))
Try this:
> ts(1:3) - structure(lag(1:3), class = "ts")
Time Series:
Start = 1
End = 2
Frequency = 1
[1] -1 -1
or
> ts(1:3) - lag(ts(1:3))
Time Series:
Sta
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of William Dunlap
> Sent: Monday, March 22, 2010 9:31 AM
> To: Martin Batholdy; r help
> Subject: Re: [R] add information above bars of a barplot()
>
> > -Original Message-
Hi all,
I am wondering if there is a way to plot symbols onto the leaves of a
dendrogram. Thanks for any help.
Wade
__
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
Hello,
Downey, Patrick wrote:
Can anyone tell me what's going on here?
x <- matrix(data=c(1,2,3,4,5),ncol=1)
x1 <- lag(x,k=1)
x
x1
x - x1
That's with x specified as a column vector, but the same thing happens when
it's a row vector.
x <- c(1,2,3,4,5)
x1 <- lag(x,k=1)
x
x1
x - x1
I'm not su
Hans W Borchers googlemail.com> writes:
>
> For an application in image processing -- using R for statistical purposes --
> I need to solve the following task:
>
> Given n (e.g. n = 100 or 200) points in the unit square, more or less
> randomly distributed. Find a rectangle of maximal area withi
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Martin Batholdy
> Sent: Monday, March 22, 2010 7:53 AM
> To: r help
> Subject: [R] add information above bars of a barplot()
>
> hi,
>
>
> I have a barplot with six clusters of
Or use the replicate function (which is basically a wrapper for lapply).
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
> project.org] O
Can anyone tell me what's going on here?
x <- matrix(data=c(1,2,3,4,5),ncol=1)
x1 <- lag(x,k=1)
x
x1
x - x1
That's with x specified as a column vector, but the same thing happens when
it's a row vector.
x <- c(1,2,3,4,5)
x1 <- lag(x,k=1)
x
x1
x - x1
When the documentation says "Vector or matrix
Hi,
Can anyone tell me how to embed R code in a C++ file.
I am actually using a mac running on the OSX 10.6.2 and the IDE Xcode
Version 3.2 and I would like to embed the basic function like geometric,
binomial, normal and hyper geometric distributions in a sample cpp file.
I heard about the l
Dear Joshua,
Thank you so much for such fast reply. Here is my thought:
I don't know if it is fair to compare means because the total samples from
each factory can be very different (like, In Factory_5 with 150 total
samples vs. Factory_9 with 70 total samples). Maybe it is more fair to
compare f
Did you back-calculate to estimate an intercept?
Alternately, I've been working on a function that takes a fitted sem
and gets predicted values given an input. Contact me off-list and
I'll send it to you.
On Mar 22, 2010, at 8:37 AM, Tryntsje Wesselius wrote:
Hi everyone,
I just conducte
-Mensaje original-
De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] En
nombre de Martin Batholdy
Enviado el: lunes, 22 de marzo de 2010 15:53
Para: r help
Asunto: [R] add information above bars of a barplot()
hi,
I have a barplot with six clusters of four bars eac
Dear Xiang,
Unequal sample size is not a problem for t-tests. If I understand
correctly, you do not want to pool your data because you believe the
variance of individual factories is heterogenous. Are you willing to
pool the means? You could calculate the variance for factories
individually and
Hi everyone,
I just conducted a structural equation model for estimating a response
model. This model should predict the probability that someone is responding
to a direct mailing. I used the sem package for this. When I have my
coefficients I want to know how well my model predicts the probabili
1 - 100 of 136 matches
Mail list logo
