Thanks Berwin, Spencer, and Gabor!!!
Berwin A Turlach wrote:
G'day all,
On Sat, 09 May 2009 08:01:40 -0700
spencerg spencer.gra...@prodsyse.com wrote:
The harmonic mean is exp(mean(logs)). Therefore, log(harmonic
mean) = mean(logs).
Does this make sense?
I think
?anova.coxph
will tell you that there's an additional parameter, test, taking
values F, Cp, or Chisq which instructs the anova method to
perform the stated test comparing the two models and spit out a p-
value (for F and Chisq at least).
example(anova.chisq) provides some examples.
Cheers,
Can anyone tell me what is skip=2, skip =7 and %in% mean here?
fromcsv=read.csv
('2_2005_top200_postdoc.csv',header=FALSE,skip=7,stringsAsFactors=FALSE)
fromreadxls = read.xls
('2_2005_top200_postdoc.xls',stringsAsFactors=FALSE,header=FALSE,skip=2)
fromreadxls = fromreadxls[fromreadxls$V7
Dear List,
First off, my deepest gratitude to the Sweave developers: this tool has
improved my quality greatly.
A question in my work I use \Sexpr{} statements scalar values and the xtable
package for all manner of tables. What I'd like to do is to use a vector
inline, rather than a whole
How do you get the design matrix R used when calculating ANOVA?
--
View this message in context:
http://www.nabble.com/how-to-get-design-matrix--tp23464638p23464638.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org
Stephen:
I am trying to use the plm package for panel econometrics. I am just
trying to get started and load my data.
Look at the Journal of Statistical Software paper that introduces the
plm package:
http://www.jstatsoft.org/v27/i02/
It seems from most of the
sample documentation that
Hello,
I've written a relatively complex interconnected population model in R.
When changing a certain parameter, the outputs end up containing NA's.
I would like to find out, in which step the model (in form of a loop)
starts to produce NA's. Does anyone know how to achieve this?
Since it
Hi all,
I wanted to let you know about our training seminar on predictive analytics
- coming May, Oct, and Nov in NYC, Stockholm, DC and other cities. This is
intensive training for marketers, managers and business professionals to
make actionable sense of customer data by predicting buying
Kyle Matoba kmmatoba at ucdavis.edu writes:
A question in my work I use \Sexpr{} statements scalar values and the xtable
package for all manner of tables. What I'd like to do is to use a vector
inline, rather than a whole separate table. Something like:
% Sweave block:
=
covmat -
Tung86 wyutung at gmail.com writes:
Can anyone tell me what is skip=2, skip =7 and %in% mean here?
fromcsv=read.csv
('2_2005_top200_postdoc.csv',header=FALSE,skip=7,stringsAsFactors=FALSE)
Did you check the docs?
skipinteger: the number of lines of the data file to skip before
Hi Kyle,
First off, my deepest gratitude to the Sweave developers: this tool has
improved my quality greatly.
A question in my work I use \Sexpr{} statements scalar values and the xtable
package for all manner of tables. What I'd like to do is to use a vector
inline, rather than a whole
Richard Chirgwin wrote:
Hi all,
I am trying to install Spatstat on OpenSUSE 11.1.
install.packages(spatstat, dependencies = TRUE)
fails on the basis of various compiler packages (full message below).
I have gcc version 4.3.2, which should include gfortran and g++ - so I'm not
Felipe Carrillo mazatlanmexico at yahoo.com writes:
I have created a MS Access table named 'PredictedValues' through the statement
below:
myDB - odbcConnectAccess(C:/Documents and Settings/Owner/Desktop/Rpond
Farming.mdb,uid=admin,pwd=)
sqlSave(myDB,PredictedValues,rownames=FALSE)
On 10/05/2009 1:01 AM, Benjamin Leutner wrote:
Hello,
I've written a relatively complex interconnected population model in R.
When changing a certain parameter, the outputs end up containing NA's.
I would like to find out, in which step the model (in form of a loop)
starts to produce NA's.
On Sat, 2009-05-09 at 20:50 +0200, Carlos J. Gil Bellosta wrote:
Hello,
I am building a package that creates a new kind of object not unlike a
dataframe. However, it is not an extension of a dataframe, as the data
themselves reside elsewhere. It only contains metadata.
I would like to be
On Sat, 2009-05-09 at 19:29 -0700, linakpl wrote:
If I was doing an ANOVA analysis how can I get the design matrix R used?
You can do ANOVA several different ways in R, and that is just the few
functions to do this that I am aware of. Showing us what function/code
you used would be really
Most common styles (e.g. APA, Harvard) include the date of access for an
electronic reference. While this may be an artifact of history, both
reviewers and editors are justified in asking authors to adhere to the
style used in a particular journal. That said, I don't see why Nature or
any
Hello,
Could someone please tell me how to find the estimate of inertia for the
second axis of an RLQ analysis using ade4? Using the example from the
ade4 package, I *suspect* that the inertia for the 2nd axis of the R
table would be 4.139332 (see below results summary from rlq1). However,
Searching the mail archives I found that using legend.position as in
p.ring.3 + opts(legend.position=top)
is a known bug. I tried doing
p.ring.3 + opts(legend.position=c(0.8, 0.2))
which works, but the legend background is transparent, i.e. I see the
plot background through the legend. Adding
On Sun, May 10, 2009 at 10:32 AM, Zeljko Vrba zv...@ifi.uio.no wrote:
Searching the mail archives I found that using legend.position as in
p.ring.3 + opts(legend.position=top)
is a known bug. I tried doing
p.ring.3 + opts(legend.position=c(0.8, 0.2))
which works, but the legend background
JC,
If each row are the counts for a 2 x 2 contingency table - so for the
ith contingency table you have counts for row 1 c(Y08[i],Z08[i]) and
row 2 (Y09[i],Z09[i]) then you could use apply:
X - cbind(vdata$Y08,vdata$X08-vdata$Y08,vdata$Y09,vdata$X09-vdata$Y98)
f.chisq - function(x){
m -
One simple adjustment is the following:
apply(matrix(x[rows,],nr=length(rows)),2,mean)
Quoting Peter Kharchenko peter.kharche...@post.harvard.edu:
Dear fellow R users,
I can't figure out how to do a simple thing properly: apply an
operation to matrix columns on a selected subset of rows.
hi,
I am searching for a nice function which computes correlations out of
a data.frame and adds asterix-signs after each correlation when they
are significant...
...or a function which just show correlations greater than x in the
output.
thanks!
Hi,
I'm not used to thinking along these lines, and wanted to ask your advice:
Suppose you have a sample of around 100, consisting of patients according to
doctors, in which patients and doctors are given a questionnaire with
categorical responses. Each patient somehow has roughly 3 doctors, or
Hello,
Could someone please tell me how to find the estimate of inertia for the
second axis of an RLQ analysis using ade4? Using the example from the
ade4 package, I *suspect* that the inertia for the 2nd axis of the R
table would be 4.139332 (see below results summary from rlq1). However,
Hi all R gurus out there,
Im a kind of newbie to kalman-filters after some research I have found that
the dlm package is the easiest to start with. So be patient if some of my
questions are too basic.
I would like to set up a beta estimation between an asset and a market index
using a
Quick question:
Which function do you use to calculate partial derivatives from a model
equation?
I've looked at deriv(), but think it gives derivatives, not partial
derivatives. Of course my equation isn't this simple, but as an example,
I'm looking for something that let's you control
Did you try the imgFFTConvolve function in the biOps package?
GreenBrower wrote:
Anybody know how to do a fast 2 dimension convolve in R? I used a
traditional method, but it's very slow, even can not stand it!
--
View this message in context:
On 5/10/2009 3:26 PM, Martin Batholdy wrote:
hi,
I am searching for a nice function which computes correlations out of a
data.frame and adds asterix-signs after each correlation when they are
significant...
...or a function which just show correlations greater than x in the output.
On Sat, 09 May 2009 12:17:57 -0400, Carl Witthoft c...@witthoft.com
wrote:
Sorry, but your professor offered me $500 NOT to do your assignments.
Actually all he wanted was your name so you could be expelled for
plagiarism!
--
Bob Doherty
__
Peter Dalgaard wrote:
Richard Chirgwin wrote:
Hi all,
I am trying to install Spatstat on OpenSUSE 11.1.
install.packages(spatstat, dependencies = TRUE)
fails on the basis of various compiler packages (full message below).
I have gcc version 4.3.2, which should include
If you want `numerical' partial derivatives, check out:
require(numDeriv)
?grad
Ravi.
Ravi Varadhan, Ph.D.
Assistant Professor,
Division of Geriatric Medicine and Gerontology
School of Medicine
Johns Hopkins University
On May 10, 2009, at 10:05 AM, Paul Heinrich Dietrich wrote:
Quick question:
Which function do you use to calculate partial derivatives from a
model
equation?
I've looked at deriv(), but think it gives derivatives, not partial
derivatives.
Your reading of the help page and the examples
Have you considered genD{numDeriv}?
If this does not answer your question, I suggest you try the
RSiteSearch package. The following will open a list of options in a
web browser, sorted by package most often found with your search term:
library(RSiteSearch)
pd -
Maybe pairs.panels(df,scale=T) from the psych library - se more here:
http://www.personality-project.org/r/html/00.psych-package.html
setting scale=T scales the cor coefficient according to their value. I
have seen an implementation with added asterix' but couldn't find it
right now.
On Sun, May
The dataset package is being loaded apparently by one of the packages that I
am using. The loading of the datasets takes a long time and I would like to
eliminate it. I thought the datasets were effectively examples so don't
understand why they would be required at all.
1) How can I determine
R 2.8.1
Windows XP
I am trying to plot the results of a coxph using plot(survfit()). The plot
should, I believe, show two lines one for survival in each of two treatment
(Drug) groups, however my plot shows only one line. What am I doing wrong?
My code is reproduced below, my figure is
Hi everyone! Thank you for the help you have been given to
me, and here I'm with another problem with my dataframes:
I have two dataframes (with much more observations), like
these:
Dataframe1
Firm Year cash
5004002002007 100
5004002002006 200
5004002002005
R 2.8.1
Windows XP
How does one plot the -log(log(survival)) from a coxph? Survfit does not seem
to be up to the task.
Thanks,
John
John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology
Baltimore VA Medical Center
10
I did not see an attached figure to compare against, perhaps it was removed
by the server. I suspect that:
plot(survfit(fit0,newdata=GVHDdata)
will produce the plot you're looking for. It will generate survival curves for
the data in GVHDdata using the estimates from the fit0 coxph model. At
On Sun, 10 May 2009, John Sorkin wrote:
R 2.8.1
Windows XP
I am trying to plot the results of a coxph using plot(survfit()). The plot
should, I believe, show two lines one for survival in each of two treatment
(Drug) groups, however my plot shows only one line. What am I doing wrong?
On Sun, 10 May 2009, John Sorkin wrote:
R 2.8.1
Windows XP
How does one plot the -log(log(survival)) from a coxph? Survfit does not seem
to be up to the task.
It depends on what you mean.
You can convert the plots to a -log(log(survival)) scale with the option
fun=cloglog. In the example
On 11/05/2009, at 9:17 AM, David A Vavra wrote:
The dataset package is being loaded apparently by one of the
packages that I
am using. The loading of the datasets takes a long time and I would
like to
eliminate it. I thought the datasets were effectively examples so
don't
understand why
?subset
?%in%
(I have gotten tired of converting dataframes that are presented in a
non-executable form, such as is supported by the dput function.
So, ... you should read those help pages and take the obvious path to
success.)
Something along the lines of:
subset(df1, Firm %in%
On May 10, 2009, at 7:33 PM, Rolf Turner wrote:
On 11/05/2009, at 9:17 AM, David A Vavra wrote:
The dataset package is being loaded apparently by one of the
packages that I
am using. The loading of the datasets takes a long time and I would
like to
eliminate it. I thought the datasets
Paul,
I suggest looking up observer agreement. The description of your study
sounds like a classical
categorical observer agreement problem. I can't give
a reference off the top of my head, but if you get
stuck, e-mail me and I'll try and find a ref to get you started.
Murray M Cooper, Ph.D.
Thanks. Perhaps something else is going on. There is a large time period
(about 20 sec.) after the message about loading the package. More
investigation, I suppose.
Thanks again,
DAV
-Original Message-
From: Rolf Turner [mailto:r.tur...@auckland.ac.nz]
Sent: Sunday, May 10, 2009 7:34
At the moment I'm just reading the large file to see how fast it goes.
Eventually, if I can get the read time down, I'll write out a processed
version. Thanks for suggesting scan(); I'll try it.
Rob
jim holtman wrote:
Since you are reading it in chunks, I assume that you are writing out each
G'day David,
On Sun, 10 May 2009 17:17:38 -0400
David A Vavra dava...@verizon.net wrote:
The dataset package is being loaded apparently by one of the packages
that I am using. The loading of the datasets takes a long time and I
would like to eliminate it. I thought the datasets were
G'day David,
On Sun, 10 May 2009 21:35:30 -0400
David A Vavra dava...@verizon.net wrote:
Thanks. Perhaps something else is going on. There is a large time
period (about 20 sec.) after the message about loading the package.
More investigation, I suppose.
What R version are you using? I do
Hi,
I'm trying to build some maps of the US by county that will have the
following characteristics:
Feature/Map
Map 1
Map2
Both
Broken out by county
Yes
Yes
Yes
Heatmaps of US Census Data for income by county
Yes
No
Yes
Thank you for suggesting other functions, I will look into them.
When I read the deriv() function, it did mention partial, but I (being a
newbie) wasn't able to get partials for a simple MNL equation. I'm sure I
did something wrong then, but here's what I tried the following and got
different
To look up observer agreement, you might consider the
RSiteSearch package. The RSiteSearch.function looks only for
matches in help pages of contributed packages.
library(RSiteSearch)
oa - RSiteSearch.function(observer agreement)
attr(oa, hits) # 4 functions matching this term
HTML(oa)
53 matches
Mail list logo