Re: [R] Help with apply

2011-05-13 Thread Berend Hasselman
M.Ribeiro wrote: Ok, I am sorry, I was building a reproducible example to post here and I made a mistake. Th second loop (in the variable k) updates ycorr. Therefore, it would be ycorr - rnorm(10) Nindiv - 10 x -matrix(sample(c(0:2),50,re = T),10) gtemp - rnorm(10) for (i in 1:Nindiv)

Re: [R] Uploading CSV file into R

2011-05-13 Thread peter dalgaard
On May 13, 2011, at 05:57 , Me wrote: Sorry if I don't know all the vocabulary. I'm taking a research methods class that uses R, but I'm a journalism major, so I'm not very good with a lot of computer jargon. The read.csv operator didn't work. I got the following error message: Error in

Re: [R] selecting data from table with timestamp

2011-05-13 Thread Pablo Rosado
Thank you Don so much, It was very helpful and it worked from the first try. Thanks for the extra comments on my script. On Wed, May 11, 2011 at 4:41 PM, MacQueen, Don macque...@llnl.gov wrote: Try something similar to this: ## unchanged full - read.table(March_15.dat,

Re: [R] Problem with constrained optimization with maxBFGS

2011-05-13 Thread peter dalgaard
On May 13, 2011, at 00:15 , Prof. John C Nash wrote: Is this a homework problem in finding the largest eigensolution of W? Not as stated. It said sum D = 1, not sum D^2 =1. With a linear constraint, it looks like weighted least squares, but then you should be minimizing, not maximizing. As

Re: [R] Embedding R's math expressions fonts for graphics generated by pdf() and postscript()

2011-05-13 Thread Prof Brian Ripley
Please study the help page for pdf (as the posting guide asked you to do before posting). This is a documented problem with fontconfig/cairo-based viewers, which evince is. On Thu, 12 May 2011, Vinh Nguyen wrote: Dear list, First, I am not writing to ask about embedding Computer Modern

Re: [R] Fw: Help with PLSR

2011-05-13 Thread Bjørn-Helge Mevik
Amit Patel amitrh...@yahoo.co.uk writes: str(FullDataListTrans) num [1:40, 1:94727] 42 40.9 65 56 61.7 ... - attr(*, dimnames)=List of 2 ..$ : chr [1:40] X X.1 X.12 X.13 ... ..$ : NULL I have also created a vector GroupingList which gives the groupnames for each respective

[R] routine for dependent correlation test with stratified random sample

2011-05-13 Thread D. Alain
Dear R-List,   I would like to have a large number of stratified random subsamples drawn from my dataframe and automatically test for correlation differences in every subsample.   Let this be my dataframe  

Re: [R] Uploading CSV file into R

2011-05-13 Thread Jannis
Did you download this data to the location given in the command? 'Downloads/UNdata_Export_20110511_130705015.csv' The way I read this eror message is that the file either does not exist or the path is speelled incorrectly! You probably have to set the working directory ( setwd() ) to one

Re: [R] R won't start keeps crashing

2011-05-13 Thread Ivan Calandra
I'm no expert, but checking the R-help list, I got the impression that saving a workspace often leads to problems. The interesting question (for me) would be to know why and how the corrupted file got corrupted, and then, what to do to avoid the problem if one really wants to save the

[R] wrong argument to TO_CHAR (RODBC)?

2011-05-13 Thread Nuno Prista
Dear all, I can't seem to run what seems (to me!) like a common SQL statement in RODBC: data - sqlQuery(a, 'SELECT TO_CHAR (DB1.DATE, '') year FROM DB1') In contrast, data - sqlQuery(a, 'SELECT TO_CHAR(DB1.DT_INICIO) year FROM DB1') works but it does not come up as I want it to

[R] to check if a group of values is randomly distributed

2011-05-13 Thread elaine kuo
Dear list, I have 603 numbers depicting range sizes of birds in Japan. I would like to learn if the 603 range sizes are randomly distributed or not, in order to check if they meet mid-domain effects. Please kindly advise if any R package or function can check the fit. Also, any more references

[R] Powerful PC to run R

2011-05-13 Thread Michael Haenlein
Dear all, I'm currently running R on my laptop -- a Lenovo Thinkpad X201 (Intel Core i7 CPU, M620, 2.67 Ghz, 8 GB RAM). The problem is that some of my calculations run for several days sometimes even weeks (mainly simulations over a large parameter space). Depending on the external conditions, my

Re: [R] to check if a group of values is randomly distributed

2011-05-13 Thread Steven Kennedy
try the fitdistr() function in the MASS package __ 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 commented, minimal,

Re: [R] Power Spectrum from STFT (e1071)?

2011-05-13 Thread Mike Marchywka
Date: Tue, 10 May 2011 11:55:34 -0700 From: ajf...@psu.edu To: r-help@r-project.org Subject: [R] Power Spectrum from STFT (e1071)? Hello. Does anyone know how to generate a power spectrum from the STFT function in package e1071? The output for this function supplies the Fourier

Re: [R] Powerful PC to run R

2011-05-13 Thread Mike Marchywka
Date: Fri, 13 May 2011 12:38:51 +0200 From: haenl...@escpeurope.eu To: r-help@r-project.org Subject: [R] Powerful PC to run R Dear all, I'm currently running R on my laptop -- a Lenovo Thinkpad X201 (Intel Core i7 CPU, M620, 2.67 Ghz, 8 GB RAM). The problem is that some of my

[R] Excel to R

2011-05-13 Thread Asan Ramzan
Hello R-help   I am trying to copy and paste large column of data from windows Excel into R base speed sheet using:   Data1-edit(data.frame())   But I can only copy and paste one cell at a time which is time consuming.   Also how can assign each column to a variable. [[alternative HTML

Re: [R] Powerful PC to run R

2011-05-13 Thread Hrishi Mittal
Have you considered using Amazon EC2? If you search on http://www.r-bloggers.com, people have written about their experiences. On Fri, May 13, 2011 at 11:38 AM, Michael Haenlein haenl...@escpeurope.euwrote: Dear all, I'm currently running R on my laptop -- a Lenovo Thinkpad X201 (Intel Core

Re: [R] Excel to R

2011-05-13 Thread Gabor Grothendieck
On Fri, May 13, 2011 at 7:02 AM, Asan Ramzan asanram...@yahoo.com wrote: Hello R-help I am trying to copy and paste large column of data from windows Excel into R base speed sheet using: Data1-edit(data.frame()) But I can only copy and paste one cell at a time which is time consuming.

Re: [R] Powerful PC to run R

2011-05-13 Thread Barry Rowlingson
On Fri, May 13, 2011 at 11:38 AM, Michael Haenlein haenl...@escpeurope.eu wrote: I'm now thinking about buying a more powerful desktop PC or laptop. Can anybody advise me on the best configuration to run R as fast as possible? I will use this PC exclusively for R so any other factors are of

[R] problem with makeSOCKcluster depending on R patch version

2011-05-13 Thread Ulrich Halekoh
Dear, I encountered a problem using the makeSOCKcluster function depending the patched version of R-2.13.0 I used. library(snow) cl - makeSOCKcluster(rep(localhost, 2)) this works fine for the R-13.0 patch (2011-04-28 r55678) but not for the patch R-13.0 patch (2011-05-10 r55826) In the

Re: [R] Can ROC be used as a metric for optimal model selection for randomForest?

2011-05-13 Thread Frank Harrell
Using anything other than deviance (or likelihood) as the objective function will result in a suboptimal model. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context:

Re: [R] Excel to R

2011-05-13 Thread Robert Baer
-- From: Asan Ramzan asanram...@yahoo.com Sent: Friday, May 13, 2011 6:02 AM To: r-help@r-project.org Subject: [R] Excel to R Hello R-help I am trying to copy and paste large column of data from windows Excel into R base speed sheet using:

Re: [R] Can ROC be used as a metric for optimal model selection for randomForest?

2011-05-13 Thread Max Kuhn
XiaoLiu, I can't see the options in bootControl you used here. Your error is consistent with leaving classProbs and summaryFunction unspecified. Please double check that you set them with classProbs = TRUE and summaryFunction = twoClassSummary before you ran. Max On Thu, May 12, 2011 at 7:04

Re: [R] Exporting interactive 3D plots with axes and labels

2011-05-13 Thread van Rijn, Peter
Snapsnot3d just creates a flat image. It is possible to export the actual interactive 3D plot and put in a PDF, just not with the axes and labels. Thanks, Peter -Original Message- From: Remko Duursma [mailto:remkoduur...@gmail.com] Sent: Thursday, May 12, 2011 7:17 PM To: van Rijn,

[R] Survival Rate Estimates

2011-05-13 Thread Terry Therneau
---begin included message -- Is there an automated way to use the survival package to generate survival rate estimates and their standard errors? To be clear, *not *the survivorship estimates (which are cumulative), but the survival *rate * estimates... --- end -- The classic method in

Re: [R] Can ROC be used as a metric for optimal model selection for randomForest?

2011-05-13 Thread Max Kuhn
Frank, It depends on how you define optimal. While I'm not a big fan of using the area under the ROC to characterize performance, there are a lot of times when likelihood measures are clearly sub-optimal in performance. Using resampled accuracy (or Kappa) instead of deviance (out-of-bag or not)

[R] Quick question: Omitting rows and cols with certain percents of missing values

2011-05-13 Thread Vickie S
Hi naive question. It is possible to get R command for omitting rows or cols with missing values present. But if i want to omit rows or cols with i.e . 20% missing values, I could´t find any package-based command, probably because it is too simple for anyone to do that manually, though not

Re: [R] Quick question: Omitting rows and cols with certain percents of missing values

2011-05-13 Thread David Winsemius
On May 13, 2011, at 9:42 AM, Vickie S wrote: Hi naive question. It is possible to get R command for omitting rows or cols with missing values present. But if i want to omit rows or cols with i.e . 20% missing values, I could´t find any package-based command, probably because it is too

Re: [R] Quick question: Omitting rows and cols with certain percents of missing values

2011-05-13 Thread John Kane
--- On Fri, 5/13/11, Vickie S is...@live.com wrote: From: Vickie S is...@live.com Subject: [R] Quick question: Omitting rows and cols with certain percents of missing values To: r-help@r-project.org Received: Friday, May 13, 2011, 9:42 AM Hi naive question. It is possible to get R

Re: [R] R won't start keeps crashing

2011-05-13 Thread Duncan Murdoch
On 11-05-13 6:03 AM, Ivan Calandra wrote: I'm no expert, but checking the R-help list, I got the impression that saving a workspace often leads to problems. The interesting question (for me) would be to know why and how the corrupted file got corrupted, and then, what to do to avoid the problem

Re: [R] Quick question: Omitting rows and cols with certain percents of missing values

2011-05-13 Thread Jorge Ivan Velez
Hi Vickie, You might try the following: # some data set.seed(123) X - matrix(rnorm(1000), ncol = 20) X[sample(1000, 100)] - NA # excluding rows with NA 20% X[!rowMeans(is.na(X)) 0.2, ] # excluding columns with NA 10% X[, !colMeans(is.na(X)) 0.1] See ?is.na, ?rowMeans and ?colMeans for more

Re: [R] Quick question: Omitting rows and cols with certain percents of missing values

2011-05-13 Thread Peter Ehlers
On 2011-05-13 06:42, Vickie S wrote: Hi naive question. It is possible to get R command for omitting rows or cols with missing values present. But if i want to omit rows or cols with i.e .20% missing values, I could´t find any package-based command, probably because it is too simple for

[R] How to store a triangular matrix

2011-05-13 Thread Estefania Ruiz Vargas
Hello, I want to create a triangular matrix and only keep the lower triangle entries without having to allocate memory for the whole matrix, is there any way I can do something like  A-matrix(data, nrow=50)  but for a triangular matrix? Thanks ERV [[alternative HTML version

Re: [R] Saving misclassified records into dataframe within a loop

2011-05-13 Thread John Dennison
I appreciate all the help and patience. My pea sized brain had yet to grasp the concept of a function return. results-function() absolutely works. i wanted it in a dataframe to most easily transfer it to a database. Ive found the Rpostgres like data frames better. Im sure a character vector would

Re: [R] How to store a triangular matrix

2011-05-13 Thread David Winsemius
On May 13, 2011, at 10:52 AM, Estefania Ruiz Vargas wrote: Hello, I want to create a triangular matrix and only keep the lower triangle entries without having to allocate memory for the whole matrix, is there any way I can do something like A-matrix(data, nrow=50) but for a triangular

[R] biplots for PCA

2011-05-13 Thread Anna Renwick
Hi all I have produced a biplot for a PCA (see attached pdf) that I ran however the names of the variables which are placed at the end of the arrows overlap and are thus unreadable. Similarly some of the numbered points overlap. I was wondering if there was a way to edit the biplot to move the

[R] Error Subscript out of bounds in ca

2011-05-13 Thread annekatrin . st
Hello, I want to do a simple correspondence analysis in R, using the ca package. First, I create and define a 2-by-2-table (table=table(var1,var2)), which looks like this: Low(var1) high(var1) Low(var2) 35 28 High(var2) 26 32 Var1 and Var2 are categorical factors.

[R] 回复: Uploading CSV file into R

2011-05-13 Thread shizhao cao
when uploading csv file to R, the file.csv should save as.csv from. xls,if u just modify suffix of files, files cannot be read . BTW,who know the problem i meet, strsplit(txt$Sentence,,|。|!|?|、) how to understand this order~ txt=read.csv(1988.csv,colClasses=character)

Re: [R] problem converting character to dates

2011-05-13 Thread Assu
Thanks Uwe, it works! I just have a problem at week 53 which gives me an NA. I guess, this week 53 has to do with the data type set at the database design and collection of data. Depends on the week/year kind of counting. -- View this message in context:

[R] More effective calculation for loop

2011-05-13 Thread Wolfgang Wu
I am getting confused in how to use R effectively. Below is the example that I would like to optimize, which is basically a loop. In its current version it takes about 20 seconds to run. I have looked at the apply function and but didn't seem to get it to work. I am sure this is relatively

[R] PLSR error

2011-05-13 Thread romzero
Hi, this is my R-Script library(pls) file - C:\\TXT\\brix.txt d - as.matrix(read.table(file, header=T, sep=,, row.names = NULL)) plsdata - data.frame(NIR=c(1:nrow(X))) plsdata$NIR - I(d[,3:603]) plsdata$Brix - d[,2] results - plsr(Brix ~ NIR, data=plsdata) after the last string i have this error

[R] issue with odfWeave running on Windows XP; question about installing packages under Linux

2011-05-13 Thread Cormac Long
Good morning R community, I have two questions (and a comment): 1) A problem with odfWeave. I have an odf document with a table that spans multiple pages. Each cell in the table is populated using \sexpr{R stuff}. This worked fine on my own machine (windows 7 box using any R2.x.y, for x=11) and

Re: [R] Uploading CSV file into R

2011-05-13 Thread skfgla...@gmail.com
Usually this kind of error indicates one or two errors. First check that the path to the input file is correct. Make sure that your input file ends with csv. Next problem could be that a spelling error is causing the problem. Make sure everything is spelled correctly in both the path

Re: [R] lm and anova

2011-05-13 Thread Sara Sjöstedt de Luna
Thanks a lot! Best, Sara -Ursprungligt meddelande- Från: istaz...@gmail.com [mailto:istaz...@gmail.com] För Ista Zahn Skickat: den 12 maj 2011 16:46 Till: Sara Sjöstedt de Luna Kopia: r-help@r-project.org; Katie Andrle (katie.and...@gmail.com) (katie.and...@gmail.com) Ämne: Re: [R] lm

[R] file size plot pdf

2011-05-13 Thread Marc Jekel
Dear R-Fans, I have been lately working on some plots in R that I save as pdf via the pdf() command. I have realized that when I open those files in Adobe and then re-save it within Adobe (save as...), the size of the pdf files decreases rapidly (e.g., from 4mb to 1mb). This can also be

Re: [R] Powerful PC to run R

2011-05-13 Thread John C Frain
Given the price of the Thinkpad X201 I would be very upset if it overheated and would contact my supplier for a replacement. I also understand that this CPU may slow down if it is overheated. You would probably get a better performance from a desktop workstation or a larger laptop with better

Re: [R] Instances of a C++ class in R

2011-05-13 Thread soeren . vogel
On 28.04.2011, at 12:18, soeren.vo...@uzh.ch wrote: On 27.04.2011, at 11:59, soeren.vo...@uzh.ch wrote: We are working on a class in C++. The files compile fine (R CMD SHLIB ...) and run in R. A bzipped tar archive with source code can be downloaded from here:

[R] help with probabilites in cph

2011-05-13 Thread viostorm
I would like to get the probabilities from cph into a vector so I can put them programatically on a graph. Does anyone know how to get a list of the probabilities? All I can seem to do is $coef which only contains coefficients. cph(formula = srv.comb ~ collaterals +

[R] RMySQL problem on Windows 7 32 bit

2011-05-13 Thread byzin
I am attempting to analyze Affymetrix exonchip data using the exon map package. I have installed this package but have had problems installing the associated programs, specifically RMySQL since no windows binary exists. I have followed the directions for compiling the package as found on -

Re: [R] biplots for PCA

2011-05-13 Thread Scott Chamberlain
Have you looked into the options for par?, which it looks like can be included in your call to biplot, perhaps adj or something else to reduce overlap. Scott On Friday, May 13, 2011 at 4:38 AM, Anna Renwick wrote: Hi all I have produced a biplot for a PCA (see attached pdf) that I ran

Re: [R] biplots for PCA

2011-05-13 Thread Kevin Wright
Did you try tweaking the value of the expand parameter? See the help page of biplot for more info. Kevin On Fri, May 13, 2011 at 4:38 AM, Anna Renwick anna.renw...@bto.org wrote: Hi all I have produced a biplot for a PCA (see attached pdf) that I ran however the names of the variables

Re: [R] More effective calculation for loop

2011-05-13 Thread John Kane
What is xts? What package? --- On Fri, 5/13/11, Wolfgang Wu wobw...@yahoo.de wrote: From: Wolfgang Wu wobw...@yahoo.de Subject: [R] More effective calculation for loop To: R-help@r-project.org Received: Friday, May 13, 2011, 5:32 AM I am getting confused in how to use R effectively. Below

Re: [R] assigning creating missing rows and values

2011-05-13 Thread Adele_Thompson
The problem with using cumsum, is that the measured output is the cumulative feed consumed throughout the day. When the animals do not eat for 30 minutes or so, it will not output a new value, but as soon as they do eat, the scale will weigh the difference and then output the cumulative feed

Re: [R] More effective calculation for loop

2011-05-13 Thread Joshua Ulrich
On Fri, May 13, 2011 at 4:32 AM, Wolfgang Wu wobw...@yahoo.de wrote: I am getting confused in how to use R effectively. Below is the example that I would like to optimize, which is basically a loop. In its current version it takes  about 20 seconds to run. I have looked at the apply function

Re: [R] help with probabilites in cph

2011-05-13 Thread David Winsemius
On May 13, 2011, at 10:49 AM, viostorm wrote: I would like to get the probabilities of what? from cph into a vector so I can put them programatically on a graph. Does anyone know how to get a list of the probabilities? require(rms) ?rms Perhaps: ?Predict Or ?survest (You do

[R] Creating dataframes

2011-05-13 Thread Woida71
I would like to create a certain number of dataframes out of one dataframe where each of the dataframes is related to a factor. This should be possible with a loop or a function, as is very clumsy to do it manually when there are quite a lot factors. For example having a dataframe called exmpl:

[R] effects package for adjusted predictions

2011-05-13 Thread Xu Jun
Dear R experts, I am trying to use John Fox's effects package to conduct some post-estimation analysis.  I think it is a great package for graphing predictions after GLM type models.  However, I am having problems figuring out how to produce an adjusted mean.  For example, using the following two

Re: [R] Creating dataframes

2011-05-13 Thread Phil Spector
The way to do what you want is to use the split function, for example, mydataframes = split(exmpl,exmpl$Site) This will return a list with all the data frames. If you don't understand why this is a better solution than creating many separate data frames, let us know what you plan to do with

Re: [R] Creating dataframes

2011-05-13 Thread Marc Schwartz
On May 13, 2011, at 10:28 AM, Woida71 wrote: I would like to create a certain number of dataframes out of one dataframe where each of the dataframes is related to a factor. This should be possible with a loop or a function, as is very clumsy to do it manually when there are quite a lot

[R] RES: How to store a triangular matrix

2011-05-13 Thread Filipe Leme Botelho
Estefania, if you do something like your.matrix - matrix(1, 50, 50) your.matrix[upper.tri(your.matrix)] - NA You will have NAs for the entire upper triangle. If you opt to do your.matrix[upper.tri(your.matrix)] - The upper triangle will be empty but it will coerce

[R] RES: Excel to R

2011-05-13 Thread Filipe Leme Botelho
---BeginMessage--- Why not using Statconn/RExcel? -Mensagem original- De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Em nome de Robert Baer Enviada em: sexta-feira, 13 de maio de 2011 09:33 Para: Asan Ramzan; r-help@r-project.org Assunto: Re: [R] Excel to R

[R] period in scale parameter of prcomp

2011-05-13 Thread Greg Hirson
Hi all, Looking at the help file for prcomp and the code in stats:::prcomp.default, the scale parameter for the default S3 method has a trailing period: it written as scale. I assume there is a reason for it (backwards compatibility with S was suggested) but I thought I'd ask. Thanks, Greg --

Re: [R] period in scale parameter of prcomp

2011-05-13 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Greg Hirson Sent: Friday, May 13, 2011 10:58 AM To: R-help@r-project.org Subject: [R] period in scale parameter of prcomp Hi all, Looking at the help file for prcomp and

[R] Installation of rattle complains about missing packages

2011-05-13 Thread Juergen Rose
Hi, if I try to install rattle, R fails to install: rpvm, rsprng, RBGL, pcalg, arulesViz, pkgDepTools, Rgraphviz . rpvm depends on pvm, rsprng on sprng which is not installed and which seems not so easy to install on my gentoo systems. ArulesViz depends on Rgraphviz, RBGL and pcalg depends on

Re: [R] period in scale parameter of prcomp

2011-05-13 Thread Greg Hirson
Bill, Thanks. So it is as we suspected. Kind regards, Greg On 5/13/11 11:19 AM, William Dunlap wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Greg Hirson Sent: Friday, May 13, 2011 10:58 AM To: R-help@r-project.org

Re: [R] file size plot pdf

2011-05-13 Thread Prof Brian Ripley
On Fri, 13 May 2011, Marc Jekel wrote: Dear R-Fans, I have been lately working on some plots in R that I save as pdf via the pdf() command. I have realized that when I open those files in Adobe and then re-save it within Adobe (save as...), the size of the pdf files decreases rapidly (e.g.,

[R] Plots: I've deleted axes, now to delete space

2011-05-13 Thread Schatzi
I am plotting 28 plots on one screen: par(mfrow=c(4,7)) for (i in 1:28) { a-seq(1,3,1) plot(a,a, ann=FALSE) } I want a main title for all the plots (I tried using main but that doesn't work). I deleted the axes, but am not sure how to delete the space. There are such large margins between plots

Re: [R] Plots: I've deleted axes, now to delete space

2011-05-13 Thread Greg Snow
Look at the help for par, specifically the section on 'mar' to set the per plot margins smaller and the section on 'oma' to leave room for the overall title. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original

Re: [R] Plots: I've deleted axes, now to delete space

2011-05-13 Thread Adele_Thompson
Easy fix. Under ?par, I don't see where I can enter an overall title. Should I add a text command or something? -Original Message- From: greg.s...@imail.org [mailto:greg.s...@imail.org] Sent: Friday, May 13, 2011 03:17 PM To: Thompson, Adele - adele_thomp...@cargill.com;

Re: [R] Plots: I've deleted axes, now to delete space

2011-05-13 Thread Duncan Murdoch
On 11-05-13 4:21 PM, adele_thomp...@cargill.com wrote: Easy fix. Under ?par, I don't see where I can enter an overall title. Should I add a text command or something? mtext() writes text in the margins; argument outer puts it in the outer margins. Duncan Murdoch -Original

Re: [R] Plots: I've deleted axes, now to delete space

2011-05-13 Thread Greg Snow
Use par with oma to leave room for the overall title, then you need another function to create the title. The mtext function is probably the easiest. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message-

Re: [R] Plots: I've deleted axes, now to delete space

2011-05-13 Thread Adele_Thompson
Thanks all for the replies. I am getting better slowly but surely. I imagine that I will get better at figuring out things as well so I don't have to post as many questions. I do lots of searches, but still cannot figure out how to do everything that I need. The new code is as such:

Re: [R] assigning creating missing rows and values

2011-05-13 Thread Bert Gunter
Adele: You are of course correct -- my earlier proposed solution is dumb (thank you for being polite and not saying this :-) ) With your more complete explanation -- which might have helped you get a quicker reply earlier from others had you given it at the beginning -- I was able to come up

[R] Coloring Continents/Regions of World Map

2011-05-13 Thread ang
Dear All, This may be a silly question, but I have tried searching through a few of the existing threads already, and haven't found any information on this. I am trying to color code entire continents/regions (Asia, North America, Africa, etc.) based on the number of relevant companies. The

[R] Tps error

2011-05-13 Thread Jennifer Gremer
Hi, I am trying to use the Tps procedure to visualize 3 variables in XYZ space, but keep getting the error messages below: Warning messages: 1: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, : GCV search gives a minimum at the endpoints of the grid search 2: In

Re: [R] assigning creating missing rows and values

2011-05-13 Thread Adele_Thompson
This code works great. I will improve my descriptions of what I want in the future. Thanks for the help. -Original Message- From: gunter.ber...@gene.com [mailto:gunter.ber...@gene.com] Sent: Friday, May 13, 2011 03:44 PM To: Thompson, Adele - adele_thomp...@cargill.com Cc:

Re: [R] RMySQL problem on Windows 7 32 bit

2011-05-13 Thread Steven Kennedy
It looks like your MySQL is missing the libmysql.lib file. Download it from MySQL then put it in the following location: C:/Program Files/MySQL/MySQL Server 5.5/lib/opt/libmysql.lib If you are still having problems, you could try using RODBC instead.

Re: [R] Can ROC be used as a metric for optimal model selection for randomForest?

2011-05-13 Thread Frank Harrell
Thanks for your note Max. Part of the picture is how predictions would be used. If they are used in a forced choice way (quite a shame because the best decision is often no decision - get more data) things are different. If there are gray zones or predicted probabilities are of interest then

[R] Russian language in R

2011-05-13 Thread lyolya
Hello, I am experiencing a problem in reading a database in Russian. The problem appears when it comes to char variables. I have already tried changing the encoding, i.e. options(encoding=UTF-8) and options(encoding=KOI8-R) but every time there appear to be something unreadable in the data

[R] L'abbe plot

2011-05-13 Thread Whitney . Melroy
I cannot seem to get a L'abbe plot to work on R. I do not understand what the X coordinates, or alternatively an object of class metabin, is supposed to mean. What is a class of metabin? Institute of Behavioral Genetics University of Colorado, Boulder whitney.mel...@colorado.edu

[R] graphs of gamma, normal fit to a histogram are about half as large as they should be

2011-05-13 Thread Benjamin Caldwell
Hello, I'm trying to compare the fit of two distributions, normal and gamma, to a histogram of my response variable. rate-mean(na.omit(rwb$post.f.crwn.length))/var(na.omit(rwb$post.f.crwn.length)) shape-rate*mean(na.omit(rwb$post.f.crwn.length)) hist((rwb$post.f.crwn.length),

[R] using glmer to fit a mixed-effects model with gamma-distributed response variable

2011-05-13 Thread Benjamin Caldwell
Sub: using glmer to fit a mixed-effects model with gamma-distributed response variable Hello, I'm currently trying to fit a mixed effects model , i.e.: burnedmodel1.2-glmer(gpost.f.crwn.length~lg.shigo.av+dbh+leaf.area+ bark.thick.bh+ht.any+ht.alive+(1|site/transect/plot), family=gaussian,

[R] How do I break my addiction to for loops!?!?

2011-05-13 Thread Justin Haynes
I know I'm not supposed to use them... but they're just so easy! I have trouble defining an appropriate function for plyr or apply! data-rnorm(144) groups1-c('a','b','c','d') groups2-c('aa','bb','cc','dd') machines-1:12

Re: [R] Russian language in R

2011-05-13 Thread Duncan Murdoch
On 13/05/2011 4:57 PM, lyolya wrote: Hello, I am experiencing a problem in reading a database in Russian. The problem appears when it comes to char variables. I have already tried changing the encoding, i.e. options(encoding=UTF-8) and options(encoding=KOI8-R) but every time there appear to

Re: [R] graphs of gamma, normal fit to a histogram are about half as large as they should be

2011-05-13 Thread Rolf Turner
On 14/05/11 10:00, Benjamin Caldwell wrote: Hello, I'm trying to compare the fit of two distributions, normal and gamma, to a histogram of my response variable. rate-mean(na.omit(rwb$post.f.crwn.length))/var(na.omit(rwb$post.f.crwn.length)) shape-rate*mean(na.omit(rwb$post.f.crwn.length))

Re: [R] R won't start keeps crashing

2011-05-13 Thread Rolf Turner
On 14/05/11 02:28, Duncan Murdoch wrote: SNIP One way to end up with a corrupted workspace is to save it with a dependency on a particular package, then try to load it when that package is unavailable. I'm going to see if R can handle that case more gracefully. SNIP Pardon me for

[R] Changing x-axis when dealing with time

2011-05-13 Thread Pablo Rosado
Hi, I am plotting data in which the x values are a timestamp. I am trying to change the x-ticks so they will be in specified hours of the day, but it always start from hour 4 of the day. And I need it to start from the beginning of the axis (at x=0) and then each tick be on the interval I

Re: [R] to check if a group of values is randomly distributed

2011-05-13 Thread Ben Bolker
elaine kuo elaine.kuo.tw at gmail.com writes: Dear list, I have 603 numbers depicting range sizes of birds in Japan. I would like to learn if the 603 range sizes are randomly distributed or not, in order to check if they meet mid-domain effects. Please kindly advise if any R package or

Re: [R] using glmer to fit a mixed-effects model with gamma-distributed response variable

2011-05-13 Thread Ben Bolker
Benjamin Caldwell btcaldwell at berkeley.edu writes: Sub: using glmer to fit a mixed-effects model with gamma-distributed response variable Hello, I'm currently trying to fit a mixed effects model , i.e.: burnedmodel1.2-glmer(gpost.f.crwn.length~lg.shigo.av+dbh+leaf.area+

Re: [R] R won't start keeps crashing

2011-05-13 Thread David Winsemius
On May 13, 2011, at 9:09 PM, Rolf Turner wrote: On 14/05/11 02:28, Duncan Murdoch wrote: SNIP One way to end up with a corrupted workspace is to save it with a dependency on a particular package, then try to load it when that package is unavailable. I'm going to see if R can handle that

[R] Adding same items together in data.frame

2011-05-13 Thread wong, honkit (Stephen)
Dear All, I am new to R. I have a 2 column data frame with more than ten thousand rows. Something like below. I want to add up all duplicated items, e.g. the three aa add up together to get a single value gene=a, value=74. How can I do that?? Thanks for help ! gene value aa 20 bb 10 cc

Re: [R] R won't start keeps crashing

2011-05-13 Thread Rolf Turner
On 14/05/11 14:09, David Winsemius wrote: On May 13, 2011, at 9:09 PM, Rolf Turner wrote: On 14/05/11 02:28, Duncan Murdoch wrote: SNIP One way to end up with a corrupted workspace is to save it with a dependency on a particular package, then try to load it when that package is

Re: [R] Adding same items together in data.frame

2011-05-13 Thread Dennis Murphy
Hi: There are many ways to do this sort of thing in R; one way is (naming your example data frame d) aggregate(value ~ gene, data = d, FUN = sum) gene value 1 aa74 2 bb10 3 cc 9 4 dd 100 5 ee55 This code line works for R-2.11.0 and later. HTH, Dennis On Fri, May

Re: [R] Adding same items together in data.frame

2011-05-13 Thread Steve Lianoglou
Hi, On Fri, May 13, 2011 at 7:06 PM, wong, honkit (Stephen) hon...@stanford.edu wrote: Dear All, I am new to R. I have a 2 column data frame with more than ten thousand rows. Something like below. I want to add up all duplicated items, e.g. the three aa add up together to get a single value