Thank you very much, Professor Ripley.
If possible, could you point me to other packages that you think I
should look at for estimating a derivative?
Best regards,
Hui
Prof Brian Ripley wrote:
On Tue, 10 May 2005, Hui Han wrote:
I wonder if the function density outputs the gaussian mixture
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/20509.html
-s.
Hui Han wrote:
Thank you very much, Professor Ripley.
If possible, could you point me to other packages that you think I
should look at for estimating a derivative?
Best regards,
Hui
Prof Brian Ripley wrote:
On Tue, 10 May 2005, Hui
Dear useRs,
I have just rendered 3 wallpapers with R logo:
ftp://ftp.ariadna.pl/pub/R-walls/
Hope you will enjoy. Regards,
--
Lukasz Komsta
Department of Medicinal Chemistry
Medical University of Lublin
6 Chodzki, 20-093 Lublin, Poland
Fax +48 81 7425165
Lapointe, Pierre wrote:
This is excellent. It works great even though I'm on Win2k and behind a
firewall.
I just have one question though, I don't use a function because I sometimes
need to add people to my distribtion list. In R, my system(paste()) string
is very long and when I try put an enter
Darren Weber wrote:
Has anyone compiled a binary for winXP, with python2.4 and R.2.1.0?
If anyone can do this soon, please advise how to pick up the binary
installer.
The Rpy homepage is http://rpy.sourceforge.net/
What about asking on the mailing list given on that page?
Uwe Ligges
Best,
Darren Weber wrote:
Dear Prof. Baron,
thankyou for your pdf paper on using R, at
http://www.psych.upenn.edu/~baron/rpsych.pdf
I am currently interested in the sphericity notes on pp. 45-47.
The equation on p. 46 suggests the following R calculation:
D - k^2 * ( mean(diag(S)) - mean(S) )^2
while
Dear List members
I am using the regsubsets function to select a few predictor variables
using Mallow's Cp:
sel.proc.regsub.full - regsubsets(CO2 ~ v + log(v) + v.max + sd.v +
tad + no.stops.km + av.stop.T + a + sd.a + a.max + d + sd.d + d.max +
RPA + P + perc.stop.T + perc.a.T + perc.d.T +
Thibaut Jombart wrote:
Hello,
I'm trying to insert R source code (functions) in an appendix of a latex
document. I guess the easiest way to do so is to use the package Sweaved
(file : Sweaved.sty) provided with the latest R version. Latex succeeds
in loading the package, but my problem comes
Dear list,
I'm making my current code more generic and would like some advise.
The basic problem is subset and the name of the column to be compared
for selection.
What I've come up with is
data(mammals)
set - bottompremolars
subset(mammals, eval(parse(file=,text=set)) 2)
This seems a bit
Fredrik Karlsson wrote:
Dear list,
I'm making my current code more generic and would like some advise.
The basic problem is subset and the name of the column to be compared
for selection.
What I've come up with is
data(mammals)
set - bottompremolars
The line above is
a) syntactically incorrect
On Wed, May 11, 2005 at 10:09:07AM +0200, Uwe Ligges wrote:
Thibaut Jombart wrote:
Hello,
I'm trying to insert R source code (functions) in an appendix of a latex
document. I guess the easiest way to do so is to use the package Sweaved
(file : Sweaved.sty) provided with the latest R
On Wed, 11 May 2005, Fredrik Karlsson wrote:
Dear list,
I'm making my current code more generic and would like some advise.
The basic problem is subset and the name of the column to be compared
for selection.
What I've come up with is
data(mammals)
set - bottompremolars
subset(mammals,
Thibaut Jombart wrote:
I'm trying to insert R source code (functions) in an appendix of a
latex document. I guess the easiest way to do so is to use the package
Sweaved (file : Sweaved.sty)
I have had very good results with the LaTeX package listings. It is
easy to use, and it does the job.
Le 11.05.2005 14:00, Michael Prager a écrit :
Thibaut Jombart wrote:
I'm trying to insert R source code (functions) in an appendix of a
latex document. I guess the easiest way to do so is to use the
package Sweaved (file : Sweaved.sty)
I have had very good results with the LaTeX package
hello,
I would like to compare indexes in order to see which clustering algorithm (k-
means, k-medoid, hierarchical clustering with euclidean and pearson distance)
is more efficient. So I want to calculate the c-index and the db in all those
cases(library cclust).
To do that I use the cclust
Hey,
i'd like to compose a clock-like looking plot composed out of two
circles, each showing the length of a period (to compare them). first,
to do so, it looked the easiest by using pie(), just puting multiple
pie-charts over each other. the problem is that once the second pie is
drawn, it
Le 11.05.2005 15:43, Lars a écrit :
Hey,
i'd like to compose a clock-like looking plot composed out of two
circles, each showing the length of a period (to compare them). first,
to do so, it looked the easiest by using pie(), just puting multiple
pie-charts over each other. the problem is that
hec.data -array(c(5,15,20,68,29,54,84,119,14,14,17,26,16,10,94,7),
dim=c(4,4),
dimnames=list(eye=c(Green,Hazel, Blue, Brown),
hair=c(Black, Brown, Red, Blond)))
#--
dfr -
Have you considered matplot:
matplot(1:3, array(1:6, dim=c(3,2)), type=b)
spencer graves
Xiao Shi wrote:
I have a matrix gene=array(rnorm(180), dim=c(18,10)) and a vector time=c(0,
0.5,2,4,6,8,12,24,48,72).
So i can get a plot with plot(time,gene[1,],type=b,col=green) for the
first row.So how can
See ?matplot. E.g.,
matplot(time, t(gene), type=b, lty=1)
Andy
From: Xiao Shi
I have a matrix gene=array(rnorm(180), dim=c(18,10)) and a
vector time=c(0,
0.5,2,4,6,8,12,24,48,72).
So i can get a plot with
plot(time,gene[1,],type=b,col=green) for the
first row.So how can i get all
Le 11.05.2005 17:26, Xiao Shi a écrit :
I have a matrix gene=array(rnorm(180), dim=c(18,10)) and a vector time=c(0,
0.5,2,4,6,8,12,24,48,72).
So i can get a plot with plot(time,gene[1,],type=b,col=green) for the
first row.So how can i get all the 18 rows in the same plot VS. the same
vector
On Wed, 11 May 2005, Neuro LeSuperHéros wrote:
Do give your real name and affiliation! What exactly are you trying to
hide?
Why don't you ask on a PageMaker list how to get it to accept valid PDF?
This is a bug in PageMaker, not in R.
I can generate PDF documents and they look good on acrobat
Dear all,
I am still trying to fit a model for habituation with lme (Dogs are
trotting on a treadmill once a week during 4 consecutive weeks).
I have maybe found a good one.
In order to validate it, I would like to split my data.
I have 28 dogs, I am thinking to build the model with 20 and use
Is there a way to implement this faster than doing it in a loop.
for (i in length(settle)-1:1) {
settle[i] = settle[i+1]/(1 + settle.pct[i+1])
}
I want to guarantee that i+1 is calculated before i
__
On Wed, 11 May 2005 13:05:58 -0400 Omar Lakkis wrote:
Is there a way to implement this faster than doing it in a loop.
for (i in length(settle)-1:1) {
settle[i] = settle[i+1]/(1 + settle.pct[i+1])
}
I want to guarantee that i+1 is calculated before i
It is, backwards, a cumulative product so you could use cumprod.
On Wed, 11 May 2005, Omar Lakkis wrote:
Is there a way to implement this faster than doing it in a loop.
for (i in length(settle)-1:1) {
settle[i] = settle[i+1]/(1 + settle.pct[i+1])
}
I want to guarantee
Omar Lakkis wrote:
Is there a way to implement this faster than doing it in a loop.
for (i in length(settle)-1:1) {
settle[i] = settle[i+1]/(1 + settle.pct[i+1])
}
You dont need a loop at all here. How so? Well, as it is written the
code in the for loop only
Hi all!
__
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
Hello,
I would like to remove the frame from wireframe, but that one only.
I am currently doing
trellis.par.set(axis.line,list(col=NA,lty=1,lwd=1))
but this has the flaw that also the tick lines of the axes disapper.
I read all the thread with the same titlee of this message appeared
Hi,
The recent thread on how to send mail from R got me thinking about a
solution to clearing the rgui console under WinXP (something I've often
wanted to do and I've seen others ask about). If you create a small vb
script, say called RClear.vbs:
Dim wsh
Set wsh = CreateObject(WScript.Shell)
If you used all 28 animals to find the model out of a group of candidate
models, I would have my reservations about this 'validation'. Any
confidence intervals you get from the final model are bound to be overly
optimistic because you haven't accounted for the degrees of freedom chewed
up during
On Wednesday 11 May 2005 12:35, Marco Chiarandini wrote:
Hello,
I would like to remove the frame from wireframe, but that one only.
I am currently doing
trellis.par.set(axis.line,list(col=NA,lty=1,lwd=1))
but this has the flaw that also the tick lines of the axes disapper.
I should add
Hello,
I'm having a error message when I try to load the som library from a
cgi. When a run my script via web the following error appears in the
apache log:
Error: package som was built for i686-pc-linux-gnu, referer: https://
I call the library function setting the path to the library
You could use the chron package. It represents date times without
using time zones so you can't have this sort of problem.
On 5/10/05, Carla Meurk [EMAIL PROTECTED] wrote:
Hi, I have a whole bunch of data, which looks like:
15/03/2003 10:20 1
15/03/2003 10:21 0
15/03/2003
Thank you so much, Suresh. I searched a lot on density among R email
archives. Should have searched using derivative.
Hui
Suresh Krishna wrote:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/20509.html
-s.
Hui Han wrote:
Thank you very much, Professor Ripley.
If possible, could you point me to
I have a data frame of daily open, high, low and settle prices. How
can I aggregate this data weekly?
The data frame has five columns, the first is the date column and the
rest are the prices.
__
R-help@stat.math.ethz.ch mailing list
Assuming dfr[day,o,h,l,c] and day like 2004-12-28:
dt - strptime(as.character(dfr$day),format=%Y-%m-%d) + 0
wk - format(dt,%Yw%U)
aggr - aggregate(list(dfr$o,dfr$h,dfr$l,dfr$c),list(wk),mean)
colnames(aggr) - etc
-Original Message-
From: Omar Lakkis [mailto:[EMAIL PROTECTED]
Sent:
In fact since you have dates and not datetimes use as.Date() instead
of strptime().
On 5/11/05, bogdan romocea wrote:
Assuming dfr[day,o,h,l,c] and day like 2004-12-28:
dt - strptime(as.character(dfr$day),format=%Y-%m-%d) + 0
wk - format(dt,%Yw%U)
aggr -
I have a list in a package environment
assign('refflat',list(),pos='package:locPkg')
to which I would like to make assignments like:
refflat[['a']] - read.table('fileA.txt')
refflat[['b']] - read.table('fileB.txt')
I am doing this to guard against a local version of refflat hanging
around, as I
With respect to calculating the epsilon index of sphericity for ANOVA,
discussed on pp. 45-47 of:
http://www.psych.upenn.edu/~baron/rpsych.pdf
It notes that epsilon is not required for a repeated measures design with
only k=2 levels, as the minimum value of epsilon (e) is given by:
e =
charles loboz [EMAIL PROTECTED] wrote:
A gstring is a string with variable names embedded and replaced by
values(converted to strings, lazy eval) before use.
I use the following function, which will take variables either from
named arguments or from the environment. It also concatenates all
Lars wrote:
Hey,
i'd like to compose a clock-like looking plot composed out of two
circles, each showing the length of a period (to compare them). first,
to do so, it looked the easiest by using pie(), just puting multiple
pie-charts over each other. the problem is that once the second pie is
Dear All,
I have some code that works in S-Plus for writing saving a graphics file to
disk :-
graphsheet(type = auto, format = WMF, file = G:\\north0l.wmf,
pages = auto, print.background = F,
orientation=landscape,
color.style=color)
plot(x,y)
dev.off()
This works fine in
Hi Guy
Try savePlot(MyPlot, bmp).
Cheers
Francisco
PS: Hang in there! In the long run the effort to move from S-Plus to R is
definitively worthy!
From: Guy Forrester [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Subject: [R] Graphics file to disk
Date: Thu, 12 May 2005 11:08:17 +1200
Dear All,
Many thanks to all who replied,
In no particular order:-
Tom Mulholland, Gunter Berton, Andrew Ward, Francisco Zagmutt, N. Olsen
all is well
Cheers
Guy
Guy J Forrester
Biometrician
Manaaki Whenua - Landcare Research
PO Box 69,
If you want to write directly to file,
?Devices
will get you more information
Guy Forrester wrote on 5/11/2005 7:08 PM:
Dear All,
I have some code that works in S-Plus for writing saving a graphics file to
disk :-
[...]
--
Michael H. Prager, Ph.D.
Population Dynamics Team
NOAA Center for Coastal
are you looking for something like:
InternalMean - mean(data1[,3])
TestValues - c(0,1,2,4,9) #should be in increasing order
TestResults - c(.95, .85, .7, NaN,0)
if (InternalMean==0) IntResult=1 else
IntResult=TestResults[which(TestValues==max(TestValues[TestValuesInternalMean]))]
-s.
Jones, Glen
oops, i meant something more like:
TestValues - c(0,1,2,4,9) #should be in increasing order
TestResults - c(.95, .85, .7, NaN,0)
if (InternalMean==0) IntResult=1 else
IntResult=TestResults[TestValues==max(TestValues[TestValues=InternalMean])]
-s.
Suresh Krishna wrote:
are you looking for
Darren is of course correct, but I hope the following, brief,
intentionally non-technical explanation will help:
Repeated measures analyses are needed only when there are three or more
measurements from a given experimental unit and the two or more
measurements are both on the right hand side of
Would it not just be easier to set the timezone to GMT for the duration of
the calculations? I don't see an OS mentioned here, but on most TZ=GMT
for the session will do it.
On Thu, 12 May 2005, Rich FitzJohn wrote:
Hi,
seq.dates() in the chron package does not allow creating sequences by
50 matches
Mail list logo