Re: [R] How to Plot With Different Marker ( ‘x ’ and ‘o’) Based on Condition in R

2010-07-06 Thread Peter Ehlers
On 2010-07-05 23:05, Gundala Viswanath wrote: Dear Expert, I have a data that looks like this: for_y_axis-c(0.49534,0.80796,0.93970,0.8) for_x_axis-c(1,2,3,4) count-c(0,33,0,4) What I want to do is to plot the graph using for_x_axis and for_y_axis but will mark each point with o if the

Re: [R] selection of optim parameters

2010-07-06 Thread Rubén Roa
-Mensaje original- De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] En nombre de Fabian Gehring Enviado el: lunes, 05 de julio de 2010 21:53 Para: r-help@r-project.org Asunto: [R] selection of optim parameters Hi all, I am trying to rebuild the results of

[R] Conditional Splitting a Vectors into Two Vectors

2010-07-06 Thread Gundala Viswanath
Suppose I have two vectors of same dimensions: x -c(0.49534,0.80796,0.93970,0.8) count -c(0,33,0,4) How can I group the vectors 'x' into two vectors: 1. Vector `grzero` that contain value in x with `count` value greater than 0 and 2. Vector `eqzero` with value in x with

Re: [R] finding subroutines

2010-07-06 Thread Peter Ehlers
On 2010-07-02 13:23, Erin Hodgess wrote: Dear R People: I have found the starma.c program in R. But now I need to find the R_setup_starma and R_free_starma subroutines as well. Where would I go about finding them, please? These are objects of class NativeSymbolInfo. You can see their list

Re: [R] Conditional Splitting a Vectors into Two Vectors

2010-07-06 Thread Peter Ehlers
On 2010-07-06 0:39, Gundala Viswanath wrote: Suppose I have two vectors of same dimensions: x-c(0.49534,0.80796,0.93970,0.8) count-c(0,33,0,4) How can I group the vectors 'x' into two vectors: 1. Vector `grzero` that contain value in x with `count` value greater than 0 and

Re: [R] Function to compute the multinomial beta function?

2010-07-06 Thread Robin Hankin
It's usually better to build vectorization in to functions: beta3- function (n1, n2, n3) exp(lgamma(n1)+lgamma(n2)+lgamma(n3)-lgamma(n1+n2+n3)) f - function(x){exp(sum(lgamma(x))-lgamma(sum(x)))} beta3(5,3,8) [1] 1.850002e-07 f(c(5,3,8)) [1] 1.850002e-07 rksh On 07/06/2010 01:54 AM,

[R] acf

2010-07-06 Thread nuncio m
Hi list, I have the following code to compute the acf of a time series acfresid - acf(residfit), where residfit is the series when I type acfresid at the prompt the follwoing is displayed Autocorrelations of series ‘residfit’, by lag 0. 0.0833 0.1667 0.2500 0. 0.4167 0.5000

[R] Error in affypdnn package

2010-07-06 Thread Hasan, Ahmed Ryadh - hasar001
Dear all, I am a PhD student working with Affymetrix HGU133atag array for analyzing the Latin square experiment. I was trying to generate gene expression index for hgu133atag array for PDNN model. While extracting the chiptype specific data structure, I got the following error-

[R] Assign Formulas to Arrays or Matrices?

2010-07-06 Thread McLovin
Hi, I am very new to R. I am hoping to create formulas and assign them to locations within an array (or matrix, if it will work). Here's a simplified example of what I'm trying to do: form.arr - array(31,5,3) for (i in seq(from=1, to=31, by=1)) { for (j in seq(from=1, to=5, by=1)) {

Re: [R] Assign Formulas to Arrays or Matrices?

2010-07-06 Thread Peter Ehlers
On 2010-07-06 1:13, McLovin wrote: Hi, I am very new to R. I am hoping to create formulas and assign them to locations within an array (or matrix, if it will work). Here's a simplified example of what I'm trying to do: form.arr- array(31,5,3) for (i in seq(from=1, to=31, by=1)) {

[R] Selection with changing number of columns

2010-07-06 Thread Kunzler, Andreas
Dear list, I'm looking for a way to select rows of a data.frame with changing number of columns (constraint) involved. Assume a data (d) structure like Var.1 Var.2 Var.3 9 2 1 2 9 5 1 2 1 I know the number of involved columns. Is there a way to generate

Re: [R] acf

2010-07-06 Thread Joris Meys
Your question is a bit confusing. acfresidfit is an object, of which we don't know the origin. with your test file, I arrive at the first correlations (but with integer headings) : residfit - read.table(fileree2_test_out.txt) acf(residfit) acfresid - acf(residfit) acfresid Autocorrelations

Re: [R] metafor and meta-analysis at arm-level

2010-07-06 Thread Angelo Franchini
Hello Wolfgang, Thank you very much for your response. When you mentionthe appropriate design matrix, do you mean by that the 'n1i, n2i, m1i, m2i, sd1i, sd2i' arguments of the rma function, or am I missing something? I read the documentation on metafor (introduction), rma/rma.uni and escalc, and

[R] multiple time series plot with dual 'y' axes

2010-07-06 Thread Jorge Hernandez
Hello. I would like to know how to generate dual 'y' axes on a multiple time series plot. I am using ts.plot() to get the multiple time series plot, but I would like a second vertical axis on the right to include another time series on a different scale. Thanks for any help. Cheers.

Re: [R] timeseries

2010-07-06 Thread Joris Meys
Difficult to guess why, but I reckon you should use ts() instead of as.ts. Otherwise set the tsp-attribute correctly. Eg : x - cumsum(1+round(rnorm(20),2)) as.ts(x) Time Series: Start = 1 End = 20 Frequency = 1 [1] 0.87 3.51 4.08 4.20 3.25 4.63 6.30 6.89 9.28 9.93 10.19 9.97 10.20

Re: [R] calculation on series with different time-steps

2010-07-06 Thread Joris Meys
Look at ?ifelse en ?abs, eg : data_frame$new_column_in_dataframe - ifelse(stage$julian_day == baro$julian_day abs(stage$time - baro$hour) = 30, stage$stage.cm - baro$pressure, NA ) Cheers Joris On Thu, Jul 1, 2010 at 10:09 PM, Jeana Lee jeana@colorado.edu wrote: Hello, I have two

[R] Sweave function

2010-07-06 Thread n.via...@libero.it
Dear list, I'm trying to generate a latex Document in which there are a lot of tables. I'm using the Sweave function in the package utils, but I'm having a lot of problems with the format. This is my code: \documentclass[a4paper]{amsbook} \title{Schema di bilancio} \begin{document} \maketitle

Re: [R] Help in the legend()

2010-07-06 Thread Dennis Murphy
Hi: You didn't specify what you meant in your original post by following that I had used didn't gave desired result (sic). I noted two potential problems: (i) the n in the plot is not what you expected it to be, and (ii) the legend didn't render in the range of values established by the graph you

Re: [R] Fast String operations in R ? Cost of String operations

2010-07-06 Thread Martin Maechler
RB == Ralf B ralf.bie...@gmail.com on Mon, 5 Jul 2010 02:33:13 -0400 writes: RB Hi experts, RB currently developing some code that checks a large amount of Strings RB for the existence of sub-strings and pattern (detecting sub-strings RB within URLs). I wonder if there is

[R] Function for gruping similar variables?

2010-07-06 Thread Timo W
Hi, I have a matrix of results of multiple 2x2 chi^2 tests, non- significant tests are marked as TRUE. Is there a function for grouping those variables in a similar way LSD.test from agricolae library does? I reviewed LSD.test's source but it's not helpful for me. This is my matrix: [,1]

Re: [R] Error in affypdnn package

2010-07-06 Thread Henrik Bengtsson
Hi, repost your question to the Bioconductor mailing list instead, and your chances for getting help will be much greater: http://bioconductor.org/docs/mailList.html /Henrik On Tue, Jul 6, 2010 at 9:06 AM, Hasan, Ahmed Ryadh - hasar001 ahmed.ha...@postgrads.unisa.edu.au wrote: Dear all, I

[R] xyplot: filtering out empty plots

2010-07-06 Thread Coen van Hasselt
Hello, I would like to know how I can filter out empty plots in xyplot, when stratifying on some variables. Example: I have a dataset in which I plot CONC ~ TIME, stratified for patient ID(1,2,..,100), FORM(1,2) and BOOST (1,2). Some patients (ID's) do not have values for all stratification

Re: [R] Assign Formulas to Arrays or Matrices?

2010-07-06 Thread Dennis Murphy
Hi: See inline... On Tue, Jul 6, 2010 at 12:13 AM, McLovin dave_dec...@hotmail.com wrote: Hi, I am very new to R. I am hoping to create formulas and assign them to locations within an array (or matrix, if it will work). Here's a simplified example of what I'm trying to do: form.arr -

Re: [R] Sweave function

2010-07-06 Thread Eik Vettorazzi
Hi, I don't know why, but your resulting table is not latex but html. (e.g. amp; is not recognized in latex). NA.string is an argument for print.xtable, not for xtable itself and is case sensitive. So echo=F,results=tex= report2=lapply(report1, function(x) { print(xtable(x),NA.string=-)}) @

[R] Could not find createData function

2010-07-06 Thread shabnam k
Hi, I am using *Maanova* package to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. From that readma object, i have to create data object by using

Re: [R] Assign Formulas to Arrays or Matrices?

2010-07-06 Thread Jens Oehlschlägel
for assigning formulas to arrays use an array of list nr form.arr[[31,5]]y ~ 1 + 2 Jens Oehlschlägel -Ursprüngliche Nachricht- Von: McLovin Gesendet: Jul 6, 2010 9:13:49 AM An: r-help@r-project.org Betreff: [R] Assign Formulas to Arrays or Matrices? Hi, I am very new to R. I am

Re: [R] Assign Formulas to Arrays or Matrices?

2010-07-06 Thread Patrick Burns
On 06/07/2010 08:13, McLovin wrote: Hi, I am very new to R. I am hoping to create formulas and assign them to locations within an array (or matrix, if it will work). Here's a simplified example of what I'm trying to do: form.arr- array(31,5,3) for (i in seq(from=1, to=31, by=1)) {

[R] Error in createData function

2010-07-06 Thread shabnam k
Hi, I am using *Maanova* package to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. From that readma object, i have to create data object by using

[R] how to simplify a table

2010-07-06 Thread Simone Gabbriellini
Hello List, as a result of a table(), I have this: 0 1 3 4 5 6 7 9 11 12 13 14 16 17 18 19 20 21 22 24 27 28 27 2 2 2 2 2 1 2 1 3 2 1 2 1 1 3 1 1 1 2 1 1 I would like to simplify it by grouping like: 0 0-5 6-10 11-15 16-20 21-25 26-30 25 84

Re: [R] how to simplify a table

2010-07-06 Thread Timo W
see: ?cut On 6 Lip, 12:49, Simone Gabbriellini simone.gabbriell...@gmail.com wrote: Hello List, as a result of a table(), I have this:  0  1   3  4  5  6  7  9 11 12 13 14 16 17 18 19 20 21 22 24 27 28 27  2  2  2  2  2  1  2  1   3   2   1   2  1   1   3  1   1   1   2  1   1 I would

[R] locfit

2010-07-06 Thread Arnab Maity
Hi, Can you provide me an example in R to estimate the density using locfit package with the help of multi dimensional explanatory variables and one dimensional dependent variable? Thank you. Arnab Kumar Maity Research Assistant Indian School of Business Hyderabad, India DISCLAIMER:\ This

Re: [R] Could not find createData function

2010-07-06 Thread Martin Morgan
Hi -- On 07/06/2010 03:21 AM, shabnam k wrote: Hi, I am using *Maanova* package to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. From that readma

Re: [R] multiple time series plot with dual 'y' axes

2010-07-06 Thread Jim Lemon
On 07/06/2010 06:45 PM, Jorge Hernandez wrote: Hello. I would like to know how to generate dual 'y' axes on a multiple time series plot. I am using ts.plot() to get the multiple time series plot, but I would like a second vertical axis on the right to include another time series on a

[R] Pseudo F statistics with index.G1

2010-07-06 Thread Kennedy
Hello, I have done some clustering with Agnes and want to calculate the pseudo F statistics using index.G1. It works for a low number of clusters but when i increase the number of clusters i eventually get the following message: Error in apply(x[cl == i, ], 2, mean) : dim(X) must have a

Re: [R] information reduction-database management question

2010-07-06 Thread Paul Chatfield
If you redefine your NAs as below to be detected as some arbitrary large number, then the code should work through. Any 5's left in your dataset can be replaced just as easily by NAs again. Not elegant, but effective. site - c(s1, s1, s1, s2,s2, s2) pref - c(1, 2, 3, 1, 2, 3) R1 - c(NA, NA, 1,

Re: [R] Linux-Windows problem

2010-07-06 Thread Bos, Roger
Uwe, I suspect I might be having a similar problem as the R code I generate in Windows editor (Tinn-R) doesn't always open properly in my Linux editor (Rkward), but I don't know how to change the default encoding in either one. In R on both machines, getOption(encoding) returns native.enc, so

Re: [R] Could not find createData function

2010-07-06 Thread Peter Ehlers
On 2010-07-06 4:21, shabnam k wrote: Hi, I am using *Maanova* package to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. From that readma object, i

Re: [R] locfit

2010-07-06 Thread David Winsemius
On Jul 6, 2010, at 7:08 AM, Arnab Maity wrote: Hi, Can you provide me an example in R to estimate the density using locfit package with the help of multi dimensional explanatory variables and one dimensional dependent variable? When I came up with a solution I posted it:

Re: [R] adding a row of names to data.frame

2010-07-06 Thread jamaas
Thanks Pete, It works fine. I should have explained, already had the data.frame with data in it, just wanted to add a string vector into the first column of the data.frame, could not figure out how to do that. Jim -- View this message in context:

Re: [R] Conditional Splitting a Vectors into Two Vectors

2010-07-06 Thread Christos Argyropoulos
One possible way is the following: x -c(0.49534,0.80796,0.93970,0.8) count -c(0,33,0,4) x[count==0] [1] 0.49534 0.93970 x[count0] [1] 0.80796 0.8 Christos Date: Tue, 6 Jul 2010 15:39:08 +0900 From: gunda...@gmail.com To: r-h...@stat.math.ethz.ch Subject: [R] Conditional

Re: [R] adding a row of names to data.frame

2010-07-06 Thread David Winsemius
On Jul 6, 2010, at 8:00 AM, jamaas wrote: Thanks Pete, It works fine. I should have explained, already had the data.frame with data in it, just wanted to add a string vector into the first column of the data.frame, could not figure out how to do that. ?cbind # generic, has a

Re: [R] A question about conducting crossed random effects in R

2010-07-06 Thread Paul Chatfield
Sounds distinctly like an assignment you've been set for which we wouldn't help. All I'll say is crossed random effects can be dealt with effectively in lmer. See that for more. -- View this message in context:

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

2010-07-06 Thread LosemindL
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| 12 14:25:14| 6 14:25:16| 21 14:25:18| 37 14:25:20|

[R] Help With ANOVA

2010-07-06 Thread Amit Patel
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, however I can not seem to re-create it in R. I have created a list (zzzanova) which contains 1)Intensity Values 2)Group Number (6 Different Groups) 3)Sample Number (54 different

Re: [R] Conditional Splitting a Vectors into Two Vectors

2010-07-06 Thread Henrique Dallazuanna
You can't try this also: split(x, count 0) On Tue, Jul 6, 2010 at 3:39 AM, Gundala Viswanath gunda...@gmail.comwrote: Suppose I have two vectors of same dimensions: x -c(0.49534,0.80796,0.93970,0.8) count -c(0,33,0,4) How can I group the vectors 'x' into two vectors: 1.

Re: [R] Selection with changing number of columns

2010-07-06 Thread Allan Engelhardt
I'm not sure your question completely makes sense, but perhaps this will help: set.seed(1) d- data.frame(matrix(floor(runif(100, max=10)), 10)) # Example data d[apply(d == 9, 1, any), ] # Select rows with 9 in any column ## Or more generally: d[ , c(1, 2, 3)] == c(2, 2, 9)

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

2010-07-06 Thread Amit Patel
Sorry i had a misprint in the appendix code in the last email 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, however I can not seem to re-create it in R. I have created a list (zzzanova) which contains 1)Intensity Values

Re: [R] Linux-Windows problem

2010-07-06 Thread Henrik Bengtsson
Sorry for asking the obvious, but have you confirmed that you are running the same version of R on both systems? /H On Tue, Jul 6, 2010 at 2:11 PM, Bos, Roger roger@rothschild.com wrote: Uwe, I suspect I might be having a similar problem as the R code I generate in Windows editor

Re: [R] Selection with changing number of columns

2010-07-06 Thread Henrique Dallazuanna
Try this: d[colSums(d == 9)] On Tue, Jul 6, 2010 at 5:33 AM, Kunzler, Andreas a.kunz...@bzaek.de wrote: Dear list, I'm looking for a way to select rows of a data.frame with changing number of columns (constraint) involved. Assume a data (d) structure like Var.1 Var.2 Var.3 9 2

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|

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 Gabor Grothendieck
On Tue, Jul 6, 2010 at 8:58 AM, LosemindL comtech@gmail.com 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

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

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,

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

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

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,

[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

[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,i102 i98) 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

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 Bralf.bie...@gmail.com To: r-help@r-project.orgr-help@r-project.org Subject: [R] Profiler for R ? Hi, is there such a thing as a profiler for

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 avile...@gmail.com wrote: Hi, (I googled for this answer but didn't find anything) I am using density plot and I want to

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

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

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 amitrh...@yahoo.co.uk wrote: Hi I needed some help with ANOVA I have a problem with My ANOVA analysis. I have a dataset with a known

[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

[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

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 lig...@statistik.tu-dortmund.de 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

[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

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

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 ehl...@ucalgary.ca wrote: On 2010-07-03 21:33, Changbin Du wrote: HI, Dear community, I am using the

[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

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

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

[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

[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

[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

[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

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

[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

[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

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

[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 ‘lme4’

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.

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

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 alex_foley_...@yahoo.com wrote: Hi, I

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

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 http://cran.r- project.org/bin/linux/ubuntu/) and/or 3) if You wish to install normal packages from source (that is IMHO what

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

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

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

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

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

[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

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

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 b1-

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

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

  1   2   >