Re: [R] Sweave and multipage lattice

2006-07-19 Thread Friedrich Leisch
> On Wed, 19 Jul 2006 11:01:02 -0500, > Sundar Dorai-Raj (SD) wrote: > Dieter Menne wrote: >> Dear R-Listeners, >> >> as the Sweave faq says: >> >> http://www.ci.tuwien.ac.at/~leisch/Sweave/FAQ.html >> >> creating several figures from one figure chunk does not work, and

Re: [R] Automating package building packages and repository uploading

2006-07-19 Thread Prof Brian Ripley
On Thu, 20 Jul 2006, Carlos J. Gil Bellosta wrote: > Dear Rusers, > > I have developed two packages for a client of mine. After new features > are added or bugs corrected, I upload them to my own web repository. I > create both source and binary versions. binary Linux packages, it seems. The la

Re: [R] Confirmação de pedido para entrar no grup o python-brasil

2006-07-19 Thread Savio Ramos
On 20 Jul 2006 02:05:32 - Yahoo! Grupos <[EMAIL PROTECTED]> wrote: > > Olá [EMAIL PROTECTED], > > Recebemos sua solicitação para entrar no grupo python-brasil > do Yahoo! Grupos, um serviço de comunidades online gratuito e > super fácil de usar. > > Este pedido expirará em 7 dias. > > PA

Re: [R] Rearrange data.

2006-07-19 Thread Gabor Grothendieck
This is FAQ 7.17: http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-does-outer_0028_0029-behave-strangely-with-my-function_003f On 7/19/06, Ritwik Sinha <[EMAIL PROTECTED]> wrote: > Hi, > > I am trying to rearrange the following data > > d.f <- data.frame(x=c(1,1,2,2), y=c(1,2,1,2), vals=c("a11",

[R] Rearrange data.

2006-07-19 Thread Ritwik Sinha
Hi, I am trying to rearrange the following data d.f <- data.frame(x=c(1,1,2,2), y=c(1,2,1,2), vals=c("a11", "a12", "a21", "a22")) to look like a table with x as the rows and y as the columns, something like y 1 2 x 1a11 a12 2a21 a22 I tried doing this

Re: [R] Sweave and multipage lattice

2006-07-19 Thread Joel kincaid
On 7/18/06, Dieter Menne <[EMAIL PROTECTED]> wrote: > Dear R-Listeners, > > as the Sweave faq says: > > http://www.ci.tuwien.ac.at/~leisch/Sweave/FAQ.html > > creating several figures from one figure chunk does not work, and for > standard graphics, a workaround is given. Now I have a multipage tre

Re: [R] get rid of error in Factor Analysis

2006-07-19 Thread Andrew Robinson
Robert, try try(). Andrew. On Wed, Jul 19, 2006 at 10:27:07AM +0200, Robert Mcfadden wrote: > Dear All, > > I wrote a program and there is a loop. At each iteration I use maximum > likelihood factor analysis (?factanal). Output of factor analysis I use > later (in this loop). Unfortunately from

[R] Automating package building packages and repository uploading

2006-07-19 Thread Carlos J. Gil Bellosta
Dear Rusers, I have developed two packages for a client of mine. After new features are added or bugs corrected, I upload them to my own web repository. I create both source and binary versions. In fact, I made an script that checks, builds, and uploads them via ftp. However, I am facing two nuis

Re: [R] error when compiling "stats" library in R-2.3.1 on Solaris x86

2006-07-19 Thread Dongseok Choi
Dear Prof. Ripley, Following your advice, I read the manual and I realized that I missed two lines in my config.site. FC="f95 -xtarget=generic64" FCFLAGS="-O -I/mounts/devel/SUNWspro/prod/include" This solved my problem. It seemed that these two lines are new since R-2.3.0. I am updat

Re: [R] Wrap a loop inside a function

2006-07-19 Thread Gabor Grothendieck
foo can be written as a mapply: t(mapply(pcm, items, score, MoreArgs = list(theta = theta))) On 7/19/06, Doran, Harold <[EMAIL PROTECTED]> wrote: > I need to wrap a loop inside a function and am having a small bit of > difficulty getting the results I need. Below is a replicable example. > > > #

Re: [R] How to find S4 generics?

2006-07-19 Thread Spencer Graves
Hi, Thomas: Thanks very much. I haven't tried it yet, but it looks very useful. Best Wishes, Spencer Graves Thomas Lumley wrote: > On Wed, 19 Jul 2006, Spencer Graves wrote: >>Am I correct then that the 'methods' function could, at least >> theoretically, be revised so methods(class=..

[R] Data from Ying, Jung and Wei (1995)

2006-07-19 Thread Mai Zhou
Dear all, I am looking for the Small Cell Lung Cancer data from the paper Ying, Z., Jung, S. H., and Wei, L. J. (1995), Survival Analysis with Median Regression Models, Journal of the American Statistical Association} {\bf 90}, 178--184. Do any one know if it exist in an R readable format some

Re: [R] How to find S4 generics?

2006-07-19 Thread Thomas Lumley
On Wed, 19 Jul 2006, Spencer Graves wrote: > Am I correct then that the 'methods' function could, at least > theoretically, be revised so methods(class=...) could identify both S3 > and S4 methods (ignoring inheritance, as it does now, I believe)? > Here is a function to find methods for a

Re: [R] error when compiling "stats" library in R-2.3.1 on Solaris x86

2006-07-19 Thread Dongseok Choi
It was given iin the config.site. CC="cc -xtarget=generic64" CFLAGS="-O -I/mounts/devel/SUNWspro/prod/include" F77="f95 -xtarget=generic64" FFLAGS="-O -I/mounts/devel/SUNWspro/prod/include" CXX="CC -xtarget=generic64" CXXFLAGS="-O -I/mounts/devel/SUNWspro/prod/include" CPPFLAGS="-I/mounts/devel/S

Re: [R] error when compiling "stats" library in R-2.3.1 on Solaris x86

2006-07-19 Thread Prof Brian Ripley
How did 'cc -xtarget=generic64' get there? AFAIK R does not know about it, so presumably you specified it for CC. You need to do the same thing for *all* the compilers, that is CC, CXX, F77 and FC. The INSTALL file asked you to read the R-admin manual: there you will find a very similar examp

[R] error when compiling "stats" library in R-2.3.1 on Solaris x86

2006-07-19 Thread Dongseok Choi
Hello, I tried to compile v2.3.1 on Solaris x86 with SUN Pro compilers. I had an error while stats libarary was being compiled and I notice that "-xtarget=generic64" was not passed to f95 while cc used it. Could you tell me how to fix this problem? f95 -PIC -O -I/mounts/devel/SUNWspr

Re: [R] how to use large data set ?

2006-07-19 Thread Dylan Beaudette
Hi, While R is generally flexible enough for just about anything you can throw at it, detailed analysis of imagery might be better accomplished in a specialized piece of software. One option might be GRASS, which would allow you to do further processing on a subset of the original data in R. C

Re: [R] voronoi tessellations

2006-07-19 Thread Duncan Murdoch
On 7/19/2006 4:21 PM, [EMAIL PROTECTED] wrote: > Okay, been working with tripack, seems the most mature package for this. Got > it to work well with their test data set - data(tritest). When i tried > random numbers to explore further, i am getting some results that don't > reconcile. > > e

Re: [R] Wrap a loop inside a function

2006-07-19 Thread Doran, Harold
Very true, the resounding echo was large. Thanks, Doug. > -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf > Of Douglas Bates > Sent: Wednesday, July 19, 2006 4:20 PM > To: Doran, Harold > Cc: r-help@stat.math.ethz.ch > Subject: Re: [R] Wrap a loop inside

Re: [R] how to use large data set ?

2006-07-19 Thread mahesh r
Hi, I would like to extend to the query posted earlier on using large data bases. I am trying to use Rgdal to mine within the remote sensing imageries. I dont have problems bring the images within the R environment. But when I try to convert the images to a data.frame I receive an warning message f

Re: [R] voronoi tessellations

2006-07-19 Thread Thomas Lumley
On Wed, 19 Jul 2006, [EMAIL PROTECTED] wrote: > Okay, been working with tripack, seems the most mature package for this. > Got it to work well with their test data set - data(tritest). When i > tried random numbers to explore further, i am getting some results that > don't reconcile. > > examp

Re: [R] Correlation Mapping

2006-07-19 Thread justin rapp
Thank you for the advice. This is what I was looking for. If I would like to label the axes, how can I do that? Specifically I would like to label the x axis (AAA, AA, A, BBB) and same with the y axis. jdr On 7/17/06, Jim Lemon <[EMAIL PROTECTED]> wrote: > > justin rapp wrote: > > > On the co

Re: [R] Wrap a loop inside a function

2006-07-19 Thread Douglas Bates
I'm sure when you read this there will be the sound of a palm slapping a forehead. :-) The value of a function is the value of the last expression evaluated in the function. If you rewrite foo as foo <- function(theta,items, score){ like.mat <- matrix(numeric(length(items) * length(theta)),

Re: [R] voronoi tessellations

2006-07-19 Thread binabina
Okay, been working with tripack, seems the most mature package for this. Got it to work well with their test data set - data(tritest). When i tried random numbers to explore further, i am getting some results that don't reconcile. example run this: library(tripack) y <- runif(100) x <- runi

Re: [R] Wrap a loop inside a function

2006-07-19 Thread Sundar Dorai-Raj
Doran, Harold wrote: > I need to wrap a loop inside a function and am having a small bit of > difficulty getting the results I need. Below is a replicable example. > > > # define functions > pcm <- function(theta,d,score){ > exp(rowSums(outer(theta,d[1:score],'-')))/ > apply(exp(apply

[R] Wrap a loop inside a function

2006-07-19 Thread Doran, Harold
I need to wrap a loop inside a function and am having a small bit of difficulty getting the results I need. Below is a replicable example. # define functions pcm <- function(theta,d,score){ exp(rowSums(outer(theta,d[1:score],'-')))/ apply(exp(apply(outer(theta,d, '-'), 1, cumsum)), 2, s

Re: [R] how to use large data set ?

2006-07-19 Thread Greg Snow
You did not say what analysis you want to do, but many common analyses can be done as special cases of regression models and you can use the biglm package to do regression models. Here is an example that worked for me to get the mean and standard deviation by day from an oracle database with over

[R] fracdiff

2006-07-19 Thread Melissa Ann Haltuch
Hi, I'm using the function fracdiff and can not figure out how to get the estimated values for sigma2 or confidence intervals for the parameter estimates. Does anyone know how to obtain these values? Thanks, Melissa __ R-help@stat.math.ethz.ch mailing

Re: [R] Fitting a distribution to peaks in histogram

2006-07-19 Thread hadley wickham
> Can you be a bit more excact? I a biologist and relatively new to R In that case, I would _strongly_ advise that you get advice from a local statistician. > I am measureing the amount of DNA in cells, and I need to know the > percentage of cells in a part of the cell cycle; that the percentage

Re: [R] mgcv::gam error message

2006-07-19 Thread Simon Wood
> Could anyone please tell me what to do to resolve this error message? - Any chance you could send me some code (and associated data) that I can cut and paste in order to chase this up? Whatever's wrong here, `gam' should probably pick it up before things get this far, but I can't tell what's

Re: [R] gamm

2006-07-19 Thread Simon Wood
Sorry for the delay replying: I was on holiday, but have foolishly come back. > I am a bit confused about gamm in mgcv. Consulting Wood (2006) or > Ruppert et al. (2003) hasn't taken away my confusion. > > In this code from the gamm help file: > > b2<-gamm(y~s(x0)+s(x1)+s(x2)+s(x3),family=

Re: [R] How to find S4 generics?

2006-07-19 Thread Prof Brian Ripley
On Wed, 19 Jul 2006, Spencer Graves wrote: > Dear Prof. Ripley: > > Thanks very much. > > Am I correct then that the 'methods' function could, at least > theoretically, be revised so methods(class=...) could identify both S3 and S4 > methods (ignoring inheritance, as it does now, I belie

Re: [R] How to find S4 generics?

2006-07-19 Thread Spencer Graves
Dear Prof. Ripley: Thanks very much. Am I correct then that the 'methods' function could, at least theoretically, be revised so methods(class=...) could identify both S3 and S4 methods (ignoring inheritance, as it does now, I believe)? I ask, because it's not always obvious

Re: [R] Reproducible Research - Examples

2006-07-19 Thread Jean lobry
> Question: Reproducible research is clearly desirable and feasible. Could > anyone provide examples (stand-alone URLs or supplementary material in > journals) that you would recommend as a models for reproducible research in > R? Here are some attempts based on RWeb, contributed by Jeff Banfi

Re: [R] plain shading (not residuals) in mosaic plot

2006-07-19 Thread Muhammad Subianto
Maybe like this: mosaic(allmorph, direction = "v", pop = FALSE, gp=gpar(fill=c(grey(0.8),grey(0.4 Best, Muhammad Subianto On 7/19/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > If you look at ?mosaic the ... argument says it gets passed to strucplot and > looking at ?strucplot we see it

Re: [R] Fitting a distribution to peaks in histogram

2006-07-19 Thread Berton Gunter
With this much data, I think it makes more sense to fit a nonparametric density estimate. ?density does this via a kernel density procedure, but RSiteSearch('nonparametric density') will find many alternatives. The ash and mclust packages are two that come to mind, but there are certainly others.

Re: [R] plain shading (not residuals) in mosaic plot

2006-07-19 Thread Muhammad Subianto
Maybe like this: mosaic(allmorph, direction = "v", pop = FALSE, gp=gpar(fill=c(grey(0.8),grey(0.4 Best, Muhammad Subianto On 7/19/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > If you look at ?mosaic the ... argument says it gets passed to strucplot and > looking at ?strucplot we see it

[R] Answer found: Re: How to write a function in a graph

2006-07-19 Thread junguo liu
Dear Thomas and Simon, Thank you for your messages. I have found the best solution based on your advice. FAQ 7.13 does provide an effective way to answer my question. Kind regards, Junguo Thomas Lumley <[EMAIL PROTECTED]> wrote: On Tue, 18 Jul 2006, junguo liu wrote: > Dear R-ers

Re: [R] Stirling numbers

2006-07-19 Thread davidr
Well, given Martin's 2nd kind function and the fact that the 1st and 2nd kind matrices are inverses, you can get at least some of what you want by: ... # fill a matrix S2 with second kind numbers and zeros > S2 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,]1000000 [2,]

Re: [R] Fitting a distribution to peaks in histogram

2006-07-19 Thread Ulrik Stervbo
On 7/19/06, hadley wickham <[EMAIL PROTECTED]> wrote: > > I would like to fit a distribution to each of the peaks in a histogram, such > > as this: > > http://photos1.blogger.com/blogger/7029/2724/1600/DU145-Bax3-Bcl-xL.png > > As a first shot, I'd try fitting a mixture of gamma distributions (say

Re: [R] trellis.focus with postscript device

2006-07-19 Thread Soukup, Mat
Thanks Deepayan. Adding the argument highlight=FALSE to each trellis.focus() call worked marvelously. Cheers, Mat -Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Wednesday, July 19, 2006 11:58 AM To: Soukup, Mat Cc: r-help@stat.math.ethz.ch Subject: Re: [R] tr

Re: [R] WLS ins systemfit question

2006-07-19 Thread Arne Henningsen
On Wednesday 19 July 2006 16:22, Benn Fine wrote: > How does one specify the weights for WLS in the > systemfit command ? > > That is, there is a weight option in lm(), but there > doesn't seem to be weight option for systemfit("WLS") WLS in systemfit means that the different _equations_ (not the

Re: [R] Fitting a distribution to peaks in histogram

2006-07-19 Thread hadley wickham
> I would like to fit a distribution to each of the peaks in a histogram, such > as this: > http://photos1.blogger.com/blogger/7029/2724/1600/DU145-Bax3-Bcl-xL.png As a first shot, I'd try fitting a mixture of gamma distributions (say 3), plus a constant term for the highest bin. You could do th

Re: [R] trellis.focus with postscript device

2006-07-19 Thread Richard M. Heiberger
It works fine as long as you give it the layout. I used xyplot(yy ~ xx | pp, groups=gg, layout=c(2,2)) This is needed to make sure that the rows and columns of the trellis that you reference are consistent with the layout that xyplot has chosen. Another way to handle this is to change the facto

Re: [R] Sweave and multipage lattice

2006-07-19 Thread Sundar Dorai-Raj
Dieter Menne wrote: > Dear R-Listeners, > > as the Sweave faq says: > > http://www.ci.tuwien.ac.at/~leisch/Sweave/FAQ.html > > creating several figures from one figure chunk does not work, and for > standard graphics, a workaround is given. Now I have a multipage trellis > plot with an a-prior

Re: [R] Random structure of nested design in lme

2006-07-19 Thread ESCHEN Rene
Although I know it's not correct, this is what I tried in lme: anova(lme(NA.1~soiltype*habitat,random=~1|destination/habitat/origin/soiltype)) # numDF denDF F-value p-value #(Intercept) 1 130 12.136195 0.0007 #soiltype 1 130 15.099792 0.0002 #habitat

Re: [R] trellis.focus with postscript device

2006-07-19 Thread Deepayan Sarkar
On 7/19/06, Soukup, Mat <[EMAIL PROTECTED]> wrote: > Hello. > > First: R 2.3.1 on Windows XP. > > I am trying to add information (sample size) to the Trellis strips which > I am successful using the trellis.focus function with the default > Windows device. However, I typically use the postscript de

[R] how to use large data set ?

2006-07-19 Thread Yohan CHOUKROUN
Hello R users, Sorry for my English, i'm French. I want to use a large dataset (3 millions of rows and 70 var) but I don't know how to do because my computer crash quickly (P4 2.8Ghz, 1Go ). I have also a bi Xeon with 2Go so I want to do computation on this computer and show the results on

[R] Fwd: Classification error rate increased by bagging - any ideas?

2006-07-19 Thread Anthony Staines
Hi, I'm analysing some anthropometric data on fifty odd skull bases. We know the gender of each skull, and we are trying to develop a predictor to identify the sex of unknown skulls. Rpart with cross-validation produces two models - one of which predicts gender for Males well, and Females poorly

Re: [R] plain shading (not residuals) in mosaic plot

2006-07-19 Thread Gabor Grothendieck
If you look at ?mosaic the ... argument says it gets passed to strucplot and looking at ?strucplot we see it accepts a gp= arg so try this (same as your plus gp= arg): cols <- c(grey(0.8),grey(0.4)) mosaic(allmorph, direction = "v", pop = FALSE, gp = list(col = cols)) On 7/19/06, Kie Zuraw <[EMAI

[R] trellis.focus with postscript device

2006-07-19 Thread Soukup, Mat
Hello. First: R 2.3.1 on Windows XP. I am trying to add information (sample size) to the Trellis strips which I am successful using the trellis.focus function with the default Windows device. However, I typically use the postscript device as I use LaTeX and \includegraphic for incorporating graph

[R] connection network - influence of site

2006-07-19 Thread Guillaume Blanchet
Dear R users, I'm trying to construct a distance matrix based on a nb object created by spdep where sites must have a larger influence in one direction then in the other. Here is an example to better illustrate what I need: Let say I have the following Gabriel connection network library(spde

Re: [R] How would you export a 3-dimensional array to an SQL database?

2006-07-19 Thread Neuro LeSuperHéros
Here's a way to flatten and unflatten your array. a <-array(1:24, 2:4) d <-dim(a) flattened <-cbind(expand.grid(1:d[1],1:d[2],1:d[3]),as.vector(a)) unflattened <-array(flattened[,4],d) identical(a,unflattened) a==unflattened NeuroRox >From: "Lapointe, Pierre" <[EMAIL PROTECTED]> >To: "'r-hel

[R] RE : Re: Progress in a loop

2006-07-19 Thread Florent Bresson
Thanks, it is exactly what I was looking for. Florent Bresson --- Gavin Simpson <[EMAIL PROTECTED]> a écrit : > On Wed, 2006-07-19 at 16:12 +0200, Florent Bresson > wrote: > > Hi, I have to use a loop to perform a quite > computer > > intensive estimation and I would like to know the > > progres

Re: [R] Problem with ordered logistic regression using polr function.

2006-07-19 Thread Thomas Lumley
On Wed, 19 Jul 2006, Debarchana Ghosh wrote: > Hi, > I'm trying to fit a ordered logistic regression. The response variable > (y) has three levels (0,1,2). > The command I've used is: > > /ordlog<-polr(y~x1+x2+x3+x4, data=finalbase, subset=heard, weight=wt, > na.action=na.omit) > / > (There are n

Re: [R] R and DDE (Dynamic Data Exchange)

2006-07-19 Thread Richard M. Heiberger
Gabor, Thank you for that information. The reason for choosing DDE is that GNU emacs doesn't speak COM. Every year or so I ask on the emacs-for-windows list, and every year I get the answer no. DDE is the technology that I have been using in ESS for sending information to the S-Plus Commands win

[R] plain shading (not residuals) in mosaic plot

2006-07-19 Thread Kie Zuraw
Hello. I've been using R for a couple of months and enjoying it a lot. This is my first post to R-help. I'm using the vcd package to make mosaic plots with labels on the tiles indicating the number of items in each cell. For example, I've made this plot: > allmorph<-structure(c(10, 26, 17, 10

Re: [R] Progress in a loop

2006-07-19 Thread Gavin Simpson
On Wed, 2006-07-19 at 16:12 +0200, Florent Bresson wrote: > Hi, I have to use a loop to perform a quite computer > intensive estimation and I would like to know the > progress of the loop during the process. I tried to > include something like > > print(paste(k,date(),sep=" : ")) > > where k is

Re: [R] Progress in a loop

2006-07-19 Thread Doran, Harold
Look for progress() in the svMisc package > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Florent Bresson > Sent: Wednesday, July 19, 2006 10:13 AM > To: r-help@stat.math.ethz.ch > Subject: [R] Progress in a loop > > Hi, I have to use a loop to per

Re: [R] Progress in a loop

2006-07-19 Thread Uwe Ligges
Florent Bresson wrote: > Hi, I have to use a loop to perform a quite computer > intensive estimation and I would like to know the > progress of the loop during the process. I tried to > include something like > > print(paste(k,date(),sep=" : ")) > > where k is the number of the iteration, but t

[R] WLS ins systemfit question

2006-07-19 Thread Benn Fine
How does one specify the weights for WLS in the systemfit command ? That is, there is a weight option in lm(), but there doesn't seem to be weight option for systemfit("WLS") Thanks! __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] Stirling numbers

2006-07-19 Thread Martin Maechler
> "Robin" == Robin Hankin <[EMAIL PROTECTED]> > on Wed, 19 Jul 2006 14:25:21 +0100 writes: Robin> Hi anyone coded up Stirling numbers in R? "Sure" ;-) Robin> [I need unsigned Stirling numbers of the first kind] but with my quick search, I can only see those for which I had "

[R] Progress in a loop

2006-07-19 Thread Florent Bresson
Hi, I have to use a loop to perform a quite computer intensive estimation and I would like to know the progress of the loop during the process. I tried to include something like print(paste(k,date(),sep=" : ")) where k is the number of the iteration, but the result appears only at the end of th

Re: [R] R and DDE (Dynamic Data Exchange)

2006-07-19 Thread Gabor Grothendieck
On 7/19/06, Philippe Grosjean <[EMAIL PROTECTED]> wrote: > Richard M. Heiberger wrote: > > I am thrilled to learn tcltk2 has DDE capability. > > It is the piece I have been needing to make ESS work directly > > with the RGUI on Windows. GNU emacs on Windows has a ddeclient, > > but no access to CO

[R] Output and Word

2006-07-19 Thread Kuhn, Max
On page 2 of the manual, Section 2 ("Requirements") has "The package also requires a utility to zip and unzip compressed files, such as unzip, Winzip or jar." With a footnote that unzip can be found at http://www.info-zip.org/ for free. You can use any utility that can zip/unzip pkzip fo

[R] Stirling numbers

2006-07-19 Thread Robin Hankin
Hi anyone coded up Stirling numbers in R? [I need unsigned Stirling numbers of the first kind] cheers Robin -- Robin Hankin Uncertainty Analyst National Oceanography Centre, Southampton European Way, Southampton SO14 3ZH, UK tel 023-8059-7743 ___

Re: [R] Fitting a distribution to peaks in histogram

2006-07-19 Thread Petr Pikal
Hi There are some packages for mass spectra processing (spectrino, caMassClass). I did not use them so I do not know how they suit your needs. However you can compute area (integrate) by these functions # uses information interactively from plot(x,y) # first it replots data between corners *re

[R] How would you export a 3-dimensional array to an SQL database?

2006-07-19 Thread Lapointe, Pierre
Hello, How would you export a 3-dimensional array to an SQL database? a<-array(1:24, 2:4) Is there an open source DB that would be more adequate for this type of operation? Is there a way to reshape/flatten a 3-dimensional array? Regards, Pierre Lapointe *

Re: [R] Test for equality of coefficients in multivariate multipleregression

2006-07-19 Thread Berwin A Turlach
G'day John, > "JF" == John Fox <[EMAIL PROTECTED]> writes: JF> Simply stacking the problems and treating the resulting JF> observations as independent will give you the correct JF> coefficients, but incorrect coefficient variances Yes, after Andrew's (off-list) answer I realised t

Re: [R] Test for equality of coefficients in multivariate multipleregression

2006-07-19 Thread John Fox
Dear Berwin, Simply stacking the problems and treating the resulting observations as independent will give you the correct coefficients, but incorrect coefficient variances and artificially zero covariances. The approach that I suggested originally -- testing a linear hypothesis using the coeffic

Re: [R] Random structure of nested design in lme

2006-07-19 Thread Doran, Harold
Can you provide an example of what you have done with lme so we might be able to evaluate the issue? > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of ESCHEN Rene > Sent: Wednesday, July 19, 2006 7:37 AM > To: r-help@stat.math.ethz.ch > Subject: [R]

[R] Random structure of nested design in lme

2006-07-19 Thread ESCHEN Rene
All, I'm trying to analyze the results of a reciprocal transplant experiment using lme(). While I get the error-term right in aov(), in lme() it appears impossible to get as expected. I would be greatful for any help. My experiment aimed to identify whether two fixed factors (habitat type and

Re: [R] Plotting lines and points on the second plot when usinggap.plot in plotrix

2006-07-19 Thread michael watson \(IAH-C\)
I've answered my own question, the x-values of the line/points I want to plot must be adjusted by the gap size Thanks Mick -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of michael watson (IAH-C) Sent: 19 July 2006 11:58 To: r-help@stat.math.ethz.ch Subjec

[R] Plotting lines and points on the second plot when using gap.plot in plotrix

2006-07-19 Thread michael watson \(IAH-C\)
Hi My question is simple - the gap.plot function in the plotrix package allows users to draw graphs that have a broken axis. However, I want to then add a line to the "second" plot, but can't. Eg: twogrp<-c(rnorm(10)+4,rnorm(10)+20) gap.plot(twogrp,rnorm(20),gap.bounds=c(8,16),gap.axis="x",xlab

Re: [R] Aligning ragged text columns

2006-07-19 Thread hadley wickham
> I find that things line up better in data.frames That's a good idea, although I was hoping there would be something in R to do it for me. I have ended up with: fwidth <- max(nchar(x[,1])) descs <- strwrap(x[,2], width=width - fwidth - 5, simplify=FALSE) output <- do.ca

Re: [R] conditional plot

2006-07-19 Thread jim holtman
does this do what you want? x <- data.frame(company=LETTERS, spec1=runif(26), spec2=runif(26)) plot(range(x$spec1), range(x$spec2), type='n') abline(v=.5, h=.5) text(x$spec1, x$spec2, labels=x$company) On 7/19/06, Manoj <[EMAIL PROTECTED]> wrote: > > Hi, > Can anyone pls help me in plotting

Re: [R] Test for equality of coefficients in multivariate multiple regression

2006-07-19 Thread Ulrich Keller
Hello and thank you for your answers, Andrew and Berwin. If I'm not mistaken, the mixed-model version of Berwin's approach would be: #My stuff: DF<-data.frame(x1=rep(c(0,1),each=50),x2=rep(c(0,1),50)) tmp<-rnorm(100) DF$y1<-tmp+DF$x1*.5+DF$x2*.3+rnorm(100,0,.5) DF$y2<-tmp+DF$x1*.5+DF$x2*.7+rnorm(

[R] Fitting a distribution to peaks in histogram

2006-07-19 Thread Ulrik Stervbo
Hello list! I would like to fit a distribution to each of the peaks in a histogram, such as this: http://photos1.blogger.com/blogger/7029/2724/1600/DU145-Bax3-Bcl-xL.png . The peaks are identified using Petr Pikal peaks function ( http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html), but a

Re: [R] Output and Word

2006-07-19 Thread Prof Brian Ripley
On Wed, 19 Jul 2006, David Hajage wrote: > thank you Greg Snow for this information ! > > But I have this message : > > > odfWeave("c:/simple.odt", "c:/essai.odt") > Setting wd > Copying c:/simple.odt > Decompressing ODF file using unzip -o > "C:\DOCUME~1\Maud\LOCALS~1\Temp\RtmpF0hdqb/sim

[R] Aligning ragged text columns

2006-07-19 Thread ken knoblauch
Just a follow-up, by giving names with variable amounts of white space, I could make the column and rownames invisible (although this is probably abusing the mechanisms, however, this is only for display purposes, I suppose). dd <- data.frame(c1 = c("Heading 1", "", "Heading 2", ""), c2 = c("T

[R] Aligning ragged text columns

2006-07-19 Thread ken knoblauch
Hi Hadley, I find that things line up better in data.frames data.frame(c1 = c("Heading 1", "", "Heading 2", ""), + c2 = c("This is some info", "about heading 1", "This is some info", "about heading ")) c1c2 1 Heading 1 This is some info 2 about heading 1 3

Re: [R] Aligning ragged text columns

2006-07-19 Thread hadley wickham
> heading1 <- "Heading1" > heading2 <- "Heading2" > > a <- matrix(c( > "Heading 1", paste("This is some info\nabout", heading1, sep=""), > "Heading 2", paste("This is some info\nabout", heading2, sep=""), > ), byrow=T, nrow=2) I wasn't so concerned about the redundancy in my exam

Re: [R] Aligning ragged text columns

2006-07-19 Thread John Wiedenhoeft
Am Mittwoch, den 19.07.2006, 10:03 +0100 schrieb hadley wickham: > Can anyone please suggest how I can print: > > a <- matrix(c( > "Heading 1", "This is some info\nabout heading 1", > "Heading 2", "This is some info\nabout heading 2", > ), byrow=T, nrow=2) > > to look like: > > Hea

Re: [R] conditional plot

2006-07-19 Thread Chuck Cleland
Manoj wrote: > Hi, > Can anyone pls help me in plotting the following data? > > The data-set contains company name, specification-1, specification-2. > > The graph would basically plot company name with specification-1 > on x-axis, & specification-2 on y-axis. > > Simply put

[R] Aligning ragged text columns

2006-07-19 Thread hadley wickham
Can anyone please suggest how I can print: a <- matrix(c( "Heading 1", "This is some info\nabout heading 1", "Heading 2", "This is some info\nabout heading 2", ), byrow=T, nrow=2) to look like: Heading 1 This is some info about heading 1 Heading 2 This is some info

Re: [R] Output and Word

2006-07-19 Thread David Hajage
thank you Greg Snow for this information ! But I have this message : > odfWeave("c:/simple.odt", "c:/essai.odt") Setting wd Copying c:/simple.odt Decompressing ODF file using unzip -o "C:\DOCUME~1\Maud\LOCALS~1\Temp\RtmpF0hdqb/simple.odt" Erreur dans odfWeave("c:/simple.odt", "c:/essai.odt

Re: [R] get rid of error in Factor Analysis

2006-07-19 Thread Prof Brian Ripley
?try On Wed, 19 Jul 2006, Robert Mcfadden wrote: > Dear All, > > I wrote a program and there is a loop. At each iteration I use maximum > likelihood factor analysis (?factanal). Output of factor analysis I use > later (in this loop). Unfortunately from time to time I get an error message > (I pa

[R] get rid of error in Factor Analysis

2006-07-19 Thread Robert Mcfadden
Dear All, I wrote a program and there is a loop. At each iteration I use maximum likelihood factor analysis (?factanal). Output of factor analysis I use later (in this loop). Unfortunately from time to time I get an error message (I paste it below) and everything is stopped. Is it possible to writ

Re: [R] A contingency table of counts by case

2006-07-19 Thread Jacques VESLOT
sorry, answered to quickly... actually it's easier using paste(): df <- df[order(df$case),] apply(combinations(9,2), 1, function(y) table(factor(do.call(paste, c(with(df[df$id %in% y, ], split(x, id)), sep="")), levels=c("00","01","10","11" -

Re: [R] R and DDE (Dynamic Data Exchange)

2006-07-19 Thread Philippe Grosjean
Gabor Grothendieck wrote: > You can access DDE via COM as in this example > which uses DDE to open an Excel file. Note that > Excel also supports COM directly and normally > one would use COM with Excel, not DDE, so you > might check if your application also supports COM. > > # opens an excel spr

Re: [R] R and DDE (Dynamic Data Exchange)

2006-07-19 Thread Philippe Grosjean
Richard M. Heiberger wrote: > I am thrilled to learn tcltk2 has DDE capability. > It is the piece I have been needing to make ESS work directly > with the RGUI on Windows. GNU emacs on Windows has a ddeclient, > but no access to COM. So if R, or tcltk2 talking in both directions to R, > has a dde

Re: [R] how can I delete rows?

2006-07-19 Thread Philipp Pagel
On Tue, Jul 18, 2006 at 12:25:47PM -0300, raul sanchez wrote: > I have the Barro-Lee data set which contains 98 countries and I want > to run the regressions only for the Latin America countries, so what > do you recomend? How can I delete all the other countries or how can > I select the c

[R] [R-pkgs] odfWeave Package

2006-07-19 Thread Kuhn, Max
The odfWeave package is now available on CRAN at http://lib.stat.cmu.edu/R/CRAN/src/contrib/Descriptions/odfWeave.html and your local mirror. The package extends Sweave to Open Document Format (ODF) text document files. Latex-style code chunks and in-line Sexpr commands can be used to embed R

Re: [R] voronoi tessellations

2006-07-19 Thread Prof Brian Ripley
On Tue, 18 Jul 2006, Don MacQueen wrote: > I'll suggest going to the CRAN packages page and doing a search for "voronoi". The problem here is that `Voronoi tessellation' is a secondary name. The concept has many names, including Dirichlet tessellation and Thiessen polygons, and Dirichlet has p

Re: [R] R & Visual Basic

2006-07-19 Thread Prof Brian Ripley
On Wed, 19 Jul 2006, Vito Ricci wrote: > Der R-UserRs, > > I need a little help, I wish to know if exists a way > to use R in Visual Basic environment, in creation of > VB applications, embedding R in VB. Is it possible? > Can anyone help me? Is there something such as books, > articles, other av

[R] Classification error rate increased by bagging - any ideas?

2006-07-19 Thread Anthony Staines
Hi, I'm analysing some anthropometric data on fifty odd skull bases. We know the gender of each skull, and we are trying to develop a predictor to identify the sex of unknown skulls. Rpart with cross-validation produces two models - one of which predicts gender for Males well, and Females poorl

[R] how can I delete rows?

2006-07-19 Thread raul sanchez
Hello, I am very new in R so I am so sorry for this question. I have the Barro-Lee data set which contains 98 countries and I want to run the regressions only for the Latin America countries, so what do you recomend? How can I delete all the other countries or how can I select the countries of

Re: [R] How to find S4 generics? (was: inames() function and lmer())

2006-07-19 Thread Prof Brian Ripley
On Tue, 18 Jul 2006, Spencer Graves wrote: > * > * "methods" * > * > You have asked an excellent question. I can provide a partial answer > below. First, however, I wish to pose a question of my own, which could > help answer your question: > > How can one obta