This is a question regarding the 'regTermTest' function in the 'survey' package.
Imagine Z as a three level factor variable, and code ZB and ZC as the two
corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z and
X, say, how do the two test specifications
Alternate title: How can I persuade my students that R is for them?
Alternate title: Can R replace SAS, SPSS or Stata for clinicians?
I am teaching introductory statistics to twelve physicians and two veterinarians
who have enrolled in a Mentored Clinical Research Training Program. My course
glmmPQL calls lme, and it is that which differs between S-PLUS and R.
These are optimization problems with multiple local maxima, and like any
complex statistical fitting problem you should not expect all programs to
give the same answer. The short answer is to believe all of them, and to
Thanks Mike.
My data has longitude and latitude coords and I used distAB {clim.pact}
then mst {ape} to calculate my minimum spanning tree. The nodes are
telecoms sites from all over Australia. My goal is to determine the
minimum cost of linking them via cabling, and I'm starting by
Thanks Mike.
My data has longitude and latitude coords and I used distAB {clim.pact}
then mst {ape} to calculate my minimum spanning tree. The nodes are
telecoms sites from all over Australia. My goal is to determine the
minimum cost of linking them via cabling, and I'm starting by
Hi Dave,
This has come up a couple of times recently.
I don't remember by heart but searching the archives should help you find
the answer. It's best to try R-SIG-MAC list archive.
Bye, ingmar
On 8/19/04 4:34 AM, David L. Van Brunt, Ph.D. [EMAIL PROTECTED]
wrote:
It's been a while since I used
On Thu, 2004-08-19 at 07:45, Jacob Wegelin wrote:
Alternate title: How can I persuade my students that R is for them?
Alternate title: Can R replace SAS, SPSS or Stata for clinicians?
I am teaching introductory statistics to twelve physicians and two veterinarians
who have enrolled in a
i want to upload different sheets from one excel file.
does anyone know of a way to do so? (i use read.csv to load the excel but
i would like to load form the sheets themselves.)
thanks
sivan
__
[EMAIL PROTECTED] mailing list
One thought -- the first course I took, I taught physicians
XLispStat. They appreciated it since it allowed them to do logistic
regression without spending large $$$.
Now, that isn't quite true any more. So might depend on how times
have changed.
Jacob Wegelin [EMAIL PROTECTED] writes:
Sivan Aldor wrote:
i want to upload different sheets from one excel file.
does anyone know of a way to do so?
Yes, the R Data Import/Export manual knows.
Uwe Ligges
(i use read.csv to load the excel but
i would like to load form the sheets themselves.)
thanks
sivan
-Ursprüngliche Nachricht-
Von: Wildi Marc, wia
Gesendet: Mittwoch, 18. August 2004 10:11
An: [EMAIL PROTECTED]
Betreff:
Hi
Does anybody know from an R-package devoted to sample selection problems
Content-Type: text/plain; charset=iso-8859-1
Content-Transfer-Encoding: 7bit
Received-SPF: none (hypatia: domain of [EMAIL PROTECTED] does not designate permitted
sender hosts)
X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch
X-Spam-Checker-Version: SpamAssassin 2.63 (2004-01-11) on
Hi,
Take a look at read.xls() from gregmisc package.
HTH,
Kevin
On Thu, 19 Aug 2004, Sivan Aldor wrote:
i want to upload different sheets from one excel file.
does anyone know of a way to do so? (i use read.csv to load the excel but
i would like to load form the sheets themselves.)
thanks
I've got the same problem with help.start(). It seems to be a problem with
the 'open'
command when executed through R's 'system' function -
ie system('open something.html')
gives an error message, whilst opening a new xterm (or console)
and entering
open something.html
works fine (here
-Original Message-
From: Nils Aschenbruck
Sent: Thursday, August 19, 2004 11:51 AM
To: [EMAIL PROTECTED]
Subject: [R] Erlang distribution
Hello,
is there a packet that supports the Erlang distribution?
I want to use this distribution for tests against empirical
Here is one way to read the mysheet1 worksheet from the
mydata.xls file:
library(RODBC)
z - odbcConnectExcel(c:\\mydata.xls)
myframe - sqlFetch(z, mysheet1)
close(z)
Your data will be in myframe. Thanks to Michael Lapsley
and Brian Ripley for the RODBC package. See the R
On Wed, 18 Aug 2004, Moises Hassan wrote:
Is there a way to get the distance between each point and the center of
the cluster it was assigned to in cluster methods such as agnes, pam,
and clara.
I don't think that it's directly available, but for pam and clara you can
easily get it from your
Hi,
I'd like to use the font Frutiger in my pdf files generated by R.
but I only manage this if I create a postscript file with
postscript() and the option family = c( LT_14361.afm,
LT_14363.afm, LT_14373.afm, LT_14375.afm ) and convert this
postscript file to pdf.
If I specify the option
On 19 Aug 2004 at 12:05, Sivan Aldor wrote:
i want to upload different sheets from one excel file.
does anyone know of a way to do so? (i use read.csv to load the excel
but i would like to load form the sheets themselves.) thanks sivan
Hi
Open Excel, go to the sheet, select rectangular
I'm baffled as to how the Lattice package achieves clipping. Would
someone mind explaining this to me?
Firstly, my attempt using just the grid package:
x-seq(0,3,by=0.3)/2.8
y-seq(0,1,by=0.1)
grid.newpage()
grid.rect(gp=gpar(fill=pink))
vp-viewport(width=0.8, height=0.8)
pushViewport(vp)
Have you read the help file? It says (R 1.9.1)
family: the font family to be used, one of 'AvantGarde',
'Bookman', 'Courier', 'Helvetica',
'Helvetica-Narrow', 'NewCenturySchoolbook', 'Palatino'
or 'Times'.
Now, it would seem if you can get away with specifying a
The following is not what I expected in sorting characters (single letters
and the same letters with preceding spaces).
Can someone enlighten me as to why the following might be a correct result
for sorting?
; x - c(LETTERS[1:3], paste( , LETTERS[1:3], sep=))
; x
[1] A B C A B C
; sort(x)
Dear all.
This is probably more PDF then R issue but anyway, maybe
somebody can help me.
I try to persuade pdf device to use eastern european (Czech)
encoding but I was not successful yet. Working on W 2000, R
1.9.0, lattice plots.
pdf(redprov.pdf, width = 5, height = 5, bg=white,
It is documented to depend on your locale. I get
sort(x)
[1] A B C A B C
in the C locale. The help page does say so:
The sort order for character vectors will depend on the collating
sequence of the locale in use: see 'Comparison'.
The default collation sequences for
Hi there,
I am using Fortran 77 g77 compiler to write Fortran subprogram to R. I
would like to change and write Fortran 90 code instead (no more
specification of matrix!). Does gcc support Fortran 90? And if Yes,
does it works in R?
Ben
__
[EMAIL
Sorting depends on the locale. For example, I get
x - c(LETTERS[1:3], paste( , LETTERS[1:3], sep=))
x
[1] A B C A B C
sort(x)
[1] A B C A B C
On Linux (Fedora Core 1) I set LANG=C, which is necessary for some
other (non-R) things to work.
-roger
[EMAIL PROTECTED] wrote:
The
I do get (R-1.9.1 on WinXPPro):
x - c(LETTERS[1:3], paste( , LETTERS[1:3], sep=))
sort(x)
[1] A B C A B C
The `right' sequence is locale dependent, and I wonder if that's the
discrepancy here.
Andy
From: [EMAIL PROTECTED]
The following is not what I expected in sorting characters
Dear Rui,
From my understanding of time-dependent covariates (not an expert but have
been working on a similar problem), it would appear that the coding of the
status column is not correct. Unless you have observed an event at each
interval you should only have status=1 for the last interval. In
Gabor Grothendieck ggrothendieck at myway.com writes:
Roger Levy rog at stanford.edu writes:
: Are there labeled break statements in R? i.e., something along the
: lines of
:
: TOPLOOP: for(i in 1:m) {
: for(j in 1:n) {
: ...
: if(condition) {
:break TOPLOOP
: }
On Thu, 19 Aug 2004, Ben Yip wrote:
I am using Fortran 77 g77 compiler to write Fortran subprogram to R. I
would like to change and write Fortran 90 code instead (no more
specification of matrix!). Does gcc support Fortran 90?
Isn't that a question for a gcc-help list? Take a look at
Thank you very much to Brian Ripley, Roger Peng, and Andy Liaw. Everyone
pointed out the same solution.
Setting
LC_COLLATE=C
did it. This default setting is indeed odd to me.
Thanks again!
Andreas Krause
__
[EMAIL PROTECTED] mailing list
You need to find out what put /usr/local/lib/libxml2.2.dylib
there and correct it. Look at PR#7063 in the R-bugs archive.
Stefano Iacus said there
Please remove your copy of libxml2 in /usr/local/lib, you no longer
need it on Panther and the one included in 10.3.4 is incompatible with
the
I'm not in your target audience but I consult with and teach to your target audience,
so I'll chime in anyway.
1. Look at the user friendliness of the alternatives you offered to R. I haven't
used SPSS recently or Stata at all but I know SAS is NOT that user friendly. Note that
my institute
I have a suggestion for the posting guide. One problem with some posts is that
they do not provide an example that can be reproduced. I think that many
people just do not know how to easily specify some data and some technical
assistance should be provided in the posting guide. If the problem
On Thu, Aug 19, 2004 at 01:50:47PM +0100, Prof Brian Ripley wrote:
On Thu, 19 Aug 2004, Ben Yip wrote:
I am using Fortran 77 g77 compiler to write Fortran subprogram to R. I
would like to change and write Fortran 90 code instead (no more
specification of matrix!). Does gcc support
[EMAIL PROTECTED] writes:
Thank you very much to Brian Ripley, Roger Peng, and Andy Liaw. Everyone
pointed out the same solution.
Setting
LC_COLLATE=C
did it. This default setting is indeed odd to me.
It's odd to everyone, except perhaps to the POSIX locale mafia... Once
upon a time all
Content-Type: text/plain; charset=iso-8859-1
Content-Transfer-Encoding: 7bit
Received-SPF: none (hypatia: domain of [EMAIL PROTECTED] does not designate permitted
sender hosts)
X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch
X-Spam-Checker-Version: SpamAssassin 2.63 (2004-01-11) on
Technical point:
Gabor Grothendieck wrote
``data(iris) will attach data set iris''
This is not true --- data(iris) put ``iris'' into the Global
environment; if one wants it attached, one has to do the follow up by
attach(iris).
cheers,
Hi,
Z is a three level factor variables and ZB and ZC are its corresponding binary dummy
variables. In a 'glm' of Y on Z and X, say, how do the two test specifications
test.terms = c(ZB:X,ZC:X) # and
test.terms = ~ ZB:X + ZC:X
in 'regTermTest' differ? I thought that both
Hi Susanne,
* Susanne Osterkamp [EMAIL PROTECTED] [040819 15:52]:
Hello,
I want to validate my data for working with CART.
I need to now how to create a randomized column of some 6800 intgers
within R.
Do you want to sample 6800 random integers? In that case try ?runif
(for uniform
On Thu, Aug 19, 2004 at 09:36:22AM -0300, Hanke, Alex wrote:
Dear Rui,
From my understanding of time-dependent covariates (not an expert but have
been working on a similar problem), it would appear that the coding of the
status column is not correct. Unless you have observed an event at each
Hi,
I've been struggling with RMySQL installation
for several days. I'm running Solaris 2.9
and MySQL 4.0. RMySQL (v0.4.3) seems to install
OK, but when I load it in R I get the message:
---
library(RMySQL)
Error in dyn.load(x, as.logical(local),
as.logical(now)) :
unable to load
A few comments.
1. What students use.
I've gotten a few undergraduate psychology students to use R to
analyze their data for independent study projects. They do fine.
I prepare the data for them and show them how to do t tests and
correlations, and how to define new variables from old ones.
On Thursday 19 August 2004 06:38, Matthew Walker wrote:
I'm baffled as to how the Lattice package achieves clipping. Would
someone mind explaining this to me?
Firstly, my attempt using just the grid package:
x-seq(0,3,by=0.3)/2.8
y-seq(0,1,by=0.1)
grid.newpage()
Dear R-Help list,
I have a problem with convertions of strings. I want to use the function
paste() to create an object name and then use that character string to
call on that object. So, for example:
dat99 - matrix(rbind(1,1,2),3,3)
no - 99
dat - paste(dat,no,sep=)
dat
[1] dat99
What should I
Use sample() in either case:
## randomly shuffle the elements of the vector `x':
rx - sample(x, replace=FALSE)
## randomly draw 6800 integers between 1 and 6800:
x - sample(6800, 6800, replace=TRUE)
Andy
From: Olaf Mersmann
Hi Susanne,
* Susanne Osterkamp [EMAIL PROTECTED] [040819 15:52]:
See ?get.
Andy
From: Gudrun Jonasdottir
Dear R-Help list,
I have a problem with convertions of strings. I want to use
the function
paste() to create an object name and then use that
character string to
call on that object. So, for example:
dat99 - matrix(rbind(1,1,2),3,3)
no -
On Thu, 2004-08-19 at 01:45, Jacob Wegelin wrote:
Alternate title: How can I persuade my students that R is for them?
Alternate title: Can R replace SAS, SPSS or Stata for clinicians?
I am teaching introductory statistics to twelve physicians and two veterinarians
who have enrolled in a
Hi Gudrun,
try to use
eval(parse(text=dat))
Best,
Dimitris
Dimitris Rizopoulos
Doctoral Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/396887
Fax: +32/16/337015
Web:
Gudrun Jonasdottir wrote:
Dear R-Help list,
I have a problem with convertions of strings. I want to use the function
paste() to create an object name and then use that character string to
call on that object. So, for example:
dat99 - matrix(rbind(1,1,2),3,3)
no - 99
dat - paste(dat,no,sep=)
dat
assign(dat, get(paste(dat, no, sep=)))
or simply
dat - get(paste(dat, no, sep=))
On Thu, 2004-08-19 at 15:15, Gudrun Jonasdottir wrote:
Dear R-Help list,
I have a problem with convertions of strings. I want to use the function
paste() to create an object name and then use that
And this is all in the FAQ.
-thomas
On Thu, 19 Aug 2004, Adaikalavan Ramasamy wrote:
assign(dat, get(paste(dat, no, sep=)))
or simply
dat - get(paste(dat, no, sep=))
On Thu, 2004-08-19 at 15:15, Gudrun Jonasdottir wrote:
Dear R-Help list,
I have a problem with
On Thu, 19 Aug 2004, Göran Broström wrote:
On Thu, Aug 19, 2004 at 09:36:22AM -0300, Hanke, Alex wrote:
Dear Rui,
From my understanding of time-dependent covariates (not an expert but have
been working on a similar problem), it would appear that the coding of the
status column is not
On Thu, 19 Aug 2004, Almirall, Daniel wrote:
This is a question regarding the 'regTermTest' function in the 'survey' package.
Imagine Z as a three level factor variable, and code ZB and ZC as the two
corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z
and X,
Dear Alex,
Thanks a lot for your help! The end is not proportional to the
covariate in my data set. I recoded the status as you suggested, but the
problem still exists. I tried to fit the model with partial data, it has
no problem when I use the first 25 rows (from the data in the form
start,
I want to run a GEE for a time series of counts. The data are daily
respiratory mortality counts and so there aren't any 'clusters' in the
longitudinal sense. Neither the gee or geese packages work. The gee one
wont run at all and the geese one produces NaNs or just runs
indefinitely depending on
Hi,
Has anyone had any joy getting the rgdal package to compile under windows?
I've been trying with MinGW using Duncan Murdoch's Rtools but to no avail.
After unpacking the tarball I put it here:
D:\rw1091\src\library\rgdalls
DESCRIPTION Rconfigure inst man
INDEX
I compiled the 1.9.1 src.rpm with the standard gnu tools and it works.
I tried compiling the 1.9.1 src.rpm with the Intel 8 C and FORTRAN
compilers and it bombs out during the testing phase:
comparing 'd-p-q-r-tests.Rout' to './d-p-q-r-tests.Rout.save' ...267c267
df = 0.5[1] Mean
Dear group,
Apologies for asking the most chomped FAQ.
I have a file with list of gene names(genes.txt):
EGF
EGFR
PTPN6
TIEG2
MAPK1
I have another object in R, I do not know the data
type for that object that looks like this
(lidnames):
RABGGTA MAPK3 TIE CYP2C19
lidnames[1:10]
Hello all.
I was wondering if there is a way to hold constant the fixed effects correlation
structure across multiple groups?
For example, I have two groups and I fit a three parameter logistic growth curve where
the fixed effects are free to vary across the groups. I'll paste in the code as a
Use the get() function:
dat99 - matrix(c(1,1,2),3,3)
no - 99
dat.name - paste(dat,no,sep=)
get(dat.name)
[,1] [,2] [,3]
[1,]111
[2,]111
[3,]222
At 04:15 PM 8/19/2004 +0200, Gudrun Jonasdottir wrote:
Dear R-Help list,
I have a problem with convertions
We are finding more and more clinical researchers interested in learning
R and are starting to teach R to clinicians. We have put some teaching
material on our site: http://biostat.mc.vanderbilt.edu
--
Frank E Harrell Jr Professor and Chair School of Medicine
S Peri wrote:
Dear group,
Apologies for asking the most chomped FAQ.
I have a file with list of gene names(genes.txt):
EGF
EGFR
PTPN6
TIEG2
MAPK1
I have another object in R, I do not know the data
type for that object that looks like this
(lidnames):
RABGGTA MAPK3 TIE CYP2C19
You wrote:
What's the difference between t.test(x, y) and pairwise.t.test()? Is
it just that the former takes two vectors, whereas the latter takes a
vector and a factor?
No. The pairwise.t.test() function (according to the help
file) does a multiplicity of t-tests, on more
Petr Pikal wrote:
Dear all.
This is probably more PDF then R issue but anyway, maybe
somebody can help me.
I try to persuade pdf device to use eastern european (Czech)
encoding but I was not successful yet. Working on W 2000, R
1.9.0, lattice plots.
pdf(redprov.pdf, width = 5, height = 5,
Hi,
Is there any IDE or any editor or any pice of code
for .vimrc is available?
Please let me know.
Thanks
SP
___
Win 1 of 4,000 free domain names from Yahoo! Enter now.
__
[EMAIL PROTECTED] mailing
From: Rolf Turner [EMAIL PROTECTED]
No. The pairwise.t.test() function (according to the help
file) does a multiplicity of t-tests, on more than two
samples, adjusting the p-value to compensate for the
multiplicity by various methods.
Thank you, Rolf, that's helpful. So pairwise.t.test() is used
On Thu, 19 Aug 2004 13:42:21 -0300 (ADT), Rolf Turner
[EMAIL PROTECTED] wrote :
You wrote:
What's the difference between t.test(x, y) and pairwise.t.test()? Is
it just that the former takes two vectors, whereas the latter takes a
vector and a factor?
No. The pairwise.t.test()
Dear List,
Is the test for proportional hazards valid when the model contains a cluster
variable? The output looks strange with the cluster variable.
My intervals are based on calendar time and the clustering variable is
related to the season the event occurs in.
From: Duncan Murdoch
On Thu, 19 Aug 2004 13:42:21 -0300 (ADT), Rolf Turner
[EMAIL PROTECTED] wrote :
You wrote:
What's the difference between t.test(x, y) and
pairwise.t.test()? Is
it just that the former takes two vectors, whereas the
latter takes a
vector and a factor?
Is there a way to list the contents of a save file,
i.e. the variable names, without load'ing from it first?
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hi,
I am running a ANOVA on a factorial design, and using TukeyHSD
for post hoc comparisons. I have 2 factors with three levels each:
Factor B
Factor A1 2 3
1
2
3
When I look at the Tukey output (on the interaction of the factors) the
The Intel compiled version also fails the below test:
### Very big and very small
umach - unlist(.Machine)[paste(double.x, c(min,max), sep='')]
xmin - umach[1]
xmax - umach[2]
tx - unique(outer(-1:1,c(.1,1e-3,1e-7)))# 7 values (out of 9)
tx - unique(sort(c(outer(umach,1+tx#
Jack Tanner [EMAIL PROTECTED] writes:
From: Rolf Turner [EMAIL PROTECTED]
No. The pairwise.t.test() function (according to the help
file) does a multiplicity of t-tests, on more than two
samples, adjusting the p-value to compensate for the
multiplicity by various methods.
Thank you,
Users,
Where I can to download rmutil package?
thanks
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Peri,
Please stop posting to both R-help and BioConductor simultaneously.
Please read the posting guide and decide which list is more appropriate.
http://www.R-project.org/posting-guide.html
On Thu, 2004-08-19 at 17:20, S Peri wrote:
Dear group,
Apologies for asking the most chomped
On Thu, 19 Aug 2004, Harris A. Jaffee wrote:
Is there a way to list the contents of a save file,
i.e. the variable names, without load'ing from it first?
I believe not. But that's not too bad:
x - system.file(package=MASS, data, hills.rda)
print(load(x , new.env()))
[1] hills
--
Brian D.
On Thu, 19 Aug 2004, Hem Nalini Morzaria-Luna wrote:
Hi,
I am running a ANOVA on a factorial design, and using TukeyHSD
for post hoc comparisons. I have 2 factors with three levels each:
Factor B
Factor A1 2 3
1
2
3
When I look at the
I am not sure what does an IDE consist besides syntax highlighting and
ability to parse lines from script into console. Nor am I familiar with
VIM but IMHO emacs with ESS is the best for R.
On Thu, 2004-08-19 at 19:13, S Peri wrote:
Hi,
Is there any IDE or any editor or any pice of code
for
Where I can to download rmutil package?
You can use the reposTools package from Bioconductor:
install.packages2(rmutil)
There's also the website:
http://popgen0146uns50.unimaas.nl/~jlindsey/rcode.html
__
[EMAIL PROTECTED] mailing list
On Thu, 19 Aug 2004, Savano wrote:
Where I can to download rmutil package?
See the FAQ 5.1.5.
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
--
Brian D. Ripley, [EMAIL PROTECTED]
Professor of Applied Statistics,
On Thu, 2004-08-19 at 14:42, Liaw, Andy wrote:
From: Duncan Murdoch
On Thu, 19 Aug 2004 13:42:21 -0300 (ADT), Rolf Turner
[EMAIL PROTECTED] wrote :
You wrote:
What's the difference between t.test(x, y) and
pairwise.t.test()? Is
it just that the former takes two
Hi!!
Has anyone used the coagulation data for statistical analysis ? I
managed to
get the data from the web but unsure of the way its supposed to read. I am
new to R so trying to gets myself familiarized with the statistical
tools using
available data. I am appending the data I got of the web.
Hi there!
Is ROracle available for Windows, please?
I found a download site, but it's really for UNIX/Linux.
Here is a thought question, please: Why do the vector elements start at
location 1 rather than zero, as C does?
Thanks in advance!
R Version 1.9.1 Windows
Sincerely,
Laura Holt
mailto:
Hi
I have a question about how to increase my memory size, could someone answer it for
me??
I am using Bioconductor in R to calculate gene expression values with mas5, dchip, and
mas4. I have only 18 samples, all from Affymetrix U133A Plus 2 arrays, which has
~54,000 genes. My machine
It is frustrating to see the labels I want in the dimensions of a list but not be able
to extract those labels into titles for plots generated from component objects. If
someone could set me straight, I would appreciate it. For your amusement, I have
provided an example of the Byzantine code I
Please read the rw-FAQ, as the posting guide asks. The answer is there.
On Thu, 19 Aug 2004, Stephanie Tsung wrote:
I have a question about how to increase my memory size, could someone
answer it for me??
I am using Bioconductor in R to calculate gene expression values with
mas5, dchip,
Laura Holt [EMAIL PROTECTED] writes:
Hi there!
Is ROracle available for Windows, please?
I found a download site, but it's really for UNIX/Linux.
As I understand it, a Windows version can easily (?) be built, but to
do so you need Oracle libraries and header files, which you can't have
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Jake,
this is a great question.
I have experience teaching R to forestry and natural resources graduate and
undergraduate students. These students are amongst the least numerically comfortable
students I've taught. I have found that the numerically comfortable among them leap
into the
Hi
Deepayan Sarkar wrote:
On Thursday 19 August 2004 06:38, Matthew Walker wrote:
I'm baffled as to how the Lattice package achieves clipping. Would
someone mind explaining this to me?
Firstly, my attempt using just the grid package:
x-seq(0,3,by=0.3)/2.8
y-seq(0,1,by=0.1)
grid.newpage()
Dear R-help,
Please can someone explain how to put a column name on an output data.frame.
##Starting with a data.frame with 3 columns (d$Year, d$NoIndiv, d$wtd_tl)
yr_ind - split (d$NoIndiv, d$Year)
yr_tl - split (d$wtd_tl, d$Year)
ann_ind - sapply (yr_ind, sum)
ann_tl - sapply (yr_tl, sum)
Professor Ripley,
These are optimization problems with multiple local maxima, and
like any complex statistical fitting problem you should not expect all
programs to give the same answer.
Examining the log-likelihoods from the anova statements, while fit1 and m1 are
similar, the maximum
Is names(d2) - c(year,av_t1) what you're looking for?
steve
-Original Message-
From: [EMAIL PROTECTED] [mailto:r-help-
[EMAIL PROTECTED] On Behalf Of Louize Hill
Sent: Thursday, August 19, 2004 7:26 PM
To: [EMAIL PROTECTED]
Subject: [R] column names in data.frame
Dear R-help,
Hello, all;
I get an unexpected result when trying to plot a probability histogram
with R1.9.1 on windows xp:
#with the following code:
x - runif(100,0,1)
hist(x)
hist(x, freq=F)
h - hist(x, freq=F)
summary(h)
#Length Class Mode
#breaks 11 -none- numeric
#counts 10
apparently not
that gives me the following error:
Error in names-.default(`*tmp*`, value = c(year, av_t1)) :
names attribute [2] must be the same length as the vector [1]
presumably because the 2nd column (av_tl) already has a column name, and
therefore seems longer?
- Original
Alec Stephenson
Department of Statistics
Macquarie University
NSW 2109, Australia
Joseph LeBouton [EMAIL PROTECTED] 08/20/04 09:56am
Hello, all;
I get an unexpected result when trying to plot a probability histogram
with R1.9.1 on windows xp:
Alec,
Thanks for your reply. I guess what I'm getting at is that I to plot
the histogram such that the HEIGHT of each bar represents the proportion
of that class in the sample. From your reply I gather that the AREA of
each bar is currently representing the proportion.
My current work-around
Hi Louize,
Sorry to lead you astray, but I misinterpreted what you were doing (I could
still be leading you astray :-(. You should be able to make a names
attribute for d2 by names(d2) - c(year, col2namewhatyouwant,
col3namewhatyouwant, ...) to the length of av_t1. When you created the
data frame
I suspect you are interpreting the column names as a column, whereas in
fact you only have a single column data frame, so using
names(d2) - whatever
will work.
Thanks,
Alec
Alec Stephenson
Department of Statistics
Macquarie University
NSW 2109,
1 - 100 of 110 matches
Mail list logo