Re: [R] Urgent help requested to modify a script

2009-10-17 Thread Jim Lemon
On 10/16/2009 11:01 PM, one2luv wrote: I am hoping someone will tak up this chalenge (I am new to R) I have inheritied an R script but need to change it. The script currently includes hardcoded file locations on lines 12,166 and 167. I need to modify this script to allow the folder to be passed

Re: [R] Generating a stochastic matrix with a specified second dominant eigenvalue

2009-10-17 Thread Jeff Newmiller
Ravi Varadhan wrote: Hi, Given a positive integer N, and a real number \lambda such that 0 \lambda 1, I would like to generate an N by N stochastic matrix (a matrix with all the rows summing to 1), such that it has the second largest eigenvalue equal to \lambda (Note: the dominant eigenvalue

Re: [R] Modify the plot from countCDFxt

2009-10-17 Thread Peter Ehlers
Michelle DePrenger-Levin wrote: Hello, I am running the plot from countCDFxt (popbio). I would like to report the y-axis as a percent instead of the log scale (e^01...). I can add an axis with axis(2, 0:1, line =2) but I'm having trouble understanding how to assign the tic marks (with 'at =').

Re: [R] how to have tkchooseDirectory resize in windows?

2009-10-17 Thread Peter Dalgaard
David Gattrell wrote: R-2.8.0 / tcltk8.5 In windows, Rgui.exe has a directory browser that can be resized, but when I call tkchooseDirectory(), it is a fixed size. In linux, when I call tkchooseDirectory() it can be resized. How do I get a windows version that I can resize? You're at the

[R] sum two columns with same value

2009-10-17 Thread Alfredo Alessandrini
Hi, I've two dataframe: snag_totale AREA snag_ha 12 1.628128 23 10.274249 34 2.778503 45 73.764307 57 12.015985 log_totale AREAlog_ha 11 22.29846 22 17.16889 33 48.80377 44 144.18996 55 70.30962 66 61.81850 77 13.24876 How can

Re: [R] attach

2009-10-17 Thread Peter Dalgaard
Dieter Menne wrote: Christophe Dutang1 wrote: I would like to know what happens on the memory side when I use attach(inputdata) Is there a second allocation of memory for inputdata? Not, it just guides the syntax. Wrong. There's a virtual copy of data plus a conversion from data frame

[R] the number of components of the mixture model and the stratfied data

2009-10-17 Thread lybaomc
Hi,all with my data,there are more than 1000 quantitative results of antibody concentrations, there may be 2 components(positive and negative), or 3 components (may be strong positive, positive, and negative), or 4-6 components. Could you tell me how to determine the number of components of

Re: [R] sum two columns with same value

2009-10-17 Thread smu
On Sat, Oct 17, 2009 at 09:36:50AM +0200, Alfredo Alessandrini wrote: Hi, I've two dataframe: snag_totale AREA snag_ha 12 1.628128 23 10.274249 34 2.778503 45 73.764307 57 12.015985 log_totale AREAlog_ha 11 22.29846 22 17.16889 33

Re: [R] ozone data

2009-10-17 Thread Liviu Andronic
On 10/17/09, Julius Tesoro jutes...@yahoo.com wrote: I am looking at the plyr package and I am intrigued at how data(i.e. ozone, baseball) is loaded without having to type data(ozone). Are they automatically loaded when i call library(plyr)? I want to do the same thing when I make my

Re: [R] ozone data

2009-10-17 Thread Liviu Andronic
On 10/17/09, Julius Tesoro jutes...@yahoo.com wrote: wow thanks but how do i load data just by typing ozone on the console? It would probably suffice to include the dataset as an object in your package, and then `data(ozone)' should bring it up as soon as your package is loaded. Look at the

Re: [R] avoiding loops in equation

2009-10-17 Thread Julius Tesoro
thank you very much --- On Sat, 10/17/09, Kenn Konstabel lebats...@gmail.com wrote: From: Kenn Konstabel lebats...@gmail.com Subject: Re: [R] avoiding loops in equation Thank God for R-help mailing list. Thanks.. To: Julius Tesoro jutes...@yahoo.com Date: Saturday, October 17, 2009, 1:51 PM

Re: [R] custom selfStart model works with getInitial but not nls

2009-10-17 Thread Douglas Bates
On Fri, Oct 16, 2009 at 7:09 PM, Michael A. Gilchrist mi...@utk.edu wrote: Hello, I'm having problems creating and using a selfStart model with nlme.  Briefly, I've defined the model, a selfStart object, and then combined them to make a selfStart.default model. If I apply getInitial to the

Re: [R] Frequencies, proportions cumulative proportions

2009-10-17 Thread Muenchen, Robert A (Bob)
David, I use CrossTable, so that was my first guess. It'll do proportions/percents by row, column or total in a 2-way table. For 1-way tables, it still tries looks like a 2-way table, unless you specify max.width=1. Then it does one column, but no cumulative proportions (see below). I appreciate

Re: [R] sum two columns with same value

2009-10-17 Thread Gabor Grothendieck
Here are a few ways: Here are a few ways: # 1. using merge m - merge(snag_totale, log_totale, all = TRUE) m[is.na(m)] - 0 with(m, data.frame(AREA, sum = snag_ha + log_ha)) # 2. Use fact that AREA = 1:7 in log_totale sum_totale - log_totale; names(sum_totale)[2] - sum

Re: [R] sum two columns with same value

2009-10-17 Thread Peter Ehlers
smu wrote: On Sat, Oct 17, 2009 at 09:36:50AM +0200, Alfredo Alessandrini wrote: Hi, I've two dataframe: snag_totale AREA snag_ha 12 1.628128 23 10.274249 34 2.778503 45 73.764307 57 12.015985 log_totale AREAlog_ha 11 22.29846 22 17.16889 33

Re: [R] custom selfStart model works with getInitial but not nls

2009-10-17 Thread Michael A. Gilchrist
On Sat, 17 Oct 2009, Douglas Bates wrote: On Fri, Oct 16, 2009 at 7:09 PM, Michael A. Gilchrist mi...@utk.edu wrote: Hello, I'm having problems creating and using a selfStart model with nlme.  Briefly, I've defined the model, a selfStart object, and then combined them to make a

Re: [R] External signal in ODE written in C (using deSolve and approx1?)

2009-10-17 Thread Thomas Petzoldt
Dear Glenn, dear list, this is just a short notice, that a new version 1.5 of package deSolve was released yesterday. It now supports the feature requested below. Details are documented in the package vignette Writing Code in Compiled Languages that comes with the package and is also

[R] repeating values in levels()

2009-10-17 Thread Donald Braman
Can someone help me understand this results? levels(as.factor(miset1$facts_convict)) [1] 1 1 2 3 4 5 6 converting to numeric and back doesn't seem to help: levels(as.factor(as.numeric(miset1$facts_convict))) [1] 1 1 2 3 4 5 6 It's messing up my ologits.  Any way to correct this?

Re: [R] generalization of tabulate()

2009-10-17 Thread Ted Harding
On 16-Oct-09 11:27:06, Gabor Grothendieck wrote: Using the generalized inner product defined in this post: https://www.stat.math.ethz.ch/pipermail/r-help/2006-July/109311.html try this: cbind(S, d = rowSums(inner(S, obs, identical))) The function inner() is defined at the above

[R] How do I access with the name of a (passed) function

2009-10-17 Thread Ajay Shah
How would I do something like this: f - function(x, g) { s - as.character(g) # THIS DOES NOT WORK sprintf(The %s of x is %.0f\n, s, g(x)) } f(c(2,3,4), median) f(c(2,3,4), mean) and get the results The median of x is 3 The mean of x is 3 -- Ajay Shah

[R] lmer function and Inverse mills ratio

2009-10-17 Thread saurav pathak
Dear R users I have two questions, I have been on this problem for last 3 months, please help First question: *How can I use the lmer function for a three level probit ( ie please help me with the command syntax)?* The second question is, *how can I then subsequently calculate the Inverse

[R] Mixing LaTeX and R Code in Loops and Functions in Sweave

2009-10-17 Thread Jacob R. Marcus
I have a question about mixing LaTeX and R code in loops and functions in Sweave. Here's my problem: I want to do something like this: = # some R code for a loop for(i in 1:10) { @ My LaTex code here would describe what I do in the loop. For instance, in this loop print out the numbers one

[R] doing a Tukey HSD post-hoc

2009-10-17 Thread Leslie J Seltzer
Hi all I have a large spreadsheet (Excel) file with many rows and columns in it. I have four treatment groups and using the lm function shows that there is a significant omnibus R with respect to several of the variables. I want to do a Tukey HSD in R. How can I do this? Leslie J. Seltzer,

[R] More polyfit problems

2009-10-17 Thread chris carleton
Hi Everyone, I'm continuing to run into trouble with polyfit. I'm using the fitting function of the form; fit - lm(y ~ poly(x,degree,raw=TRUE)) and I have found that in some cases a polynomial of certain degree can't be fit, the coefficient won't be calculated, because of a singularity. If I

[R] looking for reference that covers convergence in distribution

2009-10-17 Thread Peng Yu
I am looking for a good probability book that describes convergence in distribution. I have looked through Introduction to Probability by Charles M. Grinstead, J. Laurie Snell, but I don't find any formal description on convergence in distribution. Could somebody recommend a good book that cover

Re: [R] How do I access with the name of a (passed) function

2009-10-17 Thread Gabor Grothendieck
Try this: show.name - function(x) deparse(substitute(x)) show.name(pi) [1] pi On Sat, Oct 17, 2009 at 7:26 AM, Ajay Shah ajays...@mayin.org wrote: How would I do something like this: f - function(x, g) {  s - as.character(g)               # THIS DOES NOT WORK  sprintf(The %s of x is

Re: [R] Mixing LaTeX and R Code in Loops and Functions in Sweave

2009-10-17 Thread Gabor Grothendieck
Try this: results=latex,echo=false= for(i in 1:10) { cat(...latex code...) } @ or check out the brew package. On Fri, Oct 16, 2009 at 3:25 PM, Jacob R. Marcus jmar...@u.washington.edu wrote: I have a question about mixing LaTeX and R code in loops and functions in Sweave. Here's my

[R] Let's stay in touch on LinkedIn

2009-10-17 Thread John Cardinale
LinkedIn John Cardinale requested to add you as a connection on LinkedIn: -- Arnaud, I'd like to add you to my professional network on LinkedIn. - John Accept invitation from John Cardinale

[R] R292 on AIX53 using gcc

2009-10-17 Thread Chuck White
Hello -- I am unable to build R 2.9.2 on IBM PowerPC AIX5.3. I would appreciate any help in this matter. ===details== Machine: IBM PowerPC_POWER5 / 4 proc, 1499 MHz 64-bit / AIX 5.3.0.0 Building R 2.9.2 using gcc/g++/gfortran 4.2.4 Config.site

Re: [R] looking for reference that covers convergence in distribution

2009-10-17 Thread Ben Bolker
Peng Yu wrote: I am looking for a good probability book that describes convergence in distribution. I have looked through Introduction to Probability by Charles M. Grinstead, J. Laurie Snell, but I don't find any formal description on convergence in distribution. Could somebody recommend

Re: [R] re peating values in levels()

2009-10-17 Thread Ben Bolker
Donald Braman wrote: Can someone help me understand this results? levels(as.factor(miset1$facts_convict)) [1] 1 1 2 3 4 5 6 converting to numeric and back doesn't seem to help: levels(as.factor(as.numeric(miset1$facts_convict))) [1] 1 1 2 3 4 5 6 It's messing up my ologits.  

Re: [R] doing a Tukey HSD post-hoc

2009-10-17 Thread John Kane
?TukeyHSD --- On Fri, 10/16/09, Leslie J Seltzer lselt...@wisc.edu wrote: From: Leslie J Seltzer lselt...@wisc.edu Subject: [R] doing a Tukey HSD post-hoc To: r-help@r-project.org Received: Friday, October 16, 2009, 5:38 PM Hi all I have a large spreadsheet (Excel) file with many

Re: [R] looking for reference that covers convergence in distribution

2009-10-17 Thread Peter Dalgaard
Ben Bolker wrote: Peng Yu wrote: I am looking for a good probability book that describes convergence in distribution. I have looked through Introduction to Probability by Charles M. Grinstead, J. Laurie Snell, but I don't find any formal description on convergence in distribution. Could

[R] Easy way to `iris[,-Petal.Length]' subsetting?

2009-10-17 Thread Liviu Andronic
Dear all What is the easy way to drop a variable by using its name (and not its number)? Example: data(iris) head(iris) Sepal.Length Sepal.Width Petal.Length Petal.Width Species 1 5.1 3.5 1.4 0.2 setosa 2 4.9 3.0 1.4 0.2

Re: [R] looking for reference that covers convergence in distribution

2009-10-17 Thread Peng Yu
On Sat, Oct 17, 2009 at 3:28 PM, Peter Dalgaard p.dalga...@biostat.ku.dk wrote: Ben Bolker wrote: Peng Yu wrote: I am looking for a good probability book that describes convergence in distribution. I have looked through Introduction to Probability by Charles M. Grinstead, J. Laurie Snell,

[R] optimization problem with constraints...

2009-10-17 Thread kathie
Dear R users, I need some advises on how to use R to optimize a nonlinear function with the following constraints. f(x1,x2,x3,x4,x5,x6) s.t 0 x1 1 0 x2 1 0 x1+x2 1 -inf x3 inf -inf x4 inf 0 x5 inf 0 x6 inf Is there any built-in function or something for these constraint??

[R] Putting names on a ggplot

2009-10-17 Thread John Kane
Putting names on a ggplot Can anyone suggest what I am doing wrong here. I am plotting daily temperatures at Ottawa Ontario for 2008 broken down by months, I seperate them by lines and want to put the names of the months at the top of the chart ( with in the graphing area) Everything is working

Re: [R] Easy way to `iris[,-Petal.Length]' subsetting?

2009-10-17 Thread Gabor Grothendieck
Try this: subset(iris, select = - Petal.Length) On Sat, Oct 17, 2009 at 4:49 PM, Liviu Andronic landronim...@gmail.com wrote: Dear all What is the easy way to drop a variable by using its name (and not its number)? Example: data(iris) head(iris)  Sepal.Length Sepal.Width Petal.Length

Re: [R] Putting names on a ggplot

2009-10-17 Thread ml
hi, On Sat, Oct 17, 2009 at 02:04:43PM -0700, John Kane wrote: Putting names on a ggplot p - p + geom_text(aes(x = namposts + 2.5, y = temprange[2], label = mlabs), data = year, size = 2.5, colour='black', hjust = 0, vjust = 0) you shouldn't use aes in this case since nampost,

Re: [R] repeating values in levels()

2009-10-17 Thread Peter Ehlers
Donald Braman wrote: Can someone help me understand this results? levels(as.factor(miset1$facts_convict)) [1] 1 1 2 3 4 5 6 Don't know how you got your data that way, but I wonder if you've done str() on your data after whatever procedure you used to get to this stage. Here's one way to

[R] predict factanal scores

2009-10-17 Thread zubin
Hello, how does one use factanal to produce factor scores for a new data set? I have a factor solution estimated from historical data. I have a 'new' data set I just like to create factor scores using the prior estimated factor model. I see a reference to a package tsfa, for time series

[R] how to cluster data for use with lmer

2009-10-17 Thread saurav pathak
Dear R users My data set is e names(e) [1] yearctry discent age gender gemeduc gemhhinc ref_groupfearfail_ref knowent_ref nbgoodc_ref [11] nbstatus_ref estbbuso_ref lngdplngdpsq es_gdppcppp sq_gdppcppp estbbo_m es_gdpchg hear I have variables

Re: [R] Putting names on a ggplot

2009-10-17 Thread smu
hi, On Sat, Oct 17, 2009 at 02:04:43PM -0700, John Kane wrote: Putting names on a ggplot p - p + geom_text(aes(x = namposts + 2.5, y = temprange[2], label = mlabs), data = year, size = 2.5, colour='black', hjust = 0, vjust = 0) you shouldn't use aes in this case since nampost,

Re: [R] Mixing LaTeX and R Code in Loops and Functions in Sweave

2009-10-17 Thread Henrik Bengtsson
...or the R.rsp package. -H On Sat, Oct 17, 2009 at 10:37 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: Try this: results=latex,echo=false= for(i in 1:10) {  cat(...latex code...) } @ or check out the brew package. On Fri, Oct 16, 2009 at 3:25 PM, Jacob R. Marcus

Re: [R] RODBC sqlSave does not append the records to a DB2 table

2009-10-17 Thread Felipe Carrillo
Elaine: Try this. It works for me using my own database. Perhaps you don't have admin rights. Good luck library(RODBC) myDB - odbcConnectAccess(DB2.mdb,uid=admin,pwd=) sqlSave(myDB,se2,rownames=FALSE,append=TRUE) close(myDB) Felipe D. Carrillo Supervisory Fishery Biologist

[R] function to convert lm model to LaTeX equation

2009-10-17 Thread Ista Zahn
Dear list, I've tried several times to wrap my head around the Design library, without much success. It does some really nice things, but I'm often uncomfortable because I don't understand exactly what it's doing. Anyway, one thing I really like is the latex.ols() function, which converts an R

Re: [R] RODBC sqlSave does not append the records to a DB2 table

2009-10-17 Thread Orvalho Augusto
Hello! Before you try from R. Can you tell us what happens when you use another DB2 client when you try with those credentials? One thing to note: If you are trying to access DB2 from ODBC on Linux on DB2 versions before 9.4 there are some particular issues and better check this