On Tue, 22 Jun 2004, B. Dan Wood wrote:
Is there a package or procedure that implements restricted estimation using
glm? I need to do a poisson regression under glm with one parameter fixed at
a specific value. Just at a glance, I don't see a way to do it using glm.
Use offset(), which allows
You can a logistic non-linear regresion with a continuous response -- see
SSlogis in package stats (using nls).
However, I think
I believe I should run a logistic regression on this
is the problem. You have six groups of multivariate observations and you
want to know if and how they differ
Paul Roebuck wrote:
I have the following contrived code in package format.
On Solaris and Mac OS X, code runs just fine. On Windows,
it crashes the R environment with the Send Bug Report
dialog. I tried R 1.8.1 (Win2K) and R 1.9 (WinXP) binaries
with the same result. PCs otherwise appear properly
Does anybody know of any R functions to perform chronological clustering
as explained in:
Legendre, P., S. Dallot L. Legendre. 1985. Succession of species
within a community: chronological clustering, with applications to
marine and freshwater zooplankton. American Naturalist 125: 257-288.
Dear R People:
Is there a way to forecast with GARCH modeling as found in tseries, please?
When I use the predict command, I get an output of length 100, regardless of
what I put in the n.ahead steps.
R Version 1.9.0
Thanks in advance.
Sincerely,
Laura
mailto: [EMAIL PROTECTED]
download!
Dear Sir/Madam,
I encounter some problem duuring compilation of R 1.9.1 on AIX 5.1, after
running ./configure then I type make to compile:
# make
.
gcc -I../../src/extra/zlib -I../../src/extra/bzip2 -I../../src/extra/pcre
-I. -I../../src/include -I../../src/include -I/usr/local/include
We do have alpha/beta test periods for new releases of R, so would anyone
who did test AIX please confirm that they did succeed. (No one reported
an error, including yourself: are any AIX users interested in helping
having R available for AIX?) It is much better to have such reports
during
Monica Palaseanu-Lovejoy wrote:
Hi,
I would like to write a little automatic routine in R, but i am a too
much of a beginner for that. I will appreciate any help regarding this
particular problem.
Lets suppose I have a data.frame with j columns (from 1 to n) and i
rows (from 1 to p). I would
How about
newDF - DF
newDF[] - lapply(DF, function(x) pmax(0, pmin(x, 1)))
as in
DF - data.frame(x=rnorm(5), y=rnorm(5))
newDF - DF
newDF[] - lapply(DF, function(x) pmax(0, pmin(x, 1)))
newDF
x y
1 0.31872426 1.000
2 1. 0.000
3 1. 0.4510969
4
DeaR UseRs:
I want to put a legend in my plot. In the first line of the legend I want to
put a box filled but in the second one I would like to put a lty=2
Of course it must appear with different colors.
Thanks.
__
[EMAIL PROTECTED] mailing list
Hi,
Can someone tell me how to put labels of differente colors on a tree
dendrogram.
Thanks
--
Philippe Hupé
UMR 144 - Service Bioinformatique
Institut Curie
Laboratoire de Transfert (4ème étage)
26 rue d'Ulm
75005 Paris - France
Email : [EMAIL PROTECTED]
Tél :+33 (0)1 44 32 42 75
Hi,
- Original Message -
From: Perez Martin, Agustin [EMAIL PROTECTED]
To: lista R help (E-mail) [EMAIL PROTECTED]
Sent: Wednesday, June 23, 2004 9:19 PM
Subject: [R] legend
DeaR UseRs:
I want to put a legend in my plot. In the first line of the legend I want
to
put a box filled
Hi,
Thank you so much for the little example - this was exactly what i
was looking for. Just i didn't know what to search after in R. I will
change it for my data needs.
Thanks again,
Monica
__
[EMAIL PROTECTED] mailing list
Hi,
This works beautifully as well.
Thanks,
Monica
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Hi,
I am using optim with a L-BFGS-B method to minimize a function. As I've
understood, the way to specify a tolerance for stopping optimization is
through factr argument.
My function, is by construction, minimal when equal to 1. I wonder if there
is any way to pass this info to optim. If not,
Uwe Ligges wrote:
Monica Palaseanu-Lovejoy wrote:
Hi,
I would like to write a little automatic routine in R, but i am a too
much of a beginner for that. I will appreciate any help regarding this
particular problem.
If all columns of your data.frame are numeric:
z[z0] - 0
z[z1] - 1
For added
On Wed, 23 Jun 2004, Barry Rowlingson wrote:
Uwe Ligges wrote:
Monica Palaseanu-Lovejoy wrote:
Hi,
I would like to write a little automatic routine in R, but i am a too
much of a beginner for that. I will appreciate any help regarding this
particular problem.
If all columns
Hi -
platform powerpc-apple-darwin6.8
status
major1
minor9.0
year 2004
I am trying to deal with the output of apply(). As indicated, when
each call to 'FUN' returns a vector of length 'n', then 'apply'
returns an array of dimension 'c(n, dim(X)[MARGIN])'. However, I would
like
On Tuesday 22 June 2004 17:41, Prof Brian Ripley wrote:
data.frame(time=mytime, ts=myts)
would appear to be what you are looking for.
On Tue, 22 Jun 2004 [EMAIL PROTECTED] wrote:
I have defined a daily timeseries for the 365 days of 2003 issuing:
myts = ts(dati[,2:10],frequency=365,)
Dear R group
I am making multiple xyplot, and would like to have tick marks on bottom and left in
EACH panel, but only tick labels at the bottom and left of the whole graph.
I have browsed the internet, as well as the help page without success.
If anyone could help me find the path to the
Hi,
Why not try with the data.frame structure, wich internally yet consists in
a list:
lapply(as.data.frame(m1),function(x) which(x==1))
$V1
[1] 8 9
$V2
[1] 9 10
[...]
Eric
At 12:53 23/06/2004, Andrew Beckerman wrote:
Hi -
platform powerpc-apple-darwin6.8
status
major1
minor9.0
year
Hi,
anybody know if Rcmdr or any other GUI like this work with R/Aqua?
Thanks
Ronaldo
--
Let's organize this thing and take all the fun out of it.
--
| // | \\ [***]
| ( õ õ ) [Ronaldo Reis Júnior]
| V [UFV/DBA-Entomologia
This may have been answered before, but a nice possibility to plot the
Titanic data is offered by mosaic plots, as introduced (I think) by
John Hartigan and implemented in the VCD package
Ivailo Partchev
Jena
Marc Schwartz wrote:
On Tue, 2004-06-22 at 10:54, Patrick Lenon wrote:
Good day
Anybody use Sciviews with Linux or MacOS X? It work with these systems?
Thanks
Ronaldo
--
A hacker does for love what others would not do for money.
--
| // | \\ [***]
| ( õ õ ) [Ronaldo Reis Júnior]
| V [UFV/DBA-Entomologia
On Wednesday 23 June 2004 06:36, Ingolfsson, Olafur wrote:
Dear R group
I am making multiple xyplot, and would like to have tick marks on
bottom and left in EACH panel, but only tick labels at the bottom and
left of the whole graph. I have browsed the internet, as well as the
help page
Hi Again,
First of all thank you for all the responses to my previous query.
Your answers were very helpful and I did the job ;-). Now I hope you
can answer as quick the following (sorry I am invading you with
trivial questions):
Lets use again the following data.frame example:
DF -
Ronaldo Reis Jr wrote:
Anybody use Sciviews with Linux or MacOS X? It work with these systems?
Not yet... but we are rewritting a part of the code so that some features
become platform-independent. Please, leave us a little time to archive this.
There are tens of thousands of lines of code in
Uwe Ligges wrote:
Paul Roebuck wrote:
I have the following contrived code in package format.
On Solaris and Mac OS X, code runs just fine. On Windows,
it crashes the R environment with the Send Bug Report
dialog. I tried R 1.8.1 (Win2K) and R 1.9 (WinXP) binaries
with the same result. PCs
Uwe Ligges wrote:
Paul Roebuck wrote:
I have the following contrived code in package format.
On Solaris and Mac OS X, code runs just fine. On Windows,
it crashes the R environment with the Send Bug Report
dialog. I tried R 1.8.1 (Win2K) and R 1.9 (WinXP) binaries
with the same result. PCs
Sounds to me that you are looking for something like:
subDF - DF[3:5,]
Andy
From: Monica Palaseanu-Lovejoy
Hi Again,
First of all thank you for all the responses to my previous query.
Your answers were very helpful and I did the job ;-). Now I hope you
can answer as quick the
Perez Martin, Agustin wrote:
DeaR UseRs:
I want to put a legend in my plot. In the first line of the legend I want
to put a box filled but in the second one I would like to put a lty=2
Of course it must appear with different colors.
I think I wrote this function about a year ago for
The documentation is good for things like this. In most (all recent?)
versions of R, help.start() brings up a help page. Select An
Introduction to R, and go to Vector Indices. Also look at Arrays and
Matrices: Array Indexing. There, you will find that DF[3:5,] might be
what you want.
Hi Spencer,
Your answer is very helpful. I was wondering if i should write again
to the list to ask where i suppose to read about indices and things
liek that since it seems they are very useful in lots of things. But
you already gave me the answer.
Thanks again,
Monica
On Wed, 2004-06-23 at 13:53, Monica Palaseanu-Lovejoy wrote:
Hi Again,
First of all thank you for all the responses to my previous query.
Your answers were very helpful and I did the job ;-). Now I hope you
can answer as quick the following (sorry I am invading you with
trivial
Dear All!
How can I cope with overlapping or covered labels (covered by labels
from other data points) in plots?
Martina Renninger
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED] mailing list
On Wed, 2004-06-23 at 09:06, Martina Renninger wrote:
Dear All!
How can I cope with overlapping or covered labels (covered by labels
from other data points) in plots?
Presuming that you are using text() to identify points in a plot, you
can use the 'cex' argument (which defaults to 1) to
Dear R community
When plotting an its object (with 'plot(..)), grid lines are
automatically added which I cannot control with parameter 'tck/tcl'.
Furthermore the vertical grid lines (time axis) do not comply with the
automatically generated tickmark locations.
Can somebody tell me how to
From: Prof Brian Ripley
We do have alpha/beta test periods for new releases of R, so
would anyone
who did test AIX please confirm that they did succeed. (No
one reported
an error, including yourself: are any AIX users interested in helping
having R available for AIX?) It is much
Hello,
Yes, indeed. All you have to do it's: launch your X server (take a look
at your /Applications/Utilities folder), and then launch R from an xterm.
If you don't have your X server installed, take a look at the RAqua-FAQ
and the R-FAQ on CRAN.
If you need more help, please don't hesitate
R users:
My R is 1.8.1 in Linux. How can I call R in Perl process? And call Perl from R?
Thanks
Xiao
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hi,
I hope someone has seen this problem before. I am
trying to install the XML package from
http://cran.r-project.org/src/contrib/XML_0.95-6.tar.gz.
Using
R CMD INSTALL XML_0.95-6.tar.gz
I get the following error message:
/usr/vac/bin/xlc_r -I/usr/local/lib/R/include -DLIBXML
Hallo!
I want to fit a function. The function is e.g.:
y = c+m1*x if x0, c+m2*x if x=0
where m1, m2 and c is a parameter and x, y are
variables of a data frame.
I think using nls is appropriate. But I do not know,
how to type this formula in nls. Can anybody help?
(If there is a possibility to
You already have a few solutions but here is one more:
(z = 0 z = 1) * z + (z 1)
It uses the fact that TRUE and FALSE act as 1 and 0 in arithmetic.
Slightly shorter but trickier is:
(z*z = z) * z + (z 1)
which uses the fact that z*z = z only in the interval [0,1] .
Monica
Hi,
look at http://www.omegahat.org/RSPerl/index.html.
regards,
Arne
--
Arne Muller, Ph.D.
Toxicogenomics, Aventis Pharma
arne dot muller domain=aventis com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of XIAO LIU
Sent: 23
On Wed, 23 Jun 2004, Karl Knoblick wrote:
I want to fit a function. The function is e.g.:
y = c+m1*x if x0, c+m2*x if x=0
where m1, m2 and c is a parameter and x, y are
variables of a data frame.
I think using nls is appropriate. But I do not know,
how to type this formula in nls. Can
I am trying to batch the printing of some graphs using win.print when I leave the
printer option blank it prompts me for a printer each time and then works fine when I
put the printer name in however using the following code I get the error below. Does
the printer need to be locally installed
Do remember that you need to double backslashes in R strings ... see the
rw-FAQ. As I recall forward slashes work here too.
On Wed, 23 Jun 2004, Joy, Tom wrote:
I am trying to batch the printing of some graphs using win.print when I leave the
printer option blank it prompts me for a printer
Hi,
I was wondering if there is a way of plotting smoothed contours of PCA
loadings onto a map. At the moment I have the loadings plotted as
positive/negative (red/blue) scaled circles centred around where the
physical stations are on a map - but I was wondering if there is a nicer
way of
Prof Brian Ripley wrote:
Do remember that you need to double backslashes in R strings ... see the
rw-FAQ. As I recall forward slashes work here too.
On Wed, 23 Jun 2004, Joy, Tom wrote:
I am trying to batch the printing of some graphs using win.print when I leave the
printer option blank it
On Wed, 23 Jun 2004, Uwe Ligges wrote:
Paul Roebuck wrote:
I have the following contrived code in package format.
On Solaris and Mac OS X, code runs just fine. On Windows,
it crashes the R environment with the Send Bug Report
dialog. I tried R 1.8.1 (Win2K) and R 1.9 (WinXP) binaries
On Wed, 23 Jun 2004, XIAO LIU wrote:
My R is 1.8.1 in Linux. How can I call R in Perl process?
And call Perl from R?
This is only one direction but worth a look as an
alternative.
http://search.cpan.org/~gmpassos/Statistics-R-0.01/
--
Thanks eric... I figured this routine out as well.
cvt-function(dat){
x-as.list(rep(0,dim(dat)[2]))
for(i in 1:dim(dat)[2]){
x[[i]]-dat[,i]
x}}
# get ragged array of 1's
dat-apply(mat,2,function(x) which(x==1))
# deal with this using cvt to creat list
if(is.null(dim(dat))) dat2-dat else
Hi,
I have data like this
print(x)
ID VAL1VAL2
1 2 6
2 4 9
3 45 12
4 99 44
What I would like is data like this...
ID VAL1VAL2
1 2 6
2 4 9
3 12 45
4 44 99
So that my analysis of
Do you mean:
pmax(x[,1], x[,2]) / pmin(x[,1], x[,2])
[1] 3.00 2.25 3.75 2.25
??
Andy
From: Dan Bolser
Hi,
I have data like this
print(x)
IDVAL1VAL2
1 2 6
2 4 9
3 45 12
4 99 44
What I would like is data like this...
ID
On Wed, 23 Jun 2004, Liaw, Andy wrote:
Do you mean:
pmax(x[,1], x[,2]) / pmin(x[,1], x[,2])
Magic!
I needed a p!
Andy
From: Dan Bolser
Hi,
I have data like this
print(x)
ID VAL1VAL2
12 6
24 9
345 12
499 44
What I
I know that if I have a function that returns multiple values, I should
do return(list(foo, bar)). But what do I do on the recieving end?
fn - function(x) {
return(list(foo, bar))
}
I know that at this point I could say
values.list - fn(x)
and then access
values.list[1]
values.list[2]
But that's
On Wed, 23 Jun 2004, Jack Tanner wrote:
I know that if I have a function that returns multiple values, I should
do return(list(foo, bar)). But what do I do on the recieving end?
fn - function(x) {
return(list(foo, bar))
}
I know that at this point I could say
values.list - fn(x)
Hi R-users
Does anybody know if there are functions or packages in R for sphericity
tests like Bartlett's or Mauchly's?
Thanks!
__
Eng. Agr., M.Sc. Eduardo Dutra de Armas
__
Centro de
Here are two approaches assuming foo is zz and bar is 3.
FIRST
You could pass the return variables in the argument list and then
assign them in the caller's frame like this:
fn - function(x,y) {
assign(as.character(substitute(x)), zz, sys.frame(-1))
My $0.02:
If Jack feels that the fact that functions usually return a single list, and
one needs to access the components of the list separately, is somehow
hideous, then I'd rather suggest that R is perhaps the wrong language for
him.
To me the suggested `workarounds' are by far much more
Paul Roebuck wrote:
On Wed, 23 Jun 2004, Uwe Ligges wrote:
Paul Roebuck wrote:
I have the following contrived code in package format.
On Solaris and Mac OS X, code runs just fine. On Windows,
it crashes the R environment with the Send Bug Report
dialog. I tried R 1.8.1 (Win2K) and R 1.9 (WinXP)
Liaw, Andy [EMAIL PROTECTED] writes:
My $0.02:
If Jack feels that the fact that functions usually return a single list, and
one needs to access the components of the list separately, is somehow
hideous, then I'd rather suggest that R is perhaps the wrong language for
him.
To me the
On Wed, 23-Jun-2004 at 07:32AM -0500, Deepayan Sarkar wrote:
[...]
| 1. You can specify the limits, tick positions, and axis labels
| individually for each panel (see documentation for 'scales' in ?xyplot;
| all the relevant components - at, lab, etc - can be lists).
How difficult would it
Rolf Turner wrote:
fn - function(x) {
list(local_foo=foo, local_bar=bar)
}
OK, that works nicely. (And thanks to Paul Roebuck, who mentioned the
values.list$foo notation.)
(the ``return(...)'' is brought to you by your Department of
Redundancy Department.)
Aha! LISP reards its head. So
Dear, R experts.
Does anybody have experience with 'optim' function?
I have an error message as the following.
Error in optim(transcoefs, fn = hfdeviance, gr = hfdeviance.grad, method
= BFGS, :
initial value in vmmin is not finite
I want to make a comment when this happen.
Is there way I
Quoting Patrick Connolly [EMAIL PROTECTED]:
On Wed, 23-Jun-2004 at 07:32AM -0500, Deepayan Sarkar wrote:
[...]
| 1. You can specify the limits, tick positions, and axis labels
| individually for each panel (see documentation for 'scales' in ?xyplot;
| all the relevant components - at,
On Mon, 21-Jun-2004 at 03:50PM +0100, Barry Rowlingson wrote:
| Duncan Murdoch wrote:
|
|
| I seem to recall that S-PLUS has such a function, but I forget the
| name of it. Probably R does too, on CRAN if not in the base
| packages.
|
| objects.summary() I think it was.
|
| It always
I'm trying to better understand the nlme package and have a few questions.
1.)
Other than using various coding strategies (e.g., dummy coding, effect coding), is
there a way to identify group membership (i.e., treatment) directly? For example, the
following code will fit a two group
Could someone please tell me if it is a way to define
the number of decimals (for example 4 digits) in a
presentation of anova() results. I couldn't apply round
()function.
Thanks in advance.
Alex
__
Acabe com aquelas
Hi!
options(digits = 5)
Good look.
On Wed, 23 Jun 2004, fciclone wrote:
Could someone please tell me if it is a way to define
the number of decimals (for example 4 digits) in a
presentation of anova() results. I couldn't apply round
()function.
Thanks in advance.
Alex
I think I've found a workaround that avoids the two problems in the
replacement function approach.With the definitions of list
and [-.result shown, one can write list[a,b] on the
left side of an assignment where the right side of the assignment
evaluates to a list of the same length (or if
Does anyone have example code for a truncated normal regression using the
survival package? Or any other package? Thanks.
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Just a few more examples:
# swap a and b without explicitly creating a temporary
a - 1; b - 2
list[a,b] - list(b,a)
# get eigenvectors and eigenvalues
list[eval, evec] - eigen(cbind(1,1:3,3:1))
# get today's month, day, year
require(chron)
list[Month, Day, Year] -
I note that ?min actually contains a relevant example:
plot(x, pmin(cH, pmax(-cH, x)), type='b', main= Huber's function)
^^
It's always worth timing things. The timing tests below are the best
of several repetitions, to avoid worries about things like paging
Hi R Users,
Sorry if its out of topic. I would like to ask you about twin peaks -
normal distribution. How R can handle it, any example to explain it in
R.
Thanks,
regards,
Unung
__
[EMAIL PROTECTED] mailing list
Martina Renninger wrote:
Dear All!
How can I cope with overlapping or covered labels (covered by labels
from other data points) in plots?
This doesn't solve every such problem, but it has helped me in the past.
You will probably have to expand the xlim and ylim a bit to fit the extreme
76 matches
Mail list logo