Re: [R] ANOVA non-sphericity test and corrections (eg, Greenhouse-Geisser)

2007-06-25 Thread Peter Dalgaard
DarrenWeber wrote: I'm an experimental psychologist and when I run ANOVA analysis in SPSS, I normally ask for a test of non-sphericity (Box's M-test). I also ask for output of the corrections for non-sphericity, such as Greenhouse-Geisser and Huhn-Feldt. These tests and correction factors

[R] Random numbers from skewed distributions

2007-06-25 Thread Anup Nandialath
Dear Friends, I was wondering if there is any package to get random numbers from the Burr 10 distribution. I checked the rmutil and actuar package. Both seems to implement the Burr 12 distribution. thanks in advance Regards Anup -

Re: [R] Source code for rlogis

2007-06-25 Thread Peter Dalgaard
Anup Nandialath wrote: Dear friends, I was trying to read the source code for rlogis but ran into a roadblock. It shows [[1]] function (n, location = 0, scale = 1) .Internal(rlogis(n, location, scale)) environment: namespace:stats Is is possible to access the source code for the same.

Re: [R] ANOVA non-sphericity test and corrections (eg, Greenhouse-Geisser)

2007-06-25 Thread Simon Blomberg
On Mon, 2007-06-25 at 17:53 +1000, Simon Blomberg wrote: If you use lme, you can fit a general correlation structure to the within-subject data, and compare the fit to a model assuming uncorrelated within-subjects errors. That should tell you whether your data are ... correlated... (damn

Re: [R] How to run mathematica or c programs in R?

2007-06-25 Thread Schmitt, Corinna
Hallo, I just know a solution if you use MATLAB. Here you need the library R.matlab and there the functions writeMat and readMat. You can download the package on CRAN. Corinna ** Corinna Schmitt, Dipl.Inf.(Bioinformatik)

Re: [R] ANOVA non-sphericity test and corrections (eg, Greenhouse-Geisser)

2007-06-25 Thread Simon Blomberg
If you use lme, you can fit a general correlation structure to the within-subject data, and compare the fit to a model assuming uncorrelated within-subjects errors. That should tell you whether your data are Aren't the G-G and H-F corrections only approximate fixes? Surely it is better to work

Re: [R] adding lines to stripchart

2007-06-25 Thread Chuck Cleland
James Root wrote: I have two points of collection across 20 subjects (pre and post for each), so 20 pairs of data points. I would like to plot the actual raw data points for each subject for both pre and post and connect lines between these two points (20 in all) to depict real change between

Re: [R] R-2.5.0 compilation problem on Linux powerpc

2007-06-25 Thread Prof Brian Ripley
Is this ppc32 or ppc64? (What does uname -a say?) If the former, you might need to set CPICFLAGS and FPICFLAGS to -fPIC (rather than -fpic): please look these up in the R-admin.html file (which INSTALL points you to). For ppc64 configure should have found -fPIC. On Mon, 25 Jun 2007, Vivek

[R] Re : Half of a heatmap

2007-06-25 Thread Neil Shephard
I am trying to produce a heatmap of pairwise correlations, but since the matrix is symmetric, I only need either the upper or the lower triangle. I have scoured the web and R documentation, but I have not been able to find a way to produce such a figure. Is there a simple way to

[R] fractional calculations

2007-06-25 Thread Federico Calboli
Hi All, is there a function in R that allows me to work with fractions without transforming them to floats (or whatever) in between? Something that would calculate something like: (1/2 + 1/8) * 1/2 = 5/16 without ever transforming to 0.5 and 0.125? Best, Federico -- Federico C. F. Calboli

Re: [R] fractional calculations

2007-06-25 Thread Henrique Dallazuanna
require(MASS) ?as.fractions as.fractions(1/2+1/8) -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 Ohttp://maps.google.com/maps?f=qhl=enq=Curitiba,+Brazillayer=ie=UTF8z=18ll=-25.448315,-49.276916spn=0.002054,0.005407t=kom=1 On 25/06/07, Federico Calboli [EMAIL PROTECTED]

Re: [R] fractional calculations

2007-06-25 Thread Stefan Grosse
is there a function in R that allows me to work with fractions without transforming them to floats (or whatever) in between? You could use the ryacas (CAS) package: http://code.google.com/p/ryacas/ and/or http://cran.r-project.org/doc/vignettes/Ryacas/Ryacas.pdf e.g: library(Ryacas)

[R] gam function in the mgcv library

2007-06-25 Thread Bill Wheeler
I would like to fit a logistic regression using a smothing spline, where the spline is a piecewise cubic polynomial. Is the knots option used to define the subintervals for each piece of the cubic spline? If yes and there are k knots, then why does the coefficients field in the returned object

Re: [R] Re : Half of a heatmap

2007-06-25 Thread Rosenfeld, Jeffrey
Thank you for your reply. I have looked at LDheatmap, but it does not seem to do what I want and seems to only work well for LD data. I was looking for something that would produce a figure identical to what heatmap.2 gives me, including the proper X and Y-axis labels and a dendogram, except

[R] Peter Moser ist außer Haus.

2007-06-25 Thread peter . moser
Ich werde ab 25.06.2007 nicht im Büro sein. Ich kehre zurück am 28.06.2007. in dringenden Fällen bin ich unter 079 79 73 74 6 erreichbar __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] R-excel

2007-06-25 Thread Christophe Pallier
cf. R Data Import/Export file:///tmp/RtmpnFqONj/.R/doc/manual/R-data.htmlin the standard documentation. Christophe On 6/25/07, Erika Frigo [EMAIL PROTECTED] wrote: Good morning to everybody, I have a problem : how can I import excel files in R??? thank you very much Dr.sa. Erika Frigo

Re: [R] fractional calculations

2007-06-25 Thread Federico Calboli
Henrique Dallazuanna wrote: require(MASS) ?as.fractions as.fractions(1/2+1/8) I thought that as.fractions did transform the fractions *first* into floats and *then* found the rational approssimation (a passage I'd rather avoid): fractionspackage:MASSR

[R] Help in HMM for count data

2007-06-25 Thread Regina Verghis
Respected Sir, I am Regina M.Verghis, Final Year PG in Biostatistics. My thesis is based on a count data. I am trying to fit the data using HMM. I have the count of new infectives each month. I am not considering any other covariates. I am using the 'chidden(repeated package)'. Can you

[R] R-excel

2007-06-25 Thread Erika Frigo
Good morning to everybody, I have a problem : how can I import excel files in R??? thank you very much Dr.sa. Erika Frigo Università degli Studi di Milano Facoltà di Medicina Veterinaria Dipartimento di Scienze e Tecnologie Veterinarie per la Sicurezza Alimentare (VSA) Via Grasselli, 7

Re: [R] JRI and Axis Web Service

2007-06-25 Thread Rajarshi Guha
On Jun 24, 2007, at 7:58 PM, Yifan (Eric) Jiang wrote: I've been asked to develop a Java Axis web service to retrieve an R-script file from the client side and then using JRI jar file in order to call the R program that installed on the web service server. After reading the online

Re: [R] R-2.5.0 compilation problem on Linux powerpc

2007-06-25 Thread Prof Brian Ripley
On Mon, 25 Jun 2007, Vivek Menon wrote: uname -a gives me this:== Linux 2.6.21.1-xserve #17 SMP Thu Jun 14 19:45:57 MDT 2007 ppc64 ppc64 ppc64 GNU/Linux === Also when I configure I see the foll. output: R is now configured for

[R] transposing data.frames

2007-06-25 Thread Christoph Heibl
Hello, This must be simple... Thanks a lot - Christoph # Imagine you have a list, e.g: K - list(1:10, 2:11, 9:18) K # Transforming to dataframe... KK - as.data.frame(K) # ... one obtaines the list elements as column. KK # But I need the list elements as rows # How can I achieve this? Is

Re: [R] gam function in the mgcv library

2007-06-25 Thread Simon Wood
On Monday 25 June 2007 13:26, Bill Wheeler wrote: I would like to fit a logistic regression using a smothing spline, where the spline is a piecewise cubic polynomial. Is the knots option used to define the subintervals for each piece of the cubic spline? - if you use something like

Re: [R] R-excel

2007-06-25 Thread Stefan Grosse
http://cran.r-project.org/doc/manuals/R-data.html#Reading-Excel-spreadsheets plus there is a package xlsReadWrite that might be of your interest. Stefan Original Message Subject: [R] R-excel From: Erika Frigo [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 25.06.2007

Re: [R] R-2.5.0 compilation problem on Linux powerpc

2007-06-25 Thread Vivek Menon
uname -a gives me this:== Linux 2.6.21.1-xserve #17 SMP Thu Jun 14 19:45:57 MDT 2007 ppc64 ppc64 ppc64 GNU/Linux === Also when I configure I see the foll. output: R is now configured for powerpc64-unknown-linux-gnu Source directory: .

Re: [R] transposing data.frames

2007-06-25 Thread Stefan Grosse
Original Message Subject: [R] transposing data.frames From: Christoph Heibl [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Date: 25.06.2007 15:13 Hello, This must be simple... Thanks a lot - Christoph # Imagine you have a list, e.g: K - list(1:10, 2:11, 9:18) K #

Re: [R] transposing data.frames

2007-06-25 Thread David Barron
t(KK) will transpose your data frame On 25/06/07, Christoph Heibl [EMAIL PROTECTED] wrote: Hello, This must be simple... Thanks a lot - Christoph # Imagine you have a list, e.g: K - list(1:10, 2:11, 9:18) K # Transforming to dataframe... KK - as.data.frame(K) # ... one obtaines the

Re: [R] R-2.5.0 compilation problem on Linux powerpc

2007-06-25 Thread Vivek Menon
uname -a gives me this:== Linux 2.6.21.1-xserve #17 SMP Thu Jun 14 19:45:57 MDT 2007 ppc64 ppc64 ppc64 GNU/Linux === Also when I configure I see the foll. output: R is now configured for powerpc64-unknown-linux-gnu That might be the

Re: [R] R-2.5.0 compilation problem on Linux powerpc

2007-06-25 Thread Vivek Menon
After setting the flags to -fPIC using ./configure --with-x=no --with-lapack=/apps/lib/LAPACK/lapack_LINUX.a CPICFLAGS=-fPIC FPICFLAGS=-fPIC I still get the following errors:cc -std=gnu99 -shared -L/usr/local/lib -o grDevices.so chull.o devNull.o devPicTeX.o devPS.o devQuartz.o init.o make[5]:

Re: [R] fractional calculations

2007-06-25 Thread Prof Brian Ripley
On Mon, 25 Jun 2007, Henrique Dallazuanna wrote: require(MASS) ?as.fractions as.fractions(1/2+1/8) I think Federico wanted 1/as.fractions(2) + 1/as.fractions(8) that is avoiding computing 1/8 in float (although it is exact). You might be better off with package gmp: 1/as.bigq(2) +

[R] R behaviour related to user input (readline()) and run selection

2007-06-25 Thread Bernzweig, Bruce \(Consultant\)
When I run the below section of code I get the following error: Error in if (co == A || co[1] == O) { : missing value where TRUE/FALSE needed When I run the code in two parts where I first get the user's input then afterwards run the if/else section, there is no

Re: [R] R-2.5.0 compilation problem on Linux powerpc

2007-06-25 Thread Prof Brian Ripley
On Mon, 25 Jun 2007, Vivek Menon wrote: uname -a gives me this:== Linux 2.6.21.1-xserve #17 SMP Thu Jun 14 19:45:57 MDT 2007 ppc64 ppc64 ppc64 GNU/Linux === Also when I configure I see the foll. output: R is now configured for

Re: [R] R-2.5.0 compilation problem on Linux powerpc

2007-06-25 Thread Vivek Menon
I did not get what you replied in the previous email. I still get the same error with CPICFLAGS and FPICFLAGS set to -fPIC. Vivek On 6/25/07, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Mon, 25 Jun 2007, Vivek Menon wrote: uname -a gives me this:== Linux

[R] how to plot this?

2007-06-25 Thread Weiwei Shi
Hi, there: Suppose I have a couple of data.frames and each one has five columns (one for x-axis, two for y-axis and two for std of y's.) There is another dimensions (besides x and y) which is continuous. My question is, how to plot such series of data frames in one plot (thus, 3-dimensional plot)

[R] eps in odfWeave

2007-06-25 Thread Dieter Menne
Dear Weavers, Does someone have an example of using eps or any other vector graphics with odfWeave? It tried the example below (and commented variants) with simple.odt in the examples directory, and got an error. Dieter #--- library(odfWeave) plotInfo - getImageDefs() plotInfo$type = eps

Re: [R] R-2.5.0 compilation problem on Linux powerpc

2007-06-25 Thread Prof Brian Ripley
On Mon, 25 Jun 2007, Vivek Menon wrote: I did not get what you replied in the previous email. I still get the same error with CPICFLAGS and FPICFLAGS set to -fPIC. And a completely fresh build? Vivek On 6/25/07, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Mon, 25 Jun 2007, Vivek Menon

Re: [R] R-excel

2007-06-25 Thread Sébastien
Hi, Directly import from Excel files should be possible using the RODBC package. Yu may want to read the R data import/export manual about the RODBC package for further details. Another solution is to save each of your Excel sheets as *.csv file and import the data in R using the read.table

Re: [R] Imputing missing values in time series

2007-06-25 Thread Achim Zeileis
On Fri, 22 Jun 2007, Horace Tso wrote: Thanks to Mark and Erik for different versions of locf, also Erik's pointer to archive where I found another function due to Simon Fear. I haven't tested the zoo locf function. Just as an addition to what Marc already wrote: zoo offers at least two

Re: [R] eps in odfWeave

2007-06-25 Thread Kuhn, Max
Dieter, What is the version of odfWeave? If you use those specifications in the current version (0.5.9): plotInfo - getImageDefs() plotInfo$type = eps plotInfo$device = postscript setImageDefs(plotInfo) you will probabiliy need to set horizontal = FALSE, onefile = FALSE, and paper = special

[R] a string to enviroment or function

2007-06-25 Thread Weiwei Shi
Hi, I am wondering how to make a function Fun to make the following work: t0 - (paste(hgu133a, ENTREZID, sep=)) xx - as.list(Fun(t0)) # make it work like xx-as.list(hgu133aENTREZID) thanks, -- Weiwei Shi, Ph.D Research Scientist GeneGO, Inc. Did you always know? No, I did not. But I

[R] problem building first package

2007-06-25 Thread Richard Chandler
Hi, I am trying to build a package from source for the first time. I'm using Windows XP. After R CMD INSTALL or R CMD check I get an error message that I don't understand. I've tried to follow the instrucions provided in the R Installation and Administration .pdf and the text file that comes with

Re: [R] eps in odfWeave

2007-06-25 Thread Dieter Menne
Kuhn, Max Max.Kuhn at pfizer.com writes: Thanks, Max. What is the version of odfWeave? If you use those specifications in the current version (0.5.9): Yes, version 0.5.9 (it was further below in my message, I forgot to put the odfWeave flag before it. I had tried your horizontal...

[R] simultaneous actions of grep ???

2007-06-25 Thread Ana Patricia Martins
Hello R-users and developers, Once again, I'm asking for your help. There is other way to do the same more easily for applied simultaneous grep??? c-subset(c,!rownames(c) %in% grep(.1,rownames(c),value=T)) c-subset(c,!rownames(c) %in% grep(.5,rownames(c),value=T))

Re: [R] a string to enviroment or function

2007-06-25 Thread jim holtman
I think that you might want: t0 - (paste(hgu133a, ENTREZID, sep=)) xx - as.list(get(t0)) # make it work like xx-as.list(hgu133aENTREZID) On 6/25/07, Weiwei Shi [EMAIL PROTECTED] wrote: Hi, I am wondering how to make a function Fun to make the following work: t0 - (paste(hgu133a,

[R] manipulate a matrix

2007-06-25 Thread Jon Hak
I have read everything I can find on how to manipulate a results matrix in R and I have to admit I'm stumped. I have set up a process to extract a dataset from ArcGIS to compute a similarity index (Jaccards) in Vegan. The dataset is fairly simple, but large, and consists of rows = sample

Re: [R] Re : Half of a heatmap

2007-06-25 Thread Gavin Kelly
I am trying to produce a heatmap of pairwise correlations, but since the matrix is symmetric, I only need either the upper or the lower triangle. I have scoured the web and R documentation, but I have not been able to find a way to produce such a figure. Is there a simple way to produce

Re: [R] R behaviour related to user input (readline()) and run selection

2007-06-25 Thread jim holtman
First of all you don't appear to be returning anything from your function, so 'co' is NULL and your 'if' statement fails. You were also probably cut/paste the code in, and the readline was reading from the statements you pasted. Try sourcing in the code. On 6/25/07, Bernzweig, Bruce

Re: [R] how to plot this?

2007-06-25 Thread jim holtman
You might want to check out this link to the type of graphs that R can produce and find one you like; the code will be with it. http://addictedtor.free.fr/graphiques/allgraph.php On 6/25/07, Weiwei Shi [EMAIL PROTECTED] wrote: Hi, there: Suppose I have a couple of data.frames and each one

Re: [R] how to plot this?

2007-06-25 Thread jim holtman
You might want to check out this link to the type of graphs that R can produce and find one you like; the code will be with it. http://addictedtor.free.fr/graphiques/allgraph.php On 6/25/07, Weiwei Shi [EMAIL PROTECTED] wrote: Hi, there: Suppose I have a couple of data.frames and each one

[R] degrees of freedom in lme

2007-06-25 Thread Jean-Baptiste Ferdy
Dear all, I am starting to use the lme package (and plan to teach a course based on it next semester...). To understand what lme is doing precisely, I used balanced datasets described in Pinheiro and Bates and tried to compare the lme outputs to that of aov. Here is what I obtained:

Re: [R] simultaneous actions of grep ???

2007-06-25 Thread Marc Schwartz
On Mon, 2007-06-25 at 17:15 +0100, Ana Patricia Martins wrote: Hello R-users and developers, Once again, I'm asking for your help. There is other way to do the same more easily for applied simultaneous grep??? c-subset(c,!rownames(c) %in% grep(.1,rownames(c),value=T))

Re: [R] problem building first package

2007-06-25 Thread Prof Brian Ripley
You can ignore the message about latex (and 2.5.1 RC does not give it). hhc is part of HTML Help Workshop: you either have not installed that or not put it in your path. The settings in MkRules affect the types of help for building R, not 'R CMD INSTALL'. If you only want text help you need

[R] LanguageR pvals.fnc error message

2007-06-25 Thread Mike Ford
Hi. I get an error message about not converging when I try and use the pvals.fnc from the languageR library. The LMER analysis worked fine (See below). I am not an expert so I don't understand why the LMER worked but not the pvals.fnc Any help gratefully received. - Mike AIC BIC logLik

Re: [R] Re : Half of a heatmap

2007-06-25 Thread Brad McNeney
Just for the record, LDheatmap can display any upper-triangular matrix of measures between 0 and 1, not just LD measures. Users with their own matrix should pass it as the first argument (gdat) to the function. -b -- Brad McNeney Statistics and Actuarial Science Simon Fraser University Burnaby,

Re: [R] manipulate a matrix

2007-06-25 Thread jim holtman
Is this what you want? x [,1] [,2] [,3] [,4] [,5] [1,]16 11 16 21 [2,]27 12 17 22 [3,]38 13 18 23 [4,]49 14 19 24 [5,]5 10 15 20 25 cbind(row=as.vector(row(x)), col=as.vector(col(x)), value=as.vector(x)) row col

Re: [R] degrees of freedom in lme

2007-06-25 Thread Doran, Harold
This is such a common question that it has a an FAQ-like response from Doug Bates. Google lmer p-values and all that to find the response. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jean-Baptiste Ferdy Sent: Monday, June 25, 2007 12:26 PM

Re: [R] simultaneous actions of grep ???

2007-06-25 Thread Marc Schwartz
On Mon, 2007-06-25 at 11:58 -0500, Marc Schwartz wrote: On Mon, 2007-06-25 at 17:15 +0100, Ana Patricia Martins wrote: Hello R-users and developers, Once again, I'm asking for your help. There is other way to do the same more easily for applied simultaneous grep???

Re: [R] manipulate a matrix

2007-06-25 Thread Marc Schwartz
On Mon, 2007-06-25 at 10:58 -0400, Jon Hak wrote: I have read everything I can find on how to manipulate a results matrix in R and I have to admit I'm stumped. I have set up a process to extract a dataset from ArcGIS to compute a similarity index (Jaccards) in Vegan. The dataset is fairly

Re: [R] how to plot this?

2007-06-25 Thread hadley wickham
On 6/25/07, jim holtman [EMAIL PROTECTED] wrote: You might want to check out this link to the type of graphs that R can produce and find one you like; the code will be with it. http://addictedtor.free.fr/graphiques/allgraph.php Or for examples using the ggplot2 package:

Re: [R] problem building first package

2007-06-25 Thread Richard Chandler
Thanks so much. I was able to get it installed using: R CMD INSTALL --docs=normal removal_1.0.tar.gz However, setting --docs=txt gave me this error: ERROR: invalid --docs value `txt' The --help says it needs to be a list, but no matter, --docs=normal worked fine. Quoting Prof Brian Ripley

Re: [R] Re : Half of a heatmap

2007-06-25 Thread Rosenfeld, Jeffrey
I have correlation values between -1 and 1 and I need to keep the sign. Jeff -Original Message- From: Brad McNeney [mailto:[EMAIL PROTECTED] Sent: Mon 6/25/2007 1:01 PM To: Rosenfeld, Jeffrey Cc: Neil Shephard; r-help@stat.math.ethz.ch Subject: Re: [R] Re : Half of a heatmap Just for

Re: [R] problem building first package

2007-06-25 Thread Prof Brian Ripley
It's a buglet: --docs=txt, works. On Mon, 25 Jun 2007, Richard Chandler wrote: Thanks so much. I was able to get it installed using: R CMD INSTALL --docs=normal removal_1.0.tar.gz However, setting --docs=txt gave me this error: ERROR: invalid --docs value `txt' The --help says it needs

Re: [R] complex contrasts and logistic regression

2007-06-25 Thread Nicholas Lewin-Koh
Hi, Sorry to take so long to reply, I was travelling last week. Thanks for your suggestions. Actually in this case contrast and predict gave the same result, and what I was looking at was the correct odds from the model. What is still confusing me is the 1st part of my question, looking for a

Re: [R] eps in odfWeave

2007-06-25 Thread Dieter Menne
Dieter Menne dieter.menne at menne-biomed.de writes: library(odfWeave) plotInfo - getImageDefs() plotInfo$type = eps plotInfo$plotWidth = 4 plotInfo$plotHeight = 4 plotInfo$device = postscript plotInfo$args = list( horizontal = FALSE,onefile = FALSE, paper = special)

Re: [R] ANOVA non-sphericity test and corrections (eg, Greenhouse-Geisser)

2007-06-25 Thread John Wilkinson
Darren, Further to Peter Dalgaard's help; Take a look at the example in --- library(car) ?Anova # note upper case 'A' The example in the Anova help page following the --- ## a multivariate linear model for repeated-measures data ## See ?OBrienKaiser for a description of the data set used

Re: [R] Post-hoc tests for interactions of between- and within-subject factors

2007-06-25 Thread Richard M. Heiberger
The glht function in package multcomp requires aov objects. As you discovered, it does not work with aovlist objects. See the documentation ?MMC in the HH package from CRAN. The maiz example is similar to your example. Come back to the list if you need further help. Rich

Re: [R] degrees of freedom in lme

2007-06-25 Thread Manuel Morales
On Mon, 2007-06-25 at 13:15 -0400, Doran, Harold wrote: This is such a common question that it has a an FAQ-like response from Doug Bates. Google lmer p-values and all that to find the response. Isn't this a different question, though, since Jean-Baptiste is using nlme. Details on the

[R] correlation structure

2007-06-25 Thread Daniel O'Shea
I have been using a nlme model and wish to specify the correlation structure. My data is grouped (bas), but I have no time component or adequate spatial description beyond the grouping variable. So I chose the simplest structure and updated my original nlme model by: update(model,

[R] debug()

2007-06-25 Thread Nair, Murlidharan T
I am having problems using debug() when I am running R under linux. It works fine under the windows environment. Is there anything special that need to specify when running under linux. I would appreciate your help. Cheers../Murli [[alternative HTML version deleted]]

Re: [R] a string to enviroment or function

2007-06-25 Thread Weiwei Shi
then how to do this f1 - function(mylab){ library(mylab) ... } it seems that if you call library(hgu133a) # which is file # but library(mylab) # even you pass hgu133a as parameter, it still complains about mylab does not exist. It seems that it consider mylab as package instead of its value.

[R] How to shadow 'power' area?

2007-06-25 Thread Andrej Kastrin
Dear all, Suppose I plot two normal distributions (A and B) side by side and add vertical line which hipotheticaly represent alpha value; e.g.: x - seq(-3.5,5, length=1000) y - dnorm(x) # Plot distribution A plot(y~x, type='l',axes=F,xlab=,ylab=,lwd=2) # Plot distribution B y2 - dnorm(x-1.5)

Re: [R] How to shadow 'power' area?

2007-06-25 Thread Christos Hatzis
Look at the following link: http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=88 This should be pretty close to what you want. HTH -Christos -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Andrej Kastrin Sent: Monday, June 25, 2007

Re: [R] How to shadow 'power' area?

2007-06-25 Thread Marc Schwartz
On Mon, 2007-06-25 at 22:50 +0200, Andrej Kastrin wrote: Dear all, Suppose I plot two normal distributions (A and B) side by side and add vertical line which hipotheticaly represent alpha value; e.g.: x - seq(-3.5,5, length=1000) y - dnorm(x) # Plot distribution A plot(y~x,

Re: [R] LanguageR pvals.fnc error message

2007-06-25 Thread Ben Bolker
Mike Ford mford at csl.psychol.cam.ac.uk writes: Hi. I get an error message about not converging when I try and use the pvals.fnc from the languageR library. The LMER analysis worked fine (See below). I am not an expert so I don't understand why the LMER worked but not the pvals.fnc

[R] Problem installing MCMCpack

2007-06-25 Thread Vivek Menon
I am having a problem installing MCMCpack on a Linux ppc64 machine. Has someone had this problem??This is what I see when I do an install package: install.packages(c(MCMCpack)) trying URL ' http://www.sourcekeg.co.uk/cran/src/contrib/MCMCpack_0.8-2.tar.gz' Content type 'application/x-tar' length

[R] fft and the derivative

2007-06-25 Thread Todd Remund
Can one take f(t) and transform to F(omega) in the frequency domain using fft(), and use the properties of the fft and find the derivative of f(t)? For example, f(t) - F(omega) = f(t)^n - (i*omega)^n * F(omega) Use this and get, f(t)^n = F^(-) [ (i*omega)^n * F(omega) ] to get the nth

[R] Quiver in 1d

2007-06-25 Thread fernando espindola
Dear all R-users, I am try to plot quiver in one dimension, but there is only code for two dimension in different R web pages. Somebody know same code that can help me in plot quiver in one dimension. Any help gratefully received . Thanks for all Fernando

Re: [R] a string to enviroment or function

2007-06-25 Thread Martin Morgan
Weiwei See ?library and the character.only argument. f - function(x) library(x) f(hgu95av2) Error in library(x) : there is no package called 'x' f - function(x) library(x, character.only=TRUE) f(hgu95av2) search() [1] .GlobalEnvpackage:hgu95av2 package:stats [4]

Re: [R] a string to enviroment or function

2007-06-25 Thread Duncan Murdoch
On 25/06/2007 4:19 PM, Weiwei Shi wrote: then how to do this f1 - function(mylab){ library(mylab) ... } it seems that if you call library(hgu133a) # which is file # but library(mylab) # even you pass hgu133a as parameter, it still complains about mylab does not exist. It seems

[R] changing the position of the y label (ylab)

2007-06-25 Thread Etienne
How can I change the position of the ylab, after enlarging the margins with par(mar=...)? Here is the relevant code snippet par(mar=c(5.1,5.1,4.1,2.1)) plot(c(1979,2003),c(40,50),ylim=c(1,73),lab=c(20,10,1),pch=21,col='blue',bg='blue',axes=FALSE,xlab=Years,ylab=Onset/Withdrawl

Re: [R] How to shadow 'power' area?

2007-06-25 Thread Greg Snow
You can look at the power.examp function in the TeachingDemos package to see if it gives you the graph you want, or look at the code to see how you could modify it to give you the plot you want. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED]

Re: [R] degrees of freedom in lme

2007-06-25 Thread Doran, Harold
I don't believe the issue is different. In some respects, lme and lmer essentially accomplish the same thing for mixed linear models. lmer differs in many ways than lme in how estimates are computed and in the type of models it can handle. It is probably best to read the new vignettes in the

Re: [R] fft and the derivative

2007-06-25 Thread Ravi Varadhan
Todd, Your idea is correct for continuous Fourier transform, but I am not sure how one could apply that to fft, which corresponds to the discrete Fourier transform. For instance, what values of omega would you use for the term i*omega to get the discrete fourier transform of the derivative of

Re: [R] changing the position of the y label (ylab)

2007-06-25 Thread Etienne
Thanks (Merci) Christophe! that did it --- Christophe Bonenfant [EMAIL PROTECTED] wrote: Hi Etienne - consider to use the mtext function: par(mar=c(5.1,5.1,4.1,2.1)) plot(c(1979,2003),c(40,50),ylim=c(1,73),lab=c(20,10,1), pch=21,col='blue',bg='blue',axes=FALSE,xlab=,ylab=,font.lab=2)

Re: [R] running Rcmdr

2007-06-25 Thread Manuel
Thanks for your answer, but i think i dont do correctly my question. I need the command line to run Rmdr, like that: R Rcmdr or R loadRcmdr.R where loadRcmdr has library(Rcmdr). or something like that. I tried the last example, but when Rcmdr is executed, later it is closed. About

Re: [R] R-excel

2007-06-25 Thread Moshe Olshansky
Also try xlsReadWrite package on CRAN. --- Erika Frigo [EMAIL PROTECTED] wrote: Good morning to everybody, I have a problem : how can I import excel files in R??? thank you very much Dr.sa. Erika Frigo Università degli Studi di Milano Facoltà di Medicina Veterinaria Dipartimento

[R] Has anyone tryed out my software?

2007-06-25 Thread Robert Wilkins
Hello all, Has anyone ( who uses a Linux desktop ) tryed out my stuff I mentioned a few weeks ago? Perhaps installed it and run a couple of example programs? If you have, tell me what you think. Robert ( it's the tarball in the download section at http://code.google.com/p/vilno , discussed

[R] a-priori orthogonal contrasts

2007-06-25 Thread Keith Alan Chamberlain
Dear R-helpers, I have some stacked data with which I want to generate orthogonal contrasts - with the caveat that I determine what means comparisons are taking place a-priori. Does it matter for e.g. 'contr.SUM' that the data are stacked, and does 'contr.SUM' allow me to set the contrasts

Re: [R] R-excel

2007-06-25 Thread Erich Neuwirth
The R(D)COM server (contaiong the RExcel Excel addin) and/or the rcom package allow (among other things) to select a range in Excel and directly transfer it to R as an array or as a dataframe. It only works on Windows with Excel and R installed. More information on these packages is available at

[R] Howto fix colors in colorRampPalette?

2007-06-25 Thread Ashoka Polpitiya
Hi All, I am plotting a correlation matrix using the heatmap function. I create my color palette as follows (say, for red-green palette): colorRampPalette( c(green, black, red), space=rgb)(10) But I want to fix colors for a specific range of values, i.e. green=-1, black=0, red=+1. The purpose of

Re: [R] R-excel

2007-06-25 Thread Stephen Tucker
There are also some notes about this in the R Data Import/Export manual: http://cran.r-project.org/doc/manuals/R-data.html#Reading-Excel-spreadsheets But I've gathered the following examples from the R-help mailing list archives [in addition to the option of saving the spreadsheet as a .csv file