Re: [R] Unexpected behavior of predict and interval=confidence

2006-10-30 Thread Gavin Simpson
On Sun, 2006-10-29 at 12:54 -0500, Kevin Middleton wrote: Based on some recent r-help discussions, I have been trying out plotting confidence intervals using predict and matplot. Matplot appeared to not be plotting the linear regression when using the default column names generated by

[R] how to combine imputed data-sets from mice for classfication

2006-10-30 Thread Eleni Rapsomaniki
Dear R users I want to combine multiply imputed data-sets generated from mice to do classfication. However, I have various questions regarding the use of mice library. For example suppose I want to predict the class in this data.frame: data(nhanes) mydf=nhanes mydf$class=pos

[R] Composition like classes functions

2006-10-30 Thread Gregor Gorjanc
Hello! I am working on breed composition of animals in my project. Say I have an animal 1 with parents 2 (father) and 3 (mother). If father is of breed A and mother of breed B, then their descendant is of breed AB or 50 % of A and 50 % of B. I would like to have a general way to represent this

[R] for importing data

2006-10-30 Thread amna khan
*I am a very new user of R. I've spent several hours trying to import data, so I feel okay asking the list for help. * *I had an Excel file, then I turned it into a tab-delimited file, as instructed by directions My filename is lahore.txt I amusing the following commands for read.delim but i am

Re: [R] for importing data

2006-10-30 Thread Gabor Grothendieck
Move to the appropriate directory first (or else move your file to where you are now). For example, 1. if the file is in \ then setwd(/) 2. or getwd() which shows where you are now and then you can move your file to that spot. 3. Another possibility is: read.delim(file.choose())

Re: [R] for importing data

2006-10-30 Thread David Barron
Most likely the file lahore.txt isn't in the current working directory, which is where R will be looking for it as you don't specify a path (I'm assuming that you are working in Windows). If so, you have three options. 1) Move lahore.txt to the current working directory (you can find with

Re: [R] Multivariate regression

2006-10-30 Thread Pfaff, Bernhard Dr.
Hello Ravi, have you considered the SUR method proposed by Zellner? An implementation of it is provided in CRAN-package 'systemfit' (see ?systemfit for more information). Best, Bernhard Suppose I have a multivariate response Y (n x k) obtained at a set of predictors X (n x p). I would like to

Re: [R] Multivariate regression

2006-10-30 Thread Andris Jankevics
Also you can take a look on Partial Least Squares (PLS) regression. http://www.statsoft.com/textbook/stpls.html R-package: http://mevik.net/work/software/pls.html Andris Jankevics On Sestdiena, 28. Oktobris 2006 06:04, Ritwik Sinha wrote: You can use gee (

Re: [R] for importing data

2006-10-30 Thread Roger Bivand
On Mon, 30 Oct 2006, amna khan wrote: *I am a very new user of R. I've spent several hours trying to import data, so I feel okay asking the list for help. * *I had an Excel file, then I turned it into a tab-delimited file, as instructed by directions My filename is lahore.txt I amusing the

Re: [R] How to best divide table by table

2006-10-30 Thread Martin Maechler
But *please* do learn from Peter Dalgaard's answer that this has nothing to do with tables (Petr, I agree you have helped Serguei, but...) Using vague language is not helpful in such contexts. And hence, calling data frames 'tab1' is probably not a good idea. {Calling them 'data' is not much

[R] [R-pkgs] adehabitat version 1.5

2006-10-30 Thread Clément Calenge
Dear all, I have just uploaded to CRAN the version 1.5 of the package 'adehabitat'. Significant changes are listed below: * the package now contains a new function allowing a factorial decomposition of the Mahalanobis distances in habitat selection studies, named madifa(). * The ENFA can now

Re: [R] Debugging R's code: boxplot.stats

2006-10-30 Thread Martin Maechler
Duncan == Duncan Murdoch [EMAIL PROTECTED] on Sat, 28 Oct 2006 08:03:33 -0400 writes: [...] Duncan Not sure why boxplot.stats is in grDevices rather Duncan than graphics, but that's not really relevant to Duncan your problems. because it does computation for Graphics

Re: [R] Multivariate regression

2006-10-30 Thread Robin Hankin
Hi I discovered the other day that lm() does some of the work for you: library(mvtnorm) X - matrix(rnorm(60),ncol=3) beta - matrix(1:6,ncol=2) sig - matrix(c(1,0.7,0.7,1),2,2) Y - X %*% beta + rmvnorm(n=20,sigma=sig) lm(Y ~ X-1) Call: lm(formula = Y ~ X - 1) Coefficients: [,1]

[R] reading only some columns from a table

2006-10-30 Thread Amir Safari
Dear R users, Sometimes it is needed to read only some columns from a table, in particulare for high frequency data. How it is possible to read just some certain columns using read.table ( ). The reason could be keeping space in R and in particular accelerating in reading data when

Re: [R] reading only some columns from a table

2006-10-30 Thread Gabor Grothendieck
See the colClasses argument of read.table. e.g. read.table(myfile, header = TRUE, colClasses = c(person = NULL)) assuming you don't want the column labelled person in the header. On 10/30/06, Amir Safari [EMAIL PROTECTED] wrote: Dear R users, Sometimes it is needed to read only some

Re: [R] for importing data

2006-10-30 Thread Hans-Peter
Hi, *I had an Excel file, then I turned it into a tab-delimited file, as ... lahore-read.delim(lahore.txt) I suppose you are on windows, hence you can use the package xlsReadWrite to read the Excelfile directly. It's on the CRAN and you can download it from the RGui under the menu: Packages

[R] which duplicated rows to delete

2006-10-30 Thread Søren Merser
Hi Say I've this vector with several duplicates x-c(1,2,3,4,2,6,2,8,2,3) which(duplicated(x)) [1] 5 7 9 10 11 But what I realy want is somthing like: List({2,5,7}, {3,10}, ...) Then from each sublist I can specify which of the duplicate items to drop res-NULL for(vec in myDuplicateList)

Re: [R] reading only some columns from a table

2006-10-30 Thread Amir Safari
Thank you so much for reply. I tried to do your suggested idea. It did not work. what could be my mistake ? I got the same result with and without colClasses=c() arguement. What could be the reason? even using header=TRUE the result change slightly but not into desired one. Gabor

Re: [R] which duplicated rows to delete

2006-10-30 Thread Petr Pikal
Hi you can use apply(outer(( (1:10)[1:10%in%x]), x, ==), 1, which) to get list of duplicates. But then you will need to specify which duplicates you want to discard which can be problematic. HTH Petr On 30 Oct 2006 at 11:11, Sřren Merser wrote: From: Sřren Merser [EMAIL

Re: [R] reading only some columns from a table

2006-10-30 Thread Gabor Grothendieck
Here is an example: Lines - conid person construct + 1 1 offence + 2 1 insight read.table(textConnection(Lines), header = TRUE, + colClasses = c(person = NULL)) conid person construct 1 1 1 offence 2 2 1 insight On 10/30/06, Amir Safari [EMAIL PROTECTED] wrote:

Re: [R] which duplicated rows to delete

2006-10-30 Thread Gabor Grothendieck
Try this. The first line breaks it up into lists and the second line drops any list that is not greater than 1 in length: out - tapply(seq(x), x, function(x)x) out[sapply(out, length) 1] On 10/30/06, Søren Merser [EMAIL PROTECTED] wrote: Hi Say I've this vector with several duplicates

Re: [R] reading only some columns from a table

2006-10-30 Thread Gabor Grothendieck
Sorry, here it is fixed up. The c should have been list: Lines - conid person construct 1 1 offence 2 1 insight read.table(textConnection(Lines), header = TRUE, colClasses = list(person = NULL)) On 10/30/06, Gabor Grothendieck [EMAIL PROTECTED] wrote: Here is an example: Lines - conid

[R] [Fwd: Re: Organisation of medium/large projects with multiple analyses]

2006-10-30 Thread Mark Wardle
Daniel Elliott wrote: Mark, It sounds like your data/experiment storage and organization needs are more complicated than mine, but I'll share my methodology... Many thanks for this, and for the other replies received off-list. It is much appreciated, and confirms that with something as

Re: [R] reading only some columns from a table

2006-10-30 Thread Amir Safari
Thank you. That's the case. Gabor Grothendieck [EMAIL PROTECTED] wrote: Sorry, here it is fixed up. The c should have been list: Lines - conid person construct 1 1 offence 2 1 insight read.table(textConnection(Lines), header = TRUE, colClasses = list(person = NULL)) On 10/30/06,

[R] HELP:R for windows GUI front-end error

2006-10-30 Thread yang baohua
I use R for the first time. The install file is R-2.4.0-win32.exe, and my OS is windows XP SP2. I chose the full install, but when i open it , one Error occured:R for windows GUI front-end occurs error, You must end it Does anyone know how to solve this? Waiting for your reply. Thanks. -- Baohua

Re: [R] HELP:R for windows GUI front-end error

2006-10-30 Thread Duncan Murdoch
On 10/30/2006 6:52 AM, yang baohua wrote: I use R for the first time. The install file is R-2.4.0-win32.exe, and my OS is windows XP SP2. I chose the full install, but when i open it , one Error occured:R for windows GUI front-end occurs error, You must end it Does anyone know how to solve

[R] match lists

2006-10-30 Thread Heymans, MW
Dear list, I have this problem, please your advice. I have list A that contains two matrix elements: [[1]] a b [1,] 2 1 [2,] 3 2 [3,] 3 2 [[2]] c d [1,] 3 5 [2,] 3 1 [3,] 2 3 and list B, that also contains 2 matrices: [[1]] e f g [1,] 1 20 30 [2,] 2 40 50 [3,] 3 60

[R] nlme Error: Subscript out of bounds

2006-10-30 Thread Lisa Avery
Hello, I am new to non-linear growth modelling in R and I am trying to reproduce an analysis that was done (successfully) in S-Plus. I have a simple non-linear growth model, with no nesting. I have attempted to simplify the call as much as possible (by creating another grouped object,

[R] spreading activation

2006-10-30 Thread Kaustubh Patil
Hi, is there spreading activation function available for R? any suggestions appreiated... regards, Kaustubh - [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

[R] gui for R

2006-10-30 Thread Weiwei Shi
Hi, I met with many bugs when I used JGR and am looking for some other GUI's for Mac. I found SciViews-R by googling but it is only for windows to my best knowledge. Is there any other good gui like that, usable for Mac? Thanks, -- Weiwei Shi, Ph.D Research Scientist GeneGO, Inc. Did you

Re: [R] gui for R

2006-10-30 Thread Duncan Murdoch
On 10/30/2006 11:03 AM, Weiwei Shi wrote: Hi, I met with many bugs when I used JGR and am looking for some other GUI's for Mac. I found SciViews-R by googling but it is only for windows to my best knowledge. Is there any other good gui like that, usable for Mac? What do you find wrong with

[R] Problem setting TMPDIR on the fly

2006-10-30 Thread Benjamin Lloyd-Hughes
Hi folks, I'm having a spot of bother with tempdir(). If I set the environment variable TMPDIR in my shell prior to invoking R it works as expected: Sys.getenv(TMPDIR) TMPDIR /eurotempest/tmp/R tempdir() [1] /eurotempest/tmp/R/Rtmp0xY4XD However if I don't set it

[R] installing package help with limited permissions

2006-10-30 Thread Eva Goldwater
Hello, I administer a student (Windows based) computer lab, which has R installed. I want to permit students to install R packages, but of course they don't have write permission to R root folder. Following FAQ suggestion, I've set up a folder that they can write to, and set the R_LIB

Re: [R] Problem setting TMPDIR on the fly

2006-10-30 Thread Peter Dalgaard
Benjamin Lloyd-Hughes [EMAIL PROTECTED] writes: Hi folks, I'm having a spot of bother with tempdir(). If I set the environment variable TMPDIR in my shell prior to invoking R it works as expected: Sys.getenv(TMPDIR) TMPDIR /eurotempest/tmp/R tempdir()

[R] memory management

2006-10-30 Thread Federico Calboli
Hi All, just a quick (?) question while I wait my code runs... I'm comparing the identity of the lines of a dataframe, doing all possible pairwise comparisons. In doing so I use identical(), but that's by the way. I'm doing a (not so) quick and dirty check, and subsetting the data as

[R] (no subject)

2006-10-30 Thread Tancredi Caruso
Please, help me. I had a trouble in downloading from CRAN the package untb. Everything seems to work. I type install.packages(c(untb)) after selecting the CRAN mirror and I get a message of successful downloading. Indeed I have the package in the library and if I check with command

[R] no possible to load a package correctly dowloaded by the R prompt

2006-10-30 Thread Tancredi Caruso
Please, help me. I had a trouble in downloading from CRAN the package untb. Everything seems to work. I type install.packages(c(untb)) after selecting the CRAN mirror and I get a message of successful downloading. Indeed I have the package in the library and if I check with command

Re: [R] memory management

2006-10-30 Thread bogdan romocea
This was asked before. Collapse the data frame into a vector, e.g. v - apply(DF,1,function(x) {paste(x,collapse=_)}) then work with the values of that vector (table, unique etc). If your data frame is really large run this in a DBMS. -Original Message- From: [EMAIL PROTECTED]

Re: [R] POSIXct time zone and daylight savings issues

2006-10-30 Thread Sebastian P. Luque
Hi again, A related issue I can't quite understand is: R tt - as.POSIXct(2006-05-24, tz=EEST) R tt [1] 2006-05-24 EEST R seq(tt, length=12, by=months) [1] 2006-05-24 EEST 2006-06-24 EEST 2006-07-24 EEST 2006-08-24 EEST [5] 2006-09-24 EEST 2006-10-24 EEST 2006-11-24 EEST 2006-12-24 EEST [9]

Re: [R] no possible to load a package correctly dowloaded by the R prompt

2006-10-30 Thread Michael Kubovy
Hi, On Oct 30, 2006, at 11:43 AM, Tancredi Caruso wrote: install.packages(c(untb)) library(untb) Carico il pacchetto richiesto: partitions (I'm loading requested package:partitions) Errore: pacchetto 'partitions' non caricato (Error: package partitions not loadad) Warning messages:

[R] correlation structure in lme without random effect

2006-10-30 Thread Benjamin Tyner
I was hoping to fit along the lines of g-gl(20,5) y-runif(100) fit-lme(fixed=y~g,correlation=corAR1(0,~1|g)) But I get the error Incompatible formulas for groups in random and correlation Any help would be greatly appreciated. Ben __

[R] Random intercept-slope correlation (nlme)

2006-10-30 Thread Antonio Revilla
Dear list members, I am working with a multilevel growth curve, that in its simplest form goes like follows: Yit = Ai + Bi t + eit (the error term is assumed to follow an AR(1) autorregressive process) One major topic in my research is the convergence in the values of Y over time. Thus, I

Re: [R] correlation structure in lme without random effect

2006-10-30 Thread David Barron
You haven't specified a random equation. Try: fit-lme(fixed=y~g,random=~1|g,correlation=corAR1(0,~1|g)) On 30/10/06, Benjamin Tyner [EMAIL PROTECTED] wrote: I was hoping to fit along the lines of g-gl(20,5) y-runif(100) fit-lme(fixed=y~g,correlation=corAR1(0,~1|g)) But I get the error

Re: [R] correlation structure in lme without random effect

2006-10-30 Thread Douglas Bates
On 10/30/06, Benjamin Tyner [EMAIL PROTECTED] wrote: I was hoping to fit along the lines of g-gl(20,5) y-runif(100) fit-lme(fixed=y~g,correlation=corAR1(0,~1|g)) But I get the error Incompatible formulas for groups in random and correlation Use the gls function in the nlme package to fit

Re: [R] correlation structure in lme without random effect

2006-10-30 Thread Benjamin Tyner
Thanks! Douglas Bates wrote: On 10/30/06, Benjamin Tyner [EMAIL PROTECTED] wrote: I was hoping to fit along the lines of g-gl(20,5) y-runif(100) fit-lme(fixed=y~g,correlation=corAR1(0,~1|g)) But I get the error Incompatible formulas for groups in random and correlation Use the gls

Re: [R] POSIXct time zone and daylight savings issues

2006-10-30 Thread Gabor Grothendieck
I don't know the answer to your question but if you are using dates with no times and don't need time zones (both of which appear to be the case here) then you could use Date class and avoid the issue altogether. See the help desk article in R News 4/1 where there is a discussion of how to choose

Re: [R] Problem setting TMPDIR on the fly

2006-10-30 Thread Duncan Murdoch
On 10/30/2006 11:22 AM, Benjamin Lloyd-Hughes wrote: Hi folks, I'm having a spot of bother with tempdir(). If I set the environment variable TMPDIR in my shell prior to invoking R it works as expected: Sys.getenv(TMPDIR) TMPDIR /eurotempest/tmp/R

Re: [R] Is there any method to do spatial sampling?

2006-10-30 Thread Marshall Feldman
What kind of spatial sampling do you want to do: systematic, stratified, simple random, etc.? What will you be sampling: points, lines, polygons, etc.? Marsh Feldman URI Center for Urban Studies and Research The University of Rhode Island -Original Message- From:

[R] how to set debug breaks

2006-10-30 Thread waverley palo
Hi, I am new in R and have some frustrations as how to set debug breaks during emacs R debug. I use debug () as where or which function to debug. But during the debug, e.g., I have a for loop at the beginning of the function code and want the code execution to jump through that for loop and set

[R] read.fwf and header

2006-10-30 Thread Gregor Gorjanc
Hi! I have data (also in attached file) in the following form: num1 num2 num3 int1 fac1 fac2 cha1 cha2 Date POSIXt 11 f q 1900-01-01 1900-01-01 01:01:01 2 1.0 131.5 2 a g r z1900-01-01 01:01:01 3 1.5 1188830.5 3 b h s y 1900-01-01 1900-01-01 01:01:01 4

[R] Additional earnings

2006-10-30 Thread Clark Pray
Reply to: [EMAIL PROTECTED] Denmark software company, Gricina Software Company is looking for employees in Australia and Iceland! Demands: 21-70 years old, male(female), smart, communicative. You will work for yourself, wage 1000 USD- 3000 USD per week, work at morning time 9.00 AM- 11.00 AM.

[R] Which executable is associated with R CMD INSTALL?

2006-10-30 Thread Patrick Connolly
I'm still having trouble installing the lme4 package on RHEL 3. I've asked this list and it seems my problem is not universal. Brian Ripley indicated that the problem was with recognising the Matrix package, even though I've taken care to get the most recent versions of Matrix and lme4. It

Re: [R] read.fwf and header

2006-10-30 Thread Marc Schwartz
On Mon, 2006-10-30 at 19:51 +0100, Gregor Gorjanc wrote: Hi! I have data (also in attached file) in the following form: num1 num2 num3 int1 fac1 fac2 cha1 cha2 Date POSIXt 11 f q 1900-01-01 1900-01-01 01:01:01 2 1.0 131.5 2 a g r z1900-01-01 01:01:01

[R] psigamma derivative

2006-10-30 Thread anamaria
Hello, I am trying to find a hessian matrix that involves log(gamma(1/p)) second derivative, p being one of the parameters of the function. I am using a function deriv with the hessian=TRUE option, but psigamma is not on the list of derivative functions. I know that it is possible to use

Re: [R] read.fwf and header

2006-10-30 Thread Daniel Nordlund
Gregor, According to the help for read.fwf, sep needs to be set to a value that occurs only in the header record. I changed the spaces to commas in the header record of your example and used the following syntax and was able to read the file just fine. new.data-read.fwf(file=test.txt,

Re: [R] psigamma derivative

2006-10-30 Thread Christos Hatzis
Try ?Special For information on the special functions included in the base R distribution. Also, RSiteSearch(digamma) gives several hits on additional packages that might be useful. -Christos -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL

[R] Hmisc or Lattice plot with error bars, Depth (independent variable) on Y axis

2006-10-30 Thread Mike Bock
I am trying to create some plots of concentration versus depth. I want Depth on the Y axis and the concentrations on the X axis. I also need to plot error bars. The xYplot function in Hmisc is very nearly ideal but I have to put depth on the x-axis to get it to work. When I transpose the X

Re: [R] how to set debug breaks

2006-10-30 Thread Gavin Simpson
On Mon, 2006-10-30 at 10:37 -0800, waverley palo wrote: Hi, I am new in R and have some frustrations as how to set debug breaks during emacs R debug. I use debug () as where or which function to debug. But during the debug, e.g., I have a for loop at the beginning of the function code and

Re: [R] Hmisc or Lattice plot with error bars, Depth (independent variable) on Y axis

2006-10-30 Thread Frank E Harrell Jr
Mike Bock wrote: I am trying to create some plots of concentration versus depth. I want Depth on the Y axis and the concentrations on the X axis. I also need to plot error bars. The xYplot function in Hmisc is very nearly ideal but I have to put depth on the x-axis to get it to work.

Re: [R] read.fwf and header

2006-10-30 Thread Gregor Gorjanc
Marc Schwartz wrote: On Mon, 2006-10-30 at 19:51 +0100, Gregor Gorjanc wrote: Hi! I have data (also in attached file) in the following form: num1 num2 num3 int1 fac1 fac2 cha1 cha2 Date POSIXt 11 f q 1900-01-01 1900-01-01 01:01:01 2 1.0 131.5 2 a g r z

Re: [R] how to set debug breaks

2006-10-30 Thread Duncan Murdoch
On 10/30/2006 1:37 PM, waverley palo wrote: Hi, I am new in R and have some frustrations as how to set debug breaks during emacs R debug. I use debug () as where or which function to debug. But during the debug, e.g., I have a for loop at the beginning of the function code and want the

Re: [R] Which executable is associated with R CMD INSTALL?

2006-10-30 Thread Duncan Murdoch
On 10/30/2006 2:32 PM, Patrick Connolly wrote: I'm still having trouble installing the lme4 package on RHEL 3. I've asked this list and it seems my problem is not universal. Brian Ripley indicated that the problem was with recognising the Matrix package, even though I've taken care to get

Re: [R] match lists

2006-10-30 Thread bogdan romocea
What is it that you don't know how to do? Loop over the matrices from the 2 lists and merge them two by two, for example AB - list() ; id - 1 for (i in 1:length(A)) for (j in 1:length(B)) { AB[[id]] - merge(A[[i]],B[[j]],...) id - id + 1 } To better keep track of who's who, you may want to

[R] VGAM-vglm

2006-10-30 Thread minta
Can somebody explain to me why I get different results using multinom(package nnet) and vglm(package VGAM)? I think I understand the result I get with multinom, but what has vglm done? How can I interpret the results? (paarer: nominal variable, classes for different jobs, K26A:number of

Re: [R] installing package help with limited permissions

2006-10-30 Thread Duncan Murdoch
On 10/30/2006 11:25 AM, Eva Goldwater wrote: Hello, I administer a student (Windows based) computer lab, which has R installed. I want to permit students to install R packages, but of course they don't have write permission to R root folder. Following FAQ suggestion, I've set up a

[R] plot history

2006-10-30 Thread Rohini Mulford
___ Hi, When I create multiple graphs subsequent graphs overwrite previous graphs. How do I keep all my graphs in the current workspace (but not all on the same page) so I can scroll through them?

Re: [R] Debugging R's code: boxplot.stats

2006-10-30 Thread Matthew Walker
Hi Martin, Sorry, I did intend to put some examples in but must have forgotten. This example is easy to understand: boxplot(c(1,Inf,Inf,Inf,Inf)) has lower bound: 1, lower quartile: Inf, median: Inf, upper quartile: Inf, upper bound: Inf. The command currently errors with the obscure message

Re: [R] plot history

2006-10-30 Thread David Barron
Are you working in Windows? If so, you can turn on recording either in the History menu on the active graphics device window, or by using windows(record=TRUE) to open a new device. On 30/10/06, Rohini Mulford [EMAIL PROTECTED] wrote:

Re: [R] read.fwf and header

2006-10-30 Thread Gregor Gorjanc
Daniel Nordlund wrote: Gregor, According to the help for read.fwf, sep needs to be set to a value that occurs only in the header record. I changed the spaces to commas in the header record of your example and used the following syntax and was able to read the file just fine.

Re: [R] plot history

2006-10-30 Thread Leeds, Mark \(IED\)
I send them all to one postscript file and then bring it up in ghostview but maybe there is another way. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Rohini Mulford Sent: Monday, October 30, 2006 5:29 PM To: r-help@stat.math.ethz.ch Subject: [R] plot

[R] help_aov

2006-10-30 Thread Thanjavur Bragadeesh
Hi, I am trying to run an analysis of variance using R. in my data table x is a continuous variable lengthof 200 and p is a categorical variable also of length 200 and p is anyone of three categories 1,2 or ,3. if I run summary(aov(x~p,data=test)) I get Response: x Df Sum Sq Mean

Re: [R] plot history

2006-10-30 Thread Mihai Nica
Or maybe see ?win.graph. Before each graph use something like: win.graph(width=3,height=3, pointsize=8) You will get a sandwich of graphs. hth, Mihai Nica 170 East Griffith St. G5 Jackson, MS 39201 601-914-0361 - Original Message From: Rohini Mulford [EMAIL PROTECTED] To:

Re: [R] help_aov

2006-10-30 Thread Peter Dalgaard
Thanjavur Bragadeesh [EMAIL PROTECTED] writes: Hi, I am trying to run an analysis of variance using R. in my data table x is a continuous variable lengthof 200 and p is a categorical variable also of length 200 and p is anyone of three categories 1,2 or ,3. if I run

Re: [R] Hmisc or Lattice plot with error bars, Depth (independent variable) on Y axis

2006-10-30 Thread Deepayan Sarkar
On 10/30/06, Mike Bock [EMAIL PROTECTED] wrote: I am trying to create some plots of concentration versus depth. I want Depth on the Y axis and the concentrations on the X axis. I also need to plot error bars. The xYplot function in Hmisc is very nearly ideal but I have to put depth on

[R] as.zoo question

2006-10-30 Thread Leeds, Mark \(IED\)
this is probably a question for gabor or achim but maybe someone else can answer it in case they are not around ? Ifr you paste the code below ( again you would need the zoo and chron libraries to be installed ), you should get the error Error in Ops.POSIXt(frequency, freq) : %% not defined

Re: [R] as.zoo question

2006-10-30 Thread Leeds, Mark \(IED\)
oops, for the pasting to work, You would have to get rid of the blank line that occurs during the bayesfactor construction. Thanks. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Leeds, Mark (IED) Sent: Monday, October 30, 2006 6:15 PM To:

[R] how to plot a data.frame?

2006-10-30 Thread Weiwei Shi
hi, i have a data frame like this: 0.3 0.7 0.4 0.8 i am trying to plot this data frame and each cell is filled with different colors based on the value. Is there a function which can do this? thanks, -- Weiwei Shi, Ph.D Research Scientist GeneGO, Inc. Did you always know? No, I did not.

Re: [R] how to plot a data.frame?

2006-10-30 Thread Jeffrey Robert Spies
Yep: ?image Creates a grid of colored or gray-scale rectangles with colors corresponding to the values in z. This can be used to display three- dimensional or spatial data aka “images”. This is a generic function. Hope that helps, Jeff. On Oct 30, 2006, at 7:01 PM, Weiwei Shi wrote: hi,

Re: [R] how to plot a data.frame?

2006-10-30 Thread Waverley
try heatmap function On 10/30/06, Weiwei Shi [EMAIL PROTECTED] wrote: hi, i have a data frame like this: 0.3 0.7 0.4 0.8 i am trying to plot this data frame and each cell is filled with different colors based on the value. Is there a function which can do this? thanks, -- Weiwei

Re: [R] as.zoo question

2006-10-30 Thread Gabor Grothendieck
It should be zoo, not as.zoo. We are constructing a zoo object from scratch here, not converting a time series object of a different class to zoo. On 10/30/06, Leeds, Mark (IED) [EMAIL PROTECTED] wrote: oops, for the pasting to work, You would have to get rid of the blank line that occurs

[R] exteremely confused (simple question)

2006-10-30 Thread Chris Fonnesbeck
I thought I knew how to use data frames, but apparently not. I have created a data frame, and named the columns: df KSurv Growth Class 1 4808 0.86212 0.00669640 S 2 2430 0.98038 1.3054 S 3 2084 0.93579 0.44079000 S 4 2600 0.95394 2.0368 S

Re: [R] exteremely confused (simple question)

2006-10-30 Thread Charles Annis, P.E.
You need to tell R that you mean Growth, the column in your data.frame. Otherwise R is looking for an external criterion named Growth. Try this: df[df$Growth0.5,] That means choose all columns for which df$Growth has a row value 0.5 df[df$Growth0.5,] IDKSurvGrowth Class 1 1

[R] Problem mit den Spracheinstellungen nach Installation von R

2006-10-30 Thread Dirk Hüser
Hallo, ich habe ein Problem mit R. Ich hoffe hier auf die freundliche Hilfe der Community. Schilderung des Problems: Ich habe R 2.4.0 auf einem Rechner (1) installiert. - Am Anfang der Installation klickte ich auf Setup Deutsch. - Dann führte ich die Installation weiter mit den

[R] Confidence interval calculation in prop.test

2006-10-30 Thread Richard Johnston
The confidence interval calculation in prop.test appears to be incorrect when alternative=greater. The upper limit is always set to 1.000. Am I missing something? total=c(250,250) success=c(55,31) prop.test(success,total,alternative=greater,correct=TRUE) 2-sample test for

[R] plotting multiple groups (newbie Q)

2006-10-30 Thread Sumitrajit Dhar
Hi Folks, After loading a data set, I run the following: kSum - orderBy(~group,(summaryBy(DP_Level~F2 +group,data=kdata,FUN=c(mean,sd),na.rm=T))) kSum looks like this: kSum F2 group DP_Level.mean DP_Level.sd 1 1.0 N -1.55186475 11.022245 4 2.0 N -2.48013300

[R] Saving a function

2006-10-30 Thread patopatasfrias
Hey, I'm a real novice... how do I go about saving a function that I created eg. If I created this function: Fun-function(x1,X2) {..}, how can I store it so that the next time I use R I can load it and not have to type it out again Cheers Matt -- View this message in context:

Re: [R] [Fwd: Re: Organisation of medium/large projects with multiple analyses]

2006-10-30 Thread David Farrar
It's good to see this sort of thing discussed. For my current approach, I keep a fairly static directory for function libraries, another one for large data sets, and others for projects. I try to define tasks (probably like your analyses) within projects. There is a project

[R] [R-pkgs] rcompletion: TAB completion for the R command line

2006-10-30 Thread Deepayan Sarkar
CRAN now has a package called rcompletion that attempts to provide TAB completion for R using the GNU readline library, intended for R sessions run from a command line. From the package help page: Description: This package provides pseudo-intelligent TAB completion for a readline enabled

[R] power.z.test

2006-10-30 Thread Ethan Johnsons
Does anyone have/know how to write the power function for z tests something like below? function(a,m0,m1,n,s){ t1 = -qnorm(1-a) num = abs(m0-m1) * sqrt(n) t2 = num/s pow = pnorm(t1 + t2) } thx much ej __ R-help@stat.math.ethz.ch mailing list

[R] problem in intallasation

2006-10-30 Thread amna khan
Sir I have downloaded new release of R 2.4.0. But there is no R.2.4.0.exefile for installing it. I request you to please guid me. Regards -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. Email: [EMAIL PROTECTED] [EMAIL PROTECTED] [EMAIL PROTECTED]

[R] how to make a data file

2006-10-30 Thread amna khan
Sir after importing data from excel to R, I am not understanding how to make this data file. So that I can use it in extRemeToolkit and other packages. Thank you -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. Email: [EMAIL PROTECTED] [EMAIL PROTECTED] [EMAIL

Re: [R] plotting multiple groups (newbie Q)

2006-10-30 Thread Petr Pikal
Hi I suspect that value in kSum$group is not exactly Y (some space like Y ) or maybe there is some other variable named group elsewhere but without value Y. try to look at kSum$group or do ls(). HTH Petr On 30 Oct 2006 at 20:13, Sumitrajit Dhar wrote: To:

Re: [R] Saving a function

2006-10-30 Thread Petr Pikal
Hi 1. You can create your own package (set of functions) and load it each time you start R 2. When quitting R you can save workspace. It will create a file named .Rdata in which your function is defined. Then you can load this workspace and use your function 3. You can use some of dput,

Re: [R] read.fwf and header

2006-10-30 Thread Martin Maechler
Gregor == Gregor Gorjanc [EMAIL PROTECTED] on Mon, 30 Oct 2006 23:33:21 +0100 writes: Gregor Daniel Nordlund wrote: Gregor, According to the help for read.fwf, sep needs to be set to a value that occurs only in the header record. I changed the spaces to commas