Fred J. wrote:
Philippe Grosjean [EMAIL PROTECTED] wrote: What Fred is looking for is
local minima/maxima, also known as turning
points, or pits/peaks in a series. You can look at ?turnpoints in
pastecs library.
x - c(1:4,0:5, 4, 11)
x
[1] 1 2 3 4 0 1 2 3 4 5 4 11
Dear R-users,
I tried to export a dataframe form R to Access, that way:
library(RODBC)
channel - odbcConnectAccess(d:/test.mdb)
sqlSave(channel=channel, flore, rownames=F)
odbcClose(channel)
But I always got this error message:
Erreur dans sqlSave(channel = channel, flore, Florist) :
[RODBC]
Hi Context,
It seems to me that the changes in entropy are simply due to the binning.
The usual entropy would be the limit of the discrete entropy when the number
of bins goes to infinity. For the discrete entropy to be a useful
approximation to the entropy the number of bins should be reasonably
OK, I finally found what's wrong - date column name.
Jacques VESLOT a écrit :
Dear R-users,
I tried to export a dataframe form R to Access, that way:
library(RODBC)
channel - odbcConnectAccess(d:/test.mdb)
sqlSave(channel=channel, flore, rownames=F)
odbcClose(channel)
But I always got this
hadley wickham wrote:
Thing is, for one reason or another, the number of marbles per bag may
systematically vary with age too. However, I am not interested in the
number of marbles per bag, so I would like to group the students into 8
groups such that each group has the same total number of
Gabor Grothendieck wrote:
On 3/19/06, Dan Bolser [EMAIL PROTECTED] wrote:
Gabor Grothendieck wrote:
On 3/18/06, Dan Bolser [EMAIL PROTECTED] wrote:
Gabor Grothendieck wrote:
If you are just looking for something simple that may be good enough
then assign the largest one to group 1, the
Dear list mates,
Are {utils} dialog box functions, and winMenuAdd... functions used to change
(e.g. Console) menus, platform dependent?
I'm writing a script loaded with .First that provides first time users in a
lab course with the ability to select, load, and change between what I
called a
Are you saying that your data might look like this ?
set.seed(1) # For reproducibility only - remove this
mydf - data.frame( age=round(runif(100, min=5, max=65), digits=1),
nred=rpois(100, lambda=10),
nblue=rpois(100, lambda=5),
Hello everybody,
I am wondering if it is possible to create a gui to plot a time series
that is very big, it's an EEG signal of 20mins. What I would like to do
is plot the first 5mins, then have a button on the gui that plots the
next 5mins when pushed.
Is it possible?
Thanks in advance !
I haven't been able to find an answer on this that's direct, only
implied. In several places I have read that when people asked for
sphericity tests they were guided toward lme or mlm models. But,
there is no direct indication that the lme method is not subject to
the sphericity
http://cran.r-project.org/Rlogo.jpg
In the future all you have to do (on PC) is to right click on image you
want and press properties to get the exact link.
Jarek Tuszynski
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Erin Hodgess
Sent: Wednesday,
See demos in the tcltk package. I remember there was some example of
interactive graphics, with buttons to change graphical parameters. You can
start from there...
Antonio, Fabio Di Narzo.
2006/3/20, Gael de Lannoy [EMAIL PROTECTED]:
Hello everybody,
I am wondering if it is possible to
John:
Linear models can have different covariance structures to accommodate
certain dependencies in the data. The functions for data analysis in the
lme4 package are flexible and can be structured as such. For example,
when random intercepts only are included, this is akin to compound
symmetry,
I might be missing something but I thought that AUC was a measure for
comparing ROC curves, so there is nothing else needed to compare them. The
larger AUC is the higher correlation of 2 variables compared. No other
measures or calculations are needed.
Jarek Tuszynski
-Original Message-
Keith Chamberlain wrote:
Dear list mates,
Are {utils} dialog box functions, and winMenuAdd... functions used to change
(e.g. Console) menus, platform dependent?
Yes, Windows-only as the 'win' prefix suggests it.
I'm writing a script loaded with .First that provides first time users in a
R Users,
I was wondering if someone might point me in the right direction. I am using a
Cox model (survival package) to evaluate survival of pen-reared birds (time to
event data collected daily) and I have been trying to determine how I can
estimate a 'daily' survival rate and std error from
hello,
i need the data from hist() but i do not want the plot.
e.g.
z=hist(data)$counts #returns absolute frequency
but when i execute this command the plot occurs also. is it possible to
suppress the plot?
many thanks,
best regards gg
--
---
Le 20.03.2006 14:23, Gottfried Gruber a écrit :
hello,
i need the data from hist() but i do not want the plot.
e.g.
z=hist(data)$counts #returns absolute frequency
but when i execute this command the plot occurs also. is it possible to
suppress the plot?
many thanks,
best regards gg
From help page
If plot=TRUE, the resulting object of class histogram is plotted by
plot.histogram, before it is returned.
HTH
Petr
On 20 Mar 2006 at 14:23, Gottfried Gruber wrote:
From: Gottfried Gruber [EMAIL PROTECTED]
To: r-help
hist(data, plot=FALSE)$counts
On Mon, 2006-03-20 at 14:23 +0100, Gottfried Gruber wrote:
hello,
i need the data from hist() but i do not want the plot.
e.g.
z=hist(data)$counts #returns absolute frequency
but when i execute this command the plot occurs also. is it possible to
Look at:
?hist
z = hist(data, plot=FALSE)
Best,
Matthias
hello,
i need the data from hist() but i do not want the plot.
e.g.
z=hist(data)$counts #returns absolute frequency
but when i execute this command the plot occurs also. is it
possible to
suppress the plot?
many thanks,
Hello!
Book by Pinheiro and Bates is being repeatedly suggested in connection
with lme() from nlme package. Is this book also suggested for a
reference use with lmer in lme4 package?
--
Lep pozdrav / With regards,
Gregor Gorjanc
On 20 Mar 2006 at 14:23, Gottfried Gruber wrote:
From: Gottfried Gruber [EMAIL PROTECTED]
To: r-help r-help@stat.math.ethz.ch
Date sent: Mon, 20 Mar 2006 14:23:26 +0100
Subject:[R] hist-data without plot
hello,
i need the
Gottfried Gruber wrote:
hello,
i need the data from hist() but i do not want the plot.
e.g.
z=hist(data)$counts #returns absolute frequency
but when i execute this command the plot occurs also. is it possible to
suppress the plot?
Yes, as mentioned in ?hist
Uwe Ligges
many
STATISTICAL PRACTICE IN EPIDEMIOLOGY USING R
Tartu, Estonia, Thursday 8 - Tuesday 13 June 2006
Application deadline: 15 April 2006.
The course is aimed at epidemiologists and statisticians who wish to
use R for statistical modelling and analysis of
Well, not if you are interested in lmer() syntax. The book is good with
the conceptual issues in modeling the data, but does not include any
reference to fitting models with lmer.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Gregor Gorjanc
Sent:
Hi,
I don't blame anybody (i'm guilty too), but have we forget the section
'Responding to other post' in the posting guide.
The recent thread 'hist-data without plot' is a good example of a very
simple question, isn't it ?
The poster got responses in 33, 35, 35, 36, 37 and 40 minutes (according
Hi, trying this again - any help would be great!
Basically, I have two matrices of equal dimension, one of them should
produce something similar to a heatmap.. The 2nd matrix should be the
heights for each value on the heatmap - producing a sort of 3d surface
plot.
Viewing this seems like a
Hallo!
Does anybody know how to make R calculate confidence limits for the
Hazard Ratio based on the likelihood function when doing Cox
Regression? As the p Values and the confidence intervals of Wald
statistics differ from the p Value of the likelihood function
sometimes likelihood-p is .05
You can also refer to the following papers which have proposed Rsquares for
the linear mixed model:
Vonesh, E. F., Chinchilli, V. P., and Pu, K. W. (1996), Goodness-of-fit in
generalized nonlinear mixed-effects models, Biometrics, 52, 572-587.
Zheng, B. Y. (2000), Summarizing the goodness of
Adaikalavan Ramasamy wrote:
Are you saying that your data might look like this ?
set.seed(1) # For reproducibility only - remove this
mydf - data.frame( age=round(runif(100, min=5, max=65), digits=1),
nred=rpois(100, lambda=10),
nblue=rpois(100,
Doran, Harold wrote:
Well, not if you are interested in lmer() syntax. The book is good with
the conceptual issues in modeling the data, but does not include any
reference to fitting models with lmer.
OK, I of course suspected that lmer will not be in, but I wanted to hear
that the book has
Cashorali, Tanya wrote:
Hi, trying this again - any help would be great!
Basically, I have two matrices of equal dimension, one of them should
produce something similar to a heatmap.. The 2nd matrix should be the
heights for each value on the heatmap - producing a sort of 3d surface
plot.
Hi R people!
:-)
I am printing a table of percentages (P) and their 95% confidence
interval half-widths (CI).
I would like to
(i) create a character string of P and CI that has the plus/minus
character to the left of CI, then
(ii) use write.csv to save that file with the P, CI and the
Tuszynski, Jaroslaw W. wrote:
I might be missing something but I thought that AUC was a measure for
comparing ROC curves, so there is nothing else needed to compare them. The
larger AUC is the higher correlation of 2 variables compared. No other
measures or calculations are needed.
Jarek
Adapt the function below to suit your needs. If you really want to
plot 5 minutes at a time, round the time series to the last MM:00
times (where MM is in 5*0:11) and have idx below loop over them.
splitplot - function(x,points)
{
boundaries - c(1,points*1:floor(length(x)/points),length(x))
for
Dear R-users,
I would like to know if there is any way to minimize || y - Xb||^2 under the
constraint b'b=1.
Thank you
Giancarlo
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Dear R-users,
I would like to know if there is any way to minimize || y - Xb||^2 under the
constraint b'b=1.
Thank you
Giancarlo
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
RSiteSearch('penalized regression',restr='func')
RSiteSearch*'ridge regression', restr='func')
Or use CRAN's search capabilities in your browser. There is no need to post
such questions to the list.
Also, check out the LARS package, which might have more capabiliies and
flexibility along the
Since you did not provide an example, I'm not sure I understand your
question. I will nevertheless offer some suggestions:
1. If you have not already, please review Pinheiro and Bates (2000)
Mixed-Effects Models in S and S-Plus (Springer).
2. Also, please
I don't understand: I just ran the first example in the aov help
page, and it produced F ratios and p values.
If this does not answer your question, I suggest you read Pinheiro
and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer) if you
haven't already and try
* Adaikalavan Ramasamy [EMAIL PROTECTED] [2006-03-19 04:51:19 +]:
1) R-help is designed for and by unpaid volunteers. Therefore
sometimes RTFM without page reference is quite acceptable.
I am an unpaid volunteer maintainer of CLISP (http://clisp.cons.org).
I often answer questions with
Hi,
I'm a PhD student and I want to use the function 'daisy' from the
package 'cluster' to compute dissimilarities.
My variables are of mixed types so I use the argument 'stand' in daisy
to define the type of my variables.
I have the following error message :
Warning message:
binary
If you do not feel like answering a question, it is perfectly
fine with
me, I do not think that anyone owes me anything.
All I am asking is that if you do decide to answer, please make your
answer immediately useful, i.e., not requiring learning all
the manual by
heart (a specific manual
* Berton Gunter [EMAIL PROTECTED] [2006-03-20 09:42:49 -0800]:
If you do not feel like answering a question, it is perfectly
fine with
me, I do not think that anyone owes me anything.
All I am asking is that if you do decide to answer, please make your
answer immediately useful, i.e., not
Thanks a lot for your help!
Now I have a related problem, that I still can not solve:
I plot a t distribution using:
curve(dt(x, df = 20), xlim=c(-5,5))
The upper 5% of that distribution using:
qt(0.05, df=20, lower.tail=FALSE)
How can define the polygon to highlight the area representing
Lets say there are 10 students in the first group and denote x1 as (say)
the number of red balls for student 1 and s1 the total balls. Then I was
calculating the average the proportion ( x1/s1 + x2/s2 + ... + x10/s10 )
and you were calculating the average number of events (x1+x2
I am pretty sure I saw this mentioned in
http://cran.r-project.org/doc/manuals/R-intro.html
but I cannot recall how it was called. sorry...
how do I return (and accept) multiple return values from a function?
e.g., in lisp (multiple-value-bind (f r) (floor 10 3) (list f r))
will return list (3
[[ Please ignore the last email which was sent incomplete ]]
Lets say there are 10 students in the first group and denote x1 as (say)
the number of red balls for student 1 and s1 the total balls. Then I was
calculating the average the proportion ( x1/s1 + x2/s2 + ... + x10/s10 )
and you were
On Mon, 2006-03-20 at 13:06 -0500, Sam Steingold wrote:
I am pretty sure I saw this mentioned in
http://cran.r-project.org/doc/manuals/R-intro.html
but I cannot recall how it was called. sorry...
how do I return (and accept) multiple return values from a function?
e.g., in lisp
Dear Philippe, ( list)
I am not sure what is meant by 'floating window' but thank you for
clarifying the 'win' prefix. I had been considering different possible
senses of that prefix: win as in 'ts' window windowing functions (clearly
not), a generic 'win'dow function (abbreviated) for some GUI,
R-help,
I'm getting some unexpected behavior with subsetting a data frame
(aircraft flight data) that I can't sort out.
Here is a simplified version of my data frame and problem:
flight
FlightID TailNo FlightDate HobbsTime FlightCost Date year
1 4497 6009K NA
Gavin Simpson wrote:
On Mon, 2006-03-20 at 13:06 -0500, Sam Steingold wrote:
I am pretty sure I saw this mentioned in
http://cran.r-project.org/doc/manuals/R-intro.html
but I cannot recall how it was called. sorry...
how do I return (and accept) multiple return values from a function?
e.g.,
Gael,
Try the following to get you started:
library(tkrplot)
y - rnorm(1, 10, 2) + 5*sin( (1:1)/1000 )
tt - tktoplevel()
left - tclVar(1)
oldleft - tclVar(1)
right - tclVar(100)
f1 - function(){
lleft - as.numeric(tclvalue(left))
rright - as.numeric(tclvalue(right))
__
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
Eric Archer wrote:
R-help,
I'm getting some unexpected behavior with subsetting a data frame
(aircraft flight data) that I can't sort out.
Here is a simplified version of my data frame and problem:
flight
FlightID TailNo FlightDate HobbsTime FlightCost Date year
1
Had I just looked at flight$year 2006, I would've seen what was up.
Thanks much Peter!
Cheers,
eric
P Ehlers wrote:
[snip]
flight$year 2006 will return TRUE/FALSE, not row numbers. Try this:
errors - subset(flight, subset = year 2006)
Peter Ehlers
--
Eric Archer, Ph.D.
On 3/20/2006 10:24 AM, Smith, Phil wrote:
Hi R people!
:-)
I am printing a table of percentages (P) and their 95% confidence
interval half-widths (CI).
I would like to
(i) create a character string of P and CI that has the plus/minus
character to the left of CI, then
(ii) use
Hello, friends of mosaicplot,
I want to attach a fixed legend to mosaicplot, e.g. -3..3 independently of
the data in the graphics. This must be hidden somehow in shading/legend, but
I could not get the right handle to it.
Dieter
__
Eric Archer wrote on 20 Mar 2006 19:46:44 MET:
I'm getting some unexpected behavior with subsetting a data
frame (aircraft flight data) that I can't sort out. Here is a
simplified version of my data frame and problem:
flight
FlightID TailNo FlightDate HobbsTime FlightCost
On 3/20/06, P Ehlers [EMAIL PROTECTED] wrote:
Eric Archer wrote:
R-help,
I'm getting some unexpected behavior with subsetting a data frame
(aircraft flight data) that I can't sort out.
Here is a simplified version of my data frame and problem:
flight
FlightID TailNo
Dear R-people:
François Michonneau's method to obtain the special character plus/minus works
on Windows 2000 professional.
Many Thanks to François for his work!
Phil Smith
Centers for Disease Control and Prevention
Atlanta, GA
_
From: François MICHONNEAU [mailto:[EMAIL PROTECTED]
Hi All,
I found a solution for my question:
I have the following adjacency matrix for a directed graph:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,]00000000
[2,]00000000
[3,]10000000
[4,]
Georg Otto [EMAIL PROTECTED] writes:
Thanks a lot for your help!
Now I have a related problem, that I still can not solve:
I plot a t distribution using:
curve(dt(x, df = 20), xlim=c(-5,5))
The upper 5% of that distribution using:
qt(0.05, df=20, lower.tail=FALSE)
How can
Dear Erich and John, thankyou for providing RColorBrewer!
Are you working on options for CMYK or RGB return values?
Thanks again, Darren
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
[Berton Gunter]
[Sam Steingold]
PPS. how do I figure out the number of rows in a data.frame?
is length(attr(X,row.names)) the right way?
help.search(number of rows) immediately gets you your answer!
Hi, people. Here, I get:
Help files with alias or concept or title matching ‘number
Dan Bolser [EMAIL PROTECTED] wrote:
Gabor Grothendieck wrote:
Do you mean you want g to be in the original order of x?
No. What I mean is that I want to order x by any particular variable in
my data.frame, then group over x such that each group has roughly the
same sum.
I get the feeling
Hi there,
I use function kmeans and clara to cluster one flow cytometry
dataset. By using function plot, the clusters got from clara can be
graphed, while kmeans not. How can I get the plot of the clusters of
kmeans?
And, I hope to compare the two methods kmeans and clara, or in other
Hello,
Is it possible to pass R matrix as a parameter to an internal C procedure?
From the documentation I got the impression that only 1-dim vectors can be
passed.
Why the following wont work for me?
a-matrix(1:15,3,5)
.C(pr,as.integer(a))
void pr(int **a){
On Mon, 20 Mar 2006, Darren Weber wrote:
Dear Erich and John, thankyou for providing RColorBrewer!
Are you working on options for CMYK or RGB return values?
col2rgb() will give RGB values. convertColor() (or the colorspace
package) will give other conversions (but not CMYK)
-thomas
You need to pass it as a vector and then you can convert to an array
structure in C. Look at carray.c for one way to do this, there are others.
Browse the source code and see how different authors handle your problem.
Matt Austin
Statistician
Amgen, Inc
800 9AMGEN9 x77431
805-447-7431
From: johan Faux
Hello,
Is it possible to pass R matrix as a parameter to an
internal C procedure?
From the documentation I got the impression that only 1-dim
vectors can be passed.
Why the following wont work for me?
In short, because a matrix in R is just a vector with
On Mon, 20 Mar 2006, johan Faux wrote:
Hello,
Is it possible to pass R matrix as a parameter to an internal C procedure?
From the documentation I got the impression that only 1-dim vectors can be
passed.
Why the following wont work for me?
a-matrix(1:15,3,5)
.C(pr,as.integer(a))
Hi,
I would like to repeatedly create some simulated survival data.
Is there an R package that has a function that can do this?
I would like to simulate a clinical trial that has 2 study arms and the outcome
of interest is measured up to 30 days from randomization.
Thanks for the help.
Hi,
This is a very fundamental question. I want to export graphical results
so that I can place them in an openoffice document. I use Fedora 5. Can
anyone help? Thanks in advance.
Takatsugu Kobayashi
PhD Student
Indiana University, Dept. Geography
On Mon, 2006-03-20 at 17:33 -0500, [EMAIL PROTECTED] wrote:
Hi,
This is a very fundamental question. I want to export graphical results
so that I can place them in an openoffice document.
I use Fedora 5.
That was quick... ;-)
The best way to do this under Linux is to use the R
[EMAIL PROTECTED] wrote:
Hi,
I would like to repeatedly create some simulated survival data.
Is there an R package that has a function that can do this?
I would like to simulate a clinical trial that has 2 study arms and the
outcome of interest is measured up to 30 days from
Dear all (especially R-Win users),
I'm looking for a good way to use R functionality from within a C# application.
One way would be via (D)COM, but I would prefer a more direct solution: My
main concern is to ship all necessary files in one easy-to-install package, so
that users do not need to
Another way that works on my Windows XP system is
to use \261 .
On 3/20/06, Smith, Phil [EMAIL PROTECTED] wrote:
Dear R-people:
François Michonneau's method to obtain the special character plus/minus works
on Windows 2000 professional.
Many Thanks to François for his work!
Phil Smith
FYI:
1. That book has two authors.
2. That book is currently in 4th edition. The title of the fourth edition
is slightly different, to reflect the coverage of R. You are probably
reading an older edition (you didn't say), which does not cover R. In that
case you should consider `upgrading'.
You provided the code for LMEoptimize and lmerControl, but your
data.frame Semiconductor was not found. Therefore, I could not
replicate your problem. However, some external functions use their
inputs for temporary storage or even to return values. This was smart
programming
Gabor Grothendieck wrote:
Another way that works on my Windows XP system is to use \261 .
Please note Windows escape codes are not portable thus not recommended as
Martin Maechler pointed out in a response to a suggestion of mine:
On Tue, 14 Jun 2005, Martin Maechler wrote:
Peter == Peter
You did not provide a simple, self-contained replicable example, so I
can not say for sure. You say you have 49 observations in 7 blocks.
The issue described in the post you cite would be a problem if you had
49 blocks of size 1. I've tried to fit models like that, only to find
Hi Dr.Harrell,
thanks for the tip but I am not sure Quantile2() does what I'm need.
I need to repeatedly simulate 30 day survival data with a given sample
size and estimate a hazard ratio for each simulated data set. I want to
create a 95%CI for the hazard ratio by creating a distribution for it
Dear R Users,
I woould like to generate mixture model data.
I try two type method as followings two method, can anyone tell me which is
right? or better?
1. generate two vectors data y1 and y2 from f1 and f2 seperately, and calculate
y=alpha1*y1+alpha2*y2,
2. generate prob=unif(1), if prob
Hi there,
I use function kmeans and clara to cluster one flow
cytometry dataset. By using function plot, the clusters got
from clara can be graphed, while kmeans not. How can I
get the plot of the clusters of kmeans?
library(cluster)
data(xclara)
a1 - clara(xclara, 3)
a2 -
Keith Chamberlain wrote:
Dear Philippe, ( list)
I am not sure what is meant by 'floating window' but thank you for
clarifying the 'win' prefix. I had been considering different possible
senses of that prefix: win as in 'ts' window windowing functions (clearly
not), a generic 'win'dow
Hello,
Sometimes it is possible to copy-and-paste the R source code into your
project. Some R code is actually C code so you could create a dll or
something like that.
I did it for a Linux C++ project and it was not too dificult. You
should open MY_R_FUNCTION.R file and then you will see how R
88 matches
Mail list logo