Re: [R] R doesn't recognize R_HOME value

2012-07-25 Thread Prof Brian Ripley
On 24/07/2012 21:57, Kirk Fleming wrote: ??? Windows 7 is where I'm seeing the problem. The root problem: while I have R_HOME legitimately specified within Windows 7, and issuing 'set R_HOME' from the command line returns exactly the path I've specified, doing a Sys.getenv('R_HOME') from R

Re: [R] ERROR : cannot allocate vector of size (in MB GB)

2012-07-25 Thread Akkara, Antony (GE Energy, Non-GE)
Thank you Jim. Its working fine !. Thanks a lot. - Antony. -Original Message- From: jim holtman [mailto:jholt...@gmail.com] Sent: Tuesday, July 24, 2012 11:30 PM To: Akkara, Antony (GE Energy, Non-GE) Cc: r-help@r-project.org Subject: Re: [R] ERROR : cannot allocate vector of size (in

Re: [R] reshaping data

2012-07-25 Thread AC Del Re
Wonderful! thanks, Rui! AC On Tue, Jul 24, 2012 at 10:01 PM, Rui Barradas ruipbarra...@sapo.pt wrote: Hello, Try the following. # We are going to use this twice sl - split(long, long$id) # Remove groups with only one row l2 - lapply(sl, function(x) if(nrow(x) 1) x) l2 -

Re: [R] First value in a row

2012-07-25 Thread arun
Hi Camilo, You can either use Henrik's or mine to find it,  unlist(apply(dat1[,-(1:2)],1,function(x) tail(x[!is.na(x)],1)))  x3  x2  x1 0.6 0.3 0.1 #or you can use my functiton dat3-data.frame(NewColumn=c(unlist(lapply(dat2,function(x) tail(x[!is.na(x)],1))),NA))  

Re: [R] Collapsing a vector/data-frame based on the previous values

2012-07-25 Thread arun
Hello, Try this: #dat1 - data  diff2-c(TRUE,ifelse(diff(dat1$State)0|diff(dat1$State)0,TRUE,FALSE))   dat1[diff2,1:2]  #    Date State #1  24/07/2012 1 #6  24/07/2012    -1 #9  27/07/2012 1 #10 28/07/2012    -1 A.K. - Original Message - From: Raghuraman Ramachandran

Re: [R] Wilcoxon V = 0

2012-07-25 Thread Naomikb
This is very helpful, thank you. On Tue, Jul 24, 2012 at 10:54 PM, David Carlson [via R] ml-node+s789695n463770...@n4.nabble.com wrote: It all depends on which order you put the values. Try this wilcox.test(x=b$x,y=a$x,paired=TRUE) Either one is an extreme. A non significant value would be

Re: [R] First value in a row

2012-07-25 Thread arun
Hi Camilo, Forgot dat2: #same as in previous reply. dat2-data.frame(t(dat1[,3:5])) A.K. - Original Message - From: Camilo Mora cm...@dal.ca To: arun smartpink...@yahoo.com Cc: Henrik Singmann henrik.singm...@psychologie.uni-freiburg.de; R help r-help@r-project.org Sent: Tuesday,

Re: [R] First value in a row

2012-07-25 Thread arun
Hi, I think it was due to a row with all the NAs.  I hope it happens rarely. In those cases, you can assign NAs from looking at the list. dat2-data.frame(t(dat1[,3:5])) dat3-lapply(dat2,function(x) tail(x[!is.na(x)],1))  dat3 $X1 [1] 0.6 $X2 [1] 0.3 $X3 [1] 0.1 $X4 numeric(0)  dat3$X4-NA

[R] How to filter datetime from a dataframe

2012-07-25 Thread purushothaman
I have one month data in a dataframe.. I need to filter data based on start date and end date Ex: Start date :- 01/Jul/2012 and End date :- 10/Jul/2012 Please help me to find the solution -- View this message in context:

Re: [R] strapply and characters adjacent to the matched pattern

2012-07-25 Thread arun
HI, Tried matching with data and text using strapply-Unsuccessful.  But, you can get the result from the data alone if that helps you.  dat2-strapply(data,[^\\|],c)  list1-list(paste(dat2[[1]][1:18],collapse=),paste(dat2[[1]][19:37],collapse=))  list1 [[1]] [1] Santa Fe Gold Corp [[2]] [1]

[R] creating Pivot

2012-07-25 Thread namit
Hi Friends, I'm new to R.I have a data frame : xxx having columns colorname values R XXX 10 G YYY 4

Re: [R] First value in a row

2012-07-25 Thread arun
Hi Camilo, Glad it worked well.  You mentioned replacing zeros with 1s.  You can use the same function, replace NA by 1. dat4-ifelse(sapply(dat3,length)==0,1,dat3) A.K. - Original Message - From: Camilo Mora cm...@dal.ca To: arun smartpink...@yahoo.com Cc: Sent: Wednesday,

[R] zeroinfl problem: cannot get standard errors, hessian has NaN

2012-07-25 Thread nikkks
Hi! I have three models. In the first model, everything is fine. However, in the second and third models, I have NA's for standard errors: The hessians also have NaN's (same for m2 and m3). What should I do about it? It there a way to obtain the hessian without transforming my

Re: [R] zeroinfl problem: cannot get standard errors, hessian has NaN

2012-07-25 Thread Joshua Wiley
Hi, This is not nearly enough information. Please follow the posting guide and provide us with a reproducible example, or at the bare minimum, the code for your models. Without more details we can only wildly guess, but here are a few: ---You have more parameters than your data can support

Re: [R] zeroinfl problem: cannot get standard errors, hessian has NaN

2012-07-25 Thread Achim Zeileis
It's hard to say what exactly goes wrong without a reproducible example (see the posting guide). However, a typical source of such problems is that some regressors are scaled badly. Maybe you have one regressor that takes very large values or squares of a regressor that takes moderately large

Re: [R] Annotate forest plot 'forest.rma()' for meta-analysis with metafor package

2012-07-25 Thread Viechtbauer Wolfgang (STAT)
Just replace the fixed y coordinate values in the text() calls with something that will change appropriately with the number of studies/effects included in the plot. For example: text(c(-9.5,-8,-6,-4.5), res$k+2, c(TB+, TB-, TB+, TB-)) text(c(-8.75,-5.25), res$k+3, c(Vaccinated, Control))

Re: [R] creating Pivot

2012-07-25 Thread Jim Lemon
On 07/25/2012 02:41 PM, namit wrote: Hi Friends, I'm new to R.I have a data frame : xxx having columns colorname values R XXX 10 G YYY 4

Re: [R] zeroinfl problem: cannot get standard errors, hessian has NaN

2012-07-25 Thread nikkks
Thank you, Achim! Rescaling worked and all is fine now! -- View this message in context: http://r.789695.n4.nabble.com/zeroinfl-problem-cannot-get-standard-errors-hessian-has-NaN-tp4637715p4637732.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] R doesn't recognize R_HOME value

2012-07-25 Thread peter dalgaard
On Jul 25, 2012, at 07:58 , Prof Brian Ripley wrote: On 24/07/2012 21:57, Kirk Fleming wrote: ??? Windows 7 is where I'm seeing the problem. The root problem: while I have R_HOME legitimately specified within Windows 7, and issuing 'set R_HOME' from the command line returns exactly the

Re: [R] R doesn't recognize R_HOME value

2012-07-25 Thread Prof Brian Ripley
On 25/07/2012 09:14, peter dalgaard wrote: On Jul 25, 2012, at 07:58 , Prof Brian Ripley wrote: On 24/07/2012 21:57, Kirk Fleming wrote: ??? Windows 7 is where I'm seeing the problem. The root problem: while I have R_HOME legitimately specified within Windows 7, and issuing 'set R_HOME'

Re: [R] creating Pivot

2012-07-25 Thread Rui Barradas
Hello, Try the following. dat - read.table(text= colorname values R XXX 10 G YYY 4 Y ZZZ 5 G XXX 2 Y XXX 3 R YYY 2 Y YYY 1 R ZZZ 2 G ZZZ 3 , header=TRUE) d

Re: [R] How to filter datetime from a dataframe

2012-07-25 Thread Anders Ellern Bilgrau
Hi, Check out ?as.Date which is the key function for what you want. Objects of the class Date can be compared using logical operators. Try running the following code. As you didn't provide any data to illustrate your problem I've created some mockup data. data - data.frame(cbind(date =

Re: [R] Plm on splitted data

2012-07-25 Thread goyalprtk
Thanks all for your help... Problem was in my dataset only.. 1 of the splitted part of dataset had all the dependent variables as zero. In other words if dependent variable in all the proposed panels of a dataset is zero than plm will not work -- View this message in context:

[R] lagged variables

2012-07-25 Thread saraberta
hi guys, i have some trouble in creating lagged variables to use as external regressors. i'm trying to use lag(x) but it gives me as result the same time series (x), adding this part at the end: attr(,tsp) [1]0 23231 where do i wrong?are there other functions to be used? thanks sara

[R] error in running glm() function

2012-07-25 Thread shrutibansal
m1.logit -lm(default ~ amt.fac + age.fac + duration + chk_acct + history + purpose + sav_acct +

[R] Installing packages xslx on Ubuntu (32bit)

2012-07-25 Thread Johannes Radinger
Hi, I just recently changed my OS to Ubuntu 12.04 (32bit). Now I tried to install some packages required by my old and working scripts. Unfortunately I fail when trying to install the package xslx. Maybe it is related to the 32bit version of my R (its not possible to install a 64 bit version).

[R] Integrate: compound distribution

2012-07-25 Thread tsakk
Dear All, thanks in advance for your attention. I'm struggling with this for a few weeks now, and figured it's time to get some help. I want to integrate(f(x), lower=-Inf, upper=Inf) with f(x) = ((gamma(K+1)/(gamma(r+1)*gamma(K-r+1)))*(q(x)^r)*(((1-q(x))^(K-r))*phi(x), where phi(x) is the

[R] packages for LIML and JIVE estimators

2012-07-25 Thread Chien-Ho Wang
Hi! I am a new user for R. Recently I plan to write a program to estimate linear model with instrumental variables. I try to use Limit information maximum likelihood method and Jackknife instrumental variable method, but I do not know which R packages I need to use. Whether anybody may inform

[R] packages for LIML and JIVE estimators

2012-07-25 Thread Chien-Ho Wang
Hi! I am a new user for R. Recently I plan to write a program to estimate linear model with instrumental variables. I try to use Limit information maximum likelihood method and Jackknife instrumental variable method, but I do not know which R packages I need to use. Whether anybody may inform

Re: [R] How to filter datetime from a dataframe

2012-07-25 Thread Rui Barradas
Hello, It's not the first time you post questions to R-Help, so 1. Where is the reproducible example? 2. The form you post [part of] your data example forces us to arrange it in order to become R code. Try the following. StartDate - as.Date(01/Jul/2012, format=%d/%B/%Y) EndDate -

Re: [R] Installing packages xslx on Ubuntu (32bit)

2012-07-25 Thread Ivan Calandra
Hi Johannes, As the error message says, install the latest version of R! If you don't want to do it (I don't know why you wouldn't) I think you can also install the package from source if the old version of the package is still somewhere out there. HTH, Ivan -- Ivan CALANDRA Université de

Re: [R] lagged variables

2012-07-25 Thread Rui Barradas
Hello, You are doing nothing wrong. Follow this example. x - ts(1:5) # Seems the same with different start and end lag(x) # But it's not cbind(x, lag(x)) Are there other functions? I know of at least one. (In the end.) These two, though not lags, might also be of interess to you. ?window y -

Re: [R] Installing packages xslx on Ubuntu (32bit)

2012-07-25 Thread Pascal Oettli
Hello, I was able to install it for R version 2.15.1, after reconfiguring R to enable the support of Java (R CMD javareconf) and installing rJava. Did you try after upgrading R? Regards. Pascal Le 25/07/2012 17:47, Johannes Radinger a écrit : Hi, I just recently changed my OS to Ubuntu

Re: [R] Installing packages xslx on Ubuntu (32bit) [On R 2.14.1]

2012-07-25 Thread Prof Brian Ripley
Hmm, there is no such package: did you mean xlsx? However, your R is old and quite a few packages are not available for it. Please do as the posting guide suggests and update to R 2.15.1 (or R-patched). On 25/07/2012 09:47, Johannes Radinger wrote: Hi, I just recently changed my OS to

Re: [R] error in running glm() function

2012-07-25 Thread Rolf Turner
U, doesn't the error message tend to indicate to you that at least one of the factors in the model has only one level? That would make that factor a constant, whence it could not possibly contribute any predictive power, whence it should be eliminated from the model. Note that you *must*

Re: [R] Installing packages xslx on Ubuntu (32bit) [On R 2.14.1]

2012-07-25 Thread Johannes Radinger
Hi Hmm, there is no such package: did you mean xlsx? of course xlsx :), I upgraded to 2.15.1 and installed openjdk-6-jdk and run R CMD javareconf. Then I was successful in installing package xlsx. Thank you! /Johannes However, your R is old and quite a few packages are not available

Re: [R] How to filter datetime from a dataframe

2012-07-25 Thread purushothaman
Thanks for your reply.. The solution what you have provided, working fine as per if the dataframe has only the date in the date-column, But, in my scenario, i have date and time in the date-column. So could you please help me to filter the date and time from the dataframe. Thanks in advance,

[R] reshape - reshape 2: function cast changed?

2012-07-25 Thread Johannes Radinger
Hi, I used to use reshape and moved to reshape2 (R 2.15.1). Now I tried some of my older scripts and was surprised that my cast function wasn't working like before. What I did/want to do: 1) Melt a dataframe based on a vector specifying column names as measure.vars. Thats working so far: dfm -

Re: [R] Between-group variance from ANOVA

2012-07-25 Thread Ista Zahn
There is nothing about R in your question, hence it is not appropriate for this list. Please consult with a local statistician, or post on a stats help list such as http://stats.stackexchange.com/ On Tue, Jul 24, 2012 at 8:55 PM, tedtoal twt...@ucdavis.edu wrote: I'm trying also to understand

Re: [R] How to filter datetime from a dataframe

2012-07-25 Thread Rui Barradas
Hello, inline Em 25-07-2012 11:08, purushothaman escreveu: Thanks for your reply.. The solution what you have provided, working fine as per if the dataframe has only the date in the date-column, But, in my scenario, i have date and time in the date-column. But that is not the data example

Re: [R] lagged variables

2012-07-25 Thread Gabor Grothendieck
On Wed, Jul 25, 2012 at 4:10 AM, saraberta sara.bertape...@hotmail.it wrote: hi guys, i have some trouble in creating lagged variables to use as external regressors. i'm trying to use lag(x) but it gives me as result the same time series (x), adding this part at the end: attr(,tsp) [1]

Re: [R] strapply and characters adjacent to the matched pattern

2012-07-25 Thread Gabor Grothendieck
On Tue, Jul 24, 2012 at 5:06 PM, mdvaan mathijsdev...@gmail.com wrote: Hi, In the example below, one of the searched patterns SE is matched in the word second. I would like to ignore all matches in which the character following the match is one of [:alpha:]. How do I do this without removing

Re: [R] dump.frames and global environment

2012-07-25 Thread Jannis
Thanks again, Hendrik, for your help. Your suggestion would work. You could even wrap the save image and the dump.frames command in one function that is called when an error occours. I now, however, wrapped a dummy function around the call to source that starts my script and supplies a copy

[R] old weird error message

2012-07-25 Thread Ivan Calandra
Dear users, I have a weird questions. A friend of mine, some years ago, supposedly trying to access help files without Internet connection, got something like this error message: The Pythia is not available, please go to Delphi When I've heard about it, I found it very funny and I'd like to

[R] Plotting LDA results

2012-07-25 Thread Meffy
Dear Users! I think I still have some problems in understanding LDA and the methods of plotting the results. The case is the following: I'm having a dataset containing two classes where each datapoint has 19 dimensions. Training with lda(...) works fine, and I'm getting 19 LD coefficients. So far

Re: [R] MANOVA polynomial contrasts

2012-07-25 Thread Manzoni, GianMauro
Dear Greg Snow, thank you very much for your suggestions. However, I need an example in order to understand fully. I was told that, given the ordinal factor, I do not need to specify the contr.poly function because R does it automatically. However, I don not know if I have to add an argument into

[R] Obtain residuals from a Principal Component Analysis

2012-07-25 Thread petohtalrayn
Hi everyone, I am relatively new to R, and I need to perform the principal components analysis of a data matrix. I know that there are a bunch of methods to do it (dudi.pca, princomp, prcomp...) but I have not managed to find a method that can return the residuals obtained by retaining X

[R] bilinear time series model

2012-07-25 Thread Soham
Hello, I am interested in fitting a bilinear time series model in R after determining its order for a given data. Is there any command or function that I can use? -- View this message in context: http://r.789695.n4.nabble.com/bilinear-time-series-model-tp4637760.html Sent from the R help

[R] Boxplot graphic

2012-07-25 Thread phillen
Dear R-users! I boxplotted some data. the class of the data is numeric. There are some outliers and I would like to see their names in the graphic. So, instead that the data points of the outliers are plotted as points, I would like to have their names plotted. First of, how can I give my data

[R] network plot problems

2012-07-25 Thread Vlatka Matkovic Puljic
Hi, I wanted to create a network of drugs that are being studied together. So, I created a file as a graph. But plot of my network is not corect! It looks like to me that something else is lying behind this network links created. Could someone help me with this? Thanx!

Re: [R] Boxplot graphic

2012-07-25 Thread John Kane
PLEASE do read the posting guide. You really need to supply some sample data and the code you are using. We can make guesses about what you are doing but it is much easier to actually look at the data and code. Please use dput() ( see ?dput) to supply some sample data. John Kane Kingston ON

Re: [R] MANOVA polynomial contrasts

2012-07-25 Thread John Fox
Dear Gian, How contrasts are created by default is controlled by the contrasts option: getOption(contrasts) unordered ordered contr.treatment contr.poly So, unless you've changed this option, contr.poly() will be used to generate orthogonal polynomial contrasts for an

Re: [R] Boxplot graphic

2012-07-25 Thread Uwe Ligges
On 25.07.2012 13:26, phillen wrote: Dear R-users! I boxplotted some data. the class of the data is numeric. There are some outliers and I would like to see their names in the graphic. So, instead that the data points of the outliers are plotted as points, I would like to have their names

Re: [R] Boxplot graphic

2012-07-25 Thread S Ellison
Probably the simplest answer is to google 'label boxplot outliers in R'. I found half a dozen solutions, (one of them my own, but that doesn't quite answer your question). Adding names to a series of data can be as simple as saying names(x) - paste(x, 1:length(x)) After that, the solutions

Re: [R] old weird error message

2012-07-25 Thread Michael Weylandt
I dont think it's network connection related. Type ??? :-) Michael On Jul 25, 2012, at 6:49 AM, Ivan Calandra ivan.calan...@u-bourgogne.fr wrote: Dear users, I have a weird questions. A friend of mine, some years ago, supposedly trying to access help files without Internet connection,

Re: [R] Boxplot graphic

2012-07-25 Thread John Fox
Dear phillen, You can use the identify() command following boxplot() to identify outliers, or more simply, the Boxplot() function in the car package, which will do this for you (see the first example in ?Boxplot). I hope this helps, John John

Re: [R] Plotting LDA results

2012-07-25 Thread Uwe Ligges
On 25.07.2012 13:52, Meffy wrote: Dear Users! I think I still have some problems in understanding LDA and the methods of plotting the results. The case is the following: I'm having a dataset containing two classes where each datapoint has 19 dimensions. Training with lda(...) works fine, and

Re: [R] creating Pivot

2012-07-25 Thread John Kane
# note file name is xx. library(reshape2) dcast(xx, name ~ color, value.var = values) John Kane Kingston ON Canada -Original Message- From: saileshchowd...@gmail.com Sent: Tue, 24 Jul 2012 21:41:19 -0700 (PDT) To: r-help@r-project.org Subject: [R] creating Pivot Hi Friends,

Re: [R] Trouble with Installing R Package gsl

2012-07-25 Thread Uwe Ligges
Looking at the gsl code suggest you can set environmant variables that are used by configure, hence setting GSL_CFLAGS = -I/data1/zhucai/local/gsl/include GSL_LIBS = -L/data1/zhucai/local/gsl/lib -lgsl -lgslcblas before R CMD INSTALL should do the trick (untested). Uwe Ligges On

Re: [R] old weird error message

2012-07-25 Thread Ivan Calandra
Thanks for your answer! But I have to type 4 question marks and not 3 That's great :) Ivan -- Ivan CALANDRA Université de Bourgogne UMR CNRS/uB 6282 Biogéosciences 6 Boulevard Gabriel 21000 Dijon, FRANCE +33(0)3.80.39.63.06 ivan.calan...@u-bourgogne.fr

Re: [R] old weird error message

2012-07-25 Thread peter dalgaard
On Jul 25, 2012, at 15:10 , Michael Weylandt wrote: I dont think it's network connection related. Type ??? :-) Actually, type foo -pd Michael On Jul 25, 2012, at 6:49 AM, Ivan Calandra ivan.calan...@u-bourgogne.fr wrote: Dear users, I have a weird questions. A friend of

Re: [R] old weird error message

2012-07-25 Thread John Kerpel
Here's what I got... fooContacting Delphi...the oracle is unavailable. We apologize for any inconvenience. On Wed, Jul 25, 2012 at 8:34 AM, Ivan Calandra ivan.calan...@u-bourgogne.fr wrote: Thanks for your answer! But I have to type 4 question marks and not 3 That's great :) Ivan

Re: [R] How to do the same thing for all levels of a column?

2012-07-25 Thread John Kane
No it's actually telling it to split by the two variables (variable, value) if I understand your question correctly. The confusion is my fault. I tend to be lazy when running examples and did not rename the melt() output to something meaningful. I sometimes forget that it's not

Re: [R] old weird error message

2012-07-25 Thread Yihui Xie
see head(`?`, 20) there is another funny thing: Sys.sleep(2 + rpois(1, 2)) (in theory R can, although is unlikely to, sleep forever since it is a Poisson distribution...) Regards, Yihui -- Yihui Xie xieyi...@gmail.com Phone: 515-294-2465 Web: http://yihui.name Department of Statistics, Iowa

Re: [R] Between-group variance from ANOVA

2012-07-25 Thread S Ellison
I'm trying also to understand how to get the between-group variance out of a one-way ANOVA, but I'm beginning to think that in a sense, the variance does not exist. Emma said: *The model is response(i,j)= group(i)+ error(i,j)* Yes, if by group(i) you mean intercept +

Re: [R] Obtain residuals from a Principal Component Analysis

2012-07-25 Thread Kevin Wright
Just use 'predict' on the fitted model and subtract the predicted values from the data. Kevin Wright On Wed, Jul 25, 2012 at 5:02 AM, petohtalrayn h643...@rtrtr.com wrote: Hi everyone, I am relatively new to R, and I need to perform the principal components analysis of a data matrix. I know

Re: [R] reshape - reshape 2: function cast changed?

2012-07-25 Thread Jeff Newmiller
Not really answering your posed question, but a) reshape2 uses dcast and acast b) The cast step of melt-cast-melt can invoke plyr, and it seems like MCM is just a very inefficient way for you to get at the plyr functionality. Perhaps you should take a more direct route to your goal?

[R] read.xlsx

2012-07-25 Thread Sam Stewart
Hello list, Is anyone else having problems with read.xlsx? I used it on a project at the end of last week, and now the same code is throwing an error: dat=read.xlsx(data.xlsx,sheetIndex=1) Error in .jcall(RJavaTools, Ljava/lang/Object;, invokeMethod, cl, : org.apache.poi.POIXMLException:

[R] network plot problem

2012-07-25 Thread Vlatka Matkovic Puljic
**Hi, I wanted to create a network of drugs that are being studied together. So, I created a file as a graph. But plot of my network is not corect! It looks like to me that something else is lying behind this network links created. Could someone help me with this? Thanx!

[R] On Reproducible Code

2012-07-25 Thread David L Carlson
We often refer requesters to the Posting Guide and chide them for not reading it. Recently I had occasion to re-read the Posting Guide which is for all R lists not just R-help. The word reproducible does not appear anywhere in the guide. The closest it comes is the following suggestion: Sometimes

Re: [R] On Reproducible Code

2012-07-25 Thread Gabor Grothendieck
On Wed, Jul 25, 2012 at 11:50 AM, David L Carlson dcarl...@tamu.edu wrote: We often refer requesters to the Posting Guide and chide them for not reading it. Recently I had occasion to re-read the Posting Guide which is for all R lists not just R-help. The word reproducible does not appear

Re: [R] On Reproducible Code

2012-07-25 Thread Bert Gunter
PLEASE provide commented, minimal, self-contained, reproducible code. Whenever possible, provide a small example that can be easily loaded and run to illustrate your problem. The R function dput() should generally be used to do this. For a more complete discussion of how to

Re: [R] stepwise in svyglm???

2012-07-25 Thread Frank Harrell
None, except when the bootstrap is used correctly to fully document how well or poorly the modeling strategy worked and one is not interested in doing better or hasn't the time to do so. Cheers, Frank Bert Gunter wrote ... which begs the question: In what context is it valid? ;-) -- Bert

Re: [R] package memisc: recode examples

2012-07-25 Thread Marion Wenty
Thank you, John, for the clarification! :) I have written to the R core team pointing this out. Marion 2012/7/24 John Kane jrkrid...@inbox.com -Original Message- From: marion.we...@gmail.com Sent: Tue, 24 Jul 2012 15:48:10 +0200 To: e.vettora...@uke.de Subject: Re: [R]

Re: [R] On Reproducible Code

2012-07-25 Thread Gabor Grothendieck
On Wed, Jul 25, 2012 at 12:25 PM, Bert Gunter gunter.ber...@gene.com wrote: PLEASE provide commented, minimal, self-contained, reproducible code. Whenever possible, provide a small example that can be easily loaded and run to illustrate your problem. The R function dput()

Re: [R] translating IDL to R

2012-07-25 Thread R. Michael Weylandt
On Tue, Jul 24, 2012 at 9:31 PM, Tom Roche tom_ro...@pobox.com wrote: summary: I believe I have ported the GFED IDL example routines to R (following .sig to end of post). But there are some very loose ends, notably 2 for-loops which need replaced by more R-ful code. details: Tom Roche Mon,

[R] Memisc package: imported varibles from SPSS have got wrong measurement

2012-07-25 Thread Marion Wenty
Dear list members, I have got another problem. I imported an SPSS file with the Memisc package using the following commands: mz - spss.system.file(myspssfile.sav) mz - subset(mz,select=c( bsex,balt,xurb,dtaet,kartab,bgeb,boseit,bgeblan,xnuts2,kausb,xerwstat,

[R] Select rows based on matching conditions and logical operators

2012-07-25 Thread kborgmann
Hi, I have a dataset in which I would like to select rows based on matching conditions and return the maximum value of a variable else return one row if duplicate counts exist. My dataset looks like this: PGIDPTIDYear Visit Count 675553121 20091 0 675553121 2009

Re: [R] On Reproducible Code

2012-07-25 Thread Rui Barradas
Hello, This does not mean that the posting guide is useless. Nor that it couldn't or shouldn't be changed. I would say shouldn't because there's a clear call to reproducible code in another part of R, the man files created by package.skeleton: \examples{ ## Should be DIRECTLY

Re: [R] Select rows based on matching conditions and logical operators

2012-07-25 Thread Rui Barradas
Hello, Apart from the output order this does it. (I have changed 'df' to 'df1', 'df' is an R function, the F distribution density.) df1 - read.table(text= PGID PTID Year Visit Count 6755 53121 2009 1 0 6755 53121 2009 2 0 6755 53121 2009 3 0 6755 53122 2008 1 0 6755 53122 2008 2 0 6755

Re: [R] Select rows based on matching conditions and logical operators

2012-07-25 Thread kborgmann
Thanks! which.max did the trick -- View this message in context: http://r.789695.n4.nabble.com/Select-rows-based-on-matching-conditions-and-logical-operators-tp4637809p4637816.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] lagged variables

2012-07-25 Thread Federman, Douglas
There is a Lag function in Hmisc and I found this on StackExchange shift - function (x, shift_by) { #similar to lag function stopifnot(is.numeric(shift_by)) stopifnot(is.numeric(x)) if (length(shift_by)1) return(sapply(shift_by,shift, x=x)) out-NULL

[R] internal error -3 in R_decompress1

2012-07-25 Thread Jannis
Dear R community, i ran into a quite cryptic error message while running some non interactive calculations. The error message reads: Error in formals(FilterTSeriesSSA) : internal error -3 in R_decompress1 This FilterTseriesSSA is a function I wrote. Where could such an error originate

Re: [R] stepwise in svyglm???

2012-07-25 Thread Thomas Lumley
On Tue, Jul 24, 2012 at 2:48 PM, Diana Marcela Martinez Ruiz dianamm...@hotmail.com wrote: Hello, I want to know how to perform stepwise elimination of variables to svyglm If that's actually what you want to do (which as Frank points out, it probably isn't), you'll have to do it by hand. The

Re: [R] R functions to forecast using neural network

2012-07-25 Thread Sajeeka Nanayakkara
Hi, It's only for single hidden layer. Is there any other functions?   Sajeeka Nanayakkara Cc: R help r-help@r-project.org Sent: Wednesday, July 25, 2012 1:03 AM Subject: Re: [R] R functions to forecast using neural network Hi, Possibly ?nnet package.

Re: [R] R functions to forecast using neural network

2012-07-25 Thread Bert Gunter
SEARCH! Google on e.g. neural nets R package ... or search on Rseek.org Check the CRAN Machine Learning task view -- Bert On Wed, Jul 25, 2012 at 11:01 AM, Sajeeka Nanayakkara nsaje...@yahoo.com wrote: Hi, It's only for single hidden layer. Is there any other functions? Sajeeka

Re: [R] Obtain residuals from a Principal Component Analysis

2012-07-25 Thread Kevin Wright
I need to correct myself. I thought that predict.princomp allowed you to specify the number of principal components to use, but that is not the case. A bit more detail is in order. Suppose we have a matrix X, X - data.frame(E1= c(50, 55, 65, 50, 60, 65, 75.), E2= c(67, 71, 76, 80, 82, 89,

Re: [R] internal error -3 in R_decompress1

2012-07-25 Thread Jeff Newmiller
try options (error= recover) and examine the stack frames. --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

Re: [R] package memisc: recode examples

2012-07-25 Thread Peter Ehlers
On 2012-07-25 09:29, Marion Wenty wrote: Thank you, John, for the clarification! :) I have written to the R core team pointing this out. Marion I hope that you included the wording of an improved message. I believe that I speak English reasonably well and I disagree with John on this. I

[R] Simple question on finding duplicates

2012-07-25 Thread Jeff
I'm trying to find duplicate values in a column of a data frame. For example, dataframe (a) below has two 3's. I would like to mark each value of each row as either not being a duplicate of the one before (0), or as a duplicate (1) - for example, as in dataframe (b). In SPSS, I

[R] cforest varUsed

2012-07-25 Thread aelmore
Hi - Does anyone know if there is a way, in a fashion similar to varUsed in randomForest, to get a look at the variables used in a conditional inference forest, by tree? -- View this message in context: http://r.789695.n4.nabble.com/cforest-varUsed-tp4637779.html Sent from the R help mailing

[R] frailtyPenal function

2012-07-25 Thread britt dijkstra
Hello everyone, I am using the frailtyPenal function from the package frailtypack.I am using survival data with recurrent events.There are two different options to handle these in the frailtyPenal function.The AGrecurrent function and the frailty option. The program works if you say that

[R] hyperSpec user survey

2012-07-25 Thread Claudia Beleites
Dear all, I'm looking for users of the hyperSpec package for handling (hyper)spectral or spectroscopic data in R which I maintain. First of all, I made a few announcements concerning the further development which can be found in the hyperSpec-help mailing list an on which I hope to get user

[R] gsl lib version detected wrong during install

2012-07-25 Thread garunkel
Hello everyone, I've tried to install the R gsl package but without success because R claims there needs to be a gsl library version = 1.12 installed on the system. Actually gsl library version 1.12 is installed locally in my home directory and I set PATH and LD_LIBRARY_PATH in .bash_profile to

Re: [R] unable to run spatial lag and error models on large data

2012-07-25 Thread Roger Bivand
shish matt shishm at yahoo.com writes: Hi: First my apologies for cross-posting. A few days back I posted my queries ar R-sig-geo but did not get any response. Hence this post. Since your message never reached that list - did you check? - that isn't very surprising. Try again properly on

Re: [R] Simple question on finding duplicates

2012-07-25 Thread David L Carlson
duplicate - ifelse(c(0, a$col[-length(a$col)])==c(a$col), 1, 0) -- David L Carlson Associate Professor of Anthropology Texas AM University College Station, TX 77843-4352 -Original Message- From: r-help-boun...@r-project.org

Re: [R] Simple question on finding duplicates

2012-07-25 Thread David L Carlson
Minor correction: duplicate - ifelse(c(0, a$col[-length(a$col)])==a$col, 1, 0) --- David -Original Message- From: David L Carlson [mailto:dcarl...@tamu.edu] Sent: Wednesday, July 25, 2012 3:23 PM To: 'Jeff'; 'r-help@r-project.org' Subject: RE: [R] Simple question on finding

Re: [R] Simple question on finding duplicates

2012-07-25 Thread Bert Gunter
ummm... ?duplicates -- Bert On Wed, Jul 25, 2012 at 1:22 PM, David L Carlson dcarl...@tamu.edu wrote: duplicate - ifelse(c(0, a$col[-length(a$col)])==c(a$col), 1, 0) -- David L Carlson Associate Professor of Anthropology Texas AM University

Re: [R] Simple question on finding duplicates

2012-07-25 Thread Bert Gunter
Sorry... ?duplicated -- Bert On Wed, Jul 25, 2012 at 1:28 PM, Bert Gunter bgun...@gene.com wrote: ummm... ?duplicates -- Bert On Wed, Jul 25, 2012 at 1:22 PM, David L Carlson dcarl...@tamu.edu wrote: duplicate - ifelse(c(0, a$col[-length(a$col)])==c(a$col), 1, 0)

Re: [R] Select rows based on matching conditions and logical operators

2012-07-25 Thread William Dunlap
Rui, Your solution works, but it can be faster for large data.frames if you compute the indices of the desired rows of the input data.frame and then using one subscripting call to select the rows instead of splitting the input data.frame into a list of data.frames, extracting the desired row

  1   2   >