[R] R & C# (Mono)

2010-07-06 Thread bernde
Hello did anyone ever use C# in connection with R ? i am looking into R extension but would like to use C# instead of C or C++ i wonder whether anyone has experience in particular with Mono for doing so many thanks in advance bernd __ R-help@r-project.

Re: [R] Gray level mosaic plot with shading_Friendly

2010-07-06 Thread Achim Zeileis
On Tue, 6 Jul 2010, Michael Friendly wrote: Michael Kubovy wrote: Suppose we start with data("Titanic") mosaic(Titanic, shade = TRUE) How do I combine the dashed box contours of shading_Friendly to indicate negative residuals, with three levels of gray: dark for abs(Pearson Resid) > 4, ligh

Re: [R] plotmath vector problem; full program enclosed

2010-07-06 Thread Paul Johnson
On Tue, Jul 6, 2010 at 12:41 PM, Duncan Murdoch wrote: > On 06/07/2010 10:54 AM, Paul Johnson wrote: >> >> Here's another example of my plotmath whipping boy, the Normal >> distribution. >> > You want "as.expression(b1)", not "expression(b1)".  The latter means "the > expression consisting of the

Re: [R] Why do s appear in the list of predictor variables in logistic regression using 'step' or 'stepAIC' function?

2010-07-06 Thread Gabor Grothendieck
On Tue, Jul 6, 2010 at 11:46 PM, Kiyoshi Sasaki wrote: > > > Would anyone help me solve my problem with R, please? I am very new to R. I > am doing logistic regression analysis on the presence/absence of salamanders > using several predictor variables, as shown below. I have checked my data, >

[R] Why do s appear in the list of predictor variables in logistic regression using 'step' or 'stepAIC' function?

2010-07-06 Thread Kiyoshi Sasaki
Would anyone help me solve my problem with R, please? I am very new to R. I am doing logistic regression analysis on the presence/absence of salamanders using several predictor variables, as shown below. I have checked my data, but I didn't find any 'NA' or empty cells. When I used step() or s

Re: [R] Gray level mosaic plot with shading_Friendly

2010-07-06 Thread Michael Friendly
Michael Kubovy wrote: Suppose we start with data("Titanic") mosaic(Titanic, shade = TRUE) How do I combine the dashed box contours of shading_Friendly to indicate negative residuals, with three levels of gray: dark for abs(Pearson Resid) > 4, lighter for 4 > abs(Pearson Resid) > 2, and lighte

Re: [R] LatticeExtra Parallel

2010-07-06 Thread Ben Wilkinson
Thanks - that helps put the data on the same scale but it doesn't actually add the scale as it still shows min and max on the axis. On Sun, Jul 4, 2010 at 10:49 PM, Deepayan Sarkar wrote: > On Sun, Jul 4, 2010 at 12:59 PM, Ben Wilkinson > wrote: > > I have put together a chart of 1,000 monthly d

Re: [R] how to define a function in R

2010-07-06 Thread jd6688
Thanks Joshua for the example which has been great help me to start. Wish you the best Jason Ding On Tue, Jul 6, 2010 at 9:23 PM, Joshua Wiley-2 [via R] < ml-node+2280373-448579502-312...@n4.nabble.com > wrote: > Hello, > > As others have said, its hard to give specific advice without specific

Re: [R] How not to print '\\' as '\\'

2010-07-06 Thread Joshua Wiley
Try cat(), for instance: > cat("\\\n") \ > The extra \n is needed for a proper line break. HTH, Josh On Tue, Jul 6, 2010 at 7:01 PM, thmsfuller...@gmail.com wrote: > Hello All, > > '\\' is printed as '\\', but it is actually only one character. > Sometimes, I'd rather print it as a single '\'

Re: [R] xls.getshlib() error in xlsReadWrite package under W7

2010-07-06 Thread Gabor Grothendieck
On Tue, Jul 6, 2010 at 8:50 PM, milton ruser wrote: > Dear all, > > I installed xlsReadWrite package under Windows 7, and got error when trying > to run xls.getshlib(). > Please, see the errors below: > >> require(xlsReadWrite) > Loading required package: xlsReadWrite > xlsReadWrite version (cran

[R] Gray level mosaic plot with shading_Friendly

2010-07-06 Thread Michael Kubovy
Suppose we start with data("Titanic") mosaic(Titanic, shade = TRUE) How do I combine the dashed box contours of shading_Friendly to indicate negative residuals, with three levels of gray: dark for abs(Pearson Resid) > 4, lighter for 4 > abs(Pearson Resid) > 2, and lightest for bs(Pearson Resid)

Re: [R] How not to print '\\' as '\\'

2010-07-06 Thread Jeff Newmiller
?cat "thmsfuller...@gmail.com" wrote: >Hello All, > >'\\' is printed as '\\', but it is actually only one character. >Sometimes, I'd rather print it as a single '\'. Is there a function to >do so in R? > >> nchar('\\') >[1] 1 >> print('\\') >[1] "\\" > >-- >Tom > >__

[R] How not to print '\\' as '\\'

2010-07-06 Thread thmsfuller...@gmail.com
Hello All, '\\' is printed as '\\', but it is actually only one character. Sometimes, I'd rather print it as a single '\'. Is there a function to do so in R? > nchar('\\') [1] 1 > print('\\') [1] "\\" -- Tom __ R-help@r-project.org mailing list https

Re: [R] Stoch Prog in R

2010-07-06 Thread Ravi Varadhan
There is nothing in R that I know of. The task view on optimization shows `stoprog', but it is not available on CRAN. One obvious reason, that I can think of, for the lack of any packages, is that the word "stochastic" is too broad. One needs to know more precisely how randomness enters into

Re: [R] how to define a function in R

2010-07-06 Thread Joshua Wiley
Hello, As others have said, its hard to give specific advice without specific needs, but that's okay; I made up some examples needs and some (rather silly) code that might handle it. Depending what you need to do, it may help you get started. I tried to explicitly name all the arguments in any f

[R] xls.getshlib() error in xlsReadWrite package under W7

2010-07-06 Thread milton ruser
Dear all, I installed xlsReadWrite package under Windows 7, and got error when trying to run xls.getshlib(). Please, see the errors below: > require(xlsReadWrite) Loading required package: xlsReadWrite xlsReadWrite version (cran shlib) Copyright (C) 2010 Hans-Peter Suter, Treetron, Switzerland.

Re: [R] Adding two files into one and vlookup

2010-07-06 Thread Gabor Grothendieck
On Tue, Jul 6, 2010 at 8:26 PM, Gabor Grothendieck wrote: > On Tue, Jul 6, 2010 at 6:18 PM, raghu wrote: >> >> I have two files with dates and prices in each. The number of rows in each of >> them will differ. How do I create a new file which contains data from both >> these files? Cbind and merg

Re: [R] Adding two files into one and vlookup

2010-07-06 Thread Gabor Grothendieck
On Tue, Jul 6, 2010 at 6:18 PM, raghu wrote: > > I have two files with dates and prices in each. The number of rows in each of > them will differ. How do I create a new file which contains data from both > these files? Cbind and merge are not helpful. For cbind because the rows are > not the same

Re: [R] numerical derivative R help

2010-07-06 Thread Ravi Varadhan
Please do not do multiple postings. You do not need numerical derivative, but a refresher course in calculus. fitterma.deriv <- function(xtime) { a <- -0.09144115 b <- -0.01335756 c <- -2.368057 d <- -0.00600052 deriv <- b * exp(a + b*xtime) + d * exp(c + d * xtime) return(deriv) } require

Re: [R] PCA and Regression

2010-07-06 Thread Joris Meys
PCA components are orthogonal by definition so no, that doesn't make sense at all. Do yourself a favor and get a book on multivariate data analysis. Two books come to mind: Obviously the one of Hair and colleagues, called "multivariate data analysis" and easily found in a university library (or on

Re: [R] R2wd- how to open an existing document

2010-07-06 Thread Felipe Carrillo
Hi: I am doing exactly that but I still get the following error message: Error in if (!(tmp[["ActiveDocument"]][["Name"]] == filename)) tmp$Open(paste("path",  :   argument is of length zero   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Servic

Re: [R] Adding two files into one and vlookup

2010-07-06 Thread Erik Iverson
raghu wrote: I have two files with dates and prices in each. The number of rows in each of them will differ. How do I create a new file which contains data from both these files? Cbind and merge are not helpful. For cbind because the rows are not the same replication occurs. Also if I have simila

Re: [R] how to define a function in R

2010-07-06 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of jd6688 > Sent: Tuesday, July 06, 2010 3:49 PM > To: r-help@r-project.org > Subject: [R] how to define a function in R > > > 1. how to write a R script? > 2.How to write a SAS li

Re: [R] Adding two files into one and vlookup

2010-07-06 Thread Peter Alspach
Tena koe Why is merge() not helpful? From your description I would imagine merge(file1, file2, by='Date') would do what you require. HTH Peter Alspach > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of raghu > Sent: W

Re: [R] how to define a function in R

2010-07-06 Thread stephen sefick
I don't know this for sure, but I would probably read an introduction to R... #copy and paste into R your_fun <- function(argument){argument+1} your_fun(6) On Tue, Jul 6, 2010 at 5:48 PM, jd6688 wrote: > > 1. how to write a R script? > 2.How to write a SAS like macro/generic process to process m

[R] Adding two files into one and vlookup

2010-07-06 Thread raghu
I have two files with dates and prices in each. The number of rows in each of them will differ. How do I create a new file which contains data from both these files? Cbind and merge are not helpful. For cbind because the rows are not the same replication occurs. Also if I have similar data how do

[R] how to define a function in R

2010-07-06 Thread jd6688
1. how to write a R script? 2.How to write a SAS like macro/generic process to process multiple files by using the same funstion in R? Thanks in advance -- View this message in context: http://r.789695.n4.nabble.com/how-to-define-a-function-in-R-tp2280290p2280290.html Sent from the R help maili

[R] numerical derivative R help

2010-07-06 Thread Parminder Mankoo
I fit my CDF to sum of exponentials and now I want to take the numerical derivative of this function to obtain probability density.I will really appreciate your help reagrding the error messages I am getting which I don't understand. * * > fitterma <- function(xtime) { a <- -0.09144115 b <- -0.

Re: [R] Odd subsetting behaviour

2010-07-06 Thread math_daddy
That'll do it. Thanks very much. -- View this message in context: http://r.789695.n4.nabble.com/Odd-subsetting-behaviour-tp2279957p2280244.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://sta

Re: [R] HPDinterval question - nonlinear transformations/functions of parameters.

2010-07-06 Thread Christos Argyropoulos
Disregard the post, The answer to the question is actually pretty simple, had I done the math. For completeness here what happens. The definition of the HPD is the following: 1)  \int_{L}^{U}f(\theta|x)d\theta=1-a 2) f(U|x)=f(L|x) A one-to-one repameterization \theta = g(\phi) simultaneous

[R] grayscale wireframe??

2010-07-06 Thread Marlin Keith Cox
I need grayscale formatting for a wireframe. The only col.regions that I can find are color palettes are all colored: rainbow(n, s = 1, v = 1, start = 0, end = max(1,n - 1)/n, gamma = 1, alpha = 1) heat.colors(n, alpha = 1) terrain.colors(n, alpha = 1) topo.colors(n, alpha = 1) cm.colors(n, al

Re: [R] R2wd- how to open an existing document

2010-07-06 Thread Erik Iverson
Having no clue about the specifics here, but what about: wdGet(filename = "myDoc.doc", path = "c:/mydata") You are leaving out what "didn't work" means. What error are you receiving? Felipe Carrillo wrote: Hi: How can one open an existing word document with wdGet() I am getting an error messa

[R] R2wd- how to open an existing document

2010-07-06 Thread Felipe Carrillo
Hi: How can one open an existing word document with wdGet() I am getting an error message when trying to open it like this: wdGet(filename="myDoc.doc",path="c/mydata") Looking at Tal Galili's website example it appears that double backslashes are used but i tried it and didn't work either. Thanks

Re: [R] Visualization of coefficients

2010-07-06 Thread David Atkins
FYI, there is already a function coefplot in the arm package; for example, compare: > library(arm) Loading required package: MASS Loading required package: Matrix [snip] Attaching package: 'arm' The following object(s) are masked from 'package:coda': traceplot > data("Mroz", package =

Re: [R] Odd subsetting behaviour

2010-07-06 Thread Peter Alspach
Tena koe Compare 1:3%/%2 and 1:(3%/%2) HTH ... Peter Alspach > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of math_daddy > Sent: Wednesday, 7 July 2010 5:57 a.m. > To: r-help@r-project.org > Subject: [R] Odd subsetting b

Re: [R] Get the indices of non-zero entries of a sparse matrix in R

2010-07-06 Thread David Winsemius
On Jul 6, 2010, at 3:46 PM, G FANG wrote: Hi, I am trying to get the indices of non-zero entries of a sparse matrix in R sr d 1 1089 3772 1 2 1109 190 1 3 1109 2460 1 4 1109 3071 1 5 1109 3618 1 6 1109 38 1 I found that the following can create a sparse matrix, library(Matrix)

Re: [R] Odd subsetting behaviour

2010-07-06 Thread Erik Iverson
math_daddy wrote: Hello. I've observed some odd behaviour. Most likely, it is already known and explained somewhere, but I can't find an explanation anywhere, so I would appreciate being pointed in the right direction. The issue with the following code is self explanatory: mat <- matrix(c(c(

Re: [R] Selection with changing number of columns

2010-07-06 Thread Peter Alspach
Tena koe Andreas Assuming you mean d[d$Var.1==9 | d$Var.2==9 | d$Var.3==9 ,] not d[d$Var.1==9 | d$Var.1==9 | d$Var.1==9 ,] then d[apply(d, 1, function(x) any(x==9)),] is one way. HTH Peter Alspach > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun

Re: [R] Rcmdr installation under Unbuntu installatiion errors

2010-07-06 Thread stephen sefick
for the build dependencies sudo apt-get build-dep whatever whatever=r-base etc. hth, Stephen Sefick On Tue, Jul 6, 2010 at 3:44 PM, John Fox wrote: > Dear John, > > I see that you already got an answer suggesting that you install the Ubuntu > Rcmdr package, or install the tools required to bu

Re: [R] how to calculate summary statistics for each factor

2010-07-06 Thread Andrew Miles
You could try Summarize in the NCStats package, or aggregate in the epicalc package. Andrew Miles Department of Sociology Duke University On Jul 6, 2010, at 11:53 AM, karena wrote: I have a dataset like the following: subject class value 123110 1241

[R] Odd subsetting behaviour

2010-07-06 Thread math_daddy
Hello. I've observed some odd behaviour. Most likely, it is already known and explained somewhere, but I can't find an explanation anywhere, so I would appreciate being pointed in the right direction. The issue with the following code is self explanatory: > mat <- matrix(c(c(1,1,1),c(2,2,2)),nro

[R] numericDeriv and optim question:

2010-07-06 Thread Parminder Mankoo
*Looking at the complaints below on using numericDeriv and using optim, I don't understand how to use these functions properly. Any help will be appreciated for my example function below. Thanks.* * * *fitterma <- function(xtime) { * *a <- -0.09144115* *b <- -0.01335756* *c <- -2.368057* *d <

[R] Help needed with numericDeriv and optim functions

2010-07-06 Thread Parminder Mankoo
Hello All: I have defined the following function (fitterma as a sum of exponentials) that best fits my cumulative distribution. I am also attaching the "xtime" values that I have. I want to try two things as indicated below and am experiencing problems. Any help will be greatly appreciated. Bes

Re: [R] Rcmdr installation under Unbuntu installatiion errors

2010-07-06 Thread John Fox
Dear John, I see that you already got an answer suggesting that you install the Ubuntu Rcmdr package, or install the tools required to build R packages from source. If the packages shown below are the only ones that failed to install, then you're probably OK. Most are dependencies of dependenci

Re: [R] How to plot confidence bands for nls

2010-07-06 Thread Bert Gunter
Whoops! Seems that the "interval" argument for nls objects is not implemented. Maybe your best approach then is to bootstrap. Bert Gunter Genentech Nonclinical Biostatistics -Original Message- From: bgunter Sent: Tuesday, July 06, 2010 1:09 PM To: Claudia Penaloza; r-help@r-project.o

Re: [R] Pseudo F statistics with index.G1

2010-07-06 Thread Henrik Aldberg
I am afraid the source is not an option at the moment. But it is good to know the where the troubles come from. Thanks a lot for the help. Henrik On 6 July 2010 20:01, Allan Engelhardt wrote: > I guess that in R you have to be explicit about what you want to do. You > can't just drop them, so

Re: [R] How to plot confidence bands for nls

2010-07-06 Thread Bert Gunter
?predict.nls Bert Gunter Genentech Nonclinical Biostatistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Claudia Penaloza Sent: Tuesday, July 06, 2010 12:55 PM To: r-help@r-project.org Subject: [R] How to plot confidence b

[R] HPDinterval question - nonlinear transformations/functions of parameters.

2010-07-06 Thread Christos Argyropoulos
Hi, I have a quick question about HPDinterval (package coda). I am simulating from a trivariate multivariate normal with 3 components (A,B,C) and general (not necessarily diagonal) covariance matrix. Interest lies in describing the distribution of the function: f(A,B,C) = exp(A*B+C) Th

[R] How to plot confidence bands for nls

2010-07-06 Thread Claudia Penaloza
I adjusted an exponential regression to the following data and wish to plot confidence bands as well. Is this possible? Any help greatly appreciated. Claudia x <- c(1989,1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001, 2002,2003,2004,2005,2006,2007,2008,2009) y <- c(987,937,810,74

Re: [R] plotmath vector problem; full program enclosed

2010-07-06 Thread Charles C. Berry
On Tue, 6 Jul 2010, David Winsemius wrote: On Jul 6, 2010, at 1:41 PM, Duncan Murdoch wrote: On 06/07/2010 10:54 AM, Paul Johnson wrote: > Here's another example of my plotmath whipping boy, the Normal > distribution. > > A colleague asks for a Normal plotted above a series of axes that >

Re: [R] question about lpSolve package

2010-07-06 Thread David Winsemius
On Jul 6, 2010, at 3:32 PM, Xiaoxi Gao wrote: Hello R users, I have two quick questions while using "lpSolve" package for linear programming. (1) the result contains both characters and numbers, e.g., Success: the objective function is 40.5, but I only need the number, can I only store

Re: [R] plotmath vector problem; full program enclosed

2010-07-06 Thread Allan Engelhardt
On 06/07/10 18:51, David Winsemius wrote: Easily addressed in this case with "~" instead of "-". The value of "d" provides the minus: b1 <- substitute( mu ~ d*sigma, list(d=round(dividers[1],2)) ) Neat trick! But it gives a slightly different minus sign in the display, so perhaps simply

[R] Get the indices of non-zero entries of a sparse matrix in R

2010-07-06 Thread G FANG
Hi, I am trying to get the indices of non-zero entries of a sparse matrix in R sr d 1 1089 3772 1 2 1109 190 1 3 1109 2460 1 4 1109 3071 1 5 1109 3618 1 6 1109 38 1 I found that the following can create a sparse matrix, library(Matrix) Y <- sparseMatrix(s,r,x=d) but have not idea a

[R] question about lpSolve package

2010-07-06 Thread Xiaoxi Gao
Hello R users, I have two quick questions while using "lpSolve" package for linear programming. (1) the result contains both characters and numbers, e.g., Success: the objective function is 40.5, but I only need the number, can I only store the number? (2) How to set boundaries for variables?

[R] PCA and Regression

2010-07-06 Thread Marino Taussig De Bodonia, Agnese
Hello, I am currently analyzing responses to questionnaires about general attitudes. I have performed a PCA on my data, and have retained two Principal Components. Now I would like to use the scores of both the principal comonents in a multiple regression. I would like to know if it makes sense

Re: [R] Error message using mi() in mi package

2010-07-06 Thread Andrew Miles
On Jul 6, 2010, at 2:15 PM, Erik Iverson wrote: This looks suspiciously like a syntax problem. I would get my text editor to search for 'c14ordered' in the code. You might have missed some punctuation. -Peter Ehlers A good thought. I checked my own code (the stuff coding the data and u

Re: [R] Help With ANOVA (corrected please ignore last email)

2010-07-06 Thread Joshua Wiley
Hello, Are you saying that the -4.60517 values are supposed to be treated as missing? Unless you set them to NA in R, they will be treated as real values. This would make a huge difference. I can tell you that your formula: aov(Intensity ~ Group, data = zzzanova) is treating the variable 'Inten

Re: [R] Error message using mi() in mi package

2010-07-06 Thread Erik Iverson
This looks suspiciously like a syntax problem. I would get my text editor to search for 'c14ordered' in the code. You might have missed some punctuation. -Peter Ehlers A good thought. I checked my own code (the stuff coding the data and using the mi package) and, for good measure, the co

Re: [R] Error message using mi() in mi package

2010-07-06 Thread Andrew Miles
On Jul 6, 2010, at 1:30 PM, Peter Ehlers wrote: On 2010-07-06 10:37, Andrew Miles wrote: Hello! I get the following message when I run the mi() function from the mi package. Error while imputing variable: c3 , model: mi.polr Error in eval(expr, envir, enclos) : could not find function "c14o

Re: [R] Pseudo F statistics with index.G1

2010-07-06 Thread Allan Engelhardt
I guess that in R you have to be explicit about what you want to do. You can't just drop them, so you'll have to assign them some (other) value. Try which(table(C)==1) to give you the values you need to change and then decide what to change them to. The SAS documentation may tell you what it

Re: [R] Rcmdr installation under Unbuntu installatiion errors

2010-07-06 Thread Vojtěch Zeisek
Hello, 1) Try Rkward, it is much more better, 2) You can install Rcmdr directly via r-cran-rcmdr (generaly, search r-cran-*; You can also add special repository, see ) and/or 3) if You wish to install normal packages from source (that is IMHO what ins

Re: [R] plotmath vector problem; full program enclosed

2010-07-06 Thread David Winsemius
On Jul 6, 2010, at 1:41 PM, Duncan Murdoch wrote: On 06/07/2010 10:54 AM, Paul Johnson wrote: Here's another example of my plotmath whipping boy, the Normal distribution. A colleague asks for a Normal plotted above a series of axes that represent various other distributions (T, etc). I wan

Re: [R] lme4

2010-07-06 Thread Douglas Bates
On OS X you need to install the source package for lme4 as the binary package fails one of the tests. We have been unable to reproduce this failure under other operating systems, which makes it hard to debug. On Tue, Jul 6, 2010 at 11:09 AM, Alex Foley wrote: > Hi, > > I was trying to install l

Re: [R] plotmath vector problem; full program enclosed

2010-07-06 Thread Duncan Murdoch
On 06/07/2010 10:54 AM, Paul Johnson wrote: Here's another example of my plotmath whipping boy, the Normal distribution. A colleague asks for a Normal plotted above a series of axes that represent various other distributions (T, etc). I want to use vectors of equations in plotmath to do this, b

Re: [R] how to calculate summary statistics for each factor

2010-07-06 Thread Erik Iverson
?tapply is one way. karena wrote: I have a dataset like the following: subject class value 123110 1241 12 125112 223223 224 2 18 225 219 3233 21 324

Re: [R] Error message using mi() in mi package

2010-07-06 Thread Peter Ehlers
On 2010-07-06 10:37, Andrew Miles wrote: Hello! I get the following message when I run the mi() function from the mi package. Error while imputing variable: c3 , model: mi.polr Error in eval(expr, envir, enclos) : could not find function "c14ordered" Here's the situation: I am running R v. 2.

[R] lme4

2010-07-06 Thread Alex Foley
Hi, I was trying to install lme4 package, but got the following errors: > install.packages("lme4") Warning in install.packages("lme4") : argument 'lib' is missing: using '/Users/xx/Library/R/2.11/library' Warning message: In getDependencies(pkgs, dependencies, available, lib) : package ‘lme

Re: [R] how to save summary(lm) and anova (lm) in format?

2010-07-06 Thread Paul Johnson
There are R packages that can make nice R regression tables in LaTeX documents. I've used memisc and its good, there is also "apsrtable" and the old standby xtable. Also I use my own function "outreg", but that's just a 'not invented here' attitude. Your problem is that you need this to go into W

[R] package error

2010-07-06 Thread Tim Smith
Hi, I was trying to install the 'lme4' package on a mac, but keep getting installation errors. I'm trying to post the actual code that I use, and the output, but it keeps on getting bounced off the message board as spam (Moderator approval required). Is there something that I should be using/

[R] Differencing with auto.arima and xreg

2010-07-06 Thread Kim H
I am having some issues with differencing using auto.arima when also specifying an xreg dataframe. The xreg dataframe contains dummy variables that specify time periods that had a promotion running. When I model diff(y) with order (1,0,1), the coefficients for these dummy variables are very diff

Re: [R] information reduction-database management question

2010-07-06 Thread Paul Chatfield
I don't think the approach would change much with text. You would have to write a function which picks the 'min' or whatever that means to you with text and then it should work ok, Paul From: Brad Patrick Schneid [via R] [mailto:ml-node+2279677-1095983982-120...@n4.nabble.com] Sent: 06 J

[R] how to calculate summary statistics for each factor

2010-07-06 Thread karena
I have a dataset like the following: subject class value 123110 1241 12 125112 223223 224 2 18 225 219 3233 21 324 3 10 325 3

[R] lme4 installation

2010-07-06 Thread Tim Smith
Hi, I was trying to install the package 'lme4'. Here is the code and the sessionInfo() that I am using: > install.packages("lme4",dependencies=T) Warning in install.packages("lme4", dependencies = T) : argument 'lib' is missing: using '/Users/ts/Library/R/2.11/library' Warning message: In ge

[R] plotmath vector problem; full program enclosed

2010-07-06 Thread Paul Johnson
Here's another example of my plotmath whipping boy, the Normal distribution. A colleague asks for a Normal plotted above a series of axes that represent various other distributions (T, etc). I want to use vectors of equations in plotmath to do this, but have run into trouble. Now I've isolated t

[R] Add1 w/ coef estimates?

2010-07-06 Thread darckeen
I was wondering if there is anyway to have Add1() display the coefficient estimates for each candidate predictor along with the F test. This is for lm() btw. Thanks -- View this message in context: http://r.789695.n4.nabble.com/Add1-w-coef-estimates-tp2279662p2279662.html Sent from the R help

Re: [R] Help With ANOVA

2010-07-06 Thread Joris Meys
I still can't reproduce your example. The aov output gives me the following : > anova(aov(Intensity ~ Group, data = zzzanova)) Analysis of Variance Table Response: Intensity Df Sum Sq Mean Sq F value Pr(>F) Group 5 98.85 19.771 2.1469 0.07576 . Residuals 48 442.03 9.209 --- S

Re: [R] Linux-Windows problem

2010-07-06 Thread Duncan Murdoch
On 06/07/2010 11:54 AM, Bos, Roger wrote: Its not an R problem, but an editor problem (therefore slightly off-topic I admit). I do most of the coding in Windows and somestimes the Rkward editor (in linux) will not open the file and it tells me it cannot determine the encoding. If I try to source

[R] ctree ordering nodes

2010-07-06 Thread Paras Sharma
Hello, When using the ctree function, from library (party) what is the syntax to order the Variables in the nodes in a specific way? For example, how would I specify to make a binary come first, then a continuous variable? Also is there a way to force ctree to show variables which are

Re: [R] help with predict.lda

2010-07-06 Thread Changbin Du
Thanks all so much for your help! I went out for 2 days vacation and could not reply your guys email. Yes, the CV=False works. Thanks again! On Sun, Jul 4, 2010 at 2:47 AM, Peter Ehlers wrote: > On 2010-07-03 21:33, Changbin Du wrote: > >> HI, Dear community, >> >> I am using the linear disc

Re: [R] Break in the y-axis

2010-07-06 Thread beloitstudent
Thanks for the advice! It has worked for the most part. However, I am still coming up with an error message when placing my break line in the axis that I'm not sure what it means. If you could help me out, that would be fantastic...otherwise I might just have to see if I can add it on powerpoin

[R] Error message using mi() in mi package

2010-07-06 Thread Andrew Miles
Hello! I get the following message when I run the mi() function from the mi package. Error while imputing variable: c3 , model: mi.polr Error in eval(expr, envir, enclos) : could not find function "c14ordered" Here's the situation: I am running R v. 2.9.2 on Mac OSX v. 10.5.8. I am try

Re: [R] Profiler for R ? (HFWUtils package)

2010-07-06 Thread Hadley Wickham
And the profr package for an alternative display. Hadley On Tuesday, July 6, 2010, Uwe Ligges wrote: > or just see > > ?Rprof > > and > > ?Rprofmem > > > Uwe Ligges > > > On 06.07.2010 01:21, Jim Callahan wrote: > > Message: 21 > Date: Mon, 5 Jul 2010 02:26:29 -0400 > From: Ralf B > To: "r-help@r

[R] Rcmdr installation under Unbuntu installatiion errors

2010-07-06 Thread John Sorkin
Unbuntu 10.04 R 2.10 I am trying to install Rcmdr and receive the following messages: The downloaded packages are in ‘/tmp/RtmpzhjDZG/downloaded_packages’ Warning messages: 1: In install.packages("Rcmdr", dependencies = TRUE) : installation of package 'Rmpi' had non-zero exit status 2: In i

Re: [R] locfit

2010-07-06 Thread David Winsemius
I got a request from a reputable source (who nonetheless did not copy the list, a fact that I only noticed after responding) to post a full worked example: I just tested the bits that I had earlier posted and they seem to fit together. The OP was asking for some sort of density conditional

[R] Trouble Installing R

2010-07-06 Thread Amanda Peterson
Hi, I'm trying to install R on a Solaris 10 machine. >username -a SunOS discovery01 5.10 Generic_142900-09 sun4v sparc SUNW, T5440 I get the following error when running "make install-tests". Any suggestions? r...@discovery01:/opt/R-2.10.1> make install-tests installing specific tests mkdir /usr

Re: [R] Help With ANOVA

2010-07-06 Thread Joris Meys
We're missing the samp1 etc. in order to be able to test the code. Where did you get the other p-value? Cheers Joris On Tue, Jul 6, 2010 at 3:08 PM, Amit Patel wrote: > Hi I needed some help with ANOVA > > I have a problem with My ANOVA > analysis. I have a dataset with a known ANOVA p-value, how

Re: [R] Interpreting NB GLM output - effect sizes?

2010-07-06 Thread Achim Zeileis
On Tue, 6 Jul 2010, Anna Berthinussen wrote: Hi, I am trying to find out how to interpret the summary output from a neg bin GLM? I have 3 significant variables and I can see whether they have a positive or negative effect, but I can't work out how to calculate the magnitude of the effect on th

Re: [R] Linux-Windows problem

2010-07-06 Thread Bos, Roger
Its not an R problem, but an editor problem (therefore slightly off-topic I admit). I do most of the coding in Windows and somestimes the Rkward editor (in linux) will not open the file and it tells me it cannot determine the encoding. If I try to source that file into R (in the linux version)

Re: [R] hide ends of line in a density plot

2010-07-06 Thread Joshua Wiley
Hello, If you're just looking to 'zoom in' as it were, this should do it: plot(density(j), xlim = c(98, 102)) HTH, Josh On Tue, Jul 6, 2010 at 8:50 AM, Albert Vilella wrote: > Hi, > > (I googled for this answer but didn't find anything) > > I am using density plot and I want to trim the ends

Re: [R] Profiler for R ? (HFWUtils package)

2010-07-06 Thread Uwe Ligges
or just see ?Rprof and ?Rprofmem Uwe Ligges On 06.07.2010 01:21, Jim Callahan wrote: Message: 21 Date: Mon, 5 Jul 2010 02:26:29 -0400 From: Ralf B To: "r-help@r-project.org" Subject: [R] Profiler for R ? Hi, is there such a thing as a profiler for R that informs about a) how much process

[R] hide ends of line in a density plot

2010-07-06 Thread Albert Vilella
Hi, (I googled for this answer but didn't find anything) I am using density plot and I want to trim the ends of the line. eg: i = rnorm(100,100,2) j = subset(i,i<102 & i>98) summary(j) plot(density(j)) I only want the line to go from 98 to 102. How can I limit the line (and the axis) to the val

[R] Interpreting NB GLM output - effect sizes?

2010-07-06 Thread Anna Berthinussen
Hi, I am trying to find out how to interpret the summary output from a neg bin GLM? I have 3 significant variables and I can see whether they have a positive or negative effect, but I can't work out how to calculate the magnitude of the effect on the mean of the dependent variable. I used a log

Re: [R] Linux-Windows problem

2010-07-06 Thread Uwe Ligges
On 06.07.2010 15:24, Henrik Bengtsson wrote: Sorry for asking the obvious, but have you confirmed that you are running the same version of R on both systems? Since the question was about the Editor, I doubt the R version is related. /H On Tue, Jul 6, 2010 at 2:11 PM, Bos, Roger wrote:

Re: [R] Help With ANOVA (corrected please ignore last email)

2010-07-06 Thread Joshua Wiley
Hi Amit, When I copy in your data and run aov(Intensity ~ Group, data = zzzanova) I get neither the p-value you showed nor the one you expected. My suggestions at things to look at would be 1) Where/How did you get the expected p-value? Another statistics program (e.g., SPSS or SAS)? It help

Re: [R] Pseudo F statistics with index.G1

2010-07-06 Thread Henrik Aldberg
Thank you Allan, As i understand it the index.G1 function does not work if one of the clusters in the partition only contains one object. Is there a way to get around this in R? In SAS the PSF function seems to ignore the presence of singleton clusters. Sincerely Henrik [[alternative H

Re: [R] selection of optim parameters

2010-07-06 Thread Prof. John C Nash
Without the data and objective function, it is fairly difficult to tell what is going on. However, we can note: - no method is specified, but it would have to be L-BFGS-B as this is the only one that can handle box constraints - the fourth parameter is at a different scale, an

Re: [R] How do I plat timestamed (in seconds) data in R?

2010-07-06 Thread David Winsemius
On Jul 6, 2010, at 9:28 AM, David Winsemius wrote: On Jul 6, 2010, at 8:58 AM, LosemindL wrote: Hi all, I have an Nx2 array, where the first column contains the timestamps and the second column contains the corresponding data. second | ts | -- 14:25:00| 18 14:25:02| 14 14:25:0

Re: [R] How do I plat timestamed (in seconds) data in R?

2010-07-06 Thread Gabor Grothendieck
On Tue, Jul 6, 2010 at 8:58 AM, LosemindL wrote: > > Hi all, > > I have an Nx2 array, where the first column contains the timestamps and the > second column contains the corresponding data. > > second  | ts > | -- > 14:25:00| 18 > 14:25:02| 14 > 14:25:04| 11 > 14:25:06| 4 > 14:25:08| 24 >

Re: [R] Pseudo F statistics with index.G1

2010-07-06 Thread Allan Engelhardt
Always use set.seed in your examples. Running your code after set.seed(1) works fine but gives the error after set.seed(2). The problem in the latter case being that there is only one value 7 in C and you need two or more for the index.G1 code to make sense. Hope this helps a little Allan

Re: [R] How do I plat timestamed (in seconds) data in R?

2010-07-06 Thread David Winsemius
On Jul 6, 2010, at 8:58 AM, LosemindL wrote: Hi all, I have an Nx2 array, where the first column contains the timestamps and the second column contains the corresponding data. second | ts | -- 14:25:00| 18 14:25:02| 14 14:25:04| 11 14:25:06| 4 14:25:08| 24 14:25:10| 13 14:25:12|

  1   2   >