Re: [R] How to increase memory for R on Soliars 10 with 16GB and 64bit R

Dongseok Choi wrote: Thank you very much for your help!! Now, it runs without any problem. Is it going to be fixed in the next release? Of course, Brian Thanks again, Dongseok Dongseok Choi, Ph.D. Assistant Professor Division of Biostatistics Department of Public Health

Re: [R] How to increase memory for R on Soliars 10 with 16GB and 64bit R

Uwe Ligges wrote: Dongseok Choi wrote: Thank you very much for your help!! Now, it runs without any problem. Is it going to be fixed in the next release? Of course, Brian [hmmm, looks like some wrong shortcut has been used - and it must have been me who forgot to drink coffee before

Re: [R] texture in barplots?

On Thursday 14 July 2005 00:51, Peter Dalgaard wrote: Adrian Dusa [EMAIL PROTECTED] writes: ...snip... This comes up every now and then, and while it seems that everyone thinks fill patterns would be nice to have, I suspect that every attempt to actually implement it have gotten killed in

[R] Coxph with factors

Hello, I am fitting a coxph model with factors. I am running into problems when using 'survfit'. I am unsure how R is treating the factors when I fit, say: DATA-data.frame(time.sec,done,f.pom=factor(f.pom),po,vo) final-coxph(Surv(time.sec,done)~f.pom*vo+po,data=DATA)

Re: [R] Efficient testing for +ve definiteness

On Wed, 13 Jul 2005, Makram Talih wrote: Dear R-users, Is there a preferred method for testing whether a real symmetric matrix is positive definite? [modulo machine rounding errors.] The obvious way of computing eigenvalues via E - eigen(A, symmetric=T, only.values=T)$values and returning

Re: [R] High resolution plots

Thanks for the pointer! .wmf is far superior, I was just in the dark about the format and R's ability to produce it (An Introduction to R Device drivers does not mention it and I had obviously missed the deciding last two words in ?device 'windows') Thanks again. -Luis Gabor Grothendieck

[R] Pearson dispersion statistic

Thank you for your reply. I am aware of the good reasons not to use the deviance estimate in binomial, Poisson, and gamma families. However, for the inverse Gaussian, the choice seems to me less clear cut. So I just wanted to compare two different options. I have used the dispersion

Re: [R] Name for factor's levels with contr.sum

One way to do this generally is to make a copy of contr.sum, rename it, and set the dimnames appropriately. I think contr.treatment is misleading (it labels contrasts of two levels by just one of them), and Christoph's labels are informative but impractically long. But if you want to label

Re: [R] High resolution plots

Peter Dalgaard schrieb: What is the real issue here? Import trouble? If you're importing to Word/PowerPoint, why not use the Windows metafile? Perhaps they are too ugly compared to EPS by your taste? My reason was http://www.adobe.com/support/techdocs/328541.html I tried to plot emf files

[R] Ellipse Cosrrespondence Analysis(Internal Stability)

I'am a student of Educational Universiti of Indonesia. I take Internal Stabiliti in Correespondence Analysis of Normal Distributerd Data for my final task. I've had problem with ellipse contructing as STATISTICA 6.0 generate in Correspondence Analysis Session. I thing the ellipse in bivariate

Re: [R] High resolution plots

Knut Krueger [EMAIL PROTECTED] writes: Peter Dalgaard schrieb: What is the real issue here? Import trouble? If you're importing to Word/PowerPoint, why not use the Windows metafile? Perhaps they are too ugly compared to EPS by your taste? My reason was

Re: [R] High resolution plots

Peter Dalgaard schrieb: Knut Krueger [EMAIL PROTECTED] writes: Peter Dalgaard schrieb: What is the real issue here? Import trouble? If you're importing to Word/PowerPoint, why not use the Windows metafile? Perhaps they are too ugly compared to EPS by your taste? My reason

Re: [R] Memory question

Thank you Dr. Spencer Graves for your answer. What kind of matrices? They come form an image of about 3000x5000, and I need to generate arround 1024 matrices of the same size, they are not sparse matrices. What function can I use to, once generated one matrix, I could save into disk and

Re: [R] Name for factor's levels with contr.sum

This is exactly what I've been looking for (without success) when I was speaking about a more elegant and general solution. I agree with your argument that labels might be misleading. Nevertheless if a user is aware what contr.sum calculates, it is practical to have an annotation. Thank you

Re: [R] unexpected par('pin') behaviour

joerg == joerg van den hoff [EMAIL PROTECTED] on Wed, 13 Jul 2005 16:00:58 +0200 writes: joerg hi everybody, joerg I noticed the following: in one of my scripts 'layout' is used to joerg generate a (approx. square) grid of variable dimensions (depending on joerg no. of

Re: [R] How to increase memory for R on Soliars 10 with 16GB and 64bit R

From: Uwe Ligges Uwe Ligges wrote: Dongseok Choi wrote: Thank you very much for your help!! Now, it runs without any problem. Is it going to be fixed in the next release? Of course, Brian [hmmm, looks like some wrong shortcut has been used - and it must have been

Re: [R] High resolution plots

Knut Krueger [EMAIL PROTECTED] writes: Peter Dalgaard schrieb: Knut Krueger [EMAIL PROTECTED] writes: Peter Dalgaard schrieb: What is the real issue here? Import trouble? If you're importing to Word/PowerPoint, why not use the Windows metafile? Perhaps they are too ugly compared

[R] memory problem

I'm a beginner in R and, therefore, I don't know how serious my trouble is. After running a script: ** *t**-c(14598417794,649693)* *data**=data.frame(read.spss(C:\\Ginters\\Kalibracija\\cal_data.sav))* *Xs=**as.matrix(data[,1:2]) *

Re: [R] Memory question

On Thu, 14 Jul 2005, Kenneth Cabrera wrote: Thank you Dr. Spencer Graves for your answer. What kind of matrices? They come form an image of about 3000x5000, and I need to generate arround 1024 matrices of the same size, they are not sparse matrices. What function can I use to, once

Re: [R] memory problem

On 7/14/2005 7:19 AM, Ginters wrote: I'm a beginner in R and, therefore, I don't know how serious my trouble is. After running a script: ** *t**-c(14598417794,649693)* *data**=data.frame(read.spss(C:\\Ginters\\Kalibracija\\cal_data.sav))* *Xs=**as.matrix(data[,1:2]) *

[R] help: how to change the size of a window after it has been created

Hello, I wish to plot some figures in a window in turn, but the size of these figures is different, so how can I change the size of the window by resetting the parameters before each plotting? Thank you, Shengzhe __ R-help@stat.math.ethz.ch mailing

[R] Plotting greek symbols in plot titles, labels, etc.

Hi, Can anyone please tell me how to use the escape sequences (e.g. \\*a for alfa) for plotting Greek symbols in for example plot titles? I have read the reference manual, but can't seem to get it to work. Thanx! Roy -- The information contained in this communication and any

Re: [R] Plotting greek symbols in plot titles, labels, etc.

Take a look at '?plotmath'. -roger Roy Werkman wrote: Hi, Can anyone please tell me how to use the escape sequences (e.g. \\*a for alfa) for plotting Greek symbols in for example plot titles? I have read the reference manual, but can't seem to get it to work. Thanx! Roy --

Re: [R] help: how to change the size of a window after it has been created

Not sure you can do this. You might have to launch separate graphics windows. -roger wu sz wrote: Hello, I wish to plot some figures in a window in turn, but the size of these figures is different, so how can I change the size of the window by resetting the parameters before each

Re: [R] help: how to change the size of a window after it has been created

On Thu, 2005-07-14 at 14:38 +0200, wu sz wrote: Hello, I wish to plot some figures in a window in turn, but the size of these figures is different, so how can I change the size of the window by resetting the parameters before each plotting? Thank you, Shengzhe Other than dragging a plot

Re: [R] help: how to change the size of a window after it has been created

Marc Schwartz wrote: On Thu, 2005-07-14 at 14:38 +0200, wu sz wrote: Hello, I wish to plot some figures in a window in turn, but the size of these figures is different, so how can I change the size of the window by resetting the parameters before each plotting? Thank you, Shengzhe Other

Re: [R] Fieller's Conf Limits and EC50's

Hi, I didn't verify your formulas for Fieller's method of computing the confidence interval. A slightly simpler approach is to use the Delta method to compute the CI. It is also valid for any link function. It yields a simpler formula for the variance of EC50 (for any link function): varEC50 -

[R] Keeping memory usage low for a big script

Hi all, I've got a script that I run on several computers, some with much less ram than others, and I would like to try and keep the memory usage as low as possible. The script creates variables to store intermediate results of calculations, and then at the end writes a text file with the

Re: [R] Where's iris?

-Original Message- From: Gabor Grothendieck [mailto:[EMAIL PROTECTED] Sent: 13 July 2005 20:30 To: Ruben Roa Cc: R-help@stat.math.ethz.ch Subject: Re: [R] Where's iris? On 7/13/05, Ruben Roa [EMAIL PROTECTED] wrote: Hi: Where is the iris data set actually located in the R

Re: [R] Keeping memory usage low for a big script

On 7/14/2005 9:46 AM, Ken Termiso wrote: Hi all, I've got a script that I run on several computers, some with much less ram than others, and I would like to try and keep the memory usage as low as possible. The script creates variables to store intermediate results of calculations, and

[R] Fwd: Re: Problem installing R packages

Hi, I am trying to install the R libraries rmutil and repeated on a Mac OS X version 10.4.1 (which has the latest version of the Mac Developer tools installed) and I am having trouble compiling the libraries. The error message I receive is as follows (I have only included the error message

Re: [R] Plotting greek symbols in plot titles, labels, etc.

Type demo(plotmath) at the R prompt. At 2:45 PM +0200 7/14/05, Roy Werkman wrote: Hi, Can anyone please tell me how to use the escape sequences (e.g. \\*a for alfa) for plotting Greek symbols in for example plot titles? I have read the reference manual, but can't seem to get it to work.

[R] read.xport

I am trying to import data from a SAS XPORT file that contains 24 SAS files. When I use the read.xport procedure only about 16 data frames (components) are created. Any suggestions? * Gary A. Nelson, Ph.D

Re: [R] maps drawing

There is also the package maptools if you want or need to read ESRI shapefiles. Regards, Christoph Buser -- Christoph Buser [EMAIL PROTECTED] Seminar fuer Statistik, LEO C13 ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND

Re: [R] maps drawing

GIS-related packages in R: ade4 -- Analysis of Environmental Data : Exploratory and Euclidean methods in Environmental sciences adehabitat -- Analysis of habitat selection by animals fields -- Tools for spatial data GRASS -- Interface between GRASS 5.0 geographical information system and R

Re: [R] unexpected par('pin') behaviour

Martin Maechler wrote: joerg == joerg van den hoff [EMAIL PROTECTED] on Wed, 13 Jul 2005 16:00:58 +0200 writes: joerg hi everybody, joerg I noticed the following: in one of my scripts 'layout' is used to joerg generate a (approx. square) grid of variable dimensions

Re: [R] memory problem

On Thu, 14 Jul 2005, Duncan Murdoch wrote: On 7/14/2005 7:19 AM, Ginters wrote: Why does memory need so much (1.6 GB) space? How can I enlarge it? Is it possible to allocate a part of memory used to the hard drive? Or, is the trouble only with my script? This sounds like a problem with the

Re: [R] Fwd: Re: Problem installing R packages

What Fortran compiler (g77) are you using, obtained from where? See http://www.astro.gla.ac.uk/users/norman/note/2004/restFP/ Your lines are wrapped in ways that make parsing impossible, but my guess at the link line does not include -lcc_dynamic. See also the list archives, e.g.

[R] DSC 2005 registration

Early registration discount for DSC 2005: Directions in Statistical Computing, August 13-14, in Seattle closes on July 18. Registration is online, from the conference web page at http://depts.washington.edu/dsc2005/ Titles and abstracts for presentations are at

Re: [R] Coxph with factors

Kylie-Anne Richards kar at itga.com.au writes: I am fitting a coxph model with factors. I am running into problems when using 'survfit'. I am unsure how R is treating the factors when I fit, say: DATA-data.frame(time.sec,done,f.pom=factor(f.pom),po,vo)

Re: [R] stripchart usage and alternatives

Mike R mike.rstat at gmail.com writes: . stripchart(r~u,vertical=TRUE,pch=21) stripchart(rm~levels(uf),vertical=TRUE,pch=3,add=TRUE) the above code creates a scatter plot of nominal data are there alternatives to generate the same or similar kind of figure? How about stripplot in

Re: [R] Fwd: Re: Problem installing R packages

We used the Fortran compiler that came from the customised option when installing from the mac dmg g77 --version GNU Fortran (GCC) 3.4.2 We will look into the links you provided. Thanks, John On Jul 14 2005, Prof Brian Ripley wrote: What Fortran compiler (g77) are you using, obtained

[R] New Package: TeachingDemos

This is to announce the new package TeachingDemos now available on CRAN. This is a collection of functions that I have used to demonstrate different ideas in various classes taught. Some of the first ones I wrote just to help me visualize what was happening, they latter were usefull in teaching

[R] A statistical modeling problem

Hi, This is not an R related question and I apologize for that, but given the brain power of the R community it is hard for me to resist posting this here. I have a problem where each participant is shown a series of visual cues (displayed on a computer screen in a random order) and asked

Re: [R] read.xport

How about avoiding SAS XPORT altogether and exporting everything in the simple, clean, non-proprietary, extremely reliable, platform-independent ... etc text format (CSV, tab delimited etc)? -Original Message- From: Nelson, Gary (FWE) [mailto:[EMAIL PROTECTED] Sent: Thursday, July

Re: [R] read.xport

bogdan romocea wrote: How about avoiding SAS XPORT altogether and exporting everything in the simple, clean, non-proprietary, extremely reliable, platform-independent ... etc text format (CSV, tab delimited etc)? I hope the problem is fixed in the latest version of foreign (no version info

Re: [R] plot the number of replicates at the same point

Thank you for thinking about the problem for me. However, I have found that your method doesn't work at all. You may test the following example: x1=c(0.6,0.4,.4,.4,.2,.2,.2,0,0) x2=c(0.4,.2,.4,.6,0,.2,.4,0,.2) x1=rep(x1,4) x2=rep(x2,4) temp=data.frame(x1,x2) temp1=table(temp)

[R] Calculate of data frame

Hello, I have a dataframe which I read from a file: df - readtable(myFile); The dataframe has 4 columns: model, count, value and date where model and date are made of charactors and count and value are numbers. An example is like: modelcountvaluedate A4

[R] Error running lme.

I am trying to fit lme using R 2.1.1 under Windows 2k. I am getting the following message noted below. Any suggestions that would help me correct my error would be greatly appreciated. Thanks, John fit1lme-lme(Velocity~time,data=gate) Error in getGroups.data.frame(dataMix, groups) :

Re: [R] read.xport

I have the latest version of foreign, but it still doesn't work. I quickly tried the Hmisc package, but the same issue arose. I will delve into the Hmisc package further. Thanks. -Original Message- From: Frank E Harrell Jr [mailto:[EMAIL PROTECTED] Sent: Thursday, July 14, 2005 11:46

[R] Partek has Dunn-Sidak Multiple Test Correction. Is this the same/similar to any of R's p.adjust.methods?

The Partek package (www.partek.com) allows only two selections for Multiple Test Correction: Bonferroni and Dunn-Sidak. Can anyone suggest why Partek implemented Dunn-Sidak and not the other methods that R has? Is there any particular advantage to the Dunn-Sidak method? R knows about these

[R] East Asian language

Dear all, I just installed R 2.1.1. The installation program automatically recognized my Windows XP was using Chinese language, so now my R console displays everything in Chinese. How can I still let R console display in English without modifying my Window XP language setup? Thank you so much!

[R] Does R have ANOVA permutation tests?

Hi list, Does anybody know if R has functions to do the ANOVA permutation test? I googled and found R has the vegan package to do ANOVA like permutation test for Constrained Correspondence Analysis. But does R have a function for general ANOVA-like permutation tests? Thanks in advance!

Re: [R] Calculate of data frame

Zhang, Fan [EMAIL PROTECTED] writes: Now I want to get the total count and value for each model/date pair, like modelcountvaluedate A631.87/1/2005 A310.27/2/2005 B714.27/1/2005 B10

Re: [R] plot the number of replicates at the same point

On 7/14/05, Kerry Bush [EMAIL PROTECTED] wrote: Thank you for thinking about the problem for me. However, I have found that your method doesn't work at all. You may test the following example: x1=c(0.6,0.4,.4,.4,.2,.2,.2,0,0) x2=c(0.4,.2,.4,.6,0,.2,.4,0,.2) x1=rep(x1,4) x2=rep(x2,4)

Re: [R] Error running lme.

John Your model is not properly specified for lme. You have not included any random effects or a grouping variable. Let me assume just for sake of argument that you want to include a random effect for the intercept and for time. Your lme specification would be fm1 - lme(Velocity~time,

Re: [R] East Asian language

Nan Lin wrote: Dear all, I just installed R 2.1.1. The installation program automatically recognized my Windows XP was using Chinese language, so now my R console displays everything in Chinese. How can I still let R console display in English without modifying my Window XP language

Re: [R] plot the number of replicates at the same point

On Thu, 2005-07-14 at 12:30 -0500, Deepayan Sarkar wrote: On 7/14/05, Kerry Bush [EMAIL PROTECTED] wrote: Thank you for thinking about the problem for me. However, I have found that your method doesn't work at all. You may test the following example:

[R] Fwd: Re: East Asian language

Original message Date: Thu, 14 Jul 2005 19:50:41 +0200 From: Uwe Ligges [EMAIL PROTECTED] Subject: Re: [R] East Asian language To: Nan Lin [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch Nan Lin wrote: Dear all, I just installed R 2.1.1. The installation program automatically

Re: [R] Efficient testing for +ve definiteness

To reinforce Prof. Ripley's comment that, Knowing the determinant does not tell you if the matrix is close to non-positive definite, note that the determinant of the negative of the identity matrix, (-diag(k)), is (-1)^k; if k is even, the determinant is positive. This silly

Re: [R] Fwd: Re: East Asian language

Baoqiang Cao wrote: Original message Date: Thu, 14 Jul 2005 19:50:41 +0200 From: Uwe Ligges [EMAIL PROTECTED] Subject: Re: [R] East Asian language To: Nan Lin [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch Nan Lin wrote: Dear all, I just installed R 2.1.1. The installation

[R] Using system to run a stand alone program that requires input in Windows

Under Windows I want to run a stand alone program that takes a number of commands from the user. I've been running the program from the command line using C:\Data\ foo.exe params.txt where foo is the program is params.txt is a text file with a few lines (9) of parameters. I want to run this

Re: [R] Using system to run a stand alone program that requires input in Windows

Dr Carbon wrote: Under Windows I want to run a stand alone program that takes a number of commands from the user. I've been running the program from the command line using C:\Data\ foo.exe params.txt where foo is the program is params.txt is a text file with a few lines (9) of

Re: [R] Calculate of data frame

Thanks for your answer, Peter. Now I have a new data frame which is sorted by (model,date) pair. Since the sorting is done to charactors, I have 7/1/2005 followed by 7/10/2005 instead of 7/2/2005. How can I do the aggregation and at the same time, make the result sorted the way I want (i.e,

Re: [R] Using system to run a stand alone program that requires input in Windows

system(foo.exe params.txt) Alas that hangs just opens a window with foo.exe; params.txt is not passed to it. That's why I thought input might work. Any other ideas. I'm trying to think of a DOS program that requires input to test this. __

Re: [R] Partek has Dunn-Sidak Multiple Test Correction. Is this the same/similar to any of R's p.adjust.methods?

Earl F. Glynn [EMAIL PROTECTED] writes: The Partek package (www.partek.com) allows only two selections for Multiple Test Correction: Bonferroni and Dunn-Sidak. Can anyone suggest why Partek implemented Dunn-Sidak and not the other methods that R has? Is there any particular advantage to

Re: [R] Using system to run a stand alone program that requires input in Windows

Well, I'm confused because I tried to cook up an exaple and it appears to work as I think it should. If I have a simple script named myInput.bat: @echo off echo Enter y set /p Input= if /i %Input%==y (goto Proceed) echo y not entered exit /b :Proceed echo y entered and a text file with just a y

Re: [R] read.xport

On Thu, 14 Jul 2005, Nelson, Gary (FWE) wrote: I have the latest version of foreign, but it still doesn't work. I Are you sure: a new version was released a few hours ago? It may not answer your question, but please do give actual version numbers (as the posting guide asks). -- Brian D.

Re: [R] Using system to run a stand alone program that requires input in Windows

You want shell() not system(), I suspect, as you specified a shell command and I gather you are using Windows (not `DOS'). On Thu, 14 Jul 2005, Dr Carbon wrote: system(foo.exe params.txt) Alas that hangs just opens a window with foo.exe; params.txt is not passed to it. That's why I

[R] test for difference in the order of self generated sequence

Hello, I have an experiment in which we ask subjects to generate a list of thoughts after being exposed to a stimuli. The thoughts were then coded into two categories (e.g. A B). The objective is to show that the order in which thoughts are generated is affected by the experimental conditions.

Re: [R] read.xport

Nelson, Gary (FWE) wrote: I have the latest version of foreign, but it still doesn't work. I quickly tried the Hmisc package, but the same issue arose. I will delve into the Hmisc package further. This is disappointing as we rely on read.xport quite a bit. You can use the option with

[R] visual event programming in R

Hello, One of my students is building a visual environment for the package dprep developed by my research group at the University of Puerto Rico. The environment will be in the style of the Orange (specifically channels and tokens design) and Weka (specifically Knowledge flow

Re: [R] plot the number of replicates at the same point

On 7/14/05, Marc Schwartz (via MN) [EMAIL PROTECTED] wrote: On Thu, 2005-07-14 at 12:30 -0500, Deepayan Sarkar wrote: On 7/14/05, Kerry Bush [EMAIL PROTECTED] wrote: Thank you for thinking about the problem for me. However, I have found that your method doesn't work at all. You

Re: [R] plot the number of replicates at the same point

On Thu, 2005-07-14 at 15:08 -0500, Deepayan Sarkar wrote: On 7/14/05, Marc Schwartz (via MN) [EMAIL PROTECTED] wrote: On Thu, 2005-07-14 at 12:30 -0500, Deepayan Sarkar wrote: On 7/14/05, Kerry Bush [EMAIL PROTECTED] wrote: Thank you for thinking about the problem for me. However, I

Re: [R] read.xport

Sorry about the version. I have version 0.8-8 of foreign. I downloaded (and installed) it this afternoon. The XPORT files were created in SAS using PROC COPY. I used the code: temp-read.xport(c:/temp/int82ag.xpt). I still get only 16 file of the 24 SAS placed into the XPORT file. The XPORT

[R] integer codes of factors

U = c(b, b, b, c, d, e, e) F1 = factor( U, levels=c(a, b, c, d, e) ) as.numeric(F1) [1] 2 2 2 3 4 5 5 Here, the integer code of b in F1 is 2 K = factor( levels(F1) ) as.numeric(K) [1] 1 2 3 4 5 K [1] a b c d e Levels: a b c d e And again, the integer code of b in K is 2.

Re: [R] read.xport

Frank E Harrell Jr [EMAIL PROTECTED] writes: Nelson, Gary (FWE) wrote: I have the latest version of foreign, but it still doesn't work. I quickly tried the Hmisc package, but the same issue arose. I will delve into the Hmisc package further. This is disappointing as we rely on

[R] question from environmental statistics

Dear R users I want to knw if there is a way in which a raw dataset can be modelled by some distribution. besides the gof test is there any test involving gamma or lognormal that would fit the data. thank you -dev __ R-help@stat.math.ethz.ch mailing

Re: [R] read.xport

Peter Dalgaard [EMAIL PROTECTED] writes: Frank E Harrell Jr [EMAIL PROTECTED] writes: Nelson, Gary (FWE) wrote: I have the latest version of foreign, but it still doesn't work. I quickly tried the Hmisc package, but the same issue arose. I will delve into the Hmisc package

Re: [R] plot the number of replicates at the same point

Hi all, This is nowhere near as elegant as Deepayan's solution, but I read the spec as plotting symbols except where there were overlays. You are welcome to improve the following... Jim count.overplot-function(x,y,tol=NULL,...) { if(missing(x)) stop(Usage:

[R] Variance components from lm?

I often use simple nested random-effect models for interlaboratory data. The variance components are important things to know. Is there an R function or package that gets variance components from lm objects? Or can someone point me to a method of doing so?

Re: [R] plot the number of replicates at the same point

See ?sunflowerplot for a graphic indication of the number of replications at each point. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Variance components from lm?

Under normal circumstances, 'RSiteSearch(variance components)' would likely identify lme in package nlme and lmer in package lme4. I recommend lme, as it comes with a fairly complete set of helper functions described in Pinheiro and Bates (2000) Mixed-Effect Models in S and S-Plus

Re: [R] Calculation of group summaries

Several people suggested specific functions (by, tapply, sapply and others); thanks for not blowing off a simple question regarding how to do the following SQL in R: select year, site_id, visit_no, mean(undercut) AS meanUndercut, count(undercut) AS

Re: [R] Variance components from lm?

From: Spencer Graves Under normal circumstances, 'RSiteSearch(variance components)' would likely identify lme in package nlme and lmer in package lme4. I recommend lme, as it comes with a fairly complete set of helper functions described in Pinheiro and Bates (2000)

Re: [R] question from environmental statistics

Take a look at this document by Vito Ricci: http://cran.r-project.org/doc/contrib/Ricci-distributions-en.pdf Did you try RSiteSearch(Fit distribution) or a Google search? That will lead you to fit.dist{gnlm} and fitdistr{MASS} Cheers Francisco From: [EMAIL PROTECTED] To:

Re: [R] Calculation of group summaries

1. Try using more spaces so your code is easier to read. 2. Use data.frame to define your data frame (since the method in your post creates data frames of factors rather than the desired classes). 3. Given the appropriate function, f, a single 'by' statement rbind'ed together, as shown, will

Re: [R] Calculation of group summaries

There was an error in my code (after advising you to use data.frame rather than cbind I used it myself!). Here it is again: nsites - 6 yearList - 1999:2001 fakesub - data.frame( year = rep(yearList, nsites/length(yearList), each = 11), site_id = rep(c('site1','site2'), each =

Re: [R] test for difference in the order of self generated sequence

Have you considered glmmPQL in library(MASS)? spencer graves Kaiya Liu wrote: Hello, I have an experiment in which we ask subjects to generate a list of thoughts after being exposed to a stimuli. The thoughts were then coded into two categories (e.g. A B). The

[R] problems with nls function

Hello, I was trying to fit a generalized logistic curve and my code for testing if the function nls is working is like this ### #define the generalized logistic curve function glogit=function(a, b, c, m, t, x) { tmp = 1/(1 + t * exp(-b * (x - m)))^(1/t) model.func = a + tmp * c #define

Re: [R] integer codes of factors

What's the problem? As suggested by the help page, the numeric codes are assigned in the order the names appear in the levels argument. Consider the following example from the help page plus a minor modification: (ff - factor(substring(statistics, 1:10, 1:10),

Re: [R] problems with nls function

Sorry that I specified the x as x=c(0, 0.4, 0.7, 1, 1.4, 1.6, 1.8, 2:10) The error message is the same. Thanks. Hanna - Original Message - From: Yimeng Lu To: r-help@stat.math.ethz.ch Sent: Thursday, July 14, 2005 11:19 PM Subject: problems with nls function Hello,

Re: [R] A statistical modeling problem

Let's do the simplest things first: I'd start with normal probability plots of log(t_ij). If that all looked sensible, I'd then use lme to model log(t_ij) with subject as a random effect and (probably) que as a fixed effect. Then I'd probably try something like glmmPQL in

Re: [R] problems with nls function

The start list should only contain parameters, not x. Also your function appears to have multiple errors including reference to temp (presumably intended as tmp?) and b(x-m) which presumably should be b*(x-m). On 7/14/05, Yimeng Lu [EMAIL PROTECTED] wrote: Sorry that I specified the x as

Re: [R] Variance components from lm?

Yes, I'm very thankful for the support of Prof. Jonathan Baron and all the others who have contributed their time, creative energies and (for some) money to make R what it is today. It is beyond the budgets of most people on this planet to purchase licenses for every piece of