Re: [R] Tobit Fixed Effects

2011-09-14 Thread Daniel Malter
Not satisfactory in which sense? The survreg(Surv(Value,Censoring)~indvars+strata(id)) should/may work. For a discussion of Tobit fixed effects, see also Greene's website: http://pages.stern.nyu.edu/~wgreene/publications.htm under "Fixed Efects and Bias Due to the Incidental Parameters Problem in

Re: [R] Strucchange generating breakpoints

2011-09-14 Thread Achim Zeileis
On Wed, 14 Sep 2011, R. Michael Weylandt wrote: If by gbreakpoints, you simply mean breakpoints (I'm not familiar with the package so this is an honest point of ambiguity from me), It is really gbreakpoints() although it is not quite clear how this was accessed. gbreakpoints() is an unexporte

Re: [R] Can 'mosaic' be used with a continuous variable?

2011-09-14 Thread Achim Zeileis
On Wed, 14 Sep 2011, Mauricio Cornejo wrote: Thanks Z. Let me clarify my problem a bit further ... as I don't think I could use cut() or spine() to solve it. I have a data frame with three columns (A, B, Value).  'A' and 'B' are categorical and 'Value' is continuous and non-negative.  I'm look

[R] Creating a map for central asia countries

2011-09-14 Thread Batur
Dear all, I am trying to create a map for central asian countries(Kazakhstan, Uzbekstan, Kyrgyzstan, Turkmenstan, and Tajikstan). I tried google for shapefiles and Rdata files of central asia, but I can't seem to find them. There is only a world map and individual maps for each countries. If I use

Re: [R] R, Ubuntu, package installation with non-zero exit status

2011-09-14 Thread GeroR
I am running Ubuntu 11.04 natty and R 2.12.1 I have exactly the same message error when try to install any package. I have installed R in others linux's distros without problems... "installation of package 'any package' had non-zero exit status" It doesn't matter if i try to install it from a cran

[R] Model Selection with Phylogenetic Independent Contrasts

2011-09-14 Thread rjswift
I'm trying to select a model under PCA using independent contrasts. Since PICs need to be forced through the origin I've been using lmorigin for the original regression, but it doesn't appear that stepAIC recognizes it. I keep receiving an error message - "Error in na.fail.default(list(Phenology =

[R] Tobit Fixed Effects

2011-09-14 Thread Felipe Nunes
Hi there, I need to run a Tobit Fixed Effects in a panel data with 4500 units for 8 years. It is a huge data set, my dependent variable is left truncated at zero, the distribution is skewed and my panel is balanced. Any suggestions on how to do that in R? I tried stuff like survreg, censReg, and

Re: [R] Optimization package

2011-09-14 Thread Berend Hasselman
Diviya Smith wrote: > > Hi there, > > I have a complex math equation which does not have a closed form solution. > It is - > > y <- (p*exp(-a*d)*(1-exp((d-p)*(a-x[1]/((p-d)*(1-exp(-p*(a-x[1] > > For this equation, I have all the values except for x[1]. So I need to > solve > this probl

Re: [R] Reading Parts of Excel Files (within a sheet)

2011-09-14 Thread R. Michael Weylandt
I haven't used it, but I believe the XLConnect package allows for control on a cell-by-cell basis. Check out the extensive example given here: http://www.r-bloggers.com/xlconnect-%E2%80%93-a-platform-independent-interface-to-excel/ Hope this helps, Michael Weylandt On Wed, Sep 14, 2011 at 8:43 P

Re: [R] Image processing and analysis with R

2011-09-14 Thread Carl Witthoft
This isn't R, but you might find it quite useful: http://rsbweb.nih.gov/ij/features.html ImageJ I've written centroid, Hough Transform, Sobel Transform, etc. in R, so certainly it's possible to do image analysis. You just may find other tools better suited to your needs. Carl -- - Sen

Re: [R] portfolio, portfolio.optim function not found

2011-09-14 Thread Jean-Christophe BOUËTTÉ
portfolio.optim {tseries} means this function is in the tseries package... 2011/9/14 James Splinter : > Hello, > > After installing and loading the package "portfolio", I tried to run the > example code provided, and it would not run. > > this is the link: > http://rss.acs.unt.edu/Rdoc/library/tse

Re: [R] Loops

2011-09-14 Thread R. Michael Weylandt
Look at the combinations function in gtools and index by the rows of that output for a single loop. Pass values as the three parameters. Sorry for being terse - writing on my phone. Michael On Sep 14, 2011, at 5:36 PM, Fred wrote: > Dear forum, > > I would like to forecast e.g. with the ar

Re: [R] Loops

2011-09-14 Thread Weidong Gu
A little bit change of your code may do the work foreL<-8 b0f<-matrix(nrow=9, ncol=foreL) ct<-1 ### use this as the index of b0f for(ar.ord in 1:3){ for(ma.ord in 1:3){ b0f[ct,]<-c(predict(arima(para_qtr[1:(n-8),1],order=c(ar.ord,1,ma.ord)), n.ahead=foreL)$pred) ct<-ct+1 ### increment the counter

[R] Reading Parts of Excel Files (within a sheet)

2011-09-14 Thread Peter Lomas
Dear R Users, I have to read data from many excel spreadsheets, all which have some frustrating formatting (lots of titles, headers, etc.). I am trying to work directly from source data and the number of the spreadsheets I would have to go through make reformatting one by one a pain. I have foun

Re: [R] Color barplots with a conditional?

2011-09-14 Thread R. Michael Weylandt
You've almost got it. If your category variables are factors, then adding col=data$category to plot() will do it automatically since factors will be coerced to their integer representation. If you want a different color map, you can do arithmetic transformations on the integer codings directly a

[R] Questons about 'igraph' package

2011-09-14 Thread karena
Hi, I am using 'igraph' to make some plots. The problem I got is that I don't know how to label the nodes with gene names. My sample code: ## suppose I have 100 gene (nodes) ## --- graph <- set.vertex.attribute(graph, "color

[R] problem with predict.gam in package gam

2011-09-14 Thread Drew Tyre
I am trying to reproduce plots in Chapter 3 of Zuur et al Mixed Effects models and extensions in Ecology. For pedagogical reasons, they make a series of plots with gam(...) in package gam. I encounter errors that trace back to the predict.gam method. Below is a repeatable example using one of the e

[R] Loops

2011-09-14 Thread Fred
Dear forum, I would like to forecast e.g. with the arima-model. To figure out which model works best I am going to predict with this models. my first code: for(ar.ord in 1:3){ for(ma.ord in 1:3){ print(predict(arima(para_qtr[1:(n-8),1],order=c(ar.ord,1,ma.ord)), n.ahead=8)$pred) } } this one wo

[R] Color barplots with a conditional?

2011-09-14 Thread Allie818
I've made a barplot that has several bars. I'd like the bars to be colored according to the type of category they are in. Is there a way to put a conditional into the color parameter in barplot? i.e. if data$category[i] == 1, color the bar red Do I need to generate the coloring sequence before the

[R] portfolio, portfolio.optim function not found

2011-09-14 Thread James Splinter
Hello, After installing and loading the package "portfolio", I tried to run the example code provided, and it would not run. this is the link: http://rss.acs.unt.edu/Rdoc/library/tseries/html/portfolio.optim.html this is the example code, as found at the link: x <- rnorm(1000) dim(x) <- c(500,2

Re: [R] svytable including NAs

2011-09-14 Thread tesutton
Thanks for the replies. I'll go with redefining the factor levels for now. Thanks very much Tim -- Tim Elwell-Sutton The University of Hong Kong Office: 3906 2053 Mob: 6084 5654 From: Thomas Lumley [tlum...@uw.edu] Sent: Thursday, September 15, 2011 5:36

[R] Hints for Data Mining

2011-09-14 Thread Lorenzo Isella
Dear All, I am recycling a previous email of mine where I asked some questions about clustering mixed numerical/categorical data. This time I am more into data mining. I am given a set of known statistical indexes {s_i}, i=1,2...N for a N countries. These indexes in general are a both numerica

Re: [R] Image processing and analysis with R

2011-09-14 Thread Mikkel Grum
http://www.sciviews.org/zooimage/ might be useful. From: Jose Bustos Melo To: "r-help@r-project.org" Sent: Wednesday, September 14, 2011 4:06 PM Subject: [R] Image processing and analysis with R Hello everyone,  I'm working in a project to create an special

[R] Optimization package

2011-09-14 Thread Diviya Smith
Hi there, I have a complex math equation which does not have a closed form solution. It is - y <- (p*exp(-a*d)*(1-exp((d-p)*(a-x[1]/((p-d)*(1-exp(-p*(a-x[1] For this equation, I have all the values except for x[1]. So I need to solve this problem numerically. Can anyone suggest an optimi

Re: [R] Weibull point process

2011-09-14 Thread Rolf Turner
On 15/09/11 07:21, Torbjørn Ergon wrote: Dear list, I'm looking for a function to generate (simulate) a random Weibull point process. Can anyone help? Cheers, Torbjørn Ergon, University of Oslo Do you mean a renewal process with the inter-event times having a Weibull distribution? Should

Re: [R] svytable including NAs

2011-09-14 Thread Thomas Lumley
On Thu, Sep 15, 2011 at 1:49 AM, David Winsemius wrote: > > On Sep 14, 2011, at 6:12 AM, tesutton wrote: > >> Dear all >> >> I'm creating two way tables based on weighted survey data using svytable >> and want to include NAs as a level in the table. >> For unweighted data I use table(..., useNA="i

[R] Image processing and analysis with R

2011-09-14 Thread Jose Bustos Melo
Hello everyone,  I'm working in a project to create an special analysis about visualization in some animals and its behaviour, through checking eye movement. The point is this is very important to get data from the image as a resume of what this animal sees. I think an RGB histogram like photo

[R] Fw: Difficulty with saving .Rhistory

2011-09-14 Thread T
HTH, Hahaha, not sure why but I must have filled out T, but I thought that I filled the name out as Tommy...it truncated somewhere.  TC... "Theodore Calvin" the chopper pilot from Vietnam on Magnum PI?   ;-)  Nope, not me. Sorry all!  I fixed it.  It was that the group was also not root.  I

Re: [R] Weibull point process

2011-09-14 Thread Ken Hutchison
I don't know if this settles the matter, but if you are modeling Weibull as the Log-Intensity I.E. in a non-homegenous poisson process you can: #Where T is the length of the series desired Weibull.process<-function(T,shape,scale) { logLambda=rweibull(T,shape,scale) Lambda=exp(logLambda) Point.Pro

Re: [R] S4 vs Reference Classes

2011-09-14 Thread Steve Lianoglou
Hi, Just wanted to say that embedding a slot in your class that's an environment (as I shown earlier) will still solve your problem w/o you having to switch to Ref classes (since you've already done lots of work for your app in S4). Let's assume you have a slot `cache` that is an environment, usi

Re: [R] Open file for reading and writing with APPEND

2011-09-14 Thread Henrik Bengtsson
See details on *argument* 'open' for file(), cf. help("file"). /Henrik On Wed, Sep 14, 2011 at 12:38 PM, Noah Silverman wrote: > Hi, > > I want to keep a running file of some statistics generated by a running > process.  Occasionally the process will be restarted. > > On  restart, I'd like to r

Re: [R] Strucchange generating breakpoints

2011-09-14 Thread R. Michael Weylandt
If by gbreakpoints, you simply mean breakpoints (I'm not familiar with the package so this is an honest point of ambiguity from me), the documentation says the following: Value An object of class "breakpoints" is a list with the following elements: breakpointsthe breakpoints of the optimal partit

Re: [R] S4 vs Reference Classes

2011-09-14 Thread Joseph Park
Gentlemen: Steve, Martin & Doug: Thanks for the insightful comments regarding my query. I think that Martin and Doug have well assessed my position and both offer useful advice and have greatly improved my limited understanding of S4. Thanks! At this point, i'm well into the app

[R] Weibull point process

2011-09-14 Thread Torbjørn Ergon
Dear list, I'm looking for a function to generate (simulate) a random Weibull point process. Can anyone help? Cheers, Torbjørn Ergon, University of Oslo __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

[R] Nonlinear Regression

2011-09-14 Thread Liam Brown
I'm wondering what packages exist to implement nonlinear least squares regression in R other than 'nls'. Are there packages which implement methods to estimate the optimum values of the parameters which do not use the Gauss-Newton algorithm e.g. use Nelder Mead. In particular, I'd be interested whe

[R] Selecting from matrix and then stacking in array.

2011-09-14 Thread Kristian Lind
Hi, I'm solve a problem where I want to select every 3rd row from a matrix. I do this for all rows in a loop. This gives me "i" matrices with different number of rows I want to stack these matrices in an array and fill the blanks with zero. Does anyone have any suggestions on how to go about thi

Re: [R] draw points in R

2011-09-14 Thread Adel ESSAFI
You are kind 2011/9/14 Sarah Goslee > Hi, > > > The probleme is that when I use points function; it plots to me lines and > > not simpe points. > > That seems unlikely, but without your code it's hard to figure out > what you did. As usual, the telepathy is not working. > > Regardless, this will

[R] Open file for reading and writing with APPEND

2011-09-14 Thread Noah Silverman
Hi, I want to keep a running file of some statistics generated by a running process. Occasionally the process will be restarted. On restart, I'd like to read the last line of the log file to set some variables to their last state, and then continue to append values to that same file. I can'

Re: [R] Moving data in a dataset

2011-09-14 Thread MacQueen, Don
Spreadsheets have cells, but R does not. So you will have to be much more specific. The closest I can come would be like this example: If you want to copy the value that is in the 3rd row, 2nd column to the 4th row, 1st column, then mydata[4,1] <- mydata[3,2] But that is a copy, not a move (and

Re: [R] draw points in R

2011-09-14 Thread Sarah Goslee
Hi, > The probleme is that when I use points function; it plots to me lines and > not simpe points. That seems unlikely, but without your code it's hard to figure out what you did. As usual, the telepathy is not working. Regardless, this will create a graph with points: x <- read.table("clipboa

[R] Non compatible arguments in gvlma

2011-09-14 Thread Jon Toledo
Dear R community I have tried to run the diagnostics on my lm object with the gvlma package but it gave me an error with the following message: Non compatible arguments I post here my syntax and session info: RegreC3<-lm(Classification_Score~C3L+Age_+Dx,data=DB1)summary(RegreC3) library(gvlma)gv

Re: [R] Can 'mosaic' be used with a continuous variable?

2011-09-14 Thread Mauricio Cornejo
Thanks Z. Let me clarify my problem a bit further ... as I don't think I could use cut() or spine() to solve it. I have a data frame with three columns (A, B, Value).  'A' and 'B' are categorical and 'Value' is continuous and non-negative.  I'm looking to make a mosaic plot with A on one axis

Re: [R] Force regression line to a 1:1 relationship

2011-09-14 Thread John Fox
Dear S., On Wed, 14 Sep 2011 16:48:55 +0100 S Ellison wrote: > > From: John Fox [mailto:j...@mcmaster.ca] > > Subject: Re: [R] Force regression line to a 1:1 relationship > > If I understand correctly what you said, a plot of residuals > > against an X is not an added-variable plot. An added v

Re: [R] linear probability model

2011-09-14 Thread Achim Zeileis
On Wed, 14 Sep 2011, Adam Copella wrote: Dear All, Is there any function in R for fitting linear probability model, as my response variable is a uniformly distributed. If it is really uniformly distributed, I think there is not much to model... And in the textbooks I know that discuss the li

[R] draw points in R

2011-09-14 Thread Adel ESSAFI
Hello list I have this file 0 121289479 25 0 212599129 1 0 285254098 21 0 444889848 45 0 469197123 30 0 640007403 82 0 718215617 8 0 758534043 56 0 799706577 4 0 814441385 93 0 843545059 37 0.2 121289479 6 0.2 285254098 3 0.2 444889848 6 0.2

Re: [R] Can 'mosaic' be used with a continuous variable?

2011-09-14 Thread Achim Zeileis
On Wed, 14 Sep 2011, Mauricio Cornejo wrote: Hello, I'm wondering if the 'mosaic' plot of the vcd package (or any other function for that matter) can be used with a continuous variable that should be represented via various categorical variables.? All the documentation I've read lead me to b

[R] Can 'mosaic' be used with a continuous variable?

2011-09-14 Thread Mauricio Cornejo
Hello, I'm wondering if the 'mosaic' plot of the vcd package (or any other function for that matter) can be used with a continuous variable that should be represented via various categorical variables.  All the documentation I've read lead me to believe that it only works with counts of categor

Re: [R] Reading R Code aloud

2011-09-14 Thread Joshua Wiley
I wonder if reading C would work for those extra fussy nights? Of course if it didn't work, I'm not sure how to explain the cries of, "but foo wasn't protected from the garbage collector!" to the neighbors... On Wed, Sep 14, 2011 at 10:03 AM, andrewH wrote: > Dear Clarence-- > > LOL!  R really i

Re: [R] Reading R Code aloud

2011-09-14 Thread andrewH
Dear Clarence-- LOL! R really is good for an amazing range of things! How old are your kids? If no one else points to an intensional sample of "vocal R", maybe you could record one of your more inspired readings. It could be a service to the community of R learners -- and if it catches on, perha

[R] Reshape - novel error message and fail

2011-09-14 Thread dunner
Hello all, thanks in advance for the trouble. I can't reshape a particular large data.frame. I have used melt and reshape successfully before, alot. I am trying to go from wide to long format. I wish to separate The named variables below into 3 groups : SYST, DIAST and PULSE, and have them vary

Re: [R] normalizePath

2011-09-14 Thread Henrik Bengtsson
Looks like it originates from a .libPaths() call, cf. > .libPaths function (new) { if (!missing(new)) { new <- Sys.glob(path.expand(new)) paths <- unique(normalizePath(c(new, .Library.site, .Library), "/")) .lib.loc <<- paths[file.info(paths)$isdir %in% TRUE

[R] substitute games with randomForest::partialPlot

2011-09-14 Thread Johann Hibschman
I'm having trouble calling randomForest::partialPlot programmatically. It tries to use name of the (R) variable as the data column name. Example: library(randomForest) iris.rf <- randomForest(Species ~ ., data=iris, importance=TRUE, proximity=TRUE) partialPlot(iris.rf, iris, Sepal.Width)

[R] pdf font example

2011-09-14 Thread ivo welch
Dear R Users--- [Because I struggled with this for a while, I decided to post it into r-help for the benefit of others, where google search will pick it up. thanks to everyone who made these facilities available.] The task is to use your own opentype otf fonts in an R-created pdf file using the

[R] Confidence interval or p-value for difference in two c-statistics

2011-09-14 Thread Bonnett, Laura
Dear All, Apologies if you have a seen a question like this from me before. I am hoping that if I re-word my question more carefully someone may be able to offer more help than the last time I asked something similar. I am using R 2.9.2 and Windows XP. I am trying to determine if there is a

Re: [R] rJython matrix message

2011-09-14 Thread Ben qant
This fixed it: paste(errs, collapse="\\n") Since that fixed it, I did not try anything else. Thanks again, Ben On Wed, Sep 14, 2011 at 10:04 AM, Steve Lianoglou < mailinglist.honey...@gmail.com> wrote: > Hi, > > On Wed, Sep 14, 2011 at 12:01 PM, Ben qant wrote: > > That fixed it!! Thank you

Re: [R] Loops on data˜1

2011-09-14 Thread R. Michael Weylandt
The first thing I'd mention is that *apply are really just well implemented loops themselves so you can't expect true magic from them. However, if you really want to use them rather than just a regular loop (which seems more intuitive for a request like this), this might do what you are looking for

Re: [R] rJython matrix message

2011-09-14 Thread Steve Lianoglou
Hi, On Wed, Sep 14, 2011 at 12:01 PM, Ben qant wrote: > That fixed it!! Thank you very much! I should have thought of that. For the sake of others that might stumble on this thread -- what fixed it, exactly? Doing paste(errs, collapse="\\n"), or the writing to a tmp file thing, or? -steve > >

Re: [R] S4 vs Reference Classes

2011-09-14 Thread Martin Morgan
On 09/14/2011 06:01 AM, Joseph Park wrote: Thanks Martin. What i'm hoping to do is have a class object, with a member method that can change values of slots in the object, without having to assign values by external assignment to the object. Something like this: setClass ( "Element", representa

Re: [R] rJython matrix message

2011-09-14 Thread Ben qant
That fixed it!! Thank you very much! I should have thought of that. Thanks again, Ben On Wed, Sep 14, 2011 at 9:57 AM, Steve Lianoglou < mailinglist.honey...@gmail.com> wrote: > Hi, > > On Wed, Sep 14, 2011 at 11:44 AM, Ben qant wrote: > > Hello, > > > > I've posted something similar under a d

Re: [R] rJython matrix message

2011-09-14 Thread Steve Lianoglou
Hi, On Wed, Sep 14, 2011 at 11:44 AM, Ben qant wrote: > Hello, > > I've posted something similar under a different subject and never received a > solution. Trying again with (hopefully) a better description. > > Objective: Send a matrix of string data in an email message. The message > must have

Re: [R] I cannot load my workspace in R

2011-09-14 Thread Petr PIKAL
Hi > > Re: [R] I cannot load my workspace in R > > It would not load that file, it just stopped there. I do not believe that. Anyway, if R stop to respond to your commands you probably has broken installation. What type of file you are trying to load? Is it saved workspace? Any errors. Here

Re: [R] Force regression line to a 1:1 relationship

2011-09-14 Thread S Ellison
> From: John Fox [mailto:j...@mcmaster.ca] > Subject: Re: [R] Force regression line to a 1:1 relationship > If I understand correctly what you said, a plot of residuals > against an X is not an added-variable plot. An added variable 'Added variable plots' was the phrase used by the previous pos

[R] rJython matrix message

2011-09-14 Thread Ben qant
Hello, I've posted something similar under a different subject and never received a solution. Trying again with (hopefully) a better description. Objective: Send a matrix of string data in an email message. The message must have authentication and be sent via an R script. I'm almost there! Here

Re: [R] Force regression line to a 1:1 relationship

2011-09-14 Thread Patrick Breheny
On 09/14/2011 11:36 AM, David Winsemius wrote: On Sep 14, 2011, at 10:52 AM, Patrick Breheny wrote: The latter type of plot is called a "partial regression plot" or "added variable plot". They are discussed in any regression textbook, as well as wikipedia and probably dozens of other web sites.

Re: [R] Using R in commercial application

2011-09-14 Thread Marc Schwartz
On Sep 14, 2011, at 8:57 AM, Uwe Ligges wrote: > > > On 14.09.2011 14:36, snassimr wrote: >> >> Hi , >> >> I want to use a several R packages to generate some results and take them to >> my own application. >> >> This application would be sold . >> >> If there is some references can help me

Re: [R] Force regression line to a 1:1 relationship

2011-09-14 Thread David Winsemius
On Sep 14, 2011, at 10:52 AM, Patrick Breheny wrote: On 09/13/2011 04:27 PM, David Winsemius wrote: ... Also should be plotting against fitted() rather than regressors. ... Both types of plots (vs. fitted values and vs. regressors) are very common. The former is a default plot in R, the f

[R] Convert SAS NLMIXED code for zero-inflated gamma regression to R

2011-09-14 Thread a11msp
I'm trying to run a zero-inflated regression for a continuous response variable in R. I'm aware of a gamlss implementation, but I'd really like to try out this algorithm by Dale McLerran that is conceptually a bit more straightforward. Unfortunately, the code is in SAS and I'm not sure how to re-wr

Re: [R] I cannot load my workspace in R

2011-09-14 Thread tttcomrader
It would not load that file, it just stopped there. -- View this message in context: http://r.789695.n4.nabble.com/I-cannot-load-my-workspace-in-R-tp3812774p3813233.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org m

Re: [R] Reading large, non-tabular files

2011-09-14 Thread jim holtman
Here how long it might take in R to do. I created a file of 558MB and then read it it, found lines that had ''76095' in them and then wrote those out: > system.time(x <- readLines('tempyy')) # read in the 558MB file user system elapsed 65.910.82 67.40 > object.size(x) 63348864 bytes

Re: [R] I cannot load my workspace in R

2011-09-14 Thread Uwe Ligges
On 14.09.2011 17:02, Duncan Murdoch wrote: On 14/09/2011 10:30 AM, tttcomrader wrote: Something is wrong with the image, here is the link to it: http://s749.photobucket.com/albums/xx136/liverton11/?action=view¤t=rproblemloading.jpg I don't see any problem there. Are you saying that you can

Re: [R] Force regression line to a 1:1 relationship

2011-09-14 Thread John Fox
Dear S Ellison, If I understand correctly what you said, a plot of residuals against an X is not an added-variable plot. An added variable plot is constructed by regressing Y on all the Xs but the focal X, and regressing the focal X on all the other Xs; then the residuals from the first regress

Re: [R] R 2.13.1 console echo problem or am I wrong?

2011-09-14 Thread ONKELINX, Thierry
Have you tried adding print() statements? e.g. changing summary(x) to print(summary(x)) Posting a reproducible example would help us. Best regards, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek t

Re: [R] Force regression line to a 1:1 relationship

2011-09-14 Thread S Ellison
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Patrick Breheny > > The latter type of plot is called a "partial regression plot" > or "added variable plot". They are discussed in any > regression textbook, as well as wi

Re: [R] I cannot load my workspace in R

2011-09-14 Thread Duncan Murdoch
On 14/09/2011 10:30 AM, tttcomrader wrote: Something is wrong with the image, here is the link to it: http://s749.photobucket.com/albums/xx136/liverton11/?action=view¤t=rproblemloading.jpg I don't see any problem there. Are you saying that you can't type anything more? Duncan Murdoch

Re: [R] linear probability model

2011-09-14 Thread Adam Copella
Many thanks, but I thought the logit link function is for binary response? Regards, Adam 2011/9/14 Jean-Christophe BOUËTTÉ > Hi, > You might want to take a look at ?glm with logit link function. > JC > > 2011/9/14 Adam Copella : > > Dear All, > > > > Is there any function in R for fitting linea

[R] R 2.13.1 console echo problem or am I wrong?

2011-09-14 Thread Levent TERLEMEZ
Dear Users, I was using 2.12.1 version for arules and arulesViz, but I encountered with a problem. I wrote a small function that can also plotting with arulesViz. The problem occures while plotting command is on, console does not print other commands' results. So I put off the plot command, but

Re: [R] I cannot load my workspace in R

2011-09-14 Thread tttcomrader
Something is wrong with the image, here is the link to it: http://s749.photobucket.com/albums/xx136/liverton11/?action=view¤t=rproblemloading.jpg -- View this message in context: http://r.789695.n4.nabble.com/I-cannot-load-my-workspace-in-R-tp3812774p3813040.html Sent from the R help mailing lis

Re: [R] I cannot load my workspace in R

2011-09-14 Thread tttcomrader
The file name is "Homework #3 Hint" I changed the file name to just homework, but it still doesn't work. -- View this message in context: http://r.789695.n4.nabble.com/I-cannot-load-my-workspace-in-R-tp3812774p3812944.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] Force regression line to a 1:1 relationship

2011-09-14 Thread Patrick Breheny
On 09/13/2011 04:27 PM, David Winsemius wrote: ... Also should be plotting against fitted() rather than regressors. ... Both types of plots (vs. fitted values and vs. regressors) are very common. The former is a default plot in R, the first one that appears if you submit: fit <- lm(y~X) pl

Re: [R] linear probability model

2011-09-14 Thread Jean-Christophe BOUËTTÉ
Hi, You might want to take a look at ?glm with logit link function. JC 2011/9/14 Adam Copella : > Dear All, > > Is there any function in R for fitting linear probability model, as my > response variable is a uniformly distributed. > > Regards, > Adam > >        [[alternative HTML version deleted]]

Re: [R] Reading large, non-tabular files

2011-09-14 Thread Rainer Schuermann
That looks like a perfect job for (g)awk which is in every Linux distribution but also available for Windows. It can be called with something like system( "awk -f script.awk inputfile.txt" ) and does its job silently and very fast. 650MB should not be an issue. I'm not proficient in awk but wou

Re: [R] Reading large, non-tabular files

2011-09-14 Thread David Winsemius
On Sep 14, 2011, at 7:08 AM, Stefan McKinnon Høj-Edwards wrote: Dear R-help, I have a very large ascii data file, of which I only want to read in selected lines (e.g. on fourth of the lines); determining which lines depends on the lines content. So far, I have found two approaches for do

Re: [R] Reading large, non-tabular files

2011-09-14 Thread Gabor Grothendieck
2011/9/14 Stefan McKinnon Høj-Edwards : > Dear R-help, > > I have a very large ascii data file, of which I only want to read in selected > lines (e.g. on fourth of the lines); determining which lines depends on the > lines content. So far, I have found two approaches for doing this in R; 1) > Re

Re: [R] Using R in commercial application

2011-09-14 Thread Uwe Ligges
On 14.09.2011 14:36, snassimr wrote: Hi , I want to use a several R packages to generate some results and take them to my own application. This application would be sold . If there is some references can help me : from legal and development point of view , i.e. how to integrate R with Java

Re: [R] Deleting Rows based on Factor and Time Period

2011-09-14 Thread Mikkel Grum
alldat$yearmonth <- substr(alldat$mydate, 1, 7) From: Anna Dunietz Cc: "r-help@r-project.org" Sent: Wednesday, September 14, 2011 3:21 AM Subject: Re: [R] Deleting Rows based on Factor and Time Period Mikkel - thank you so much! That's a great start! I could

Re: [R] I cannot load my workspace in R

2011-09-14 Thread Uwe Ligges
On 14.09.2011 15:06, tttcomrader wrote: I was trying to open up the workspace that you help me write last night, but for some reasons, when I tried to load the saved file R reads "load("C:\\Users\\Documents\\Statistical Learning\\HW1 #3 hint") and then it basically just stopped. I changed the

Re: [R] svytable including NAs

2011-09-14 Thread David Winsemius
On Sep 14, 2011, at 6:12 AM, tesutton wrote: Dear all I'm creating two way tables based on weighted survey data using svytable and want to include NAs as a level in the table. For unweighted data I use table(..., useNA="ifany"). Is there any equivalent for svytable? If not, is there a so

Re: [R] I cannot load my workspace in R

2011-09-14 Thread Richard M. Heiberger
remove the comment from the filename argument and put it outside the statement. "load("C:\\Users\\Documents\\Statistical Learning\\HW1 #3 hint") "load("C:\\Users\\Documents\\Statistical Learning\\HW1") #3 hint On Wed, Sep 14, 2011 at 9:06 AM, tttcomrader wrote: > I was trying to open up the

[R] nls panel estimation

2011-09-14 Thread Tanker Jedi
Hello, I wanted to estimate an equation. First it seemed an easy task, but than I could not solve that the estimate function sums only a part of the data. The equation is (12) from this article http://irps.ucsd.edu/assets/022/8789.pdf Sincerely, József Kelemen _

Re: [R] Question regarding dnorm()

2011-09-14 Thread Vincy Pyne
Dear Sirs, Thanks a lot for your explanation. This was such a hugh conceptual error from my end. I never realized probability and density are two different things. I used to feel I have strated understanding stats a bit. This explanation has changed everything again. Thanks a lot again Mr Ellis

Re: [R] S4 vs Reference Classes

2011-09-14 Thread Joseph Park
Thanks Martin. What i'm hoping to do is have a class object, with a member method that can change values of slots in the object, without having to assign values by external assignment to the object. Something like this: setClass ( "Element", representation ( x = "numeric", y

[R] Using R in commercial application

2011-09-14 Thread snassimr
Hi , I want to use a several R packages to generate some results and take them to my own application. This application would be sold . If there is some references can help me : from legal and development point of view , i.e. how to integrate R with Java or .Net application Thanks -- View this

Re: [R] envfit vector labels with ordiplot3d

2011-09-14 Thread Briony
Hi Gavin, It might help to push the label away from the arrow head a little bit: ordilabel(pl$arrows * 1.075) This seems to do a pretty good job. will push the label 7.5% away from the arrow head. It seems that that is the best we can currently do with vegan. When you say that this is not wha

[R] I cannot load my workspace in R

2011-09-14 Thread tttcomrader
I was trying to open up the workspace that you help me write last night, but for some reasons, when I tried to load the saved file R reads "load("C:\\Users\\Documents\\Statistical Learning\\HW1 #3 hint") and then it basically just stopped. I changed the file to a .Rdata file but the problem still

Re: [R] Reading large, non-tabular files

2011-09-14 Thread jim holtman
What is overkill about reading in a 650MB text file if you have the space? You are going to have to process one way or another. I would use 'readLines' to read it in, and then 'grepl' to determine which lines I want to keep and then delete the rest, and then write the new file out. At this point

Re: [R] normalizePath

2011-09-14 Thread Duncan Murdoch
On 14/09/2011 9:06 AM, boyang zhe wrote: Hi, I update R from 2.10 to 2.13. Then I find soma problems when I start R. Warning infos as below: # Warning message: In normalizePath(c(new, .Library.site, .Library), "/") : path[1]="": No such file or directory R version 2.13.1 (2011-07-08)

[R] normalizePath

2011-09-14 Thread boyang zhe
Hi, I update R from 2.10 to 2.13. Then I find soma problems when I start R. Warning infos as below: # Warning message: In normalizePath(c(new, .Library.site, .Library), "/") : path[1]="": No such file or directory R version 2.13.1 (2011-07-08) Copyright (C) 2011 The R Foundation for Stat

Re: [R] Changelog maintenance

2011-09-14 Thread Duncan Murdoch
On 14/09/2011 5:49 AM, Vikram Bahure wrote: Hi, We are using r-forge to create a package. My query is that how do I keep the change log updated. I created the change log by doing the following: svn log -v>./inst/changelog Also if I do svn add ./inst/changelog and svn commit -m "Adding changelog

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-14 Thread Timothy Bates
>> %<% would extend the vocabulary established by %in%, and work in the same >> situations. >> e.g. >> # only set “flynnEffect" for people known to be under 12, not for people >> where age is NA >> twinData[twinData$Age %<% 12, "flynnEffect"] = FALSE >> Addressing Duncan's point about returni

[R] linear probability model

2011-09-14 Thread Adam Copella
Dear All, Is there any function in R for fitting linear probability model, as my response variable is a uniformly distributed. Regards, Adam [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

  1   2   >