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 http://holmwood.

[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 line

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

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 wrote: > Try this: > > <>= > 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 > wrote: >> I have a question about mixing LaT

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,

[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_group""fearfail_ref" "knowent_ref" "nbgoodc_ref" [11] "nbstatus_ref" "estbbuso_ref" "lngdp""lngdpsq" "es_gdppcppp" "sq_gdppcppp" "estbbo_m" "es_g

[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 fact

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

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] 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 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 Petal.Width Species

[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

[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 co

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 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, but I

[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 se

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 somebo

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

2009-10-17 Thread John Kane
?TukeyHSD --- On Fri, 10/16/09, Leslie J Seltzer wrote: > From: Leslie J Seltzer > 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 rows and > columns in it. I

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"

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 recomm

[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 changes

[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 http://www.linkedin.com/e/qlt6CtWzi7sEoE_As_C0_wMfE

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

2009-10-17 Thread Gabor Grothendieck
Try this: <>= 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 wrote: > 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: > >

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 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 %.0f\n", s

[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 th

[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] 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, P

[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 o

[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 Mills

[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

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 a

[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

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 availab

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 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. I

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] 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" sum_totale$sum[snag_

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] 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 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 selfSta

Re: [R] avoiding loops in equation

2009-10-17 Thread Julius Tesoro
thank you very much --- On Sat, 10/17/09, Kenn Konstabel wrote: > From: Kenn Konstabel > Subject: Re: [R] avoiding loops in equation Thank God for R-help mailing list. Thanks.. > To: "Julius Tesoro" > Date: Saturday, October 17, 2009, 1:51 PM > a3 <- sapply(acc, function(x) > pex(x,pga,std),

Re: [R] ozone data

2009-10-17 Thread Liviu Andronic
On 10/17/09, Julius Tesoro 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 source code of any pa

Re: [R] ozone data

2009-10-17 Thread Liviu Andronic
On 10/17/09, Julius Tesoro 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 package. How does on

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

[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 t

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] 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 c

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 merc

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] 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 eigenvalu