Thanks a ton for your help. Now I'm totally
an SDI person.
Cheers,
SI
FYI
ISO-2002-JP and us all believe in \, which is our
currency symbol. We use it in place of a backslash.
For some good reasons, ancient programmers
assigned the same code in US-ASCII fo
I left out the brackets in my last email but the problem
(a reappears after have been substituted out) still remains:
> z <- substitute( function(){a+1}, list(a=quote(b)) )
> z
function() {
b + 1
}
> eval(z)
function(){a+1}
---
Date: Wed, 17 Mar 2004 20:10:43 -0500 (EST)
From: Gabor
Hi,
I'm having trouble with glmmPQL.
I'm fitting a 2 level random intercept model, with 90,000 cases and about 330 groups.
I'm unable to get any results on the full data set. I can get it to work if I sample
down to about 30,000 cases. But for models with N's much larger than that I get the
Hello,
Is their any function in R to calculate multivariate
generalization of AIC or SBC?
Ideally, I am looking for a function which could generate a
tabular results set of AIC & SBC values for different lag length.
I think a combination of glm & AIC function could do the
Consider the following example:
# substitute a with b in the indicated function. Seems to work.
> z <- substitute( function()a+1, list(a=quote(b)) )
> z
function() b + 1
# z is an object of class call so use eval
# to turn it into an object of class expression; however,
# when z is evaluated,
SI wrote:
> Hi,
>
> thanks for your reply. But I guess I misfired a
> question.
>
> I want to save the size of _GUI_ window, NOT console
> window. GUI window is the window with "RGui"
One option to consider is to start R in "sdi" mode. My shortcut for
R-1.8.1 on WinXP is...
"C:\Program File
Hi,
thanks for your reply. But I guess I misfired a
question.
I want to save the size of _GUI_ window, NOT console
window. GUI window is the window with "RGui"
writtenon the left of title bar. Console window
is inside of the GUI window with "R Console" on
its t
Hi Marc,
As you guessed my current point is like your second example. I actually
have three groups of three, so the normal box plot drws them at positions,
1:9 by the time I've placed them where I'd like they plot at positions
(roughly) 1:8, leaving me with a large right hand space that I'd like
Cristian Pattaro asked:
> is there a automatic method to obtain Odds Ratio estimates and their
> confidence intervals from a GLM Logistic model?
something like:
oddsr<-function(x) {
require(MASS)
round(exp(cbind(Odds.Ratio=coef(x), confint(x))),3)
}
oddsr(glm(y ~ x, family=binomial))
Da
> From: Thomas Lumley [mailto:[EMAIL PROTECTED]
>
> On Wed, 17 Mar 2004, Liaw, Andy wrote:
>
> > Sorry for being dense: So What is the likely source of the
> problem, then?
> > Is it the compiler, the OS, R itself, or some combination
> of these? Any
> > suggestion on how to resolve this?
>
On Wed, 2004-03-17 at 16:32, Mark Van De Vyver wrote:
> Hi Marc,
> Thanks for the response. Apologies for the lack of detail... I have used
> 'boxwex', 'at' and ylim, to place things nicely, except I know have some
> spare room at the end of the plot, which I would like to get rid of Your
> su
?dump
?savehistory
Icabalceta, Jorge L. wrote:
Dear friends,
In R there is a file | save to file option which takes a "picture" of
the text in your screen and saves it to a text file. Is there any way to do
the same with a function?
__
[EMAIL PROTECT
Hi,
I am having a problem to understand why as.data.frame method doesn't
dispatch properly on my class:
> setClass("Foo", "character")
[1] "Foo"
> as.data.frame(list(foo=new("Foo", .Data="a")))
Error in as.data.frame.default(x[[i]], optional = TRUE) :
can't coerce Foo into a data.frame
I was e
On Wed, 17 Mar 2004, Liaw, Andy wrote:
> Sorry for being dense: So What is the likely source of the problem, then?
> Is it the compiler, the OS, R itself, or some combination of these? Any
> suggestion on how to resolve this?
In a sense it is R itself, making assumptions that hold on many, but
Dear friends,
In R there is a file | save to file option which takes a "picture" of
the text in your screen and saves it to a text file. Is there any way to do
the same with a function? I need this because I run a long program and when
I get the results and try to use file | save to file, the p
Prof Brian Ripley <[EMAIL PROTECTED]> writes:
> I think NA+0 should be special, though. All of mine give
>
> > NA+NaN
> [1] NA
> > NaN+NA
> [1] NA
>
> That's with four separate compilers on Sparc-Solaris as well as Windows
> and Linux i686.
>
> I have long wondered why it works
I get the
Hi Marc,
Thanks for the response. Apologies for the lack of detail... I have used
'boxwex', 'at' and ylim, to place things nicely, except I know have some
spare room at the end of the plot, which I would like to get rid of Your
suggestion to use 'add' per the example is one I had not thought o
I think the easiest way is using na.omit function which omits pairs of data
Look at the following example:
> x <- c(0,1,2,3,4,5)
> y <- c(NA,2, 3, 5, NA, 6)
> data.1 <- data.frame(x,y)
> data.1
x y
1 0 NA
2 1 2
3 2 3
4 3 5
5 4 NA
6 5 6
> new.data <- na.omit(data.1)
> new.data
x y
2 1 2
3
What is it "Package bundles" and 1.9.0 (alpha) ?
'Rcmd check Bundle' aborts with:
Error: cannot open file '.../Bundle.Rcheck/Package/R/Package' for reading
In ".../Bundle.Rcheck/" the packages of the bundle are "merged" into a single
"package". ('Rcmd INSTALL' does also "merge" the packages of th
"Marc R. Feldesman" <[EMAIL PROTECTED]> writes:
> At 12:28 PM 3/17/2004, Ed L Cashin wrote:
>
> I hope this response isn't indicative of the speed with which gmane
> posts messages. I think this entire thread was more than 7 or 8
> months ago, possibly longer.
Your hopes are fulfilled. The reas
On Wed, 2004-03-17 at 14:31, Mark Van De Vyver wrote:
> Hi,
> Could anyone hint at how to set the x axis plot range in boxplot/bxp? The
> docs are clear that usr and xlim are not passed through/down
> TIA
> Mark
It is not entirely clear what you are looking to do here.
Are you attempting to pl
Dear R colleagues:
Am trying to fit a simple NL model to determine Economical Optimum Nitrogen
Rates.
The segmented (quadratic + plateau) model only works with some y's, in some
cases I get a "singular gradient" error.
I'll appreciate any ideas in how to solve the singular gradient error.
Thank
Sorry for being dense: So What is the likely source of the problem, then?
Is it the compiler, the OS, R itself, or some combination of these? Any
suggestion on how to resolve this?
Best,
Andy
> From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
>
> I think NA+0 should be special, though. All
At 12:28 PM 3/17/2004, Ed L Cashin wrote:
I hope this response isn't indicative of the speed with which gmane posts
messages. I think this entire thread was more than 7 or 8 months ago,
possibly longer.
>"Roger D. Peng" <[EMAIL PROTECTED]> writes:
>
>> Marc R. Feldesman wrote:
>>
>>> I agree w
On 17 Mar 2004 at 11:50, Cristian Pattaro wrote:
> Dear all,
>
> is there a automatic method to obtain Odds Ratio estimates and their
> confidence intervals from a GLM Logistic model?
>
The parameters in the linear predictor in logistic regression
are log odds ratios (assuming the default cont
Hi
Prof Brian Ripley wrote:
On Wed, 17 Mar 2004, Pisut Tempatarachoke wrote:
Hi all,
Could anyone help answer the following questions, please?
(I'm using R-1.8.1-4 on Fedora Core 1 and very new to R)
(i) Is it possible to specify a font family (e.g. courier or helvetica)
when graphing?
`when
"Roger D. Peng" <[EMAIL PROTECTED]> writes:
> Marc R. Feldesman wrote:
>
>> I agree with you on the "flood" of messages lately. Often this
>> flood accompanies a new release, but this flood has continued
>> unabated for longer than I would have imagined. The good news is
>> that R is becoming mo
Hi,
Could anyone hint at how to set the x axis plot range in boxplot/bxp? The
docs are clear that usr and xlim are not passed through/down
TIA
Mark
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read th
I think NA+0 should be special, though. All of mine give
> NA+NaN
[1] NA
> NaN+NA
[1] NA
That's with four separate compilers on Sparc-Solaris as well as Windows
and Linux i686.
I have long wondered why it works
On Wed, 17 Mar 2004, Liaw, Andy wrote:
> > From: Thomas Lumley [mailto:[EMAI
Prof Brian Ripley <[EMAIL PROTECTED]> writes:
> > tt<-table(zz%/%10)
> > n <- names(tt)
> > names(tt) <- paste(n,0,"-",n,9,sep="")
> > tt
> > data.frame(count=c(tt))
> >
> > Beware that empty groups are silently zapped, though.
>
> FWIW, table(factor(zz%/%10, levels=0:9)) avoids that
I knew, b
not.ok <- is.na(y)
y.new <- y[! not.ok]
x.new <- x[! not.ok]
should work.
jan
On Wed, 17 Mar 2004, klea lambrou wrote:
>
>hello R-users.I really need your help on these one.i have two vectors
>x and y of equal length.y contains missing values,so i need to remove
>them.i know how to
Do you not have an editor?
Sd2Rd is designed for rough conversion, especially for Ssgml files, not
for accurate and complete conversion. It is purely for convenience.
In this case, prompt() in S-PLUS will never create the file you quote. We
don't actually know exactly what the allowed format
On Wed, 17 Mar 2004, C Peroni wrote:
> Dear All,
>
> I am trying to visualize the surface of a bivariate kernel density
> estimate.
> I have a vector of bivariate observations(x,y), and a function which
> computes the kernel density estimate z corresponding to each
> observation. I cannot gener
On Wed, 2004-03-17 at 11:42, klea lambrou wrote:
>hello R-users.I really need your help on these one.i have two vectors
>x and y of equal length.y contains missing values,so i need to remove
>them.i know how to do that for y but i also need the corresponding x
>value to be removed t
> From: Thomas Lumley [mailto:[EMAIL PROTECTED]
>
> On Wed, 17 Mar 2004, Liaw, Andy wrote:
>
> > Dear R-devel,
> >
> > Has anyone seen this problem? We tried building R-1.8.1
> (and R-1.9.0 alpha
> > 2004-03-17) on an Irix 6.5 box using
> >
> > ./configure CC="cc -64" F77="f77 -64" --with-tclt
Assuming x contains the data, taking the solutions so far and
adding some minor improvements gives:
groups <- x %/% 10
lev <- min(groups):max(groups)
lab <- factor( paste( lev, "0-", lev, "9", sep = "" ) )
groups <- factor( groups, lev = lev, lab = lab )
tab <- table( groups,
Good point
Sean
-Original Message-
From: [EMAIL PROTECTED]
To: '[EMAIL PROTECTED] '
Sent: 3/17/2004 11:51 AM
Subject: RE: [R] projection pursuit
Yes, xgobi does projection pursuit. Though
I'm not so sure about 'projection pursuit density estimation'.
Not that I know what that is.
C Peroni wrote:
I am trying to visualize the surface of a bivariate kernel density
estimate.
I have a vector of bivariate observations(x,y), and a function which
computes the kernel density estimate z corresponding to each
observation. I cannot generate the (x,y) data in the ascending order
neede
On Wed, 17 Mar 2004, Liaw, Andy wrote:
> Dear R-devel,
>
> Has anyone seen this problem? We tried building R-1.8.1 (and R-1.9.0 alpha
> 2004-03-17) on an Irix 6.5 box using
>
> ./configure CC="cc -64" F77="f77 -64" --with-tcltk=no --enable-R-shlib
>
> make check failed because NA + 0 gave NaN ins
miss <- is.na(y)
y.good <- y[!miss]
x.good <- x[!miss]
HTH,
Andy
> From: klea lambrou
>
>hello R-users.I really need your help on these one.i have
> two vectors
>x and y of equal length.y contains missing values,so i
> need to remove
>them.i know how to do that for y but i also n
hello R-users.I really need your help on these one.i have two vectors
x and y of equal length.y contains missing values,so i need to remove
them.i know how to do that for y but i also need the corresponding x
value to be removed too.i cannot find or at least think of a command
which
See ?interp in the `akima' package. Also, I believe the `sm' package does
have function that produce the surface plot of 2D kernel density surface.
HTH,
Andy
> From: C Peroni
>
> Dear All,
>
> I am trying to visualize the surface of a bivariate kernel density
> estimate.
> I have a vector of
On Wed, 17 Mar 2004 16:20:06 +1100
Jim Lemon <[EMAIL PROTECTED]> wrote:
> Haven't there been one or two people who asked list members to submit
> peer-reviewed, published papers in which R was used for the analysis? I
> certainly have sent an email or two like this.
>
> Jim
>
Just to be provoc
Using R 1.8.1 on Solaris
Developing an Splus and R package from the same files.
Creating an alias for multiple functions to the same file in .sgml:
my.function
print.my.function
using R CMD Sd2Rd would convert to:
\name{my.function}
\alias{my.function}
\name{print.my.function}
\alias{p
Dear All,
I am trying to visualize the surface of a bivariate kernel density
estimate.
I have a vector of bivariate observations(x,y), and a function which
computes the kernel density estimate z corresponding to each
observation. I cannot generate the (x,y) data in the ascending order
needed by
Dear All,
I am trying to visualize the surface of a bivariate kernel density
estimate.
I have a vector of bivariate observations(x,y), and a function which
computes the kernel density estimate z corresponding to each
observation. I cannot generate the (x,y) data in the ascending order
needed by
Dear All,
I am trying to visualize the surface of a bivariate kernel density
estimate.
I have a vector of bivariate observations(x,y), and a function which
computes the kernel density estimate z corresponding to each
observation.
I cannot generate the (x,y) data in the ascending order needed by
Also, how many books explicitly discuss R -- and especially how many
mention R (and not S-Plus) in the title? spencer graves
Jim Lemon wrote:
Haven't there been one or two people who asked list members to submit
peer-reviewed, published papers in which R was used for the analysis? I
certainly
Yes, xgobi does projection pursuit. Though
I'm not so sure about 'projection pursuit density estimation'.
Not that I know what that is.
-Original Message-
From: Davis, Sean (NIH/NHGRI) [mailto:[EMAIL PROTECTED]
Sent: Wednesday, March 17, 2004 10:44 AM
To: '[EMAIL PROTECTED] '; '[EMAIL
Dear R-devel,
Has anyone seen this problem? We tried building R-1.8.1 (and R-1.9.0 alpha
2004-03-17) on an Irix 6.5 box using
./configure CC="cc -64" F77="f77 -64" --with-tcltk=no --enable-R-shlib
make check failed because NA + 0 gave NaN instead of NA. I've tried both
32- and 64-bit build, wi
On 17 Mar 2004, Peter Dalgaard wrote:
> Kai Hendry <[EMAIL PROTECTED]> writes:
>
>
> > 40-49 2
> > 50-59 15
> > 60-69 20
> > 70-79 19
> > 80-89 12
> > 90-99 2
>
> Here's another solution for this 10-year age group thing:
>
> tt<-table(zz%/%10)
> n <- names(tt)
> names(tt) <- paste(
Hi,
I am new at R, and I need some help with some basic stuff.
I want to do clustering of data. It works perfectly, but now I want to
export my file form R to my Linux directory. The problem is that the file
is generated only with the cluster numbers, but it does not carry over the
identifiers
Luis,
If I am not mistaken, the xgobi package (which requires xgobi executable)
does interactive projection pursuit. Others can correct me if I'm wrong and
there may be other packages, as well.
Sean
-Original Message-
From: [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Sent: 3/17/2004 9:59 AM
Kai Hendry <[EMAIL PROTECTED]> writes:
> 40-49 2
> 50-59 15
> 60-69 20
> 70-79 19
> 80-89 12
> 90-99 2
Here's another solution for this 10-year age group thing:
tt<-table(zz%/%10)
n <- names(tt)
names(tt) <- paste(n,0,"-",n,9,sep="")
tt
data.frame(count=c(tt))
Beware that empty gro
On Wed, Mar 17, 2004 at 04:55:19PM +0200, Kai Hendry wrote:
> I am trying to construct a frequency table. I guess this should be done with
> table. Or perhaps factor and split. Or prop.table. cut? findInterval? Argh!
> I got this far:
> > table(cut(zz$x9, brk))
>
> (40,50] (50,60] (60,70] (7
On Wed, 2004-03-17 at 08:55, Kai Hendry wrote:
> This must be FAQ, but I can't find it in archives or with a site search.
>
> I am trying to construct a frequency table. I guess this should be done with
> table. Or perhaps factor and split. Or prop.table. cut? findInterval? Argh!
>
> Please corre
Kai Hendry <[EMAIL PROTECTED]> writes:
> This must be FAQ, but I can't find it in archives or with a site search.
>
> I am trying to construct a frequency table. I guess this should be done with
> table. Or perhaps factor and split. Or prop.table. cut? findInterval? Argh!
>
> Please correct me i
Kai Hendry wrote:
This must be FAQ, but I can't find it in archives or with a site search.
I am trying to construct a frequency table. I guess this should be done with
table. Or perhaps factor and split. Or prop.table. cut? findInterval? Argh!
Please correct me if what I am looking for is not cal
?data.frame
data.frame( table(cut(x, seq(0, 1, by=0.1))) )
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of Kai Hendry
> Sent: 17 March 2004 14:55
> To: [EMAIL PROTECTED]
> Subject: [R] Frequency table
>
>
> This must be FAQ, but I can't find it in ar
I guess you want something like:
table(cut(zz$x9, c(-Inf, seq(40, 90, by=10), Inf)))
HTH,
Andy
> From: Kai Hendry
>
> This must be FAQ, but I can't find it in archives or with a
> site search.
>
> I am trying to construct a frequency table. I guess this
> should be done with
> table. Or perh
I like the previous suggestion of counting the number of unique e-mails in
the archive.
Another useful thing would be to count and plot the growth of number of R
(and Bioconductor) packages over the years. I know not all packages are
created equal.
Regards, Adai.
> -Original Message-
> F
Dear helpers
Does R have a package that performs projection pursuit density estimation? Or
anyone knows code in Matlab or C for example to do this?
Thank you all
Luis
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r
This must be FAQ, but I can't find it in archives or with a site search.
I am trying to construct a frequency table. I guess this should be done with
table. Or perhaps factor and split. Or prop.table. cut? findInterval? Argh!
Please correct me if what I am looking for is not called a "frequency t
Sophos Plc MailMonitor for Domino/D R1.0(4.003c)
Server: BC1
---
Your email contained infected attachment(s). For advice consult your
system administrator.
---
Dear R-list users,
I'm new to principal components and factor analysis.
I thought this method can be very useful for me to find relationships
between several variables (which I know there is, only don't know which
variables exactly and what kind of relation), so as a structure
detection method.
N
On Tue, 16 Mar 2004 13:41:38 -0800, Spencer Graves
<[EMAIL PROTECTED]> wrote :
> Some web sites have "hit counters". It should be possible to get
>a counts of the numbers of times different parts of R are downloaded.
>Do the CRAN web sites include any such?
The logs aren't currently ava
Carlos, what you see are the lattice default settings;
use trellis.par.get/set to get and change them:
library(gstat) # requires lattice
data(meuse)
trellis.par.set("background", list(col="#ff"))
trellis.par.set("plot.symbol", list(cex=1,col="#00",font=1,pch="+"))
plot(variogram(zinc~1,~x+y
On Wed, 17 Mar 2004, Prof Brian Ripley wrote:
> Do you somehow have bash set as your shell for batch
> operations? On that platform you get sh or ksh by default,
> and bash is not even installed by default.
I didn't know that bash was set as the shell for batch
operations and it wouldn't have oc
hello,
try:
tt=tktoplevel()
#...
tkwm.state(tt,"withdrawn")
# the tt window disapperas but is not destroyed
#...
tkwm.state(tt,"normal")
# the tt window appears again.
hope this help,
DLC
Hi,
At the moment, I am working with the "tcltk" library.
I wonder if it is possible to create a window, destro
On Wed, 17 Mar 2004, Pisut Tempatarachoke wrote:
> Hi all,
>
> Could anyone help answer the following questions, please?
> (I'm using R-1.8.1-4 on Fedora Core 1 and very new to R)
>
> (i) Is it possible to specify a font family (e.g. courier or helvetica)
> when graphing?
`when graphing' means
This sounds like a good case for mixture regression, for which there's Fritz
Leisch's `flexmix' package. However, I don't think flexmix has facility for
testing whether the mixture has one vs. two components. Others on the list
surely would know more than I do.
HTH,
Andy
> From: Rob Knell
>
>
Hi all,
Could anyone help answer the following questions, please?
(I'm using R-1.8.1-4 on Fedora Core 1 and very new to R)
(i) Is it possible to specify a font family (e.g. courier or helvetica)
when graphing?
(ii) How can I make ticks point inwards on all four sides of a plot? Is
it possible
Hi,
At the moment, I am working with the "tcltk" library.
I wonder if it is possible to create a window, destroy
it and open it a new time.
For example
tt <- tktoplevel()
#... instructions
tkdestroy(tt)
#... other instructions
A this point, I would like to open tt from the tt
object. Is it po
mgcv 1.0 (package providing gams etc) will be released with R 1.9.0.
(R 1.8.x compatible versions can be found at:
http://www.stats.gla.ac.uk/~simon/simon/mgcv.html)
There are quite a few changes from mgcv 0.9: hence this message.
The main new features are:
* A generalized additive mixed modellin
lr.mod <- glm( y ~ x + w, family=binomial )
exp( summary( lr.mod )$coef[,1:2] %*% rbind( c(1,1,1), 1.96*c(0,-1,-1) )
)
should do the job.
Pack it in a function if you like, see e.g. the
(so far) undocumented function ci.lin in:
http://www.biostat.ku.dk/~bxc/R/ci.lin.R
(depends on)
http://www.bio
Hello people
I am doing some thinking about how to analyse data on dimorphic animals
- where different individuals of the same species have rather different
morphology. An example of this is that some male beetles have large
horns and small wings, and rely on beating the other guys up to get
Haven't there been one or two people who asked list members to submit
peer-reviewed, published papers in which R was used for the analysis? I
certainly have sent an email or two like this.
Jim
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.et
Thanks for the advise.
I think I've found the solution, the problem was that the file had been
written from a C program using a structure containing floats and int, if
I rewrite the C code not using a structure but independent floats and
ints the data then is read smoothly into R. I still can't
Angel Lopez wrote:
Hi,
I have a binary file with mixture of short integers (2 bytes) and floats
(4 bytes).
To get the data into R I use readBin reading each value one at a time
(n=1) but it is giving me headaches, two questions:
-In the "what" option is there any difference between "integer" and
Dear all,
is there a automatic method to obtain Odds Ratio estimates and their
confidence intervals from a GLM Logistic model?
Thanks
Cristian
=
Cristian Pattaro
=
Unit of Epidemiology & Medical Statistics
D
Hi,
I have a binary file with mixture of short integers (2 bytes) and floats
(4 bytes).
To get the data into R I use readBin reading each value one at a time
(n=1) but it is giving me headaches, two questions:
-In the "what" option is there any difference between "integer" and "int"?
-To read the
Dear all,
Is there any package in R can do the GEE2? Thanks.
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
82 matches
Mail list logo
