On Wed, 02 May 2007 11:37:22 +1000
Kate Stark [EMAIL PROTECTED] wrote:
[...] My model is:
fit - glm(Mature ~ Season * Size - 1, family = binomial, data=dat)
where Mature is a binary response, 0 for immature, 1 for mature. There
are 3 Seasons.
I would use:
fit - glm(Mature
Thanks for your code! It is not exactly what I really want - but it is my
fault, because my description was wrong...
It is not sim but rhater exp(rgamma(...)) that should not exceed 500. So
I tried to modify your code but it doesn't really work. sim.test returns just
1 value and not 999.
I've computed a loglinear model on a categorical dataset. I would like to
test whether an interaction can be dropped by comparing the log-likelihoods
from two models(the model with the interaction vs. the model without).
Since R does not immediately print the log-likelihood when I use the glm
Hello,
I've got a matrix (mail end) with the colnames x, y, z. In this matrix
are different measurements. x and y are risign coordinates.
With the following line I got the median value of z for all x AND y witch
are the same (not every measurment in my list hast the same number of
x and y
Please do not send everything twice: you are using R-help in both the To:
and Cc: fields.
I disagree with Ted: it _is_ much easier to create a generator for this
purpose.
Consider
rtgamma - function(n, ..., tr = log(500))
{
p - pgamma(tr, ...)
qgamma(p*runif(n), ...)
}
as
l'm trying to find the survdiff package/function but it seems i cant get it
how do l instal it if its not there
thanks
-
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing
I think you need to learn about deviances, which R does print.
Log-likelihoods are only defined up to additive constants. In this case
the conventional constant differs if you view this as a Poisson or as a
product-multinomial log-linear model, and R gives you the log-likelihood
for a Poisson
The version 2.5.0 has left Alpha status long time ago and its final
version has been released so please try the new version.
Inman, Brant A. M.D. wrote:
This email is intended to highlight 2 problems that I encountered
running R 2.5.0 alpha on a Windows XP machine.
#1 - Open script error
On 02-May-07 07:45:48, Prof Brian Ripley wrote:
Please do not send everything twice: you are using R-help in both the
To:
and Cc: fields.
I disagree with Ted: it _is_ much easier to create a generator for this
purpose.
Consider
rtgamma - function(n, ..., tr = log(500))
{
p
hello,
I have a problem with a R program I don't understand my errors
my program looks like this
for(i in LE) {
+for(j in LEC[[i]]) {
+for(k in LR) {
+ donRep[[k]] - subset(don2, Id_Essai == 1006961
Id_Cara == j Id_Rep == k Id_Geno != 65125,
You're right. I do need to learn more. I never learned null/residual
deviance. I know the deviance is equivalent to an anova decompostion. But
I've never dealt with it seperated like this.
I understand deviance as the difference between two model's log-likelihood
difference between them
--- raymond chiruka [EMAIL PROTECTED] wrote:
l'm trying to find the survdiff package/function
but it seems i cant get it how do l instal it if its
not there
thanks
I think you want the survival package.
__
R-help@stat.math.ethz.ch mailing list
hello,
I need your help for this example
for(k in LR) {
+ donGeno[[k]] - as.numeric(levels(factor(subset(don2, Id_Essai == 1006961
Id_Cara == LC[1] Id_Rep == k, select = Id_Geno)[,1])))
+ print(donGeno[[k]])}
[1] 65125 65126 65127 65128 65129 65130 65131 65132 65133 65134 65135 65136
65137
help.search(survdiff)
tells you it is in package survival. Most probably this package is not
missing, nevertheless if it is, you can install it in a standard way
using install.packages() or from the menu (Windows).
Petr
raymond chiruka napsal(a):
l'm trying to find the survdiff
elyakhlifi mustapha napsal(a):
hello,
I have a problem with a R program I don't understand my errors
my program looks like this
for(i in LE) {
+for(j in LEC[[i]]) {
+for(k in LR) {
+ donRep[[k]] - subset(don2, Id_Essai == 1006961
Id_Cara ==
(Ted Harding) wrote:
This thread prompts me to ask about something I've
been pondering for a while, as to whether there's an
implementation somewhere ticked away in the R resources.
So far, people have been responding to the original query
in terms of increasing the numbers of points, and
Alan E. Davis wrote:
I have been trying to visualize times of lowest tides, month by month.
I have tide predictions with times either in unix time or a text
format, and heights in feet or meters. I had been able to derive the
clock times of each prediction. I would now like to graph this
Alternatively generate the log-likelihood using the sum(dpois(y,
fitted(model), log = TRUE))
Regards
Ross Darnell
Doxastic wrote:
You're right. I do need to learn more. I never learned null/residual
deviance. I know the deviance is equivalent to an anova decompostion.
But I've
On 01/05/2007 11:34 AM, Prof Brian Ripley wrote:
On Tue, 1 May 2007, Duncan Murdoch wrote:
On 01/05/2007 12:29 AM, Graham Wideman wrote:
[...]
Refman p122: Environments consist of a frame, or collection of named
objects, and a pointer to an enclosing environment.
Is the or here explaining
Hello,
I've just carried out my first good-looking model using glmmPQL, and
the output makes perfect sense in terms of how it fits with our
hypothesis and the graphical representation of the data. However,
please could you clarify whether my degrees of freedom are
appropriate?
I had 106
Thanks. I used this and it gave me the same result as the logLik function.
The reason I ask is the SAS output gives me a loglik = 1089. R gives me
-298.09583. Both for my reduced model. For the saturated (or complex)
model, SAS gives me an loglik = 1143. R gives me -298.1993. The problem
hello,
I need help because I don't understand the syntaxe else how can I write it
for example I writed a script to cut missings values and I have errors
if(na==length(C)){
+ pos=match(0,match(donGeno[[na-1]],donGeno[[na]],nomatch=0))
+ for(k in 1:(na-1)) {
+
On 02/05/2007 8:11 AM, elyakhlifi mustapha wrote:
hello,
I need help because I don't understand the syntaxe else how can I write it
for example I writed a script to cut missings values and I have errors
if(na==length(C)){
+
Hi
[EMAIL PROTECTED] napsal dne 02.05.2007 14:11:51:
hello,
I need help because I don't understand the syntaxe else how can I
write it
for example I writed a script to cut missings values and I have errors
if(na==length(C)){
+ pos=match(0,match(donGeno[[na-1]],donGeno[[na]],nomatch=0))
Use ave. Also its easier to use read.table to read it in and then
convert it to a matrix if that's what you want.
Input - 29 4.5 1.505713
29 4.6 1.580402
29 4.7 1.656875
29 4.8 1.735054
30 0 0
30 0.1 0.00096108
30 0.2 0.00323831
29 4.5 1.495148
29 4.6 1.568961
29 4.7 1.644467
30 0 0
30 0.1
Hello,
I need a data set on repeated measures and recurrent events data which had been
registered on the same subject. This is very important to my thesis work.
Thanks a lot,
Luis Guillermo Díaz Monroy
Profesor Asociado
Universidad Nacional de Colombia
Hi Kate,
try looking at the package 'drc' on CRAN and in particular look at the
example in the help page for the dataset 'daphnids' (?daphnids).
You can obtain arbitrary ED values with approximate standard errors
using the function 'ED'.
Christian
I have some comment text taken from a SAS data file.
It is stored in two vectors and is difficult to read.
I would like to simply concatentate the individual
entries and end up with a character vector that give
me one line of text per comment.
I cannot see how to do this, yet it must be very
I work with Windows and use R version 2.4.1. I am JUST starting to learn this
program...
I get this warning message 'NAs introduced by coercion' while trying to
build a distance matrix (to be analyzed with NMDS later) from a 336 x 100
data matrix. The original matrix has lots of zeros and no
Hi,
I'm trying to evaluate a Monte Carlo p-value (using truncgof package) on
a left truncated sample.
From an empirical sample I've estimated a generalized pareto
distribution parameters (xi, beta, threshold) (I've used fExtremes pkg).
I'm in doubt on what of the following command is the most
Dear all,
After upgrading to R 2.5.0 on my Intel-based Mac I find that there
are only a few binary packages available. I've installed most of the
packages I need from source, but am having trouble compiling RGtk2.
Will Mac binary packages be available for R 2.5.0?
Thanks,
Ista Zahn
Apologies, I made a mistake with my maths. The degrees of freedom look
correct, assuming they are the denominator and that glmms work this
way. I just under-estimated the number of data points I had. Sorry.
Charlotte
On 02/05/07, Charlotte Burn [EMAIL PROTECTED] wrote:
Hello,
I've just
On 5/2/07, Jim Lemon [EMAIL PROTECTED] wrote:
Alan E. Davis wrote:
I have been trying to visualize times of lowest tides, month by month.
I have tide predictions with times either in unix time or a text
format, and heights in feet or meters. I had been able to derive the
clock times of
Is paste() what you're looking for?
Andy
From: John Kane
I have some comment text taken from a SAS data file.
It is stored in two vectors and is difficult to read.
I would like to simply concatentate the individual
entries and end up with a character vector that give
me one line of text
Kate Stark wrote:
Hi,
I was wondering if someone could please help me. I am doing a logistic
regression to compare size at maturity between 3 seasons. My model is:
fit - glm(Mature ~ Season * Size - 1, family = binomial, data=dat)
where Mature is a binary response, 0 for immature, 1 for
slomascolo wrote:
I work with Windows and use R version 2.4.1. I am JUST starting to learn this
program...
I get this warning message 'NAs introduced by coercion' while trying to
build a distance matrix (to be analyzed with NMDS later) from a 336 x 100
data matrix. The original matrix has
You may have a confusion between 64 bit and 32 bit versions. (I'm not
certain; I haven't studied your transcript in complete detail, but
rather just noticed a couple of things.)
There are a lot of references in your transcript to 64 bit libraries,
e.g. gcc -I/usr/lib64 ... but I think that
Hello and thank you for this first response
I'll rephrase what I mean by more general: in a case where I have only 2
levels in groups, but in which I have more than 2 levels for variable s, and in
a case where the second bar is not necessarily higher than the first.
Thanks in advance.
Best
It was suggested that the 'NAs introduced by coercion' message might be
warning me that my data are not what they should be. I checked this using
str(PeaksMatrix), as suggested, and the data seem to be what I thought they
were:
'data.frame': 335 obs. of 127 variables:
$ Code : Factor w/
At 07:30 AM 5/2/2007, Doxastic wrote:
Thanks. I used this and it gave me the same result as the logLik function.
The reason I ask is the SAS output gives me a loglik = 1089. R gives me
-298.09583. Both for my reduced model. For the saturated (or complex)
model, SAS gives me an loglik = 1143.
dear R wizards:
I am trying to replace subset() with my own version that first checks
that each name in the select statement has a corresponding name in the
data set. preferably, it would have the same syntax and semantics as
subset() otherwise.
alas, subset works in interesting ways:
Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having
to install all the packages again?
Thanks
Jason
Dr. Iasonas Lamprianou
Department of Education
The University of Manchester
Oxford Road, Manchester M13 9PL, UK
Tel. 0044 161 275 3485
[EMAIL PROTECTED]
Hello,
Thanks to the author for writing knnFinder. I seem to have problem
with a data set found here (http://www.stat.purdue.edu/~sguha/random/
mydata.csv)
This is the R-Code
ma-read.csv(~/mydata.csv)
ma-as.matrix(ma)
rownames(ma)-NULL
colnames(ma)-NULL
#ma is a 390x2 matrix
Hi
I am trying to use RODBC in R installed on Windows to
access MySQL database (on a linux box).
I set up a DSN and specified this DSN in R as follows
library(RODBC);
channel - odbcConnect(mysqldsn);
RODB Connection 5
Details:
case=nochange
PORT=3306
Although this seems to connect
ivo welch said the following on 5/2/2007 8:13 AM:
dear R wizards:
I am trying to replace subset() with my own version that first checks
that each name in the select statement has a corresponding name in the
data set. preferably, it would have the same syntax and semantics as
subset()
Silvia Lomascolo wrote:
It was suggested that the 'NAs introduced by coercion' message might be
warning me that my data are not what they should be. I checked this using
str(PeaksMatrix), as suggested, and the data seem to be what I thought they
were:
'data.frame': 335 obs. of 127
Hi all R positive,
does anyone know how to refer R in article?
thanks
tomas
__
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
and
Dear all,
I just realized that I forgotten to write some kind of final email for
this thread and to thank you for your help.
It seems that the recommeneded procedure in such circumstances has three
steps:
1) readLines()
2) select the desired lines
3) strsplit()
Thanks Ferdinand, Jim, and Paul!
Hi all,
I have written a R-script under Windows using choose.dir. Now, I have
seen that this function is missing at MacOS. Does anybody know an
alternative?
Antje
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Thank you for your answer,
But I think it doesn't work, because the function bargraph.CI (I use it
because barplot doesn't have a standard error option) use two internally
function(one for mean and the other for standard error) which use something
like tapply (which sort the factors by alphabet).
citation()
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Tomas Mikoviny
Sent: Wednesday, May 02, 2007 11:44 AM
To: r-help@stat.math.ethz.ch
Subject: [R] reference in article
Hi all R positive,
does anyone know how to refer R in article?
try
citation()
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
Hi
I have a dataframe which has 3 columns of numeric data
A,B,C each of which has been obtained independent of
the other.
We are trying to find out, which of A or B cause C
i.e. We are hypothesising that C is the effect and
either A or B, not both is the cause.
i.e. A causes C and this
Beautiful. thank you. /ivo
Try:
methods(subset)
which will point you to subset.data.frame. The latter code will answer
your questions.
HTH,
--sundar
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
On Wed, 2 May 2007, Tomas Mikoviny wrote:
Hi all R positive,
does anyone know how to refer R in article?
Every time you start R it says (in part)
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages
Hi all,
I have written a R-script under Windows using choose.dir. Now, I have
seen that this function is missing at MacOS. Does anybody know an
alternative?
Antje
-
Kennt man wirklich jeden über 3 Ecken? Die Antworten gibt's bei Yahoo! Clever.
Luis Guillermo Díaz Monroy
Profesor Asociado
Universidad Nacional de Colombia---BeginMessage---
Hello,
I need a data set on repeated measures and recurrent events data which had been
registered on the same subject. This is very important to my thesis work.
Thanks a lot,
Luis Guillermo Díaz
I typically wouldn't use strsplit but would reread it using read.table
and textConnection as in:
http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg84752.html
On 5/2/07, Roland Rau [EMAIL PROTECTED] wrote:
Dear all,
I just realized that I forgotten to write some kind of final email for
I don't have a Mac but perhaps this works:
library(tcltk)
tclvalue(tkchooseDirectory())
On 5/2/07, Antje [EMAIL PROTECTED] wrote:
Hi all,
I have written a R-script under Windows using choose.dir. Now, I have
seen that this function is missing at MacOS. Does anybody know an
alternative?
Romain,
Try this before calling bargraph.CI():
Fst2$REGION - reorder(Fst2$REGION, Fst2$MOY_FST, mean)
This will reorder the factor levels for REGION, based upon the mean
value of MOY_FST for each region. Thus, your REGIONS should then be
sorted by increasing mean values in the plot.
For
Hi
This is not a homework assignment :)
Me and my manager are trying to understand the problem better. In the
meanwhile, we thought we would post the problem on this forum to seek some
input from statisticians who possibly do this kind of analyses everyday and
hence are possibly more proficient
I am sorry if I misinterpreted this but this questions looks/looked very
like a very basic one about statistical inference... (have a look at the
Manuals/contributed documentation e.g. Faraway)
Lalitha Viswanath wrote:
Hi
This is not a homework assignment :)
Me and my manager are trying to
Tomas Mikoviny wrote:
Hi all R positive,
does anyone know how to refer R in article?
thanks
tomas
Like this (well, almost! This is from R 2.4.1 - 2007 on the most recent
version)
citation()
To cite R in publications use:
R Development Core Team (2006). R: A language and
Hi,
I have to plot three Ziph distributions for three languages where the
x value represents the rank of a given word and the y value represents
the relative frequency of this word in the corpus. Is there some way
so that I can plot all three distributions on a single scatterplot,
preferably with
How about making your homeworks yourselfes?
lalitha viswanath wrote:
Hi
I have a dataframe which has 3 columns of numeric data
A,B,C each of which has been obtained independent of
the other.
We are trying to find out, which of A or B cause C
i.e. We are hypothesising that C is the effect
Hi,
I have to plot three Ziph distributions for three languages where the
x value represents the rank of a given word and the y value represents
the relative frequency of this word in the corpus. Is there some way
so that I can plot all three distributions on a single scatterplot,
preferably with
Iasonas Lamprianou said the following on 5/2/2007 8:25 AM:
Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without
having to install all the packages again?
Thanks
Jason
You may find the following link relevant.
Thank you. For Windows it works, I'll check whether it also works for Mac :)
Gabor Grothendieck schrieb:
I don't have a Mac but perhaps this works:
library(tcltk)
tclvalue(tkchooseDirectory())
On 5/2/07, Antje [EMAIL PROTECTED] wrote:
Hi all,
I have written a R-script under Windows
Lalitha Viswanath wrote:
We are trying to find out, which of A or B cause C
i.e. We are hypothesising that C is the effect and
either A or B, not both is the cause.
(...)
I would greatly appreciate any inputs on the best
statistcal approach to tackle this problem.
I am thinking that we
On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Iasonas Lamprianou said the following on 5/2/2007 8:25 AM:
Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without
having to install all the packages again?
Thanks
Jason
You may find the following link relevant.
Hi
Thanks for your input. I stand corrected.
The causation is not linear.
We wish to find out which is the cause and under what circumstances. i.e. at
what points along a scale for C, is A the cause and when does B become the
cause, if at all.
As a crude analyses, we assumed that the above is not
Your title and your posting do not say the same thing.
Assuming you want all three distributions on one
scatter plot does this help?
aa - 1:10
bb - 11:2
cc - bb^2
dd - c(3,4,7,9,11,32,11,14,5,9)
plot(aa,cc, col=red)
points(aa,bb, col=blue)
points(aa,dd, col=green)
Also in plotting it is a
Lisa: Can we expect to see R used [exclusively, I would hope] during
your demonstration? Learning how data mining models work: the inputs,
the outputs, and the nature of the predictive mechanism only makes
sense for me if I can follow/retrace your steps on my systems. Thank
you.
Brian J. Koch
Hello,
I have just upgraded from R-2.4.1 to R-2.5.0 for Windows.
I had installed a large number of add-on packages under 2.4.1. Is there an easy
way to install (or load, if that's the easier way) those packages under 2.5.0,
without having to install each package by hand?
Thanks,
Jeff
Hello Andy,
I had the same problem in Feisty and now I fixed it.
I updated to Java 6 through the repositories from Feisty.
Then I choosed java 6 with sudo update-alternatives --config java.
After that I started R with sudo R.
If you had already installed JGR, but it doesn't run you have to
At 01:41 PM 5/2/2007, you wrote:
On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Iasonas Lamprianou said the following on 5/2/2007 8:25 AM:
Hi I am using R version 2.4.1. How can I upgrade to version 2.5
without having to install all the packages again?
Thanks
Jason
On 02/05/2007 5:12 PM, Jeffrey Wood wrote:
Hello,
I have just upgraded from R-2.4.1 to R-2.5.0 for Windows.
I had installed a large number of add-on packages under 2.4.1. Is there an
easy way to install (or load, if that's the easier way) those packages under
2.5.0, without having
Robert A LaBudde said the following on 5/2/2007 2:39 PM:
At 01:41 PM 5/2/2007, you wrote:
On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Iasonas Lamprianou said the following on 5/2/2007 8:25 AM:
Hi I am using R version 2.4.1. How can I upgrade to version 2.5
without having to
Rename the library subdirectory in your new installation to something
like library.new. Copy the library subdirectory in the old
distribution to the new installation. Copy the contents of the
renamed library.new to the the directory copied from the old
distribution. You should then start your
Please note: in R you can specify (some of the) graphics parameters as the
appropriate length vectors. So your plot example below can also be done as,
for example:
plot(
rep.int(aa,3),c(cc,bb,dd),col=rep(c(red,blue,green),e=length(aa)))
However, this doesn't seem to fit the posted request, where
Hi All,
I found mvfft in R and fft2 in Matlab give different result
and can't figure out why. My example is:
In R:
matrix(c(1,4,2,20), nrow=2)
[,1] [,2]
[1,]12
[2,]4 20
mvfft(matrix(c(1,4,2,20), nrow=2))
[,1] [,2]
[1,] 5+0i 22+0i
[2,] -3+0i -18+0i
In Matlab:
Hi all,
Is it possible to generate the same sample number of times in R? In
SAS, using the option seed it is possible to reproduce exactly the
same sample. Is there any such feature in R which I can use?
For further clarity,
for (i in 1:2)
{
samp = sample(1:1000,100,replace = FALSE)
On 02-May-07 23:25:18, Santanu Pramanik wrote:
Hi all,
Is it possible to generate the same sample number of times in R?
In SAS, using the option seed it is possible to reproduce exactly
the same sample. Is there any such feature in R which I can use?
For further clarity,
for (i in
On 5/2/07, Robert A LaBudde [EMAIL PROTECTED] wrote:
At 01:41 PM 5/2/2007, you wrote:
On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Iasonas Lamprianou said the following on 5/2/2007 8:25 AM:
Hi I am using R version 2.4.1. How can I upgrade to version 2.5
without having to
On 5/2/07, Jeffrey Wood [EMAIL PROTECTED] wrote:
Hello,
I have just upgraded from R-2.4.1 to R-2.5.0 for Windows.
I had installed a large number of add-on packages under 2.4.1. Is there an
easy way to install (or load, if that's the easier way) those packages under
2.5.0, without
Li Li said the following on 5/2/2007 4:06 PM:
Hi All,
I found mvfft in R and fft2 in Matlab give different result
and can't figure out why. My example is:
In R:
matrix(c(1,4,2,20), nrow=2)
[,1] [,2]
[1,]12
[2,]4 20
mvfft(matrix(c(1,4,2,20), nrow=2))
[,1]
Hi,Dear R users, I have a file text nommed chif which contains 16 lines and 4
columns in the disc dur. I have a difficulty to read this file in R console I
have used the following command
chif - read.table(c:/chif.txt, header=T, sep= )
I have obtained from R console:undefined file! Can you
Hi all
I'm updating a package for submission to CRAN, and I am getting
an error message I never have seen before, and can't find out
what it implies or means (and hence, what I need to do to fix
it).
After running R CMD check, I get this:
snip
* checking foreign function calls ... OK
* checking
?set.seed
On Wed, 02 May 2007 18:25:18 -0500, Santanu Pramanik
[EMAIL PROTECTED] wrote:
Hi all,
Is it possible to generate the same sample number of times in R? In
SAS, using the option seed it is possible to reproduce exactly the
same sample. Is there any such feature in R which I can
We don't have any of the basic details the posting guide asks for!
If this is R 2.5.0, rcompgen is a recommended package and ought to be
installed. If it is not, then .check_Rd_refs cannot cross-check
references against all the standard and recommended packages, and that
looks like the cause
Peter,
On 3 May 2007 at 10:43, Peter Dunn wrote:
| After running R CMD check, I get this:
|
| snip
| * checking foreign function calls ... OK
| * checking R code for possible problems ... OK
| * checking Rd files ... OK
| * checking Rd cross-references ... WARNING
| Error in
Dirk
On 3 May 2007 at 10:43, Peter Dunn wrote:
| After running R CMD check, I get this:
|
| snip
| * checking foreign function calls ... OK
| * checking R code for possible problems ... OK
| * checking Rd files ... OK
| * checking Rd cross-references ... WARNING
| Error in
Hello all!
I am once again analyzing patient survival data with chronic liver disease.
The severity of the liver disease is given by a number which is continuously
variable. I have referred to this number as meld--model for end stage
liver disease--which is the result of a mathematical
Peter,
On 3 May 2007 at 12:23, Peter Dunn wrote:
| Are you by chance on Debian (or Ubuntu)?
|
| Bingo.
|
| The short answer seems to be that I do need rcompgen to be
| installed, and my Debian system didn't do that. Easy fixed though.
Yes, sorry. I realized relatively late that we needed
Thanks for both replies.
Then I found the ifft2 from Matlab gives different result from fft( ,
inverse=T) from R.
An example:
in R:
temp - matrix(c(1,4,2, 20), nrow=2)
fft(temp)
[,1] [,2]
[1,] 27+0i -17+0i
[2,] -21+0i 15+0i
fft(temp,inverse=T)
[,1] [,2]
[1,] 27+0i -17+0i
I should clarify. I can generate plots for each category individually but
not for all three on the same chart.
Greg
-Original Message-
From: Gregory Pierce [mailto:[EMAIL PROTECTED]
Sent: Wednesday, May 02, 2007 10:21 PM
To: 'r-help@stat.math.ethz.ch'
Subject: Survival
Discrete Fourier transforms can be normalized in different ways.
Some apply the whole normalization to the forward transform, some to
the reverse transform, some apply the square root to each, and some
don't normalize at all (in which case the reverse of the forward
transform will need
? layout()
? par()
E.g.,
layout(matrix(c(1,2,3),1,3,byrow=TRUE) #3 plots side-by-side
Then use plot() three times to generate each of your graphs.
At 11:14 PM 5/2/2007, Greg wrote:
I should clarify. I can generate plots for each category individually but
not for all three on the same chart.
Li Li said the following on 5/2/2007 7:53 PM:
Thanks for both replies.
Then I found the ifft2 from Matlab gives different result from fft( ,
inverse=T) from R.
An example:
in R:
temp - matrix(c(1,4,2, 20), nrow=2)
fft(temp)
[,1] [,2]
[1,] 27+0i -17+0i
[2,] -21+0i 15+0i
Have a look at
?lines
?points
?plot - the option add, add=TRUE
You will have to specify proper limits on both axes, otherwise you may
only see parts of the graphs. Some functions in the survival library
also allow stratified analyses in which case the plots account for
different strata.
Petr
1 - 100 of 110 matches
Mail list logo