Hi,
Can anyone help me with following? I want to write a program which does
following: for a 2dim surface generate random data, for certain points
on that surface, for which a certain correlation length exists.
Thanx,
Roy
--
The information contained in this communication and any
On Tue, 2005-10-18 at 13:13 -0700, Thomas Lumley wrote:
Both for this and for working out why Tcl/Tk isn't installed you might be
better off trying R-sig-mac rather than r-help
This is a very common piece of advice. However, this is not what you
would imagine if you read the description of
On Wed, 19 Oct 2005, Jari Oksanen wrote:
On Tue, 2005-10-18 at 13:13 -0700, Thomas Lumley wrote:
Both for this and for working out why Tcl/Tk isn't installed you might be
better off trying R-sig-mac rather than r-help
This is a very common piece of advice. However, this is not what you
I was trying to install R on a unix server. Because of a firewall i can't
install biocLite for working with Bioconductor.
With windows it wasn't a problem. I used the option '--internet2' to bypass the
firewall.
I don't have any idea, how to do it with unix.
I tried to set my proxy
Hi There,
I tried to find a function in {waveslim} or {wavethresh} in order to
recunstruct the decomposed signals. As far as I found there is no function in
{waveslim} to recunstruct decomposed data. The function wr{wavethresh}
recunstructs the results of wd function. Apart from its
From ?download.file
These environment variables must be set before the download code
is first used: they cannot be altered later by calling
'Sys.putenv'.
Also, Sys.putenv(http_proxy=http...:8080) sets that to the value
in and not the value of a variable named that: for the
Hi There, I tried to find a function in {waveslim} or {wavethresh} in order to
reconstruct the decomposed signals. As far as I found there is no function in
{waveslim} to reconstruct decomposed data. The function wr{wavethresh}
reconstructs the results of wd function. Apart from its
Guten Tag
Die Krankenkassen-Praemien 2006 sind bereit zum Vergleichen!
Gerne informieren wir Sie ueber die aktuellen Sparmoeglichkeiten fuer das Jahr
2006.
VERGLEICHEN SIE JETZT
Unter http://vergleich.krankenversicherung.ch vergleichen Sie Ihre Praemien
fuer das Jahr 2006.
PRAEMIEN DES
Dear R-list,
I have a problem to open my R workspace.
When I try to open my file .Rdata with double-clik on windows explore I
get the error like this:
Error in load(name, envir = .GlobalEnv) : error reading from connection
and on windows error:
Fatal error: unable to restore saved data in
Hi,
can you tel me how can I make a Forrest Plot with R?
It is possible and easy or are there a more practical free software available?
Than you
Mic
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
Hi everybody:
I've a output interpolation matrix files, and my question is if are there
anybody that could tell me how access to particular values of this matrix
from labels.
Thnaks in advance.
--
Antoni Viúdez Mora
Dept. Dinámica de
Hi
Your .Rdata file is probably corrupted. Unless you have a working
copy of it elsewhere or sources of your data together with history
of your commands you are probably in deep trouble .Rdata is a
binary format and it is not recommended to safe and reliable saving
of your work as you have
Hi All,
Does R have the function for FIGARCH?
SUMANTA BASAK.
---
This e-mail may contain confidential and/or privileged infor...{{dropped}}
Dennis
Try
TEST[-3]
[[1]]
[1] A1 A2
[[2]]
[1] B1 B2
for removing more than one element from the list (say 2 3) --
TEST[-c(2,3)]
[[1]]
[1] A1 A2
HTH
John
Dennis Fisher wrote---
Colleagues,
I have created a list in the following manner:
TEST- list(c(A1, A2), c(B1, B2),
BDR == Prof Brian Ripley [EMAIL PROTECTED]
on Wed, 19 Oct 2005 08:58:20 +0100 (BST) writes:
BDR On Wed, 19 Oct 2005, Jari Oksanen wrote:
On Tue, 2005-10-18 at 13:13 -0700, Thomas Lumley wrote:
Both for this and for working out why Tcl/Tk isn't
installed you might be
Hello,
2005/10/19, Michela Ballardini [EMAIL PROTECTED]:
Hi,
can you tel me how can I make a Forrest Plot with R?
It is possible and easy or are there a more practical free software available?
Maybe this will help you to find an answer:
Your .Rdata file is probably corrupted.
I will investigate. Thanks for you info.
Because of my file .RData very large about 75MB.
Best wishes, Muhammad Subianto
On this day 19/10/2005 12:21 PM, Petr Pikal wrote:
Hi
Your .Rdata file is probably corrupted. Unless you have a working
copy
Hi !
I try to compare some models obtained from glmmPQL.
model1 -
glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4
+I(freq8_4^2), random=~1|num, binomial);
model2 -
glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4
, random=~1|num, binomial);
Hi,
I discover that when I save a workspace containing a genefilter (pkg
from Bioconductor) object I cannot open no more after. I have to
restore the .RData file from a backup to be able to start R again.
I didn't upgrade to Version 2.2 but I'm not sure that it will solve the
problem.
Did
Hello to everybody,
I'd like to submit a problem I'm dealing with, and I can't get an answer to
by myself.
I have to test if my data come from a specific probability distribution, of
which I know the analytic form both of the p.d.f. and the c.d.f.
Namely, it is the hypoexponential distribution,
Dear All,
I am wonderng how to test whether the data follows poisson distribution.
Thank you so much!
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
I'm trying to run a nls on a subset of a data.frame.
In the subset, one observation is NA. So I drop the
observation but when I ask for :
sm - nls(machin$revcum ~
Lc.singh(machin$popcum,p), start=list(p=c(2,3)))
I get :
Erreur dans parse(file, n, text, prompt) : syntax
error in ~
If I put
On 10/19/2005 8:37 AM, David Ruau wrote:
Hi,
I discover that when I save a workspace containing a genefilter (pkg
from Bioconductor) object I cannot open no more after. I have to
restore the .RData file from a backup to be able to start R again.
I didn't upgrade to Version 2.2 but I'm not
Le 19.10.2005 13:56, Thomas Schönhoff a écrit :
Hello,
2005/10/19, Michela Ballardini [EMAIL PROTECTED]:
Hi,
can you tel me how can I make a Forrest Plot with R?
It is possible and easy or are there a more practical free software available?
Maybe this will help you to find an answer:
If you're really interested in reading the source for functions and aren't
interested in tracking down various methods (possibly hidden in namespaces) at
the R prompt, I think it's much easier to download the R source code from CRAN
and go through the original source files.
-roger
Wolfrum, Ed
Hi
Does anyone know where is the package: npmc (Nonparametric Multiple
Comparisons).
I found the reference on R Site Search, but not the package itself on
CRAN as suggested.
Thanks
Mauricio
__
R-help@stat.math.ethz.ch mailing list
Carlos Mauricio Cardeal Mendes wrote:
Hi
Does anyone know where is the package: npmc (Nonparametric Multiple
Comparisons).
I found the reference on R Site Search, but not the package itself on
CRAN as suggested.
The packages is ORPHANED and removed from the CRAN main repository.
You
Le 19.10.2005 15:11, Romain Francois a écrit :
Le 19.10.2005 13:56, Thomas Schönhoff a écrit :
Hello,
2005/10/19, Michela Ballardini [EMAIL PROTECTED]:
i,
can you tel me how can I make a Forrest Plot with R?
It is possible and easy or are there a more practical free software
available?
I try to compare some models obtained from glmmPQL.
model1 -
glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^
2)+freq8_4
+I(freq8_4^2), random=~1|num, binomial);
model2 -
glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^
2)+freq8_4
, random=~1|num, binomial);
BertG == Berton Gunter [EMAIL PROTECTED]
on Wed, 12 Oct 2005 10:47:40 -0700 writes:
BertG I hope that this job posting does not offend
BertG anyone. My sincere apologies if it is inappropriate
BertG -- blame me, not my company.
since you ask for it .. ;-)
BertG -- Bert
Thanks a lot Professor Ripley, it works that way.
I will upgrade to v 2.2 one day when I have time and see if it still
happen.
David
On Oct 19, 2005, at 15:58, Prof Brian Ripley wrote:
On Wed, 19 Oct 2005, Duncan Murdoch wrote:
On 10/19/2005 8:37 AM, David Ruau wrote:
Hi,
I discover that
To R-project help
I work with a clinical dataset including about 27,000 patients and I wanted to
do some model diagnostics on cox regression models with frailty effects (for
every patient). But in this case the cox.zph function did not work. I always
got the following error message:
Fehler
On Wed, 19 Oct 2005, Romain Francois wrote:
Pretty interresting. You just pointed to a good candidate for r graph
gallery.
There is an example on the R home page, and has been for some time. I have
an improved version based on grid code by Paul Murrell, but it's not on a
nearby computer. The
Hi
I have a matrix u, for which diag() gives an error:
u - structure(c(5.42334674128216, -2.31319389204264, -5.83059042218476,
-1.64112369640695, -2.31319389212801, 3.22737617646609,
1.85200668021569,
-0.57102273078531, -5.83059042231881, 1.85200668008156,
11.9488923894962,
maybe you should keep your NA but deal with the na.action option of the
nls() function ?
HTH.
Florence.
On 10/19/05, Florent Bresson [EMAIL PROTECTED] wrote:
I'm trying to run a nls on a subset of a data.frame.
In the subset, one observation is NA. So I drop the
observation but when I ask
Jose == José Ernesto Jardim [EMAIL PROTECTED]
on Tue, 18 Oct 2005 15:30:59 +0100 writes:
Jose Dennis Fisher wrote:
Colleagues,
I have created a list in the following manner:
TEST- list(c(A1, A2), c(B1, B2), c(C1, C2))
I now want to delete one element
It has to do with the NAs in dimnames, but I don't know why this is a
problem. If you either got rid of dimnames or assigned actual names, this
works fine.
Ravi.
-Original Message-
From: [EMAIL PROTECTED] [mailto:r-help-
[EMAIL PROTECTED] On Behalf Of Robin Hankin
Sent: Wednesday,
On 10/19/2005 10:55 AM, Robin Hankin wrote:
Hi
I have a matrix u, for which diag() gives an error:
u - structure(c(5.42334674128216, -2.31319389204264, -5.83059042218476,
-1.64112369640695, -2.31319389212801, 3.22737617646609,
1.85200668021569,
-0.57102273078531, -5.83059042231881,
On 10/19/2005 10:55 AM, Robin Hankin wrote:
Hi
I have a matrix u, for which diag() gives an error:
u - structure(c(5.42334674128216, -2.31319389204264, -5.83059042218476,
-1.64112369640695, -2.31319389212801, 3.22737617646609,
1.85200668021569,
-0.57102273078531,
Hi,
I am hoping some one can help with this.
I am using nlme to fit a random coefficients model. It ran for hours before
returning
Error: Singularity in backsolve at level 0, block 1
The model is
On Wed, 19 Oct 2005, Duncan Murdoch wrote:
On 10/19/2005 10:55 AM, Robin Hankin wrote:
Hi
I have a matrix u, for which diag() gives an error:
u - structure(c(5.42334674128216, -2.31319389204264, -5.83059042218476,
-1.64112369640695, -2.31319389212801, 3.22737617646609,
1.85200668021569,
On 10/19/2005 11:37 AM, Prof Brian Ripley wrote:
On Wed, 19 Oct 2005, Duncan Murdoch wrote:
On 10/19/2005 10:55 AM, Robin Hankin wrote:
Hi
I have a matrix u, for which diag() gives an error:
u - structure(c(5.42334674128216, -2.31319389204264, -5.83059042218476,
-1.64112369640695,
On 10/18/05, Kiermeier, Andreas (PIRSA - SARDI)
[EMAIL PROTECTED] wrote:
Dear all,
back in 2002 Martin Henry H. Stevens wrote
(https://stat.ethz.ch/pipermail/r-help/2002-May/019851.html)
How do I control the text in strips? Specifically, I want to put in the
ranges generated in shingle(x)
Hi,
2005/10/19, Wensui Liu [EMAIL PROTECTED]:
Dear All,
I am wonderng how to test whether the data follows poisson distribution.
Thank you so much!
Did you notice the PDF on distribution tests using R by Vito Ricci,
its found at CRAN in the docs contrib section, called FITTING
To be pedantic (I'm feeling cranky today):
One can never test whether the data follow [data is plural] a Poisson
distribution -- only whether there is sufficient evidence to cast that
assumption into doubt. Perhaps a better shorthand is whether the data are
consistent with Poisonness . This
Hi,
I have a problem about graphics. I would like to plot two graphs: a barplot
and curve. Here is the code:
barplot(dpois(0:45,20),xlim=c(0,45),names=0:45)
curve(dnorm(x,20,sqrt(20)),from=0,to=45,add=T)
Both graphs are drawn in the same figure, however the scale in both graphs
dooes not
On Wed, 19 Oct 2005 10:18:02 -0700 Berton Gunter wrote:
To be pedantic (I'm feeling cranky today):
One can never test whether the data follow [data is plural] a
Poisson distribution -- only whether there is sufficient evidence to
cast that assumption into doubt. Perhaps a better shorthand
Le 19 Octobre 2005 11:09, Martin Maechler a écrit :
[...]
Lists can have 'dim' attributes and hence be treated as arrays;
For me, this is an amazing feature that I discovered almost by accident (I
tried it an it worked)! This creates a sort of three-dimensional object ---
much like an array
So, is there another package to substitute those functions described on
ORPHANED npmc package ?
Regards,
Mauricio
Brazil
Uwe Ligges escreveu:
Carlos Mauricio Cardeal Mendes wrote:
Hi
Does anyone know where is the package: npmc (Nonparametric Multiple
Comparisons).
I found the
Hi,
I make a model:
m.null - glm(y~1)
m.comp - glm(y~x1+x2+x3)
I try make a anova like this
anova(m.null,m.comp,test=F)
The result:
Erro em anova.glmlist(c(list(object), dotargs), dispersion = dispersion, :
models were not all fitted to the same size of dataset
The problem is that I
I am looking to create what I would call a simple variation on the
boxplot. What I would like to do is to be able to plot the upper and
lower confidence limits as the box and the 10th and 90th percentile as
the whiskers. What I have done is write the code to create a dataframe,
the columns of
Hi
Mike Bock wrote:
I am looking to create what I would call a simple variation on the
boxplot. What I would like to do is to be able to plot the upper and
lower confidence limits as the box and the 10th and 90th percentile as
the whiskers. What I have done is write the code to create a
hi
It seems that, in some cases, saving a .RData
containing objects from some packages and later trying
to start R (or load the .RData) before loading the
package could cause this error. See
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/46910.html
If the .RData contains
The first example on this page is pretty close to what I want, but the
x-axis would be a category (like location), and the true range
anotations are left out. I don't; have a place to post but can e-mail a
pdf of my excel version by request.
Using base graphics you can use the bxp function. This is what
boxplot calls to do the actual plotting. Look at the return value
for boxplot in the help for boxplot to see the form of the data
that needs to be passed to bxp (also look at the help for bxp)
You could write a simple function that
On 10/19/05, Mike Bock [EMAIL PROTECTED] wrote:
I am looking to create what I would call a simple variation on the
boxplot. What I would like to do is to be able to plot the upper and
lower confidence limits as the box and the 10th and 90th percentile as
the whiskers. What I have done is
Le 19.10.2005 19:27, Rui Cardoso a écrit :
Hi,
I have a problem about graphics. I would like to plot two graphs: a barplot
and curve. Here is the code:
barplot(dpois(0:45,20),xlim=c(0,45),names=0:45)
curve(dnorm(x,20,sqrt(20)),from=0,to=45,add=T)
Both graphs are drawn in the same figure,
On Wed, 2005-10-19 at 18:27 +0100, Rui Cardoso wrote:
Hi,
I have a problem about graphics. I would like to plot two graphs: a barplot
and curve. Here is the code:
barplot(dpois(0:45,20),xlim=c(0,45),names=0:45)
curve(dnorm(x,20,sqrt(20)),from=0,to=45,add=T)
Both graphs are drawn
I am wondering if someone would have any suggestion about my issue?
I have the following code:
wgts-aggregate(subset(lendata,select=c(Length)),list(lendata$Cruise,len
data$Station,lendata$Region,lendata$Total),mean)
wgts-wgts[order(wgts$Group.3,wgts$Group.1,wgts$Group.1),]
Dear R users,
I am trying to find out what the function survReg does exactly with the
Weights parameters. I looked in Terry Therneau documentations and other
places and couldn't find anything. I tried to make an analogy with weighted
OLS and assumed that the scale parameter Sigma in the
On 19/ott/05, at 13:36, Martin Maechler wrote:
BDR == Prof Brian Ripley [EMAIL PROTECTED]
on Wed, 19 Oct 2005 08:58:20 +0100 (BST) writes:
BDR On Wed, 19 Oct 2005, Jari Oksanen wrote:
On Tue, 2005-10-18 at 13:13 -0700, Thomas Lumley wrote:
Both for this and for working out why
Hi all,
I wanted to run the hclust (or any other clustering algorithm) on a distance
matrix. I have formed the distance matrix as:
distmat:
ab c de
a0.000.961.601.601.68
b0.960.000.961.802.64
c 1.600.96
Nelson, Gary (FWE) [EMAIL PROTECTED] writes:
I am wondering if someone would have any suggestion about my issue?
I have the following code:
wgts-aggregate(subset(lendata,select=c(Length)),list(lendata$Cruise,len
data$Station,lendata$Region,lendata$Total),mean)
Hello,
I have a follow-up from Jens's question and Professor Ripley's response.
Jens wants to do quadratic optimization with 2 constraints:
# I need two constraints:
# 1. each element in par needs to be between 0 and 1
# 2. sum(par)=1, i.e. the elements in par need to sum to 1
how
Have you tried as.dist(distmat)?
Christian
On Wed, 19 Oct 2005, Tim Smith wrote:
Hi all,
I wanted to run the hclust (or any other clustering algorithm) on a distance
matrix. I have formed the distance matrix as:
distmat:
ab c de
a0.00
Dear alls,
Could someone tell me how to select a subset of string observations (e.g.
females in a sex column) with sqlQuery in the RODBC library?
Indeed, I'm trying to select a subset of observations on my access database
with:
female-sqlQuery(mychannel,SELECT Micromammiferes.sex
FROM
Dear R-Users,
how to include error bars within lattice?
How should the panel = function(x,y,...){
looks like?
Does panel.arrows works here as well?
I appreciate any help on this.
Regards,
Mario AT
[[alternative HTML version deleted]]
__
Dr. Williams,
I ran across your inquiry on one of the R-help mailing lists regarding
digital filter design and implementation. I found no response to your
email in the archives and was wondering if you were able to find anything.
Thanks,
Israel
--
Israel Christie, Ph.D.
Email: [EMAIL
Jérôme Lemaître wrote:
Dear alls,
Could someone tell me how to select a subset of string observations (e.g.
females in a sex column) with sqlQuery in the RODBC library?
Indeed, I'm trying to select a subset of observations on my access database
with:
female-sqlQuery(mychannel,SELECT
(sorry if a duplicate pops through ...)
Jérôme Lemaître wrote:
Dear alls,
Could someone tell me how to select a subset of string observations (e.g.
females in a sex column) with sqlQuery in the RODBC library?
Indeed, I'm trying to select a subset of observations on my access database
?llines, lsegments and the like can be used in the panel functions to draw
any sort of error bar that you can compute from the x,y,... data of the
panel.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific
The syntax error is that you have unescaped quotes inside quotes.
You also do not need a semicolon, nor to refer to columns in this
table.column form. Try
'select sex from Micromammiferes where sex=females'
(I suspect you do not need quotes, but keep forgetting the quirks of
various DBMSs.)
Sasha,
This is R, which has survreg with weights (no caps). You can find out for
sure from the source code, but the help file says these are `observation
weights'. I take that to mean weight=2 says `I have two cases like this'
in forming the likelihood. The C code says `case weights', the
mantelhaen.test() gives the exact conditional p-value (for independence) and
confidence intervals (CIs)for the common odds ratio for a stratified 2x2 table.
The epitools package by Tomas Aragon (available via CRAN) contains functions
which use fisher.test() to calculate mid-p exact p-values and
I'm new to R and have searched for help and consulted the the pdf
manuals, but I can't seem to figure out how to plot more than one
series on the same graph.
I've tried using multiple par(new=TRUE) statements such as
plot(series1, ci.type=line, col=red, lwd=2, ci.lty=0, ci.col=red)
par(new=TRUE)
[EMAIL PROTECTED] wrote:
I'm new to R and have searched for help and consulted the the pdf
manuals, but I can't seem to figure out how to plot more than one
series on the same graph.
I've tried using multiple par(new=TRUE) statements such as
plot(series1, ci.type=line, col=red, lwd=2,
Hello,
I have read through the manuals and can't seem to find an answer.
I have a categorical, character variable that has hundreds of values. I want
to group the existing values of this variable into a new, derived (categorical)
variable by applying conditions to the values in the data.
For
Dear Prof Ripley,
I tried what you suggested and it worked perfectly well!
I also appreciated your corrections about my SQL formulation. However, it
seems that quotes are indeed needed to get string observations (but my
access software is in french...).
Finally, do you think that it might be
I suggest you use the recode function in car package to do your job.
=== 2005-10-20 08:09:08 您在来信中写道:===
Hello,
I have read through the manuals and can't seem to find an answer.
I have a categorical, character variable that has hundreds of values. I want
to group the
79 matches
Mail list logo