[R] Marginal Effects for Beta Regression

2011-08-23 Thread Chris_h
Hi, Im just doing some Beta-Regressions with the betareg package. My question is now, is there a possibility to calculate the marginal effects with the betareg package or is there another package which can handle marginal effects on regression output for the beta class? I try to calculate the

Re: [R] Package Installation in osx

2011-08-23 Thread Prof Brian Ripley
This was garbled en route, but (a) Your R is far too old: please update as the posting guide asked of you. (You seem to be using a pre-release of 2.11.0.) (b) You need to set DISPLAY, as the message says. It usually is set at a Mac OS X console, so perhaps you need to ask on R-sig-mac

Re: [R] How to keep the command 'R CMD INSTALL -l' running, even when some error is encountered?

2011-08-23 Thread Prof Brian Ripley
On Mon, 22 Aug 2011, thmsfuller...@gmail.com wrote: Hi, 'R CMD INSTALL -l' will stop if some error is encountered. I don't find in the manual an option to keep the command running. Is there such an option? No. Why would you want to install a package with errors? If you mean that R CMD

Re: [R] Sequence

2011-08-23 Thread Petr PIKAL
Hi as far as I understand your question it seems to me that round(5*(1.4^(0:10))) gives you your sequence and cumsum(round(5*(1.4^(0:10 gives you summary sequence. Regards Petr I definitely used too much lines s of code because I still don't know how to do some staff. but it

[R] obtaining p-values for lm.ridge() coefficients (package 'MASS')

2011-08-23 Thread Liviu Andronic
Dear all I'm familiarising myself with Ridge Regressions in R and the following is bugging me: How does one get p-values for the coefficients obtained from MASS::lm.ridge() output (for a given lambda)? Consider the example below (adapted from PRA [1]): require(MASS) data(longley) gr -

Re: [R] Time series and ggplot2

2011-08-23 Thread ONKELINX, Thierry
Without a reproducilbe example, it is hard to give you good advise... -Oorspronkelijk bericht- Van: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Namens ashz Verzonden: maandag 22 augustus 2011 23:53 Aan: r-help@r-project.org Onderwerp: [R] Time series and

[R] Glmnet lambda value choice

2011-08-23 Thread Noah Silverman
Hi, When using the glmnet() function of the package glmnet, A series of coefficients is returned for a list of descending lambda values. I am unable to locate anything in the documentation that explains HOW this choice of lambda series is made. (There is documentation about how to choose my

Re: [R] Linear Regression with 2 grouping variables

2011-08-23 Thread Michael Dewey
At 02:15 23/08/2011, Nathan Miller wrote: Hi all, See comment in-line I have a data set that looks a bit like this. feed1 RFU Site Vial Time lnRFU 1 811 10 10.702075 2 4752111 20 10.768927 3 4290511 30 10.66674 4 4686711 40

[R] Survey data: Quantile Regressions

2011-08-23 Thread Anupam
I searched archives for how to do quantile regressions with complex survey data, and there was nothing that could be helpful to a first time user. I am looking for equivalent of the functions in quantreg package. A vignette and/or examples will be very helpful. Is there someone on this list who

Re: [R] likert scale analysis with R

2011-08-23 Thread Iasonas Lamprianou
Dear colleagues, I would like to run a cluster analysis on a number of variables. They are Likert Scale (0 to 10), but they also have a Don't know' option at the end of the scale. Apparently, with the 'Don't Know' option in place, they cannot be considered to be linear or ordinal. How can

Re: [R] likert scale analysis with R

2011-08-23 Thread Jim Lemon
On 08/23/2011 08:39 PM, Iasonas Lamprianou wrote: Dear colleagues, I would like to run a cluster analysis on a number of variables. They are Likert Scale (0 to 10), but they also have a Don't know' option at the end of the scale. Apparently, with the 'Don't Know' option in place, they

Re: [R] Glmnet lambda value choice

2011-08-23 Thread eugen pircalabelu
 Checked the authors' paper on page 8 it explains how the \lambda sequence is chosen. I guess this is what you were looking for. http://www.stanford.edu/~hastie/Papers/glmnet.pdf   Eugen Pircalabelu (0032)471 842 140 (0040)727 839 293 From: Noah Silverman

[R] Change Variable Labels in Quantile Plot

2011-08-23 Thread Kitty Lee
Hi Roger, Maybe I'm missing a clue. Here's an example: x-rnorm(50) x1-rnorm(50) y-x+x1+rnorm(50) out-rq(y~x+x1, tau=1:9/10) plot(summary(out)) plot.out-plot(summary(out)) #I change the variable names dimnames(plot.out)-list(c(intercept, sex, inc)) #Look at the plot #it's ONE simple

Re: [R] likert scale analysis with R

2011-08-23 Thread Iasonas Lamprianou
Thank you for the hint. But this means that I will still have missing responses in my data and this defeats my purpose to analyse the data using cluster analysis.   Dr. Iasonas Lamprianou Department of Social and Political Sciences University of Cyprus From:

[R] Fisher transformation for pooling estimates - p values 1

2011-08-23 Thread lifty . gere
Dear all, I am using Fisher r to z transformation for pooling partial correlation estimates over multiple imputed data (number of imputations = 200). The number of observations in my data is 190. Unfortunately, when i calculate p values for the pooled estimates, some of them are p 1 (ranging

[R] Testing Specific Hypothesis

2011-08-23 Thread Anna Dunietz
Hi All! I am interested in testing whether the means for the data I am investigating are equal to a specific value - let's say 0.01. I have already run a one-way ANOVA and know that the differences in the means are not significant, so now I want to know what values the means take on. otestme is

Re: [R] likert scale analysis with R

2011-08-23 Thread Jim Lemon
On 08/23/2011 09:07 PM, Iasonas Lamprianou wrote: Thank you for the hint. But this means that I will still have missing responses in my data and this defeats my purpose to analyse the data using cluster analysis. If it is really necessary to avoid missing values, and there is an indifference

[R] Matrix:::qr.qy and signature(qr = sparseQR, y = dgCMatrix)

2011-08-23 Thread Felix Lamp
sessionInfo() R version 2.13.1 (2011-07-08) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) locale: [1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1]

[R] NA in last row while reading xlsx file

2011-08-23 Thread ashz
Hi, I am using this script to read a xlsx file to a data frame: library(xlsx) File - file.path(d:, car , car95-99.xlsx) B_car - read.xlsx(File, raw_data) Car2x - data.frame(month = B_car$Date,Ch = B_car$Ch.des, lat=B_car$Latitude) The last row in the data.frame is always NA, how can I remove

[R] Correlation discrepancy

2011-08-23 Thread Vincy Pyne
Dear R list, I have one very elementary question regrading correlation between two variables. x = c(44,46,46,47,45,43,45,44) y = c(44,43,41,41,46,48,44,43) cov(x, y) [1] -2.428571 However, if I try to calculate the covariance using the formula as covariance = sum((x-mean(x))*(y-mean(y)))/8

Re: [R] Correlation discrepancy

2011-08-23 Thread Dimitris Rizopoulos
well, you don't have the correct denominator, i.e., n-1, with n denoting the sample size. Have a look at the *Details* section of the online help file for cov(), and try also sum((x-mean(x))*(y-mean(y)))/7 cov(x, y) I hope it helps. Best, Dimitris On 8/23/2011 1:18 PM, Vincy Pyne wrote:

Re: [R] Correlation discrepancy

2011-08-23 Thread Ted Harding
In addition, something has gone wrong, Vincy, with your data x,y between evaluating cov(x,y) and evaluating your explicit formula. If I repeat your commands: x = c(44,46,46,47,45,43,45,44) y = c(44,43,41,41,46,48,44,43) cov(x, y) # [1] -2.428571 sum((x-mean(x))*(y-mean(y)))/8 # [1]

[R] R system command does not work with objects/variables

2011-08-23 Thread syrvn
Hello! with the R system command I would like to call a perl script which needs an input directory and an output directory in form of a path. When I put in the path directly it works. The script line looks as follows: system(perl '../path1' '../path2' '../path3') If I store the path in

[R] Dummy variable regression

2011-08-23 Thread Marcin P?�ciennik
Dear list members, I want to apply AR(1)-GARCH(1,1) model in order to conduct a test of structural shifts in conditional correlations which I previously estimated. To be more exact, first, I estimate the conditional correlations using the DCC-GARCH model. Now I want to check whether these

Re: [R] R system command does not work with objects/variables

2011-08-23 Thread Prof Brian Ripley
See ?paste or use system2(perl, shQuote(c(p1, p2, p3))) On Tue, 23 Aug 2011, syrvn wrote: Hello! with the R system command I would like to call a perl script which needs an input directory and an output directory in form of a path. When I put in the path directly it works. The script line

Re: [R] Counting Elements Conditionally

2011-08-23 Thread Jean V Adams
There is probably a more elegant way to do this, but this worked for me: # a function to identify the first occurrence of a run of values first - function(x) { l - length(x) c(1, 1-(x[-1]==x[-l])) } # identify the first occurrence of a run of Rev values df$frst - first(df$Rev) #

Re: [R] NA in last row while reading xlsx file

2011-08-23 Thread Jean V Adams
ashz wrote on 08/23/2011 03:25:57 AM: Hi, I am using this script to read a xlsx file to a data frame: library(xlsx) File - file.path(d:, car , car95-99.xlsx) B_car - read.xlsx(File, raw_data) Car2x - data.frame(month = B_car$Date,Ch = B_car$Ch.des, lat=B_car$Latitude) The last row

Re: [R] R system command does not work with objects/variables

2011-08-23 Thread syrvn
Hi, it works great with the paste command! thanks a lot! Best, syrvn -- View this message in context: http://r.789695.n4.nabble.com/R-system-command-does-not-work-with-objects-variables-tp3762544p3762583.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Testing Specific Hypothesis

2011-08-23 Thread Jean V Adams
So, you are looking for confidence intervals for each factor level? You can use the predict() function to do that. fit - aov(values ~ ind, data=otestme) newdat - data.frame(ind=factor(levels(otestme$ind))) cbind(newdat, predict(fit, newdata=newdat, interval=confidence)) Jean Anna Dunietz wrote

Re: [R] Correlation discrepancy

2011-08-23 Thread Vincy Pyne
Dear Mr. Dimitris and Mr Harding, thanks a lot for your guidance. It will be interesting to find out how the Excel deals with this formula. I will try it. Thanks again. Regards Ashok --- On Tue, 8/23/11, ted.hard...@wlandres.net ted.hard...@wlandres.net wrote: From: ted.hard...@wlandres.net

Re: [R] Correlation discrepancy

2011-08-23 Thread Vincy Pyne
Dear Mr Dimitris and Mr Harding, by mistake I have typed my colleagues name (i.e. Ashok) while thanking you. Please excuse me for that. Regards Vincy --- On Tue, 8/23/11, ted.hard...@wlandres.net ted.hard...@wlandres.net wrote: From: ted.hard...@wlandres.net ted.hard...@wlandres.net Subject:

Re: [R] likert scale analysis with R

2011-08-23 Thread Iasonas Lamprianou
Hmmm, maybe I will try this as an excercise, thank you Jim   Dr. Iasonas Lamprianou Department of Social and Political Sciences University of Cyprus From: Jim Lemon j...@bitwrit.com.au To: Iasonas Lamprianou lampria...@yahoo.com Cc: r-help@r-project.org

Re: [R] NA in last row while reading xlsx file

2011-08-23 Thread David Winsemius
On Aug 23, 2011, at 8:21 AM, Jean V Adams wrote: ashz wrote on 08/23/2011 03:25:57 AM: Hi, I am using this script to read a xlsx file to a data frame: library(xlsx) File - file.path(d:, car , car95-99.xlsx) B_car - read.xlsx(File, raw_data) Car2x - data.frame(month = B_car$Date,Ch =

Re: [R] R system command does not work with objects/variables

2011-08-23 Thread Johann Hibschman
syrvn ment...@gmx.net writes: If I store the path in a variable/object and call the perl script again it does not run and I don't know how to overcome that issue. p1 - ../path1 p2 - ../path2 p3 - ../path3 system(perl p1 p2 p3) You want something like: system(paste(perl, p1, p2, p3))

Re: [R] Access Violation error after loading csv file

2011-08-23 Thread Uwe Ligges
Can you please make an example available that reproduces the error, e.g. some code that produces a similar csv file like yours? Uwe Ligges On 23.08.2011 06:29, Mitu De wrote: Hello, I am using Windows 7 and R version 2.13.1. I was trying to read a csv file which has about 1.5 million

Re: [R] Multiple regression in R - unstandardised coefficients a

2011-08-23 Thread Ista Zahn
On Tue, Aug 23, 2011 at 7:54 AM, JC Matthews j.c.matth...@bristol.ac.uk wrote: Thankyou for your replies, you've answered my question and given me more to think on.  I guess it is unwise to draw any conclusions from the standardised results for these reasons. No, by all means try to draw

Re: [R] Reading DESCRIPTION files to create dependency diagram

2011-08-23 Thread Rainer M Krug
2011/8/22 Uwe Ligges lig...@statistik.tu-dortmund.de On 22.08.2011 18:43, Rainer M Krug wrote: Hi I want to create a dependence diagram of a subset of the packages on CRAN and would therefore like to read the DEACRIPTION files into a list. The list should be as follow for each package:

Re: [R] obtaining p-values for lm.ridge() coefficients (package 'MASS')

2011-08-23 Thread Michael Friendly
On 8/23/2011 3:35 AM, Liviu Andronic wrote: [snip] But how does one obtain the customary 'lm' summary information for the model above? I tried supplying the chosen lambda to Design::ols() using its 'penalty' argument, but strangely the results differ. See below. require(Design)

[R] R.oo modify an object inside another classes method

2011-08-23 Thread Ben qant
Can someone show me how to modify one (R.oo) class's object inside another (R.oo) class's method? Is that possible with the R.oo package? A quick example or reference to an example would be outstanding... Thanks, Ben [[alternative HTML version deleted]]

Re: [R] Linear Regression with 2 grouping variables

2011-08-23 Thread Nathan Miller
Thanks Dennis! Worked perfectly. I keep forgetting that plyr can split data based on multiple subsetting variables. Thanks so much, Nate On Mon, Aug 22, 2011 at 10:12 PM, Dennis Murphy djmu...@gmail.com wrote: Hi: You're kind of on the right track, but there is no conditioning formula in

Re: [R] Correlation discrepancy

2011-08-23 Thread Kohta Ishikawa
Divide by 8 leads biased estimator of covariance. R cov function calculates unbiased estimator(divide by (sample size)-1). Regards, Kohta -- View this message in context: http://r.789695.n4.nabble.com/Correlation-discrepancy-tp3762457p3762491.html Sent from the R help mailing list archive at

[R] Passing parameters

2011-08-23 Thread aelmore
Hi, I've got an R script that I'm trying to turn into a ArcGis script tool so that I can run it from ModelBuilder in Arc. Arc isn't giving me any errors when I run the model holding the current tool, but the run time for the R script is 0 seconds. I don't know if the parameters aren't being

Re: [R] Different Estimated values between R and Excel

2011-08-23 Thread Bart Joosen
I you had posted your code which gave the results, we would have seen that you switched your variables. Bart -- View this message in context: http://r.789695.n4.nabble.com/Different-Estimated-values-between-R-and-Excel-tp3762508p3762738.html Sent from the R help mailing list archive at

[R] Different Estimated values between R and Excel

2011-08-23 Thread Komine
Hi, I used simple linear regression with the R software and EXCEL on the same data. Although , I find the same R2=0.84, I find different estimated values (intercept and slope). For the R software (slope =0.0009, Intercept = -0.1478), for EXCEL (slope =927.7, Intercept = 154,41). When I use the

Re: [R] Multiple regression in R - unstandardised coefficients a

2011-08-23 Thread JC Matthews
Thankyou for your replies, you've answered my question and given me more to think on. I guess it is unwise to draw any conclusions from the standardised results for these reasons. James. --On 22 August 2011 17:30 +0100 ted.hard...@wlandres.net wrote: On 22-Aug-11 15:37:40, JC Matthews

[R] duncan test (kruskal post-hoc)

2011-08-23 Thread gaiarrido
Hi, i´ve read all that i could about the post-hoc test for kruskall-Wallis, now i know that i need to run Duncan test, but i´m not able to find in any package or script. I s anybody knows where can i find it? Thanks - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología,

[R] [r] histogram with mean for every break

2011-08-23 Thread Francesco Nutini
Dear R-users, I need to produce a histogram where for every breaks there are the mean of the data. I tried tu use the function hist(x, break=20 ... ) but this return the numerosity for every breaks, not the mean. Any hint? Thanks in advance, francesco

[R] Mapping strings to integers

2011-08-23 Thread Ali.Abbas
Dear All, I am quite a newbie to R. Trying to learn it these days. Sorry for asking such a basic question, but could you kindly tell me how to map unique string values to integers in a data frame? I have a graph which has, most of its, vertices' attributes in a string format. I would like to

[R] P values for vglm(zibinomial) function in VGAM

2011-08-23 Thread suuz
Hi , I know this question has been asked twice in the past but to my knowldege, it still hasn't been solved. I am doing a zero inflated binomial model using the VGAM package, I need to obtain p values for my Tvalues in the vglm output. code is as follows

Re: [R] Different Estimated values between R and Excel

2011-08-23 Thread Komine
Thanks, I see the problem. R did the opposite of Excel. I invert the position of variables in my code and the result is now correct. -- View this message in context: http://r.789695.n4.nabble.com/Different-Estimated-values-between-R-and-Excel-tp3762508p3762870.html Sent from the R help

[R] spectral analysis

2011-08-23 Thread cassie jones
Hi all, I am trying to do spectral analysis of a time series data. But I am not sure how to do it. Can anyone tell me if there are any package I can use to do the analysis using fast fourier transform? Thanks in advance. Cassie [[alternative HTML version deleted]]

[R] Odp: Mapping strings to integers

2011-08-23 Thread Petr PIKAL
Hi Dear All, I am quite a newbie to R. Trying to learn it these days. Sorry for asking such a basic question, but could you kindly tell me how to map unique string values to integers in a data frame? I have a graph which has, most of its, vertices' attributes in a string format. I would

[R] Odp: spectral analysis

2011-08-23 Thread Petr PIKAL
Hi Hi all, I am trying to do spectral analysis of a time series data. But I am not sure how to do it. Can anyone tell me if there are any package I can use to do the analysis using fast fourier transform? Maybe ?fft or ?spectrum Regards Petr Thanks in advance. Cassie

Re: [R] [r] histogram with mean for every break

2011-08-23 Thread R. Michael Weylandt
I do not believe your code (minimal as it is) would work: the correct argument is breaks. More generally, do you really mean to say that hist(x, breaks = 20) immediately returns the bin counts? It doesn't on my machine and if you knew how to get the counts, you should be able to get the midpoints

[R] pMCMC and HPD in MCMCglmm

2011-08-23 Thread m.fen...@libero.it
Dear R users, I’d like to pose aquestion about pMCMC and HDP. I have performed a mixed logistic regression by MCMCglmm (a very good package) obtaining the following results: Iterations = 250001:799901 Thinning interval = 100 Sample size = 5500 DIC: 10.17416 G-structure: ~ID_an post.mean

[R] survival analysis of EEG data

2011-08-23 Thread Michael Simon
Hi everybody, I'm fairly new to survival analysis with R and have some questions how to apply and interpret the coxph and related functions: I have time-dependent covariates with several measurements per subject with constant delta t. The covariates change in each time step. I fitted the

Re: [R] Different Estimated values between R and Excel

2011-08-23 Thread David Winsemius
On Aug 23, 2011, at 7:48 AM, Komine wrote: Hi, I used simple linear regression with the R software and EXCEL on the same data. Although , I find the same R2=0.84, I find different estimated values (intercept and slope). For the R software (slope =0.0009, Intercept = -0.1478), for EXCEL

[R] How to assign vector value as object name

2011-08-23 Thread J. Augusiak
Dear list, I have a table with entries for 20 animals - x and y coordinates to analyze movement with the package adehabitat. The package does all the necessary analyses, but I need to create an object of class ltraj for each animal first. For this kind of object I need to define xy

[R] subsetting a list of matrices

2011-08-23 Thread Lara Poplarski
Hi all, I have an object that looks (roughly) like the following: l - list(a = matrix(rnorm(9), 3), b = matrix(rnorm(9), 3), c = matrix(rnorm(9), 3)) l$a[3,] - sample(c(Message 1, Message 2, Message 3)) l$b[3,] - sample(c(Message 1, Message 2, Message 3)) l$c[3,] - sample(c(Message 1, Message

Re: [R] spectral analysis

2011-08-23 Thread cassie jones
Is there anything in R similar to spectogram command in matlab? On Tue, Aug 23, 2011 at 10:20 AM, Petr PIKAL petr.pi...@precheza.cz wrote: Hi Hi all, I am trying to do spectral analysis of a time series data. But I am not sure how to do it. Can anyone tell me if there are any package

Re: [R] subsetting a list of matrices

2011-08-23 Thread Henrique Dallazuanna
Try this: subset(as.data.frame(do.call(rbind, lapply(l, [, , 1))), row3 == Message 1) On Tue, Aug 23, 2011 at 1:28 PM, Lara Poplarski larapoplar...@gmail.com wrote: Hi all, I have an object that looks (roughly) like the following: l - list(a = matrix(rnorm(9), 3), b = matrix(rnorm(9), 3), c

Re: [R] subsetting a list of matrices

2011-08-23 Thread Justin Haynes
His is better, but you can also use a for loop... out-data.frame(rows=1:3) for(i in 1:3){ if(l[[i]][3]=='Message 1') { out$V1[i]-l[[i]][1] } else { out$V1[i]-NA } } but shouldn't if your list is very long On Tue, Aug 23, 2011 at 9:35 AM, Henrique Dallazuanna www...@gmail.comwrote:

Re: [R] Efficiency of random and fixed effects estimator

2011-08-23 Thread Ben Bolker
Daniel Malter daniel at umd.edu writes: id-rep(c(1:100),each=2) obs-rep(c(0:1),100) d-rep(sample(c(-1,1),100,replace=T),each=2) base.happy-rep(rnorm(100),each=2) happy-base.happy+1.5*d*obs+rnorm(200) data-data.frame(id,obs,d,happy) I am statistically confused tonight. When the

[R] qnorm?

2011-08-23 Thread Jim Silverton
Hi everyone, I have the following problem. I have some small p-values but when I use qnorm(1-4e-30) I get an error. Is there anyway to get around this? -- Thanks, Jim. [[alternative HTML version deleted]] __ R-help@r-project.org mailing

Re: [R] qnorm?

2011-08-23 Thread Peter Langfelder
How about qnorm(4e-30, lower.tail = FALSE) ? You cannot subtract 4e-30 from 1 and expect to get something else than 1 (exactly). Peter On Tue, Aug 23, 2011 at 10:21 AM, Jim Silverton jim.silver...@gmail.com wrote: Hi everyone, I have the following problem. I have some small p-values but when

Re: [R] qnorm?

2011-08-23 Thread Marc Schwartz
On Aug 23, 2011, at 12:21 PM, Jim Silverton wrote: Hi everyone, I have the following problem. I have some small p-values but when I use qnorm(1-4e-30) I get an error. Is there anyway to get around this? Here is a hint: qnorm(1 - 4e-30) [1] Inf qnorm(1) [1] Inf # See ?all.equal

[R] dummy variables from factors

2011-08-23 Thread Noah Silverman
Hi, Looking at a large data set with many factors. I would like to expand each factor variable into multiple new variables for each level. (0,1) coding. My first though was just to code a big nasty loop, to take each level and cbind a column onto my data set. But, that seems painful. There

Re: [R] dummy variables from factors

2011-08-23 Thread Marc Schwartz
On Aug 23, 2011, at 12:34 PM, Noah Silverman wrote: Hi, Looking at a large data set with many factors. I would like to expand each factor variable into multiple new variables for each level. (0,1) coding. My first though was just to code a big nasty loop, to take each level and

Re: [R] dummy variables from factors

2011-08-23 Thread Joshua Wiley
On Tue, Aug 23, 2011 at 10:34 AM, Noah Silverman noahsilver...@ucla.edu wrote: Hi, Looking at a large data set with many factors. I would like to expand each factor variable into multiple new variables for each level. (0,1) coding. My first though was just to code a big nasty loop, to

Re: [R] Odp: Mapping strings to integers

2011-08-23 Thread Ali.Abbas
Thank you very much. I have managed to do the mappings. And sorry for not mentioning the data type. It was a data frame. -- View this message in context: http://r.789695.n4.nabble.com/Mapping-strings-to-integers-tp3762828p3763317.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] [r] histogram with mean for every break

2011-08-23 Thread Francesco Nutini
fantastic Michael! The code V = sort(unique(ave(x,cut(x,HistOut$breaks works good. I got what I need! that is 20 mean values. I calculate the first value in excel, it's the same (79.61429). HistOut= hist(data$ DMP_m3.jaso..10, freq=F, breaks =20, xlim=c(0,3000),col=grey26,

Re: [R] interpreting interactions in a model

2011-08-23 Thread Greg Snow
The Predict.Plot function in the TeachingDemos package can help you visualize interactions. It will work best if Month is treated as a continuous variable. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original

[R] Replacing NAs in one variable with values of another variable

2011-08-23 Thread StellathePug
Hello everyone, I am trying to figure out a way of replacing missing observations in one of the variables of a data frame by values of another variable. For example, assume my data is X X -as.data.frame(matrix(c(9, 6, 1, 3, 9, NA, NA,NA,NA,NA, 6, 4, 3,NA, NA, NA, 5, 4, 1, 3),

Re: [R] Replacing NAs in one variable with values of another variable

2011-08-23 Thread Ista Zahn
Hi, On Tue, Aug 23, 2011 at 12:29 PM, StellathePug ritacarre...@hotmail.com wrote: Hello everyone, I am trying to figure out a way of replacing missing observations in one of the variables of a data frame by values of another variable. For example, assume my data is X X

Re: [R] Replacing NAs in one variable with values of another variable

2011-08-23 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of Ista Zahn Sent: Tuesday, August 23, 2011 11:06 AM To: StellathePug Cc: r-help@r-project.org Subject: Re: [R] Replacing NAs in one variable with values of another variable

[R] Setting Windows environment variables via a batch file

2011-08-23 Thread Dennis Fisher
R 2.13 Vista Colleagues, I have encountered a problem with Windows environment variables that I don't understand. My R code is designed to execute another program that uses Intel Fortran. Intel Fortran apparently has some complicated issues regarding environment variables. As a result,

Re: [R] pMCMC and HPD in MCMCglmm

2011-08-23 Thread Ben Bolker
m.fenati at libero.it m.fenati at libero.it writes: Dear R users, I’d like to pose aquestion about pMCMC and HDP. I have performed a mixed logistic regression by MCMCglmm (a very good package) obtaining the following results: [snip] post.mean l-95% CI u-95% CIeff.samp ID_an

Re: [R] dummy variables from factors

2011-08-23 Thread Noah Silverman
@Josh Thanks! That's exactly what I need. @ Marc. I want to do this manually because I want to do many things with the data beyond R's built in functions. -- Noah Silverman UCLA Department of Statistics 8117 Math Sciences Building Los Angeles, CA 90095 On Aug 23, 2011, at 10:43 AM, Joshua

Re: [R] Setting Windows environment variables via a batch file

2011-08-23 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of Dennis Fisher Sent: Tuesday, August 23, 2011 11:40 AM To: r-h...@stat.math.ethz.ch Subject: [R] Setting Windows environment variables via a batch file R 2.13 Vista

Re: [R] Using the ConText editor?

2011-08-23 Thread Patrick Connolly
On Tue, 23-Aug-2011 at 12:22AM -0400, Jim Holtman wrote: | if you are running on Windows, the install AutoHotKey and this will allow you to setup some macros to accomplish this. this is what notepad++ does and it is easybto extend. | | Sent from my iPad | | On Aug 22, 2011, at 9:07, Robert

Re: [R] spectral analysis

2011-08-23 Thread Ben Bolker
cassie jones cassiejones26 at gmail.com writes: Is there anything in R similar to spectogram command in matlab? Don't know what it does, perhaps you could tell us, but: have you *looked* at example(spectrum) ... ? Ben Bolker __

[R] Bug or feature? sum(c(a, b, c)) != (a + b + c)

2011-08-23 Thread Daniel Lai
Greetings all, I'm porting an algorithm from MATLAB to R, and noticed some minor discrepancies in small decimal values using rowSums and colSums which are exacerbated after heavy iteration and log space transformation. This was rather perplexing as both programs claimed and appeared to use

Re: [R] Change Variable Labels in Quantile Plot

2011-08-23 Thread RKoenker
Hi Roger, Maybe I'm missing a clue. Here's an example: x-rnorm(50) x1-rnorm(50) y-x+x1+rnorm(50) out-rq(y~x+x1, tau=1:9/10) plot(summary(out)) plot.out-plot(summary(out)) #I change the variable names dimnames(plot.out)-list(c(intercept, sex, inc)) Why not just change the names of the

[R] Ideas for tricky subset of XTS object

2011-08-23 Thread Noah Silverman
Hi, I have a very large XTS object. It is about 600,000 entries over 1 year time. I would like to subset a specific piece, by number of days not a specific date. The way I do it now is awkward. Would love to find a way to do this easier to generate a new object just containing the days I

Re: [R] Bug or feature? sum(c(a, b, c)) != (a + b + c)

2011-08-23 Thread R. Michael Weylandt
Not directly related to what you said below, but might I suggest that for numerical work all.equal() might be a little more robust in a computationally heavy implementation. x = c(0.812672,0.916541,0.797810) #dont' call variables c -- just a bad idea y = x[1]+x[2]+x[3] sum(x) ==y [1] FALSE

Re: [R] Bug or feature? sum(c(a, b, c)) != (a + b + c)

2011-08-23 Thread Clint Bowman
?all.equal -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

[R] Opening package manual from within R

2011-08-23 Thread Tyler Rinker
Simple question but searching rseek did not yield the results I wanted. Question: Is there a way to open a help manual for a package from within R. For instance I would like to type a function in r for the tm package and R would open that PDF as seen here:

Re: [R] Opening package manual from within R

2011-08-23 Thread Bert Gunter
After loading the package, does help.start() do what you want? -- Bert On Tue, Aug 23, 2011 at 2:32 PM, Tyler Rinker tyler_rin...@hotmail.com wrote: Simple question but searching rseek did not yield the results I wanted. Question:  Is there a way to open a help manual for a package from

[R] ddply - how to transform df column in place

2011-08-23 Thread jjap
Dear R-users, I am trying to get the plyr syntax right, without much success. Given: d- data.frame(cbind(x=1,y=seq(20100801,20100830,1))) names(d)-c(first, daterep) d2-d # I can convert the daterep column in place the classic way: d$daterep-as.Date(strptime(d$daterep, format=%Y%m%d)) # How to

Re: [R] Opening package manual from within R

2011-08-23 Thread Tyler Rinker
I don't think help.start is what I'm looking for but I may be doing it wrong. I tried: library(tm) help.start(tm) This may be inappropriate as it returns: library(tm) help.start(tm) Error in if (update) make.packages.html(temp = TRUE) : argument is not interpretable as logical Just

Re: [R] Bug or feature? sum(c(a, b, c)) != (a + b + c)

2011-08-23 Thread Barry Rowlingson
On Tue, Aug 23, 2011 at 8:17 PM, Daniel Lai dan...@bccrc.ca wrote: Greetings all, I'm porting an algorithm from MATLAB to R, and noticed some minor discrepancies in small decimal values using rowSums and colSums which are exacerbated after heavy iteration and log space transformation. This

Re: [R] ddply - how to transform df column in place

2011-08-23 Thread Ista Zahn
Hi Jean, On Tue, Aug 23, 2011 at 6:16 PM, jjap jean.plamon...@fpinnovations.ca wrote: Dear R-users, I am trying to get the plyr syntax right, without much success. Given: d- data.frame(cbind(x=1,y=seq(20100801,20100830,1))) names(d)-c(first, daterep) d2-d # I can convert the daterep

Re: [R] ddply - how to transform df column in place

2011-08-23 Thread Justin Haynes
Jean, Ista is right, but: In your function you are asking as.Date to convert the whole data.frame df rather than just your daterep column. out-ddply(d2, .(daterep), function(df) as.Date(strptime(df$daterep,format='%Y%m%d'))) str(out) 'data.frame':30 obs. of 2 variables: $ daterep: num

Re: [R] ddply - how to transform df column in place

2011-08-23 Thread David Winsemius
On Aug 23, 2011, at 6:16 PM, jjap wrote: Dear R-users, I am trying to get the plyr syntax right, without much success. Given: d- data.frame(cbind(x=1,y=seq(20100801,20100830,1))) names(d)-c(first, daterep) d2-d # I can convert the daterep column in place the classic way:

Re: [R] ddply - how to transform df column in place

2011-08-23 Thread David Winsemius
On Aug 23, 2011, at 6:38 PM, David Winsemius wrote: On Aug 23, 2011, at 6:16 PM, jjap wrote: Dear R-users, I am trying to get the plyr syntax right, without much success. Given: d- data.frame(cbind(x=1,y=seq(20100801,20100830,1))) names(d)-c(first, daterep) d2-d # I can convert the

Re: [R] Opening package manual from within R

2011-08-23 Thread Tyler Rinker
David, For me, on a windows machine, help(package=) results in a summary window opening in R rather than the detailed help manual that is available through CRAN. Others suggested help.start() which takes me to a CRAN library site but requires me to still click with the mouse to locate the

Re: [R] Opening package manual from within R

2011-08-23 Thread David Winsemius
Try: help(package=tm) (You do not need library(). ) -- David. On Aug 23, 2011, at 6:17 PM, Tyler Rinker wrote: I don't think help.start is what I'm looking for but I may be doing it wrong. I tried: library(tm) help.start(tm) This may be inappropriate as it returns: library(tm)

Re: [R] Bug or feature? sum(c(a, b, c)) != (a + b + c)

2011-08-23 Thread Daniel Lai
Hi Barry, Shock and horror indeed, addition is not _associative_, at least for floating point arithmetic [1]. Flipping the order of the operands seems to be the explanation of the discrepancies between R and MATLAB as you suggest: In R: a = 0.812672 b = 0.916541 c = 0.797810

Re: [R] Bug or feature? sum(c(a, b, c)) != (a + b + c)

2011-08-23 Thread Gabor Grothendieck
On Tue, Aug 23, 2011 at 6:56 PM, Daniel Lai dan...@bccrc.ca wrote: Hi Barry, Shock and horror indeed, addition is not _associative_, at least for floating point arithmetic [1]. Flipping the order of the operands seems to be the explanation of the discrepancies between R and MATLAB as you

Re: [R] Where are the ticks on grid.xaxis?

2011-08-23 Thread Paul Murrell
Hi Are you looking for the axTicks() function (or possibly the pretty() function) ? Paul On 19/08/2011 3:51 a.m., David A. Johnston wrote: Hi R list, I like the default ticks that are set up using grid.xaxis() or grid.yaxis() with no arguments. Finding good values for the 'at' argument is

  1   2   >