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

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

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

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 this too.

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

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] Fitting a distribution to peaks in histogram

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

[R] Stirling numbers

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

[R] Output and Word

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

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

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 COM. So if

[R] Progress in a loop

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 the

Re: [R] Stirling numbers

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 2nd kind :

[R] WLS ins systemfit question

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

Re: [R] Progress in a loop

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 the result

Re: [R] Progress in a loop

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 perform a

Re: [R] Progress in a loop

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 the

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

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,

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

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

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

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 no NA's in

[R] RE : Re: Progress in a loop

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 progress of the

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

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] connection network - influence of site

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

[R] trellis.focus with postscript device

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

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

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 [EMAIL

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

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

[R] how to use large data set ?

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

Re: [R] trellis.focus with postscript device

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 device as

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

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] Sweave and multipage lattice

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-priori

Re: [R] trellis.focus with postscript device

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

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

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 this

Re: [R] WLS ins systemfit question

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] trellis.focus with postscript device

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]

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

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 3),

Re: [R] Stirling numbers

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,]

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

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] plain shading (not residuals) in mosaic plot

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

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

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] Reproducible Research - Examples

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

Re: [R] How to find S4 generics?

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

Re: [R] How to find S4 generics?

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 believe)?

Re: [R] gamm

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:

Re: [R] mgcv::gam error message

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] Fitting a distribution to peaks in histogram

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

[R] fracdiff

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] how to use large data set ?

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] Wrap a loop inside a function

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,

Re: [R] Wrap a loop inside a function

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],'-')))/

Re: [R] voronoi tessellations

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 -

Re: [R] Wrap a loop inside a function

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] Correlation Mapping

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 cover of

Re: [R] voronoi tessellations

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. example

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

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

Re: [R] Wrap a loop inside a function

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 a

Re: [R] voronoi tessellations

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.

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

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.

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

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

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

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

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

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

Re: [R] How to find S4 generics?

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

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

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

Re: [R] How to find S4 generics?

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=...)

Re: [R] Wrap a loop inside a function

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] error when compiling stats library in R-2.3.1 on Solaris x86

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 updating

[R] Automating package building packages and repository uploading

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

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

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

Re: [R] Sweave and multipage lattice

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 trellis

[R] Rearrange data.

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 funny -

Re: [R] Rearrange data.

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, a12,

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

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. PARA ENTRAR

Re: [R] Plot fit of a generic function

Yes, that will do. Thank you Gabor! Gabor Grothendieck wrote: You could plot y vs. fitted(y.lm) where y.lm is the output of lm or plot both y and fitted(y.lm) against x on the same chart. On 7/18/06, Gregor Gorjanc [EMAIL PROTECTED] wrote: Hello! Say I have a function, which creates a

[R] R Visual Basic

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 available on the web? Thanks in advance. Regards. Vito

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

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 obtain a

[R] how can I delete rows?

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

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

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

Re: [R] R Visual Basic

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

Re: [R] voronoi tessellations

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

[R] [R-pkgs] odfWeave Package

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] how can I delete rows?

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

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

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

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

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

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

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

[R] get rid of error in Factor Analysis

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

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

?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 paste

Re: [R] Output and Word

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) :

[R] Aligning ragged text columns

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] conditional plot

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 - company

Re: [R] Aligning ragged text columns

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: Heading 1 This is

Re: [R] Aligning ragged text columns

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 example, but how it looks -

[R] Aligning ragged text columns

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 Heading 2 This

[R] Aligning ragged text columns

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(This is

Re: [R] Output and Word

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/simple.odt

[R] Fitting a distribution to peaks in histogram

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

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

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)

Re: [R] conditional plot

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 the

Re: [R] Aligning ragged text columns

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 -

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

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=X

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

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

[R] Random structure of nested design in lme

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] Random structure of nested design in lme

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]