I have everything in the data frame now. The calculations I need to
preform for each round are slightly different. How do I specify to make
one calculation for one round and one for the other? (One round is
called "C", and the other is "F")
Thomas Levine
On Thu, 2008-01-24 at 18:23 -0500, Thomas
This is a text part of the message.
It is shown for the users of old-style e-mail clients
[[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 http
On 1/25/08, Spilak,Jacqueline [Edm] <[EMAIL PROTECTED]> wrote:
>
> -Original Message-
> From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
> Sent: Thursday, January 24, 2008 9:56 AM
> To: Spilak,Jacqueline [Edm]
> Cc: r-help@r-project.org
> Subject: *SPAM* Re: [R] Bar width and labels
On 26/01/2008, at 10:54 AM, Carson Farmer wrote:
> Dear List,
>
> I am attempting to perform a harmonic analysis on a time series of
> snow
> depth, in which the annual curve is essentially asymmetric (i.e. snow
> accumulates slowly over time, and the subsequent melt occurs
> relatively
> rap
Hi, The following code, from Angelo Canty article on line
"Resampling Methods in R: the boot Package, 2002", works fine for Angelo Canty
using R 2.6.0 on Windows XP.
It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after
installing the S-PLUS bootstrap
I think you can do this fairly easily with the reshape
package. Is this basically what you want?
x <- "Plot1 Species1Abundance1
Plot1 Species2Abundance2
Plot2 Species1Abundance1
Plot2 Species3Abundance3 "
xx <- read.table(textConnection(x), header=TRUE,
as.is=TRUE); xx
names
Take a look at ?reshape - it is used to convert between 'long' and 'wide'
formats of data tables (frames).
-Christos
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of
> [EMAIL PROTECTED]
> Sent: Friday, January 25, 2008 4:30 PM
> To: r-help@r-proj
Dear List,
I am attempting to perform a harmonic analysis on a time series of snow
depth, in which the annual curve is essentially asymmetric (i.e. snow
accumulates slowly over time, and the subsequent melt occurs relatively
rapidly). I am trying to fit a curve to the data, however, the actual
Hello.
Does anyone know of an existing function that takes a list in the form of:
Plot1 Species1Abundance1
Plot1 Species2Abundance2
Plot2 Species1Abundance1
Plot2 Species3Abundance3
.
.
.
PlotN SpeciesNAbundanceN
and translates into a matrix in the form of
Species1
Hi, I'm trying to set up AR(1) as a correlation structure in modeling some
data (attached file data.txt in text format) with lme, but have trouble
getting it to work.
Incent, Correctness, and Oppor are 3 categorical variables, Beta is a
response variable, and Time is an equally-spaced variable wit
Hi Jarek,
an alternative approach is to provide more precise starting values!
It pays off to realise that it's possible to find quite good guesses
for some of the parameters in your model function:
t ~ tr+(ts-tr)/((1+(a*h)^n)^(1-(1/n)))
The parameters tr and ts correspond to the re
YAHOO.Shortcuts.hasSensitiveText = false; YAHOO.Shortcuts.sensitivityType =
[]; YAHOO.Shortcuts.doUlt = false; YAHOO.Shortcuts.location = "us";
YAHOO.Shortcuts.lang = "us"; YAHOO.Shortcuts.document_id = 0;
YAHOO.Shortcuts.document_type = ""; YAHOO.Shortcuts.document_title = "";
YAHOO.Shortcu
On Jan 25, 2008 2:04 PM, Domenico Vistocco <[EMAIL PROTECTED]> wrote:
> You could use ggplot2:
>
> library(ggplot2)
> x=rnorm(100)
> y=rnorm(100)
> df=melt(data.frame(x,y))
> #for "vertical" histogram
> ggplot(data=df,aes(x=value))+geom_histogram()+facet_grid(.~variable)
> #for "horizontal" histogr
You could use ggplot2:
library(ggplot2)
x=rnorm(100)
y=rnorm(100)
df=melt(data.frame(x,y))
#for "vertical" histogram
ggplot(data=df,aes(x=value))+geom_histogram()+facet_grid(.~variable)
#for "horizontal" histogram
ggplot(data=df,aes(x=value))+geom_histogram()+facet_grid(.~variable) +
coord_flip
Hello, I would like to make my R script more interactive. Well, I want the
script to work like this:
I run the script from R Console and it offers me some options (press 1 for
something, or press 2 for something else). After pressing 1 or 2, do
corresponding procedure.
I am not able to find out
Benilton Carvalho wrote:
> I'm afraid others might have some other recommendations regarding your
> message... but the first one is that the bioconductor mailing list is
> more appropriate to this...
1. We all three should note that we should use sensible subject lines.
2. You, Benilton, should
Prateek Singh wrote:
> Hi,
>
> I got command to execute R code from Mark.
>
> but it was showing error..
> Error in library(hopach) : there is no package called 'hopach'
>
> Does it mean this library hopach
hopach is a package.
> is not installed. But I have installed complete Bioconductor
I'm afraid others might have some other recommendations regarding your
message... but the first one is that the bioconductor mailing list is
more appropriate to this...
Apparently you didn't install all the BioConductor packages (and I'm
glad you didn't :-) ).
Just install the package you
Hi,
I got command to execute R code from Mark.
but it was showing error..
Error in library(hopach) : there is no package called 'hopach'
Does it mean this library hopach is not installed. But I have
installed complete Bioconductor package.
I am confused.
can you help me in this
thanks
Regards,
Hi,
one possibility is to start R in a terminal window by entering
R
then R should be started.
now simply enter:
source(exmaple.r)
and everything should be fine.
Hope this helps,
Roland
Prateek Singh wrote:
> Hi,
>
> I have 'R' code in file. saved it as exmaple.r
> here is the code...
>
>
You might look at the method histbackback (back-to-back histograms) in
package Hmisc. At the time I looked at this package about a year ago,
I could not get it working on my Mac, so I don't actually know what
the results look like, but I have created them before in XLispStat and
found them
http://cran.r-project.org/src/contrib/Descriptions/popbio.html
On Jan 25, 2008, at 10:51 AM, Michelle DePrenger-Levin wrote:
> Hello,
>
>
>
> I'm in need of running a PVA and want to load {popbio}. However, I
> can't find
> it. I can't see where or how to load it. I have added repositories and
http://cran.r-project.org/src/contrib/Descriptions/popbio.html
On 25/01/2008, Michelle DePrenger-Levin <[EMAIL PROTECTED]> wrote:
> Hello,
>
>
>
> I'm in need of running a PVA and want to load {popbio}. However, I can't find
> it. I can't see where or how to load it. I have added repositories and
Hello,
I'm in need of running a PVA and want to load {popbio}. However, I can't find
it. I can't see where or how to load it. I have added repositories and
changed the CRAN mirror site. I see {polynom} and {popgen} but no {popbio}
between.
What am I doing wrong?
Thanks,
Michelle
marinercapital.com> writes:
>
> I need to write a C# program to create and display Cusum chart from any
> of the packages,
> spc, qcc or strucchange.
>
> 1-The data resides in a MS SQL Database. The C# program will handle
> obtaining the data for the requisite types of samples.
It is easy t
Rick DeShon msu.edu> writes:
..
> When I use nlsList code that is virtually identical to the example in
> the doc file I get the following error.
..
> "Error in nlsList.formula(RT ~ SSasymp(trial, Asym, R0, lrc),
> data = df.group, :
> Data must be a groupedData object if formula d
On 1/25/2008 12:55 PM, Remund, Todd wrote:
> I have had a strange problem with rgl. I am currently switching from a
> desktop computer to a workstation. When I get a very large data set and
> run spectro3D, the computations run then the window with the 3d plot
> appears it doesn't show the da
Or you can use lattice's xyplot(..., type = c("p","g"))
?panel.xyplot provides details
-- Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of John Kane
Sent: Friday, January 25, 2008 9:49 AM
To: John Lande; r-help@r-project.org
Subject
Hi All.
I'm trying to run nlsList an getting an error that makes no sense to
me. I have accuracy and reaction time data over many trials for each
person (id)
When I use nlsList code that is virtually identical to the example in
the doc file I get the following error. I've tried everything I cou
I have had a strange problem with rgl. I am currently switching from a
desktop computer to a workstation. When I get a very large data set and
run spectro3D, the computations run then the window with the 3d plot
appears it doesn't show the data inside and immediately shuts down R.
The computers h
?grid
--- John Lande <[EMAIL PROTECTED]> wrote:
> dear all,
>
> I have a very simple question but I could not figure
> out.
>
> I need to make plots with grid in the background.
>
> something like I old retrive like this
>
> a=runif(100)*10
> b=runif(100)*10
> plot(a,b, pch=20, xlim=c(0, rou
Hi All,
Thanjuvar wrote:
>> model2<-lm(lavi~age+sex+age*race+diabetes+hypertension, data=tb1)
David wrote:
>> in the second equation you are only including the interaction of
>> age*race,
>> the main effect of age, but not the main effect of race which is what
>> came out significant
>> in y
-Original Message-
From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
Sent: Thursday, January 24, 2008 9:56 AM
To: Spilak,Jacqueline [Edm]
Cc: r-help@r-project.org
Subject: *SPAM* Re: [R] Bar width and labels in barchart
On 1/23/08, Spilak,Jacqueline [Edm] <[EMAIL PROTECTED]> wrote
Hi,
I have 'R' code in file. saved it as exmaple.r
here is the code...
library("hopach")
GSE <- read.table("gene_expression_data",sep="\t",header=TRUE,row.names=TRUE)
gene.dist <- distancematrix(t(GSE),d="euclid")
gene.hobj <- hopach(t(GSE), dmat=gene.dist, mss="med")
labelstosil(gene.hobj$cl
On Wed, 23 Jan 2008, Daniel Oberski wrote:
> Hello
>
> Does anybody happen to know if it is possible to use the survey package to
> estimate a covariance matrix from a complex survey?
Yes. svyvar() in the survey package estimates a population covariance
matrix.
It doesn't give standard errors f
On 1/25/2008 11:39 AM, Karin Lagesen wrote:
> I have a data frame containing columns which are factors. I use this
> to make boxplots for the data, with one box per factor. I would now
> like to get at the data in the data frame which corresponds to the
> outliers. I have so far found the $out, whi
With thanks to all who responded, especially Charles Berry and
Don MacQueen, my student has everything working. The final
solution was:
1. Install X11 and X11SDK (or make sure they are installed).
2. There is a button for starting X11 within the R GUI.
3. Sys.setenv( PATH=paste(Sys.getenv("PATH")
I have a data frame containing columns which are factors. I use this
to make boxplots for the data, with one box per factor. I would now
like to get at the data in the data frame which corresponds to the
outliers. I have so far found the $out, which gives "the values of any
data points which lie b
Hi, I'm trying to set up AR(1) as a correlation structure in modeling some
data (attached file data.txt in text format) with lme, but have trouble
getting it to work.
Incent, Correctness, and Oppor are 3 categorical variables, Beta is a
response variable, and Time is an equally-spaced variable wit
What about reading the last four lines of this email?
G.
On Fri, Jan 25, 2008 at 04:45:04PM +0100, [EMAIL PROTECTED] wrote:
> How can I join two lists? I have q1 and q2 and I want to merge them. I
[...]
> __
> R-help@r-project.org mailing list
> http
Try c(q1,q2)
Nael
On Jan 25, 2008 4:45 PM, <[EMAIL PROTECTED]> wrote:
> How can I join two lists? I have q1 and q2 and I want to merge them. I
> have tried to use the comand merge, but not work. Any solutions? Thanks!
>
> > q1
> $Input1
> 7.84615384615385
> 0.5
>
> $Input2
> 8.92307
How can I join two lists? I have q1 and q2 and I want to merge them. I
have tried to use the comand merge, but not work. Any solutions? Thanks!
> q1
$Input1
7.84615384615385
0.5
$Input2
8.92307692307692
-3.2
$Input3
4.53846153846154
-5
> q2
$Input1
7.
Dear Svempa,
See ?legend.
In this case, you could have discovered the legend function via
apropos("legend") or help.search("legend").
I hope this helps,
John
John Fox, Professor
Department of Sociology
McMaster University
Hamilton, Ontario, Canada L8S 4M4
905-5
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
On Fri, 2008-01-25 at 14:21 +, David Winsemius wrote:
> "Paul Sweeting" <[EMAIL PROTECTED]> wrote in
> news:[EMAIL PROTECTED]:
> The most straightforward approach would be to use glm() with
> family="poisson", a.k.a Poisson regression. Don't ask
On 1/25/2008 9:46 AM, Luis Ridao Cruz wrote:
> R-help,
>
> Sorry for this question (I guess it has been addressed before but
> I could not find it in the archives)
> but how can I see a function code
> when the following comes up:
>
>> svymean
> function (x, design, na.rm = FALSE, ...)
> {
>
hits=-2.6 tests�YES_00
X-USF-Spam-Flag: NO
On Fri, 2008-01-25 at 14:46 +, Luis Ridao Cruz wrote:
> R-help,
>
> Sorry for this question (I guess it has been addressed before but
> I could not find it in the archives)
> but how can I see a function code
> when the following comes up:
>
> > sv
Hi,
This IS the code of the function. It is a generic function, the actual method
used depends on the class of the object.
You will find methods for this function by doing:
R> methods("svymean")
[1] svymean.survey.design* svymean.survey.design2* svymean.svyrep.design*
[4] svymean.twophase*
a*b is not an interaction term. An interaction term looks like a:b.
y~age*race means
y~age+race+age:race.
For
TB> model1<-lm(lavi~age+sex+race+diabetes+hypertension, data=tb1)
TB> model2<-lm(lavi~age+sex+age*race+diabetes+hypertension, data=tb1)
both main effects are present in model 2 (below) a
How can I test whether a plot I'm about to call will be within a "layout" or
not? I need to fix some par options based on that ...
Thanks, Joh
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting g
From: Bill Vorias
> Is there an Interaction function in the Random Forest
> algorithm in the current version of R?
No (at least not in the randomForest package). Prof. Breiman regarded
that feature as "experimental", and said he wouldn't mind not having it
in the R version. I don't know if Pr
R-help,
Sorry for this question (I guess it has been addressed before but
I could not find it in the archives)
but how can I see a function code
when the following comes up:
> svymean
function (x, design, na.rm = FALSE, ...)
{
.svycheck(design)
UseMethod("svymean", design)
}
Thanks i
I need to write a C# program to create and display Cusum chart from any
of the packages,
spc, qcc or strucchange.
Issues:
1-The data resides in a MS SQL Database. The C# program will handle
obtaining the data for the requisite types of samples.
Assistance needed on:
1-How can I call the cusum
"Paul Sweeting" <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> Hi
>
> I am trying to carry out some maximum likelihood estimation and I'm
> not making much headway, and I'm hoping that someone will be able to
> point me in the right direction.
>
> I am modelling mortality statistics. O
Tirthadeep wrote:
>
> Dear Users,
>
> In case of ridge logistic regression, i want to calculate the optimum
> penalty using aic and bic criteria. Here is the sample code:
>
> fit <- lrm(RES ~CAT01+NUM01+NUM02+CAT02+CAT03+CAT04+NUM03+CAT05+CAT06+NUM04+
> CAT07+CAT08+NUM05+NUM06, dat
Visit this site for free classifieds http://www.a2z-4-sale.com, have nice week
end
__
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
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
See ?grid though read the Note section of that page, which says that
more fine tuning can be achieved through the use of abline directly.
Also, note that abline is vectorized so you don't need the loop. Here is
you example, adjusted to avoid the large
On 1/25/2008 8:14 AM, John Lande wrote:
> dear all,
>
> I have a very simple question but I could not figure out.
>
> I need to make plots with grid in the background.
>
> something like I old retrive like this
>
> a=runif(100)*10
> b=runif(100)*10
> plot(a,b, pch=20, xlim=c(0, round(max(a))),
dear all,
I have a very simple question but I could not figure out.
I need to make plots with grid in the background.
something like I old retrive like this
a=runif(100)*10
b=runif(100)*10
plot(a,b, pch=20, xlim=c(0, round(max(a))), ylim=c(0, round(max(b
vs=seq(0, max(a), 0.5)
for(i in 1:le
Hi,
I have been trying to run R programs on the XGrid but without any
sucess...it would be great help if somebody could give me an
indication/steps as to how to get this running...am using a Mac OSX Leopard
and the xgrid is on a tiger cluster...
thanks very much in advance!
[[alternative
Hi
[EMAIL PROTECTED] napsal dne 24.01.2008 23:21:05:
> I need to group some subjects according to two
> variables (var1, var2)
>
> I have a data frame that looks something like this:
>
> Subject var1 var2
> 1 400150
> 2 320100
> 3 580600
> 4
Hi
I am trying to carry out some maximum likelihood estimation and I'm not
making much headway, and I'm hoping that someone will be able to point me in
the right direction.
I am modelling mortality statistics. One way to do this is to model the
mortality rate (or, more accurately, log of the mor
Dear R-Users,
I am working on an Hierarchical Bayes model and tried to replace the inner
for-loop (which loops over a list with n.observations elements) with truely
vectorized code (where I calculated everything based on ONE dataset over all
respondents).
However, when comparing the performance
Hi
I am carrying out some logit regressions and want to (a) make sure I'm
taking the right approach and (b) work out how to carry out some additional
analysis. So, to carry out a logit regression where the dependent variable
is a factor db, I use something like:
res1_l <- glm(formula = db ~ y1 +
I have one graph with two linear plots, one where the points are marked with
a triangle and the others marked with a dot. How do I label them in a legend
for example, comparable to what you automatically get in Excel (sorry for
the comparison people), so you know which type of point corresponds to
Sergey Goriatchev wrote:
> Dear members of R forum,
>
> Say I have a list:
>
> L <- list(1:3, 1:3, 1:3)
>
> that I want to turn into a matrix.
>
> I wonder why if I do:
>
> do.call(cbind, L)
>
> I get the matrix I want, but if I do
>
> cbind(L)
>
> I get something different from what I want. Why is
Thanjavur,
I'm new to R, so it is possible I'm interpreting you syntax
incorrectly, but it looks like in the second equation you are only
including the interaction of age*race, the main effect of age, but
not the main effect of race which is what came out significant in your
first model.
In effec
On Jan 25, 2008 11:27 AM, Sergey Goriatchev <[EMAIL PROTECTED]> wrote:
> Dear members of R forum,
>
> Say I have a list:
>
> L <- list(1:3, 1:3, 1:3)
>
> that I want to turn into a matrix.
>
> I wonder why if I do:
>
> do.call(cbind, L)
>
> I get the matrix I want, but if I do
>
> cbind(L)
>
> I ge
Use Rprof on your code to see where time is being spend and then focus
on those areas for improvement.
On Jan 25, 2008 2:21 AM, Juliet Hannah <[EMAIL PROTECTED]> wrote:
> Dear R Programmers,
> I am trying to run a Poisson regression on all pairs of variables in a data
> set and
> obtain the permut
Perhaps there are other opinions, but from my point
of view R is spectacularly similar in Linux and Windows
(and will be virtually identical in Linux and Unix). As
long as you are doing "R" things, you should be in home
territory.
I regularly transfer R objects between 64-bit Linux and
Windows, a
Dear members of R forum,
Say I have a list:
L <- list(1:3, 1:3, 1:3)
that I want to turn into a matrix.
I wonder why if I do:
do.call(cbind, L)
I get the matrix I want, but if I do
cbind(L)
I get something different from what I want. Why is that? How does
do.call() actually work?
I've read
Berwin A Turlach wrote:
> Of course, you could always do something like:
>
>
>> ilapply <- function(...) invisible(lapply(...))
>>
>
> ## perhaps better:
> ## ilapply <- function(X, FUN, ...) invisible(lapply(X, FUN, ...))
>
>
>> ilapply(p, function(ps) {x11(); plot(ps)})
>>
lapply
[EMAIL PROTECTED] wrote:
>
>Hi,
>
>
>i'm very new to R, so sorry for what i'm sure is a very basic question. I'm
>producing a boxplot with the data below:
>
>df3<-data.frame(
>
> x=c(10,11,115,12,13,14,16,17,18,21,22,23,24,26,27,28,29,3,30,32,33,34,35,4,4
>1,45,5,50,52,5
Asche auf mein Haupt.
Yours,
Mag. Ferri Leberl
Am Donnerstag, den 24.01.2008, 17:12 + schrieb Prof Brian Ripley:
> I suspect you have ignored
>
>'If you are using an old version of R and think it does not work
>properly, upgrade to the latest version and try that, before posting.'
>
Hello R-ill world!
My problem : I obtain a segmentation fault when passing a character
argument to a C function, compiled in a shared object and loaded by
dyn.load.
1. As manuals (and overall "Writing R extensions ") don't seem to
mention it or I failed to find the info, could someone explain
Berwin A Turlach wrote:
> G'day Rainer,
>
> On Fri, 25 Jan 2008 10:34:32 +0200
> Rainer M Krug <[EMAIL PROTECTED]> wrote:
>
> [...]
>
>> > p <- data.frame(runif(10), runif(10), runif(10))
>> > lapply( p, function(ps) {x11(); plot(ps)} )
>>
>> which results in three graphs and a printout:
> [.
Juan Pablo Fededa wrote:
> Dear Contributors:
>
> I have two vectors x and z, and I want to display the histograms of both
> vectors in the same graph, x in red bars, z in blue bars.
> If you have any clue on how to do that, I will be very glad to hear it!!
> Thanks in advance again,
>
Hi Jua
G'day Rainer,
On Fri, 25 Jan 2008 10:34:32 +0200
Rainer M Krug <[EMAIL PROTECTED]> wrote:
[...]
> > p <- data.frame(runif(10), runif(10), runif(10))
> > lapply( p, function(ps) {x11(); plot(ps)} )
>
> which results in three graphs and a printout:
[...]
> How can I avoid this printout without
Hello,
I want to compute OR from a simulation,
I simulate outcome-y and binary covariate, and I repeat it 1000 times.
I built a logistic model, and want to compute OR for every sample
At the end I'm looking for the mean(OR) of the 1000 repetitions.
Can I use cc (from epicalc)?
example: if yall is t
Lightning Protection Technology
Micro Power engineering (Pvt.) Ltd is one of the world's leading
companies in Lightning Protection using the most up to date
technology available. We offer a complete range of Lightning
Protection equipment to produce safe and reliable Lightning and Surge
Protect
Tobias Teichert wrote:
> Hi Everybody,
>
> I need to install the R-package 'foreign'. However, when I try the
> standard install.packages('foreign') I get the following error
> message:
>
> " package 'foreign' is not available in: install.packages("foreign") "
>
> The funny thing is, that I am usi
Hi
I am looking for a function which is executing some code repeatedly,
exactly like lapply() and sappl() are doing, but do not return anything
as I am not interested in the return values.
An example would be:
> p <- data.frame(runif(10), runif(10), runif(10))
> lapply( p, function(ps) {x11(
Dear Users,
In case of ridge logistic regression, i want to calculate the optimum
penalty using aic and bic criteria. Here is the sample code:
fit <- lrm(RES ~CAT01+NUM01+NUM02+CAT02+CAT03+CAT04+NUM03+CAT05+CAT06+NUM04+
CAT07+CAT08+NUM05+NUM06, data = train.data, x = TRUE, y =
TR
Dear Users,
In case of ridge logistic regression, i want to calculate the optimum
penalty using aic and bic criteria. Here is the sample code:
fit <- lrm(RES ~CAT01+NUM01+NUM02+CAT02+CAT03+CAT04+NUM03+CAT05+CAT06+NUM04+
CAT07+CAT08+NUM05+NUM06, data = train.data, x = TRUE, y = T
83 matches
Mail list logo
