On Tue, 22 Nov 2005, S. Sun wrote:
I have a question about the glm.
Not really: your question is about understanding logistic regressions.
When the events of an observation is 0,
the logit function on it is Inf. I wonder how the glm solve it.
Note that logit(0) = -Inf whereas logit(1) =
On Tue, 22 Nov 2005, Stephen wrote:
I have been asked to provide Weibull parameters from a paper using
Kaplan Meir survival analysis.
This is something I am not familiar with.
The survival analysis in R works nicely and is the same as commercial
software (only the graphs are superior in
Hi useRs!
I am in search of any R programer / user to go on with the package R2HTML.
I did change of job 3 months ago and I could only realize I have
absolutely no time anymore for that (even so little) task.
R2HTML package consists in a set of HTML.* functions that could be
used to export R
On Fri, 18 Nov 2005, Ulrike Grömping wrote:
From: Duncan Temple Lang [EMAIL PROTECTED]
[...]
For the record, the problem you are experiencing is that
you are not using the same argument names in your
method as the function $ uses. Hence there is a mismatch.
Use
setMethod($,
hi, hi all,
you can consult this link :
http://tolstoy.newcastle.edu.au/R/help/05/11/15455.html
hope this help,
cheers,
P.BADY
Selon David C. James [EMAIL PROTECTED]:
Dear R-Help,
I'm hoping to find a Multinomial Nested Logit package in R. It would
be great to find something
Dear Claus,
I had not exactly the same pb, but close enough to propose you to take a
look at the options fin and fig in the par() function.
Hope this helps.
Florence.
On 11/19/05, Claus Atzenbeck [EMAIL PROTECTED] wrote:
Hi,
I produce a series of diagrams with R in order to include them
Dear All,
I am trying to apply the function 'cox.zph' of the library survival, but I
receive this error message:
not found the object 'residuals.coxph'.
I have re-installed the library 'survival' without any change and also a search
with RSiteSearch was unsuccessful..
Any suggestion?
TIA
Thanks
Regards
Stephen
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Tue 22/11/2005 10:15
To: Stephen
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Weibull and survival
On Tue, 22 Nov 2005, Stephen wrote:
I have been asked to provide Weibull
Hi there again,
I have a problem with the parse-command.
First of all, I show you in a simplified way, what I am trying to do and what
R answers:
test
[1] u.g$par1, u.g$par2
mode(test)
[1] character
ausdruck - parse(text = test)
Error in parse(file, n, text, prompt) :
On Tue, 22 Nov 2005, Florence Combes wrote:
I had not exactly the same pb, but close enough to propose you to take a
look at the options fin and fig in the par() function.
Thanks for your answer. In the meanwhile I also found another solution.
I produce diagrams as usual and use pdfcrop to
Hi
I am conducting 2-sample Kolmogorov Smirnov tests for my Masters project to
determine if two independant tree populations have the same size-class
distribution or not. The trees have been placed into size-class categories
based on their basal diameters. Once I started running the stats on my
Another question on cointegration...
Is it possible to insert in the model dummy variables restricted in the
cointegration space?
Many thanks,
Carlo
On Nov 21, 2005 01:23 PM, Pfaff, Bernhard Dr.
[EMAIL PROTECTED] wrote:
Thanks a lot!
I have another question on cointegration, so I will go
Hello Antje,
Is there anything I have overlooked? Can't I put special characters
like $ or , into this test-character?
If I understand ?parse correctly, you are supposed to pass a valid R
expression to it and
u.g$par1, u.g$par2
is not as the comma is not a valid delimiter here.
Lisa Dowson wrote:
Hi
I am conducting 2-sample Kolmogorov Smirnov tests for my Masters project to
determine if two independant tree populations have the same size-class
distribution or not. The trees have been placed into size-class categories
based on their basal diameters. Once I started
Dear all,
I am looking for a way to identify the row and column names of all elements
within a matrix which fulfill some specified condition.
Example: I have a correlation matrix of the following form:
Y1Y2
X10.4 0.3
X20.6 0.1
Suppose, I want to know which elements
?which
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Petr Pikal wrote:
Hi
did you tried something like
write.table( tab, file.txt, sep=\t, row.names=F)
which writes to tab separated file?
Petr thanks, but I do not want a tab delimited file. I need spaces
between columns.
write.table( tab, file.txt, sep=, row.names=F)
Can it do what you
On Tue, 22 Nov 2005, giovanni parrinello wrote:
I am trying to apply the function 'cox.zph' of the library survival, but
I receive this error message: not found the object 'residuals.coxph'.
Please show us exactly what you did and use traceback() to show where the
message came from.
On 22-Nov-05 [EMAIL PROTECTED] wrote:
Dear all,
I am looking for a way to identify the row and column names of all
elements
within a matrix which fulfill some specified condition.
Example: I have a correlation matrix of the following form:
Y1Y2
X10.4 0.3
X20.6
hi all
i would like to know if anyone has a reference on how one would place
the bands on the pp plot.
i want to test whether or not a certain data set comes from a particular
distribution (not normal).
i've already plotted F(X(j)) vs j/(n+1) where F(x) is the cum dist
function, X(j) is the
Hi
can you please explain why do you need it? What do you want to do with
the exported file?
I wonder what type of software can not accept any reasonable delimiter
and requires fwf files.
The only workaround I can imagine is to transfer all columns to
character and
Dear all,
I am using the package spatstat, for analysis of inter and intra species
spatial patterns.
PPPData is may ppp object, and I am using formulas like:
Js- Jcross(PPPData,Aun,Qfa)
plot(Js, sqrt(./pi)~r, ylab=L(r), main=L function Aun,Qfa)
Js- Jest(PPPData)
plot(Js, sqrt(./pi)~r,
Is qq.plot in package 'car' of use to you? I think that it
requires your distribution to be one of those available in R.
Peter
Clark Allan wrote:
hi all
i would like to know if anyone has a reference on how one would place
the bands on the pp plot.
i want to test whether or not a certain
Hi
can you please explain why do you need it? What do you want to do with the
exported file?
I wonder what type of software can not accept any reasonable delimiter and
requires fwf files.
The only workaround I can imagine is to transfer all columns to character and
add leading spaces
I sent out a question a few days ago asking help with cv.glm and I
didn't get any responses. Is it my question? I'd appreciate any
response just to see my post makes it out there.
What I'm looking to do is to get a column of predicted responses from
cv.glm (leave-one-out). This produces a
Dear list,
I have a question concerning the above mentioned methods in the pls
package with respect to the loadings matrix produced by the call. In
some work I am doing I have found that the values produced are nearly of
the same magnitude but of opposite sign. When I use the example data
the distribution is not a standard one
P Ehlers wrote:
Is qq.plot in package 'car' of use to you? I think that it
requires your distribution to be one of those available in R.
Peter
Clark Allan wrote:
hi all
i would like to know if anyone has a reference on how one would place
On 11/22/05, Gorjanc Gregor [EMAIL PROTECTED] wrote:
Hi
can you please explain why do you need it? What do you want to do with the
exported file?
I wonder what type of software can not accept any reasonable delimiter and
requires fwf files.
The only workaround I can imagine is
Hello Carlo,
no, this is not possible, per se. In case of a structural break in terms of
a one-time level shift, you might want to consider the function cajolst().
Another possibility would be to run a regression with the dummy-variables
and use the fitted values for your data matrix x. That is,
On Tue, 2005-11-22 at 10:48 +0100, Claus Atzenbeck wrote:
On Tue, 22 Nov 2005, Florence Combes wrote:
I had not exactly the same pb, but close enough to propose you to take a
look at the options fin and fig in the par() function.
Thanks for your answer. In the meanwhile I also found
On Tue, 22 Nov 2005, Marc Schwartz wrote:
Without your code, it is hard to know exactly what you need to resolve
the problem or why your manipulation of the margins did not work as you
want, but I am guessing that the following might be helpful.
The reason why the manipulation of margins does
I tried to get old libraries (libreadline4) installed without success.
If the configure script checks for old versions of the readline library can't
we just comment the corresponding lines?
This check is in there to check for some really old readline versions
that didn't allow the style of
Why don't you just use 'sprintf' and set the format up to what you want:
x - head(iris)
cat(sprintf(%6.1f%6.1f%6.1f%6.1f %-20s\n, x[,1], x[,2], x[,3], x[,4],
x[,5]), sep='')
5.1 3.5 1.4 0.2 setosa
4.9 3.0 1.4 0.2 setosa
4.7 3.2 1.3 0.2 setosa
4.6 3.1 1.5 0.2 setosa
5.0 3.6 1.4 0.2 setosa
5.4
Clark Allan [EMAIL PROTECTED] writes:
the distribution is not a standard one
As long as it has a CDF and a density it doesn't seem to matter if you
just follow the naming conventions:
qfoo - qnorm
dfoo - dnorm
z - rnorm(50)
qq.plot(z,foo)
The code for qq.plot.default appears quite
Hi
Excellent diagnosis we have zero survival times
Unfortunately no MASS4 p. 380 - I can't locate - help appreciated.
Thanks
S.
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Tue 22/11/2005 10:15
To: Stephen
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R]
Pedro Cordeiro Estrela [EMAIL PROTECTED] writes:
I tried to get old libraries (libreadline4) installed without success.
If the configure script checks for old versions of the readline library
can't we just comment the corresponding lines?
No, because the code in question still needs to
Having finally updated from R 1.91 to R 2.2.0 with my installation of a
new computer, I discovered that something has changed drastically about
the way code for contributed packages is stored when installed in a local
version of R. In the 1.* versions it was easy for me to go in and modify
On Tue, 22 Nov 2005, giovanni parrinello wrote:
Dear All,
I am trying to apply the function 'cox.zph' of the library survival, but I
receive this error message:
not found the object 'residuals.coxph'.
I have re-installed the library 'survival' without any change and also a
search with
Re: Uwe's advice.
A quote from the late Ellis Ott (a U.S. quality control statistician of the
1970's and 80's):
Look at your data -- and think!
Cheers,
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific
Hello,
I need some help in finding the 'right'(maybe better 'the best') solution for
my cluster analysis calculations. There are so many ways to do cluster
analysis, I'm a litte bit puzzled. I've read one should choose the solution
after what makes the best sense in regard to the content of
On 22 Nov 2005, at 17:35, Brian S Cade wrote:
Having finally updated from R 1.91 to R 2.2.0 with my installation of a
new computer, I discovered that something has changed drastically about
the way code for contributed packages is stored when installed in a
local
version of R. In the 1.*
Dear R users,
My apologies for a simple question for which I suspect there
is a simple answer that I have yet to find. I'd like to plot
panels in lattice with different graphical parameters for the
axes. For example, the code
x-rnorm(100)
y-rnorm(100)
z-c(rep(1,50), rep(2,50))
Brian S Cade wrote:
Having finally updated from R 1.91 to R 2.2.0 with my installation of a
new computer, I discovered that something has changed drastically about
the way code for contributed packages is stored when installed in a local
version of R. In the 1.* versions it was easy for me
On Tue, 2005-11-22 at 08:35 -0700, Brian S Cade wrote:
Having finally updated from R 1.91 to R 2.2.0 with my installation of a
new computer, I discovered that something has changed drastically about
the way code for contributed packages is stored when installed in a local
version of R. In
dear deepayan all...
let me explain you my full problem clearly..so that you could clearly suggest
me in right direction
my matrix (a) looks like as below...I´ve done clusutering in this matrix
I made it through the first chapter (Background) of Bayesian Data
Analysis by Gelman, et al. Conceptually, the Bayesian approach
appeals to me and my curiosity is piqued. However, the next chapter
was much too terse. The math is daunting. Where can I find a gentle
introduction? Or which math
On 22 Nov 2005, [EMAIL PROTECTED] wrote:
Due to lazy loading (see a corresponding R News article) code is
stored in databases for some time now, if you are looking into
binary packages.
In order to modify/fix something, download the source version of a
package, change the code, and INSTALL
On Tue, 22 Nov 2005, Pedro Cordeiro Estrela wrote:
I tried to get old libraries (libreadline4) installed without success.
Why do you need that? The configure test works correctly with the current
readline 5.0, which also needs ncurses.
What is `R.2.2'? There is no such thing, and your line
Please read ?princomp to understand more about loadings and their
arbitrary signs.
On Tue, 22 Nov 2005, Roy Little wrote:
Dear list,
I have a question concerning the above mentioned methods in the pls
package with respect to the loadings matrix produced by the call. In
some work I am doing
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to
announce 2-day R/S-plus Fundamentals and Programming
Techniques in New York: www.xlsolutions-corp.com/Rfund.htm
New York, December 19-20, 2005
Reserve your seat now at the early bird rates! Payment due AFTER
the class
Course
It's hard for me to resist dipping my oar into this...
Variable labels are so generally useful, both in documenting a
dataset (what was 'pctx723' again?) and in producing readable
output and graphs that it is a shame they are not provided in
base R. If they were (and were used in print and plot
Does that mean, then, that the help page's x$name just accidentally coincides
with the arguments x and name that the $-function requires in S4?
Regards, Ulrike
-- Original Message ---
From: Prof Brian Ripley [EMAIL PROTECTED]
To: Ulrike Grömping [EMAIL PROTECTED]
Cc: Duncan
On Tue, 2005-11-22 at 08:53 -0400, Roy Little wrote:
Dear list,
I have a question concerning the above mentioned methods in the pls
package with respect to the loadings matrix produced by the call. In
some work I am doing I have found that the values produced are nearly of
the same
On 11/22/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Dear R users,
My apologies for a simple question for which I suspect there
is a simple answer that I have yet to find. I'd like to plot
panels in lattice with different graphical parameters for the
axes. For example, the code
Hello all,
Does anyone know how to make R graph the x-axis in descending order? I am all
set in creating a vertical bar chart using either the plot or the xychart (in
lattice), but I need the x-axis in descending order. I have tried resorting
the data but apparently R automatically graphs the
On 11/22/05, peter eric [EMAIL PROTECTED] wrote:
dear deepayan all...
let me explain you my full problem clearly..so that you could clearly
suggest me in right direction
my matrix (a) looks like as below...I´ve done clusutering in this matrix
a-c(1:10)
plot(a,ylim=c(10,1))
On Tue, 2005-22-11 at 15:03 -0500, Anthony Pitts wrote:
Hello all,
Does anyone know how to make R graph the x-axis in descending order? I am all
set in creating a vertical bar chart using either the plot or the xychart (in
lattice), but I need the x-axis in
I've read through many postings about principle component analysis in
the R-help archives, but haven't been able to piece together the
information I need. I'd like to recreate an SPSS-like experience of
factor analysis using R. Here's what SPSS produces:
1. Scatterplots of all possible
I am just starting with R and have the following problem: given a matrix
of ones and zeroes, say
mdim=4
m-matrix(round(runif(mdim^mdim)),mdim,mdim)
how to construct the matrix of those rows of m, whose elements sum to 2.
Contrary to the random matrix above, the actual matrix always has at
least
I am running a VERY LARGE regression (many factors, many rows of data) using LM.
I think I have my memory set as high as possible. ( I ran memory.limit(size
= 4000))
Is there anything I can do? ( FYI, I think I have removed all data I am
not using, and I think I have only the data
On Tue, 2005-11-22 at 21:55 +0100, Serguei Kaniovski wrote:
mdim=4
m-matrix(round(runif(mdim^mdim)),mdim,mdim)
One approach:
m
[,1] [,2] [,3] [,4]
[1,]0010
[2,]0110
[3,]0011
[4,]0000
idx - apply(m, 1, sum)
m[idx,]
I've made a package for which R CMD build isn't producing very
satisfactory results. I'll get to the details in a moment.
I wonder if it would make sense to have my own makefiles (which already
exist and are doing quite a lot) produce the .tar.gz file ordinarily
produced by R CMD build. As far
On 22 Nov 2005, [EMAIL PROTECTED] wrote:
I am just starting with R and have the following problem: given a
matrix of ones and zeroes, say
mdim=4
m-matrix(round(runif(mdim^mdim)),mdim,mdim)
how to construct the matrix of those rows of m, whose elements sum
to 2. Contrary to the random
Standard caution: == should not be used for comparison of floating point
numbers. See ?all.equal and ?identical. See also the Green Book (Chambers:
PROGRAMMING WITH DATA). Sometimes what you think are integers may be floats,
too, as numerous postings to this list have shown.
-- Bert Gunter
___
There is also rowSums (and colSums) which does the same thing without
needing apply:
m[rowSums(m) == 2, ]
Subject to Bert Gunter's caveat about comparing floating point numbers
for equality, of course.
--
Hello,
I'm a grad student in the Intelligent Transportation Systems lab at Portland
State Univ. in Portland, OR, USA. I'm trying to learn the basics of R to run a
negative binomial in the near future, and so I ran a test regression on roadway
crash data obtained from Statistical and Econometric
___
I don't know the text in question, but it looks like a difference in
parametrisation -- the R output is just the inverse of that from the
book.
--
Hong Ooi
Senior Research Analyst, IAG Limited
388 George
Dear Joshua,
You might take a look at Lancaster, An Introduction to Modern Bayesian
Econometrics.
I hope this helps,
John
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario
Canada L8S 4M4
905-525-9140x23604
http://socserv.mcmaster.ca/jfox
Dear Anthony,
You can use xlim in plot() to reverse the x-axis. For example,
x - 1:10
y - x
plot(x, y, xlim=rev(range(x)))
(I'm pretty sure, by the way, that this question has been asked before.)
I hope this helps,
John
John Fox
Department of Sociology
Dear List
I am new to R and to the list and will try best as I can be clear and
concise. My apologies if anything I write contravenes the posting code
on this list. I would also like to say I have run through most of the
material on the R website before writing this email however, I am
stuck.
Not sure about reading the first few lines, but summary(x) might be what you're
after.
Try using the menu to save the workspace. File Save Workspace. You can then
load the workspace at the start of your next session and continue working.
Murray
-Original Message-
From: Ronnie
Hi,
I'm running a repeated measure ANOVA to test the effects of the
within-subjects
factors congr, isi and their interaction, on the variable latencies
aov(latencies ~ congr*isi + Error(subj/(congr*isi)),data=dats)
congr has 4 levels, while isi has 3 levels.
I would like to check for possible
I am currently using the coxph function from R on a microarray dataset.
The curve fit is performed using each gene's expression measurements
separately. I then get a p-value for each gene, from the Coxph.summary
method. Can someone please explain how this is computed and the
significance
Dear R-hep,
I'm trying to get the first derivative of a smooth from a gam
model like:
model-gam(y~s(x,bs=cr, k=5)+z) and need the derivative: ds(x)/dx. Since
coef(model) give me all the parameters, including the parameters of the
basis, I just need the 1st derivative of the basis s(x).1,
Ronnie, try
head?
tail?
and, I also fine useful after loading
str?
On 11/22/05, Ronnie Babigumira [EMAIL PROTECTED] wrote:
Dear List
I am new to R and to the list and will try best as I can be clear and
concise. My apologies if anything I write contravenes the posting code
on this list. I
If you are desperate, save the file as tab delimited and then use vi or
some stream editor to convert the tabs to spaces. If you set the tab
stop wide enough you should be able to guarantee uniform columns.
Gorjanc Gregor wrote:
Petr Pikal wrote:
Hi
did you tried something like
write.table(
76 matches
Mail list logo