[R] Call R from Perl

2009-05-06 Thread Daniel Fernandez
Hi, I am trying to call R in my perl program but is not working. It does not give an error but my program gets stuck indefinitely. I am using the following command to call R from Perl `R BATCH regressExpr.spl temp`; Any suggestions? Best, Daniel F. [[alternative HTML version

Re: [R] Heatmap without levelplot

2009-05-06 Thread Antje
Hi Uwe, thanks a lot for your answer! And thanks a lot to all others helping me with this issue! Uwe Ligges schrieb: Antje wrote: Hi Uwe, I tried to explain my problem with the given example. I don't see any documentation which tells me that the length of col.regions should be one less

[R] Duplicating meta-regression results from PROC MIXED with lmer

2009-05-06 Thread Brant Inman
R-experts: In 2002, Hans Van Houwelingen et al. published a tutorial on how to do meta-regression in Statistics in Medicine. They used the classic BCG dataset of Colditz to demonstrate correct methodology and computed the results using PROC MIXED in SAS. In trying to duplicate the results

[R] I need your help.

2009-05-06 Thread 풍이
Hi, I need your help, so I send letter to you. I have a problem about plsr in pls package. I want to show how classfied or related each ohter samples, so I tried to use plsr and biplot. But, I failed. Because, I had to change data type of my sample. Unfortunately, I didn't know

[R] Rcompression and Java Deflator

2009-05-06 Thread Stabler, Ben
(this may be a duplicate post since I attached a file to the previous try...sorry about that) Below are the first few lines of a zlib compressed byte array written from Java with the Deflator class. readBin(row_1,raw(),1000) [1] 4c 45 50 e2 49 d5 86 bc 48 a1 32 5d 49 9d f5 90 48 e0

[R] subscript out of bounds

2009-05-06 Thread PDXRugger
Consider - if (BinNumber==1) Loc_Prob=Bin1Main.data[findInterval(Dev_Size,Bin1Main.data[,1])+1,2] . .. ... if (BinNumber==10) Loc_Prob=Bin10Main.data[findInterval(Dev_Size,Bin10Main.data[,1])+1,2] BinNumber is just referencing 1 of 10 sets of data so: where BinNumber = 1 and

[R] help

2009-05-06 Thread aliasghar talaei
I noticed that both in the fArma and fracdiff packages it is not possible to predict an object fitted with an ARFIMA model. Moreover, the ox link is no more available. However, is it possible to find the coefficient d of differentiation of a time series, and then predict the time series by apply

[R] Determining the contribution of individual variables to LOF

2009-05-06 Thread Dr D.W. Murray
Hi I have been using the LOFactor in the dprep package to identify outliers. I need to communicate to an end user, which variables have contributed the most to an individual record being identified as an outlier. Does such functionality exist in the LOF arena? I was hoping that it would

Re: [R] illustrating a report

2009-05-06 Thread Wydou
Philipp, Sweave is exactly what I was looking for. If you are familiar with LaTeX, then it's the tool of choice. Thank you for the tip! Stephan -- View this message in context: http://www.nabble.com/illustrating-a-report-tp21328500p23400798.html Sent from the R help mailing list archive at

[R] threat analysis

2009-05-06 Thread ms.com
dear all i am in great need of threat analysis technique, could anyone give me information about this method in R, which package does the Threat Analysis (in animals) thanking you Ms.com _ Hotmail® has a new way to see what's

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread viktoras d.
well, I am less proficient in R comparing with other tools/languages. Therefore my biased opinion is - it is possible in R, but it may be easier if you use other tools, especially if you have to build a user-friendly GUI. The most accessible (although limited to MS Windows only) method would

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Zeljko Vrba
Sorry for reply to the wrong person, I lost the original email. Farrel Buchinsky wrote: Is R an appropriate tool for data manipulation and data reshaping and data organizing? I think so but someone who recently joined our group thinks not. The new recruit believes that python or another

Re: [R] A question about using “by” in GAM model fitting of interaction between smooth terms and factor

2009-05-06 Thread Simon Wood
The problem here is that the help page you are looking at appears to be from an earlier version of `mgcv' than you are using (it's from a version that did not support factor `by' variables). Take a look at ?gam.models for the version that you are actually using. The reason that your models

[R] QUADRATIC LINK FUNCTIONS FOR MLE ESTIMATE OF NON-STATIONARY GEV FITS

2009-05-06 Thread Alberto Bordallo
Hi All, I am a newcomer to S/R. Could you please let me know how to model quadratic trends for the mul/sigl link functions when fitting non-stationary GEV distributions using the ismev package? Thanks [[alternative HTML version deleted]] __

Re: [R] Heatmap without levelplot

2009-05-06 Thread Uwe Ligges
Antje wrote: Hi Uwe, thanks a lot for your answer! And thanks a lot to all others helping me with this issue! Uwe Ligges schrieb: Antje wrote: Hi Uwe, I tried to explain my problem with the given example. I don't see any documentation which tells me that the length of col.regions

Re: [R] R 2.9 and XML

2009-05-06 Thread Uwe Ligges
Ian Kennedy wrote: After I installed R 2.9.0 my XML package cannot load; I get an error This application has failed to start because iconv.dll was not found. Re-installing the application may fix this problem. I already re-installed both but the problem persists. Does anyone know what is

Re: [R] Problem with update.packages on repo with a single package

2009-05-06 Thread Uwe Ligges
Is this R-2.9.0? If so, please submit to R-bugs so that the bugfix does not get lost. Otherwise, please try R-2.9.0 patched and submit a patch against that version. Thank you very much! Uwe Ligges Michael DiPersio wrote: I've set up a CRAN-style repository to distribute a single package.

Re: [R] Problem with update.packages on repo with a single package

2009-05-06 Thread Uwe Ligges
Uwe Ligges wrote: Is this R-2.9.0? If so, please submit to R-bugs so that the bugfix does not get lost. Otherwise, please try R-2.9.0 patched and submit a patch against that version. Whoops, no! Already fixed this morning by Brian Ripley. So, pelase don't send! Thanks again, Uwe

Re: [R] big databases

2009-05-06 Thread Uwe Ligges
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and tell us a bit about your machine, OS, R version, the am,ount of RAM used by your workspace. Pseudo Phobic wrote: Dear all, I have a .dta database which is about 400 MB. I cannot open it though I have no

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Emmanuel Charpentier
Le mercredi 06 mai 2009 à 00:22 -0400, Farrel Buchinsky a écrit : Is R an appropriate tool for data manipulation and data reshaping and data organizing? [ Large Snip ! ... ] Depends on what you have to do. I've done what can be more or less termed data management with almost uncountable tools

Re: [R] how to modify a function in a R package calling other (invisible?) functions

2009-05-06 Thread Uwe Ligges
Krusty the Klown wrote: Hey hey kids! Kid Klown, it is polite to sign messages with your real name and affiliation, if not clear from the mailing address. Otherwise people might be offended and do not reply. I'm facing this rough problem: I have to modify a function in a R package

Re: [R] View cluster results of large data sets

2009-05-06 Thread Uwe Ligges
Wade Wall wrote: Hi all, I am trying to run a cluster analysis with a relatively large data set (545 samples) and obviously am having trouble viewing the results using standard plot() method. Is there a way to plot the results of a large cluster analysis in a window with a scroll bar. Any

Re: [R] Quantile Regression for Longitudinal Data. Warning message: In rq.fit.sfn

2009-05-06 Thread Lola Gadea
Dear Dimitris, I have exactly the same problem than you, Do you get some solution? Thanks, Lola Lola Gadea Profesora titular de Economía Aplicada/Lecturer in Applied Economics Universidad de Zaragoza/University of Zaragoza (Spain) lga...@unizar.es

Re: [R] problems with rq.fit.sfn

2009-05-06 Thread Lola Gadea
Dear Dimitris, I have exactly the same problem than you...Do you have some any solution? Thanks, Lola Lola Gadea Profesora titular de Economía Aplicada/Lecturer in Applied Economics Universidad de Zaragoza/University of Zaragoza (Spain) lga...@unizar.es

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Philipp Pagel
On Wed, May 06, 2009 at 12:22:45AM -0400, Farrel Buchinsky wrote: Is R an appropriate tool for data manipulation and data reshaping and data organizing? I think so but someone who recently joined our group thinks not. The new recruit believes that python or another language is a far better

Re: [R] Bristol mirror GPG problem ubuntu repository

2009-05-06 Thread Daniel Brewer
Anyone know how to find the details of the mirror maintainers then? because I can't find it. Dan Uwe Ligges wrote: Daniel Brewer wrote: Hello, I am getting a GPG error with the ubuntu repository at the bristol UK mirror. When my source.list has this line: deb

Re: [R] Bristol mirror GPG problem ubuntu repository

2009-05-06 Thread Uwe Ligges
Daniel Brewer wrote: Hello, I am getting a GPG error with the ubuntu repository at the bristol UK mirror. When my source.list has this line: deb http://www.stats.bris.ac.uk/R/bin/linux/ubuntu/ intrepid/ On an apt-get update you get this: W: GPG error: http://www.stats.bris.ac.uk intrepid/

Re: [R] View cluster results of large data sets

2009-05-06 Thread Jim Lemon
Wade Wall wrote: Hi all, I am trying to run a cluster analysis with a relatively large data set (545 samples) and obviously am having trouble viewing the results using standard plot() method. Is there a way to plot the results of a large cluster analysis in a window with a scroll bar. Any

Re: [R] Plotting pairs of bars

2009-05-06 Thread Steve Murray
Thanks for the reply - the 'beside' argument certainly looks useful, although I'm still not getting the output I'd hoped for. By doing: barplot(combine86[,1:2], beside = TRUE, las = 1, xlab=rownames(combine86)) ...I get all the bars for the 'Sim Mean' column plotted on the left side of the

Re: [R] quick square root axes

2009-05-06 Thread Jim Lemon
Markus Loecher wrote: Dear R users, while I enjoy the built-in log argument to the plot() function, I wished it would be as easy to create more general custom transformed axes such as sqrt(), logit, etc... for example, instead of plot(x=exp(rnorm(10)), y=(1:10)^4, log = xy), sth. along the

Re: [R] Rcompression and Java Deflator

2009-05-06 Thread Andrew Dunn
Hi Ben, I have successfully used rcompression with archives created in Java. I have not used the Deflater to do it but rather the GZipOutputStream. (http://java.sun.com/j2se/1.5.0/docs/api/java/util/zip/GZIPOutputStream.html) in combination with a java based tar archive generator

Re: [R] Plotting pairs of bars

2009-05-06 Thread Jim Lemon
Steve Murray wrote: Thanks for the reply - the 'beside' argument certainly looks useful, although I'm still not getting the output I'd hoped for. By doing: barplot(combine86[,1:2], beside = TRUE, las = 1, xlab=rownames(combine86)) ...I get all the bars for the 'Sim Mean' column plotted on the

Re: [R] Plotting pairs of bars

2009-05-06 Thread Steve Murray
Jim and all, Thanks for the suggestion, however, I get the following error: barplot(t(combine86[,1:2], beside = TRUE, las = 1)) Error in t(combine86[, 1:2], beside = TRUE, las = 1) : unused argument(s) (beside = TRUE, las = 1) I've looked up ?t and cannot see any extra arguments that I

[R] Can't find package for decision plot anymore

2009-05-06 Thread bartjoosen
Hi, About a year ago I found some references to a package which takes a dataframe and derived rules from it. If I remember correctly, there were examples shown from the Iris data. You could choose between a graphical representation with balloons with choises, or a text version: if x y then z.

Re: [R] Plotting pairs of bars

2009-05-06 Thread Ben Bolker
Steve Murray-3 wrote: Jim and all, Thanks for the suggestion, however, I get the following error: barplot(t(combine86[,1:2], beside = TRUE, las = 1)) Error in t(combine86[, 1:2], beside = TRUE, las = 1) : unused argument(s) (beside = TRUE, las = 1) I've looked up ?t and

Re: [R] Cox Proportional Hazard with missing covariate data

2009-05-06 Thread Terry Therneau
I concur with your other reply - you have left-trunctated data. Left-truncated: Any parts that failed before 2001 (start of study) are not included in the study. Left-censored: Parts that failed before the study start are included, but with an indeterminate failure date ( 20001).

Re: [R] Call R from Perl

2009-05-06 Thread Kevin E. Thorpe
Daniel Fernandez wrote: Hi, I am trying to call R in my perl program but is not working. It does not give an error but my program gets stuck indefinitely. I am using the following command to call R from Perl `R BATCH regressExpr.spl temp`; Any suggestions? Best, Daniel F. Does `R CMD

Re: [R] Call R from Perl

2009-05-06 Thread Daniel Fernandez
I'm not sure but my guess is that part isn't working... is that the way to call R from perl? Sent from my iPhone On May 6, 2009, at 8:09 AM, Kevin E. Thorpe kevin.tho...@utoronto.ca wrote: Daniel Fernandez wrote: Hi, I am trying to call R in my perl program but is not working. It does

Re: [R] Call R from Perl

2009-05-06 Thread Kevin E. Thorpe
Daniel Fernandez wrote: I'm not sure but my guess is that part isn't working... is that the way to call R from perl? It's the way I would call R in batch from the command line. In your call you are missing the CMD part. Sent from my iPhone On May 6, 2009, at 8:09 AM, Kevin E. Thorpe

[R] model fitting using by(): how to get fitted values?

2009-05-06 Thread xavier . chardon
Hi all, I'm doing nonlinear regressions on data with several factors. I want to fit say a logistic curve with different parameter values for each factor level. So I'm doing something like: tmp - by( myData, list(myFactor1, myFactor2), function(x) nls(...) ) It works fine. However, I could not

Re: [R] Create Pie chart from .csv file

2009-05-06 Thread DonkeyRhubarb
Im more than happy to use a barplot, but at this stage time is just so short I need something to show. I can find examples of plotting a line of numbers that has been manually typed into R, but Im missing something to do it from a CSV file. The name of my csv file is OutputFromDB. Would you

[R] sigma in glmer (lme4)

2009-05-06 Thread Eric Elguero
dear R-users, I am trying to understand what is the sigma parameter returner by glmer I thought it was (an estimate of) the sigma parameter defined by Mc Cullagh Nelder (e.g. p 126 of 2nd edition) but I ran some simulations and it seems that this is something else. I simulated data

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Frank E Harrell Jr
I second what Zeljko wrote. In addition, see the data manipulation section in Chapter 4 of http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RS/sintro.pdf Frank Zeljko Vrba wrote: Sorry for reply to the wrong person, I lost the original email. Farrel Buchinsky wrote: Is R an appropriate

Re: [R] model fitting using by(): how to get fitted values?

2009-05-06 Thread Duncan Murdoch
On 06/05/2009 8:22 AM, xavier.char...@free.fr wrote: Hi all, I'm doing nonlinear regressions on data with several factors. I want to fit say a logistic curve with different parameter values for each factor level. So I'm doing something like: tmp - by( myData, list(myFactor1, myFactor2),

[R] DSC 2009 update, early registration ends May 15

2009-05-06 Thread Peter Dalgaard
As previously announced, the sixth international workshop on Directions in Statistical Computing (DSC2009) will take place at the Center for Health and Society, University of Copenhagen, Denmark, 13th-14th of July 2009. For details, see http://www.R-project.org/dsc-2009. - The paper selection

Re: [R] Plotting pairs of bars

2009-05-06 Thread Jim Lemon
Steve Murray wrote: Jim and all, Thanks for the suggestion, however, I get the following error: barplot(t(combine86[,1:2], beside = TRUE, las = 1)) Error in t(combine86[, 1:2], beside = TRUE, las = 1) : unused argument(s) (beside = TRUE, las = 1) I've looked up ?t and cannot see

[R] X and Y axis labels on a barplot

2009-05-06 Thread Steve Murray
Dear all, I have produced a barplot and wish to alter the axes a little. In place of the variable names appearing on the x-axis, I'd like to have the numbers 1986 to 1995. I have tried using the argument xlim=c(1986,1995) in the barplot command but receive: Error in plot.window(xlim, ylim,

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread hadley wickham
Take a look at plyr and reshape packages (http://had.co.nz/), I have a hunch that they would have saved me a lot of headache had I found out about them earlier :) As the author of these two packages, I'm admittedly biased, but I think R is unparalleled for data preparation, manipulation, and

[R] Summary help

2009-05-06 Thread mathallan
Hi, I have fittet a gamma model, and is wondering if I can read the shape and the scale direct from the summary Estimate Std. Errort valuePr(|t|) (Intercept) 1.612e+00 4.735e-02 34.052 2e-16 *** myvalue 3.564e-02

[R] PGF Device

2009-05-06 Thread Lasse Bombien
Hi all, I saw a thread from 2007 about the possible implementation of a PGF (http://sourceforge.net/projects/pgf/) graphics device. Does anyone know if something came of it? I would very much like if that was the case. Regards, Lasse __

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread John Kane
--- On Wed, 5/6/09, Farrel Buchinsky fjb...@gmail.com wrote: Is R an appropriate tool for data manipulation and data reshaping and data organizing? I think so but someone who recently joined our group thinks not. I only do small scale projects and am by no means a programmer. Isn't Perl

[R] tapply changing order of factor levels?

2009-05-06 Thread Chirantan Kundu
Hi, Does tapply change the order when applied on a factor? Below is the code I tried. mylevels-c(IN0020020155,IN0019800021,IN0020020064) mydata-c(IN0020020155,IN0019800021,IN0020020064,IN0020020155,IN0019800021,IN0019800021,IN0020020064,IN0020020064,IN0019800021)

Re: [R] X and Y axis labels on a barplot

2009-05-06 Thread Steve Murray
A very reasonable request! Sorry for not doing this initially, but please find below the data I am trying to plot: total_sums sums86 sums87 sums88 sums89 sums90 sums91 sums92 Sim_1986 17722203 16875889 18626582 18428415 17611182 17290016 16819289 X198615276602

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Laura Arsanto
I used R for my master thesis (with big effort, anyway) and now I find difficult to use R in my daily work, becasue it has really serious problems with datasets of big dimension, both in the data manipulation step and in the analysis step. But I really would love to use it, as I like its

Re: [R] tapply changing order of factor levels?

2009-05-06 Thread jim holtman
The result of 'tapply' is just a named vector and the names are in alphabetical order. If you want them printed in a different order, then you have to specify it. Since you have the order in 'mylevels', this will work: str(mysummary) int [1:3(1d)] 4 3 2 - attr(*, dimnames)=List of 1 ..$ :

Re: [R] tapply changing order of factor levels?

2009-05-06 Thread Chirantan Kundu
Thanks, Jim. It works! On Wed, May 6, 2009 at 7:26 PM, jim holtman jholt...@gmail.com wrote: The result of 'tapply' is just a named vector and the names are in alphabetical order. If you want them printed in a different order, then you have to specify it. Since you have the order in

Re: [R] X and Y axis labels on a barplot

2009-05-06 Thread Steve Murray
dput(total_sums) structure(c(17722202.6898231, 15276602.215475, 16875888.5155229, 14086271.625756, 18626581.9628846, 15387747.481166, 18428414.8535184, 15560882.404998, 17611181.5207881, 14905453.195546, 17290016.3934661, 14939493.120707, 16819288.8227961, 13979000.614402, 17657959.3656573,

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Simon Pickett
My institute uses SAS religiously, I am the only R heathen. I have resisted learning to use SAS because I dont see the point after years of using R and I like being able to do everything using one program. However, my colleagues maintain that SAS is better for programming without really ever

Re: [R] tapply changing order of factor levels?

2009-05-06 Thread Alain Guillet
Hi, I don't believe the problem is related to tapply. I would say it is because of the factor. In fact, the order of a factor is given by the alphanumerical order of his levels. You can see it with levels(myfactor). I you want to change the order, redefine the levels of myfactor with the

Re: [R] X and Y axis labels on a barplot

2009-05-06 Thread jim holtman
This will put the years on the plot: barplot(total_sums, beside=TRUE, col=(2:3), las=2, names.arg=1986:1995) On Wed, May 6, 2009 at 10:03 AM, Steve Murray smurray...@hotmail.comwrote: dput(total_sums) structure(c(17722202.6898231, 15276602.215475, 16875888.5155229, 14086271.625756,

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread C.H.
There is a book on data manipulation using R. Data manipulation with R. http://www.springer.com/statistics/computational/book/978-0-387-74730-9 It highlighted how comprehensive the data manipulation capabilities of R can be. Regards, CH On Wed, May 6, 2009 at 10:01 PM, Simon Pickett

Re: [R] PGF Device

2009-05-06 Thread Dieter Menne
Lasse Bombien lasse at phonetik.uni-muenchen.de writes: I saw a thread from 2007 about the possible implementation of a PGF (http://sourceforge.net/projects/pgf/) graphics device. Does anyone know if something came of it? I am always using pgf to make drawings together with Sweave. See

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Paul Emberson
I also use the approach Philipp describes below. I use Python and shell scripts for processing thousands of input files and getting all the data into one tidy csv table. From that point onwards it's R all the way (often with the reshape package). Paul Philipp Pagel wrote: On Wed, May 06,

Re: [R] X and Y axis labels on a barplot

2009-05-06 Thread Steve Murray
Thanks Jim, that's great. Based on the information in the previous messages, is it possible to change the y-axis as I'd hoped? Thanks again, Steve _ [[elided Hotmail spam]] __

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Simon Pickett
Thanks Laura, I deal with huge data sets and have to do alot of fancy juggling of data to get the job done in R. I have recently been granted access to a cluster at a university which means 64 bit machines with 8gb of memory, which could prove to be a saviour hopefully. Simon. -

Re: [R] tapply changing order of factor levels?

2009-05-06 Thread Chirantan Kundu
Hi Alain, I tried levels(myfactor) as you suggested. levels(myfactor) [1] IN0020020155 IN0019800021 IN0020020064 The order is preserved, no alphanumerical sorting done here. Regards. On Wed, May 6, 2009 at 7:35 PM, Alain Guillet alain.guil...@uclouvain.bewrote: Hi, I don't believe the

Re: [R] Summary help

2009-05-06 Thread Dieter Menne
mathallan mathanmath at gmail.com writes: Hi, I have fittet a gamma model, and is wondering if I can read the shape and the scale direct from the summary Estimate Std. Errort valuePr(|t|) (Intercept) 1.612e+00 4.735e-02 34.052

Re: [R] tapply changing order of factor levels?

2009-05-06 Thread Chirantan Kundu
Hi Ellison, Thanks for pointing it out. Bad miss on my part. Regards. On Wed, May 6, 2009 at 7:55 PM, S Ellison s.elli...@lgc.co.uk wrote: The problem is in the code. When you say mysummary-tapply(myfactor,mydata,length) mysummary you have used mydata as a factor and myfactor as the

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Mike Lawrence
I work in cognitive science where we collect one or more data files per participant in an experiment then merge those files to perform subsequent analyses. Sometimes some files are in wide format and others are in long format, necessitating reshaping. I've found R entirely satisfactory for this.*

[R] bug with RODBC ?

2009-05-06 Thread JP_RENAUD_(DTCB)
Dear R help, I've attach a dbf file (and my prog. below) from which I did not succeed in geting the right values for my variable densite... using RODBC. It returns 10 instead of 100 , without error messages !!! is it a bug with RODBC ? I've succeed in geting the right values using the

Re: [R] PGF Device

2009-05-06 Thread Lasse Bombien
Thanks, that's a good idea and I might use it for some purposes. But want I would like is a complete graphics device, which creates PGF output in the same way e.g. dev.copy2pdf() would. Lasse Am Mittwoch, den 06.05.2009, 14:15 + schrieb Dieter Menne: Lasse Bombien lasse at

[R] Une glacière campingaz offerte pour vos d éparts en weekend

2009-05-06 Thread Laura de Welcome Office
WELCOME OFFICE n°1 DU DISCOUNT AUX ENTREPRISES Si vous ne visualisez pas correctement ce message, cliquez ici :http://trc1.emv2.com/HM?a=A9X7Cq,jKYmd8XLLjqxobfPiTA Pour être certain de recevoir notre newsletter, nous vous conseillons d'ajouter servicecli...@welcomeoffice.emv1.net à votre carnet

Re: [R] [R-sig-ME] Duplicating meta-regression results from PROC MIXED with lmer

2009-05-06 Thread Brant Inman
Rolf, I actually don't believe that this is a SAS vs R issue since I have 3 sources that report the same results. I know that STATA, SAS and the mima function from R can all be used to give the correct results. The question is related more to how I can get similar results with lmer.

[R] Need help building combined model

2009-05-06 Thread E Cell
I have two kind of variables. One kind where the values are continious and ranges from -1 to 5, and the other is boolean (i.e. true or false). For example: V1   V2  V3 V4 V5  V6 0.1  0.4   1    0   -1   3.7 0.4  0.1   2   -1   0    3.7 and V1  V2  V3 0 1    0 0     0    1 I want to

[R] Estimating parameters in a nonlinear model with observations

2009-05-06 Thread mathr
Hi, I have a problem. I want to estimate some parameters in a function. I already have an empirical function (made from 100 observations), which I want to estimate the parameters from. The function is f(x) = 1-((a+1)b^x)/(a+b^x) f(x) in [0,1], x in [0,1]. I want to estimate a and b. I tried

Re: [R] tapply changing order of factor levels?

2009-05-06 Thread S Ellison
The problem is in the code. When you say mysummary-tapply(myfactor,mydata,length) mysummary you have used mydata as a factor and myfactor as the data. tapply has (correctly) used the ordered labels in the grouping factor (mydata) to label its output. If you did what you probably intended:

Re: [R] tapply changing order of factor levels?

2009-05-06 Thread Bert Gunter
Jim's advice is patently false. Please read ?tapply for correct details. Counterexample: y - rnorm(6) x - factor(rep(factor(letters[1:3],lev = letters[3:1]),2)) x [1] a b c a b c Levels: c b a tapply(y,x,mean) c b a 0.4545897 -1.0544782 0.4682773 Bert Gunter

[R] Age Calculation in years, months and days

2009-05-06 Thread utkarshsinghal
Hi All, Is there any way in R to calculate the difference between two dates in years, months days. date1 = as.Date(1993-11-23) date2 = as.Date(2009-04-15) I want my output to be:15 years 4 months 27 days Any hint would be helpful. Regards Utkarsh [[alternative HTML

Re: [R] PGF Device

2009-05-06 Thread Gabor Grothendieck
There is pgfSweave on rforge: http://www.rforge.net/pgfSweave/ On Wed, May 6, 2009 at 10:43 AM, Lasse Bombien la...@phonetik.uni-muenchen.de wrote: Thanks, that's a good idea and I might use it for some purposes. But want I would like is a complete graphics device, which creates PGF output in

Re: [R] tapply changing order of factor levels?

2009-05-06 Thread Alain Guillet
Hi, I meant that your problem occured because the levels of mylevels are not ordered whereas tapply uses the ordered levels for printing. If you order them (look under), you can see the results of the tapply has the same order as the levels of myfactor

Re: [R] PGF Device

2009-05-06 Thread baptiste auguie
I think the pgfSweave project on R-forge is working on this (as far as i know it currently relies on eps2pgf) http://r-forge.r-project.org/R/?group_id=331 HTH, Baptiste On 6 May 2009, at 15:37, Lasse Bombien wrote: Hi all, I saw a thread from 2007 about the possible implementation of a

Re: [R] smoothing spline in package gam

2009-05-06 Thread willow1980
Sorry, I did not notice you were using GAM package. Most R users are using Simon Wood's MGCV package. I recommend you to use it. I have never used GAM package, so I cannot make further comments. Good luck! 楊 詩韻 wrote: dear all, i have a little question, but it make me torment long

[R] General Question about Graphics

2009-05-06 Thread Steve_Friedman
Greetings I have encountered a situation with regards to plotting barcharts with associated error bars. My search for clues on how to accomplish this turned up some interesting information. Basically, I found that including error bars with barplots is not desirable and hence there appears that

Re: [R] A question about using “by” in GAM model fitting of interaction between smooth terms and factor

2009-05-06 Thread willow1980
Dear Simon, Thank you so much! Actually, it seems that Crawley's R book adopted the information from that earlier version and he discussed by in the context of that version. I will take a practice according to your suggestion. Thanks again! Jianghua Simon Wood-4 wrote: The problem here is

[R] changing variables in a dataset

2009-05-06 Thread schung
Hi! I just quasi-learned how to use R and I'm trying to change the feed from horsebean, linseed, etc to 1, 2, 3, etc The dataset I'm practicing on is pre-loaded in R data(chickwts) chickwts Any help is much appreciated!! Thank you!! -- View this message in context:

[R] Need help building combined model

2009-05-06 Thread E Cell
I have two kind of variables. One kind where the values are continious and ranges from -1 to 5, and the other is boolean (i.e. true or false). For example: V1   V2  V3 V4 V5  V6 0.1  0.4   1    0   -1   3.7 0.4  0.1   2   -1   0    3.7 and V1  V2  V3 0 1    0 0     0    1 I want to

[R] Houston, TX Users Group

2009-05-06 Thread Robert Sanford
I'm looking for a Users Group in or near Houston, TX. Many thanks! rjsjr __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] General Question about Graphics

2009-05-06 Thread Dieter Menne
Steve_Friedman at nps.gov writes: Basically, I found that including error bars with barplots is not desirable and hence there appears that there is no function to do this. Can someone offer suggestions on how to do this simple procedure http://markmail.org/message/3f6jaztq2cpezrkl or

Re: [R] changing variables in a dataset

2009-05-06 Thread Federman, Douglas
you could use chickwts$feed- as.numeric(levels(chickwts$feed))[as.integer(chickwts$feed)] I got this by searching archives of old messages. I can't recall why the simpler as.numeric(chickwts$feed) is not an appropriate all around solution. -- I like nonsense, it wakes up the brain

Re: [R] General Question about Graphics

2009-05-06 Thread Mike Lawrence
Check out ggplot2: http://had.co.nz/ggplot2/ Particularly the book: http://had.co.nz/ggplot2/book For error bars on bar plots specifically: http://had.co.nz/ggplot2/geom_errorbar.html Mike On Wed, May 6, 2009 at 12:34 PM, steve_fried...@nps.gov wrote: Greetings I have encountered a

Re: [R] problem with ggplot2 boxplot, groups and facets

2009-05-06 Thread stephen sefick
The code is real, yes. But I can not copy and paste it right out of the email into an R session and see what you are talking about. Look at ?dput as a way of posting data to the list. Stephen Sefick On Wed, May 6, 2009 at 1:52 AM, Zeljko Vrba zv...@ifi.uio.no wrote: On Tue, May 05, 2009 at

Re: [R] Houston, TX Users Group

2009-05-06 Thread hadley wickham
Hi Robert, I'm organising one - sign up to the mailing list, http://groups.google.com/group/houston-r. I'm hoping to organise our first meeting this summer. Hadley On Wed, May 6, 2009 at 10:15 AM, Robert Sanford wob...@gmail.com wrote: I'm looking for a Users Group in or near Houston, TX.

Re: [R] Create Pie chart from .csv file

2009-05-06 Thread DonkeyRhubarb
Ive found out a way around my problem. I was trying to plaot a histogram of strings, but I had to change it into integers. I ran an sql query on the original DB that I got the CSV file from and used COUNT to get the number of each unique item in a given column. I then used these numbers to

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Greg Snow
In my opinion, no statisticians toolbox should contain only 1 tool (even if it is as amazing a tool as R). Learning the different tools helps you appreciate when each are the most appropriate to use and learn different ways of looking at problems. There are some tasks that I (it could easily

Re: [R] Create Pie chart from .csv file

2009-05-06 Thread Gabor Grothendieck
The sqldf package allows you to apply SQL statements to R data frames: http://sqldf.googlecode.com On Wed, May 6, 2009 at 12:08 PM, DonkeyRhubarb michaeli...@gmail.com wrote: Ive found out a way around my problem. I was trying to plaot a histogram of strings, but I had to change it into

[R] [Fwd: loading SPOT file]

2009-05-06 Thread Barbara Cegielska
---BeginMessage--- I have problem with loading SPOT files. Sample of my operations in R object fromfile=TRUE fext=spot rawboxplot=TRUE flagi=Flags f-function(x) as.numeric(x$Flags -99) targets-readTargets(ApoAITargets.txt) if (fromfile){ + files-dir(pattern=fext) + } else { +

Re: [R] [Fwd: loading SPOT file]

2009-05-06 Thread Martin Morgan
Hi Barbara -- ask on the Bioconductor mailing list. http://bioconductor.org/docs/mailList.html Include the output of the command sessionInfo() after the failure. Be sure to update to a current version of R / Bioconductor, for the latter following 'Update installed Bioconductor packages' at

Re: [R] Summary help

2009-05-06 Thread mathallan
To glm is glm(log(mydata)~log(max_data)*as.factor(grp),family=Gamma(link=log)) And I was wondering if you can read the scale and shape from summary There a quite a few gamma models around, so you should tell us more. glmXXX? lmer? Dieter __

[R] 'RG' looks like a pre-2.4.0 S4 object: please recreate it

2009-05-06 Thread Barbara Cegielska
I would like to load ApoAI.RData. During the operation of reading this data an error occurs. There is also a problem with STF file. library (limma) load(ApoAI.RData) Warning message: 'RG' looks like a pre-2.4.0 S4 object: please recreate it objects() [1] RG names(RG) [1] R G Rb

[R] Positioning a legend via X and Y coordinates

2009-05-06 Thread Steve Murray
Dear R Users, I'm able to display a legend using the following code: legend(topright, c(Simulation, Observation), fill=2:3, bty=n) However, this causes the legend to be positioned too close to the bars in my barplot. I'd like to move the legend up slightly. I have been trying to determine

  1   2   >