At 03:42 PM 26/10/2005, you wrote:
We can double check this with yacas which is a free
computer algebra system that supports arbitrary
precision math:
In N(Exp(-1000))
Out 0.5075958898e-434
Or in Axiom, which is open source too:
numeric exp(-1000)
0.50759588975494567653E-434 Type:
Hi,I want to replace some elements of a dataframe a, using 2 index vectors
(x and y) and a vector of replacement (z) of the same length of x and
y.
I've tried a[cbind(x,y)]-z, but it doesn't work.How can I do it in a simply
way?Thanks,
Is
TAB3 - as.data.frame(with(TAB1, tapply(Number, Name, sum)))
colnames(TAB3) - c(Name_singular, Sum)
what you are looking for?
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Michael Graber
Sent: 25 October 2005 09:45 PM
To: R-Mailingliste
Subject:
On Wed, 26 Oct 2005, Marco Venanzi wrote:
Hi,I want to replace some elements of a dataframe aJ, using 2 index
vectors (x and y) and a vector of replacement (z) of the same length
of x and y.
I've tried a[cbind(x,y)]-z, but it doesn't work.How can I do it in a
simply way?Thanks,
Only
Ajit Chakrapani Warrier wrote:
pp - rpoispp(.01, win=owin(c(0,100),c(0,100)))
But how can I export the generated point pattern to an external file so
that I could use it as input for some network simulation
programs ?
You can get the coordinates of a point-pattern object
Hi,
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of lamack lamack
Dear all, someone could send me a introductory reference
about Survival
Frailty Models???
what about this article:
@article{vaupel79,
author={Vaupel, James W. and Manton, Kenneth G. and Stallard,
Dear R,
I´m really new to R, so it could be that my
question is one of those read the manual ones.
But I did that and found nothing.
Problem: I need to get a matrix of datapoints.
Each datapoint has to contain tree attributes.
trials like:
make.LDmatrix - function(nrbases){
b1 -
Hi all!
Does anyone knows if it exists a symbolic math package in R, that allows to
compute derivatives, integrals, etc.?
Does exist a freeware version of Maple?
Cheers,
Marco
[[alternative HTML version deleted]]
__
Dear all,
I use the COXPH function in the SURVIVAL package to estimate survival models
with a random effect.
I would like to know if there is a measure of model fit and complexity
analogous to the Akaike Information Criterion implemented in COXPH, in order
to compare together models including
Hi,
-Original Message-
[mailto:[EMAIL PROTECTED] On Behalf Of
Giannitrapani, Marco GSUK-GSSC
Subject: [R] symbolic math
Does anyone knows if it exists a symbolic math package in
R, that allows to compute derivatives, integrals, etc.?
Have a look at:
?D
Does exist a
Hi,
I have been using the partial.cor function in Rcmdr but I was wondering if
there is any easy way to get statistical significance tests (two tailed) along
with the partial correlation coefficients?
Simon Pickett
Simon Pickett
Centre for Ecology and Conservation Biology
University of Exeter
On Wednesday 26 October 2005 10:55, Giannitrapani, Marco GSUK-GSSC wrote:
Hi all!
Does anyone knows if it exists a symbolic math package in R, that allows
to compute derivatives, integrals, etc.?
Does exist a freeware version of Maple?
Cheers,
Marco
I've just learned yesterday the use
Dear All,
I have a cluster object based on a
dissimilarity matrix from about 1,100
cases and wish to know whether anyone
can think of any tips to display some
form of graphical output which would
give some sense of the similarity
between the cases.
A standard form of dendrogram would be
Good day to all,
Just for info if this can be of any use to R maintainers the update
package through the France (Paris) CRAN mirror does not seem to work.
France (Lyon) works fine though.
The Paris mirror used to work fine a while ago (can't remember exactly
when !).
Here is a copy of today's
-Original Message-
From: [EMAIL PROTECTED] [SMTP:[EMAIL PROTECTED] On Behalf Of Giannitrapani,
Marco GSUK-GSSC
Sent: Wednesday, October 26, 2005 7:56 AM
To: r-help@stat.math.ethz.ch
Subject: [R] symbolic math
Hi all!
Does anyone knows if it exists a symbolic math
Yet another option is axiom
In a debian-type linux:
apt-get install axiom
On Wed, 26 Oct 2005, Ruben Roa wrote:
-Original Message-
From: [EMAIL PROTECTED] [SMTP:[EMAIL PROTECTED] On Behalf Of Giannitrapani,
Marco GSUK-GSSC
Sent: Wednesday, October 26, 2005 7:56 AM
To:
Dear Simon,
The population partial correlation rho[12|3...p] is 0 when the regression
coefficient beta[2] for x[2] from the regression of x[1] on x[2] ... X[p] is
0. Thus, the usual t-test for a regression coefficient also tests that the
partial correlation is 0.
Now, the sample partial
Sorry for posting a possibly stupid question.
I am using aov() for calculating ANOVA as follows:
summary(aov(depVar~factor1*factor2+Error(subject/(factor1+factor2)),
data=anovaAllData))
and usually all works fine.
How exactly the call has to look if I want to use lm() directly to
obtain the
Hi all,
I've just installed R on my Mac PC.
The base demo run fine.. But I'm encountering several pb with some packages
(installed using CRAN binaries using the menu)
- Lattice:
demo(lattice,package='lattice')
demo(lattice)
~~~
Type Return to start :
require(grid)
[1]
Dear All,
How can I store hist() or any plot output graphics in
PNG format? I tried with bitmap() but not getting the
result. Appreciate your help on fixing this.
Thanks in Advance.
Regrads,
Ezhil
__
R-help@stat.math.ethz.ch mailing list
On 10/26/05 8:40 AM, A Ezhil [EMAIL PROTECTED] wrote:
Dear All,
How can I store hist() or any plot output graphics in
PNG format? I tried with bitmap() but not getting the
result. Appreciate your help on fixing this.
help.search('png')
Or
?png
Gets you the answer.
png('file.png')
Hi everone !
I know its not a epidemiologic or statistic list, but since I´m using R
to solve related problems ... here it goes:
This is a code to solve the rate problem, but the real question is:
what´s the correct way to calculate the period rate (3 years period put
together)
The third
?png
For example:
png(c:/myplot.png)
plot(rnorm(100))
dev.off()
A Ezhil wrote:
Dear All,
How can I store hist() or any plot output graphics in
PNG format? I tried with bitmap() but not getting the
result. Appreciate your help on fixing this.
Thanks in Advance.
Regrads,
Ezhil
Le 26.10.2005 14:40, A Ezhil a écrit :
Dear All,
How can I store hist() or any plot output graphics in
PNG format? I tried with bitmap() but not getting the
result. Appreciate your help on fixing this.
Thanks in Advance.
Regrads,
Ezhil
Hi,
You should try png then.
?png
BTW, There is a
r-help@stat.math.ethz.ch on Wednesday, October 26, 2005 at 6:00 AM -0500 wrote:
Ronaldo,
Try Harold's suggestion. The df still won't agree, because lmer (at least in
its current version) just puts an upper bound on the df. But that should be OK,
because all those t tests are approximations
I am trying to make horizontal violin plots. I have tried both vioplot
and simple.violinplot, but both of them seem to not be willing to take
the horizontal option. Is this correct, or am I just bungling it
somehow?
For instance, for vioplot (from the example shown, with the horizontal
Ronaldo,
According to Douglas Bates's paper in 'R' News, It would seem that the
correct model for nested split plot random effects with lmer , in your
example ,with x2 nested within x1, would be --
lmer(y~x1 + x2 +(1|x1)+(1|x1:x2))
Try it with your model any see how it compares with your aov
In addition to the response below, Doug Bates has talked about this on
this list previously. I did
RSiteSearch('bates degrees of freedom lmer')
The first one that came up has Doug's response to this question as well
Harold
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL
Hello everyone,
I wanted to know if it was possible to change the color of certain leaves in a
hclust object in order to make my graph more readable. I know I can color
certain groups but I would like to enter a vector telling the plot function
which leaves to color in which color.
Thanks in
Problem: I need to get a matrix of datapoints.
Each datapoint has to contain tree attributes.
Have you tried:
m - matrix(list(), nrbases, nrbases)
? You would then index it using m[[i,j]].A list is a basic
vector, so you can make a matrix with it, just as you can with a
vector
Hello all.
I have several time series of several variables.
One of them is the result, and others are believed to be the factors,
influencing the result.
There is a correlation between the result and an each factor.
Those factors, in turn, are correlated to each other.
I need to build a table,
Dear all,
I am using the spatstat package, in particular the Jest, Jdot and Jcross
functions.
When I plot the results using
J - Jest (SpatData)
plot.fv (J, main=Recruits)
there are 4 different lines (different colours and shapes):
black _
green
red _ _ _ _
blue _._._
Could
Le 26.10.2005 14:52, Karin Lagesen a écrit :
I am trying to make horizontal violin plots. I have tried both vioplot
and simple.violinplot, but both of them seem to not be willing to take
the horizontal option. Is this correct, or am I just bungling it
somehow?
For instance, for vioplot (from the
Hello. I am having a problem setting up a self-starting function for
use in nonlinear regression (and eventually in the mixed model version).
The function is a non-rectangular hyperbola - called NRhyperbola -
which is used for fitting leaf photosynthetic rate to light intensity.
It has one
Jorge,
If you'll send me details of the error messages, I'll see what I can do to
help. I notice in your posting there is a missing ')' in the next to last line
of model, but no error message; I don't suppose that could have anything to do
with it?
(send the details to my work email:
hy all,
I'm under linux,
I have now a big file (45 méga bytes), this file is a line of header defining
the names of the columns, the rows are only integers and there is 75 columns
(and 239096 rows included the header).
Before working with this big file i was working with another one (1.7 méga
Dear all,
I have ran a simulation in R.
This simulation was running about at least two days.
Here is below the result some part of my code about time result.
I don't understand about
Start time: Mon Oct 24, 2005 at 04:23:01 PM
Finish time: Wed Oct 26, 2005 at 03:26:19 PM
Run time:
Just a followup. I suppose you meant something like this:
library(odesolve)
y - c(10, 20, 10, 20)
parms - matrix(c(0.05, 0.1, 0.2, 0.05, 0.1, 0.2), nc=3, byrow=T)
model - function(times, y, parms) {
P - y[1:2]
V - y[3:4]
beta - parms[,1]
mu - parms[,2]
r - parms[,3]
dPdT -
Le 26.10.2005 15:17, william ritchie a écrit :
Hello everyone,
I wanted to know if it was possible to change the color of certain leaves in a
hclust object in order to make my graph more readable. I know I can color
certain groups but I would like to enter a vector telling the plot function
Hi All,
I am running R 2.2.0 on Mac OS 10.4.2, dual G5 processors with 8 Gig
RAM.
I am running a simulation with lsoda that requires ~378 s to complete
one set of time intervals. I need to optimize the parameters, and so
need to considerably speed up the simulation.
I have tried to figure
Hi
I installed R 2.2.0 from source and want to use install.packages but it
doesn't work.
http_proxy is set to http://proxy.sun.ac.za:3128
but it still can't connect to the repository.
The mirror is available, I can connect to it via the internet.
Any help welcome,
Rainer
--
NEW TELEPHONE
Hank,
I don't understand why you think memory is the problem, here. I'd try
writing my model in C or Fortran (there is an example in the odesolve
package). That speeds things up a lot, and is what I do with slow
systems.
R. Woodrow Setzer, Jr.
National Center for Computational Toxicology
US
Dear all,
I am an beginner with R and I have a question concerning apply(). My problem
is as follows:
I have data on four variables (x1,x2,x3,x4) with missing values for x1 and
x2.
data[1:9,]
x1 x2 x3 x4
1 NA NA 10 1
2 NA NA 8 3
3 NA NA 13 7
4 9 9 9 9
5 11 14 20 10
6 14 14 14
Bill Shipley [EMAIL PROTECTED] writes:
Hello. I am having a problem setting up a self-starting function for
use in nonlinear regression (and eventually in the mixed model version).
The function is a non-rectangular hyperbola - called NRhyperbola -
which is used for fitting leaf
Following on my posting of this morning, concerning a problem that I am
having constructing a self-starting function for use with nls (and
eventually with nlsList and nlme), the following is the self-starting
function called NRhyperbola:
NRhyperbola
function (Irr,theta,Am,alpha,Rd)
{
# Am
Hello,
I am trying to make the switch from SAS, and I have a fairly elemental
problem with syntax using the nlme package for analyzing mixed models.
There was a previous question on this topic posted to this list, so I
apologize for redundancy, but I didn't understand the advice given to
Hello list,
I would like to know if it is possible to get output of help(foo) with
Sweave.
If I insert ?foo in a chunk, the console give me the help but no output
is obtained in the tex file.
Perhaps something to do with options ?
Thanks in advance !
--
Stéphane DRAY ([EMAIL PROTECTED] )
Hello,
When using the plot(df) command, which scatter plots every variable of the
df against every other variable, I would like to add text labels, but cannot
find how to do this in this case so that every little sub-plot in the window
has text labels on the points...
Any help would be much
Le 26 Octobre 2005 12:48, Stéphane Dray a écrit :
Hello list,
I would like to know if it is possible to get output of help(foo) with
Sweave.
If I insert ?foo in a chunk, the console give me the help but no output
is obtained in the tex file.
Perhaps something to do with options ?
Thanks in
Those are obviously days, not seconds. A simple test would have
answered your question:
test - strptime(20051026 15:26:19,format=%Y%m%d %H:%M:%S) -
strptime(20051024 16:23:01,format=%Y%m%d %H:%M:%S)
class(test)
test
cat(test,\n)
If you prefer you can use difftime for conversion:
difftime
As usual, please read the docs (?plot.data.frame)! There you will find a
link to pairs() among whose arguments is panel for writing your own panel
function. In this panel function you can use low level plotting and text
functions like points() and text() to customize what is plotted. Look at the
1) You do need to load the packages before doing this. This is not
supposed to be needed.
2) So these are bugs in the packages, so please report them to the package
maintainers (and the posting guide suggest you contact them before
posting).
On Wed, 26 Oct 2005, Naji wrote:
Hi all,
I've
I am using the following code to get estimates of the variances of the
levels of factor 1 (factor1).
lme(sp1~factor1-1,data=z.a.1,random=pdBlocked(list(pdIdent(~factro1-1),pdIdent(~factor2-1))),
weights=varIdent(~1|factor1))
In this study factor 1 is crossed with factor 2. is this correct?
Dear R user,
I have a question on using R to analyze data with repeated measurements. I
have 2 species with several strains (12) per species, each of which has
been measured twice with for a given trait. No particular covariance, just
two measures. Now I want to analyze the data with an ANOVA
The word secs appears in Run time: 1.960625 secs because you put
it there in your cat() statement. It has nothing to do with the
number itself.
Simply try typing
end.time - begin.time
at the prompt, and see what you get.
Then see
?difftime
for more information. Example
Is there an R Package for Population Projections?
Orlando J. Rodriguez
The Center for Population Research
University of Connecticut
Unit 2068
344 Mansfield Road
Storrs, CT 06269-2068
(860)486-9269
(860)486-6356 (f)
http://popcenter.uconn.edu
__
--- Rodriguez, Orlando [EMAIL PROTECTED] wrote:
Is there an R Package for Population Projections?
Orlando J. Rodriguez
The Center for Population Research
University of Connecticut
Unit 2068
344 Mansfield Road
Storrs, CT 06269-2068
(860)486-9269
(860)486-6356 (f)
Martin,
For my simulations (same hardware with less memory).
- Optimize the code : using Narray, avoiding loops
- Installed distributed computations (PVM RPVM tools are very easy to use,
many thanks to Simon Urbanek for his help)
Hope this will help
Best
Naji
Le 26/10/05 16:43, « Martin Henry
As you have not yet received an **authoritative** answer, I'll attempt a
nonauthoritative one. Caveat emptor.
1. If strain is a fixed effect and you have but one exemplar of each strain
in species (which I assume is fixed) that is measured twice, measurement
error is all you have and a
Quick question, as I attempt to learn R. For post-hoc tests
1) Is there an easy function that will take, say the results of
tukeyHSD and create a grouping table. e.g., if I have treatments 1, 2,
and 3, with 1 and 2 being statistically the same and 3 being different
from both
Group
Depends what you want to do. Forming population projection matrices and
doing the eigenanalysis to work out population growth rates, stable stage
distributions, elasticities etc. is simple enough. Do you really need a
package to do that? See ?eigen. What would you like to see in a population
Jarrett Byrnes wrote:
Quick question, as I attempt to learn R. For post-hoc tests
1) Is there an easy function that will take, say the results of
tukeyHSD and create a grouping table. e.g., if I have treatments 1, 2,
and 3, with 1 and 2 being statistically the same and 3 being
Indeed, the following works as well
On Oct 26, 2005, at 5:23 PM, P Ehlers wrote:
fm1 - aov(breaks ~ wool*tension, data = warpbreaks)
TukeyHSD(fm1, c(wool,tension, wool:tension))
However, when working with my own dataset, I get the following errors.
I have some inkling this may be due to a
Dear List,
Is there a way to extract variance components from lmeObjects or
summary.lme objects without using intervals()? For my purposes I don't
need the confidence intervals which I'm obtaining using parametric
bootstrap.
Thanks,
Mike
[[alternative HTML version deleted]]
Looking at the errors your code produces, it looks like you need to make
Dock and Slip factors.
dock_2004_data$Dockf-factor(dock_2004_data$Dock)
dock_2004_data$Slipf-factor(dock_2004_data$Slip)
rich.aov - aov(X.open ~ Dockf*Slipf, data=dock_2004_data)
TukeyHSD(rich.aov, c(Dockf, Slipf))
Indeed, that does it. Odd. I guess as Slip was a number, it needed to
be categorized. Interesting...
Not if only I can figure out how to get the TukeyHSD grouping table
On Oct 26, 2005, at 7:35 PM, Mark Lyman wrote:
Looking at the errors your code produces, it looks like you need to
R-Help,
I am trying to reduce at data set to rows where a specified value occurs
in a specified value. Below is a screen capture of Rgui:
nsmalldata -read.csv(c:\\DATA\\UNITY\\\PASS0_DOWNFADE\\nsmall.csv)
nsmalldata
BOARDNUMBER SESSIONID MATRIXID ARRAYPOINT Temperature PS1 PS2
Hello Everybody,
I am running R 2.2.0 with Windows XP
i am trying to fit nonlinear differential equation to data sets which looks
like this:
Week N C
0 1 1
1 5 6
2 6.2 12.2
3 59 71.2
4 39 110.2
5 38 148.2
6 44 192.2
7 20.4 212.6
8 19.4 232
9 34.2 266.2
10 35.4 301.6
and i need to fit these data
One of the most useful functions:
?subset
as in
nsmSubData - subset(nsmalldata, SESSIONID==7757513)
On 10/26/05, Bill Hunsicker [EMAIL PROTECTED] wrote:
R-Help,
I am trying to reduce at data set to rows where a specified value occurs
in a specified value. Below is a screen capture of
Sorry for the repost, but I've really been looking, and can't find any
syntax direction on this issue...
Just browsing the documentation, and searching the list came up short... I
have some unbalanced data and was wondering if, in a 0 v 1
classification forest, some combo of these options might
See
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/40898.html
On 10/27/05, David L. Van Brunt, Ph.D. [EMAIL PROTECTED] wrote:
Sorry for the repost, but I've really been looking, and can't find any
syntax direction on this issue...
Just browsing the documentation, and searching the list came
Host (fixed)
Sire (random)
Dam nested within Sire (random)
Host * Sire (random)
Host * Dam within Sire (random)
So without the interactions I have:
hogmodel = lme(gain ~ host, random = ~1|sire/dam)
If I understand correctly, that sire/dam term gives me both
Sire and Dam within Sire
?VarCorr
At 12:02 PM 27/10/2005, you wrote:
Dear List,
Is there a way to extract variance components from lmeObjects or
summary.lme objects without using intervals()? For my purposes I don't
need the confidence intervals which I'm obtaining using parametric
bootstrap.
Thanks,
Mike
73 matches
Mail list logo