Re: [R] quantreg speed

2014-11-16 Thread Roger
You could try method = "pin". Sent from my iPhone > On Nov 16, 2014, at 1:40 AM, Yunqi Zhang wrote: > > Hi William, > > Thank you very much for your reply. > > I did a subsampling to reduce the number of samples to ~1.8 million. It > seems to work fine except for 99th percentile (p-values f

Re: [R] quantile fixed effects with weights

2009-06-18 Thread roger
On Jun 18, 2009, at 5:01 AM, alessia matano wrote: Dear all, I 'm implementing the koenker procedure for quantile fixed effects. Who he? Wuz 'dat? I would like also to apply weights to the procedure, so that to give more weight to the observation that better represent my original sample (

Re: [R] Oja median

2008-09-18 Thread roger
If you had followed the thread of the link that Richard reported you will see an implementation for the general d-dimensional version. Of course this isn't very speedy in higher dimensions, but that is of the nature of the beast, I'm afraid. On Sep 18, 2008, at 4:44 AM, <[EMAIL PROTECTED]>

[R] glm problem

2009-02-14 Thread Roger
Hello, I have a weird problem here. What I want to do is that I need to draw 1000 samples from a matrix, and use glm on them. when I used this command, it runs without the problem > qdata.glm = glm(X258 ~ ., family = binomial, data = q2data[sample(dim(q2data)[1], 1000), ]) but if I drew the sa

Re: [R] glm problem

2009-02-15 Thread Roger
hy this is? Thanks! On Sat, Feb 14, 2009 at 8:09 AM, Duncan Murdoch wrote: > On 13/02/2009 11:59 PM, Roger wrote: > >> Hello, >> >> I have a weird problem here. What I want to do is that I need to draw >> 1000 >> samples from a matrix, and use glm on them. &

[R] matrix transpose

2009-02-15 Thread Roger
This is supposed to be an easy operation, but R 2.8.1 on Mac OS X gives me a lot trouble. the t() function simply does not work properly. What I mean is it works sometimes but does not work at the most of the time, even with the same matrix. this is an example taken from R help > a <- matrix(1:3

[R] feature extraction on time series data

2009-02-15 Thread Roger
Hi, This is a practical question and I am sure there are many statisticians can give me a hand. I have 500 time series data (500 rows), each row contains 100 intervals, i.e., on each row, I have X1, X2, . X100. I am trying to reduce the dimension of this input because the data at the end of

[R] transfer function in R

2009-03-15 Thread Roger
I want to forecast a time series Y using a model that includes previous values of Y and an exogenous time series X using a transfer function. The standard procedure as described in Box and Jenkins and numerous other references is to first fit an ARIMA model to X. Use the ARIMA model to computer res

[R] The R Journal, Volume 9, Issue 2

2018-01-29 Thread Roger Bivand
other supplementary matter from landing pages. Many thanks to all contributors - especially reviewers and authors. Best wishes, Roger -- Roger Bivand Department of Economics, Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; e-mail: roger.biv

Re: [R] Package 'data.table' in version R-3.5.0 not successfully being installed

2018-04-27 Thread Bos, Roger
es all the work for you. Thanks, Roger -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Fox, John Sent: Friday, April 27, 2018 9:21 AM To: Akhilesh Singh Cc: r-help mailing list Subject: Re: [R] Package 'data.table' in version R-3.5.0 not succ

Re: [R] Scaling - does it get any better results than not scaling?

2018-07-17 Thread Roger Koenker
In certain fields this sort of standardization has become customary based on some sort of (misguided) notion that it induces “normality.” For example, in anthropometric studies based on the international Demographic and Health Surveys (DHS) childrens’ heights are often transformed to Z-scores pr

Re: [R] I cannot run R.EXE or RSCRIPT.EXE

2017-05-05 Thread Bos, Roger
Dominik, Folders with spaces can be very tricky sometimes. In my case single quotes didn't work, but double quotes did (see my example below). Adding the full path to your PATH variable with eliminate the need to specify the full path, making it much easier. Thanks, Roger Micr

[R] The R Journal, Volume 9, Issue 1

2017-07-01 Thread Roger Bivand
Dear all, The latest issue of The R Journal is now available at: https://journal.r-project.org/archive/2017-1/. Many thanks to all contributors - especially reviewers and authors. Regards, Bettina -- Roger Bivand Department of Economics, Norwegian School of Economics, Helleveien 30, N-5045

[R] Lattice Histogram Scaling

2017-08-14 Thread Roger Koenker
to think is less silly. Many thanks, Roger require(REBayes) # Source for the velo data require(lattice) x <- velo x <- x[!(x == 0)] bandwidths <- (1:4)*10 lp <- histogram(~x|bw, nint = 100, border = grey(.9), col = grey(.9), type = "density", panel = funct

Re: [R] Lattice Histogram Scaling

2017-08-15 Thread Roger Koenker
My apologies, the data can now be found at: url <- "http://www.econ.uiuc.edu/~roger/research/ebayes/velo.d"; x <- scan(url,skip = 1) If I could get each of the histograms to mimic what is produced by hist(x, 100, freq = FALSE) I’ve experimented with xlim, ylim, without succes

Re: [R] dissolve tiny polygons to others with unionSpatialPolygons{maptools}

2016-06-06 Thread Roger Bivand
o dissolve slivers gets a fair number of hits. There is always also the upstream question of where the slivers came from, and whether the resolution is not in the earlier process - generate a map without slivers that says exactly what you mean, rather than fudging it afterwards. Roger

[R] problems reading XML type file from ishares website

2016-07-28 Thread Bos, Roger
: Style line 14 and Style Thanks in advance, Roger This message and any attachments are for the intended recipient's use only. This message may contain confidential, proprietary or legally privileged information. No right to confidential or privileged treatment of this message is waived

Re: [R] problems reading XML type file from ishares website

2016-07-29 Thread Bos, Roger
e) xml_data <- xmlToList(xmlfile) datadate <- unlist(xml_data[["Worksheet"]][[1]])[1] d <- list() for (ii in 10:length(xml_data[["Worksheet"]])) { rd <- unlist(xml_data[["Worksheet"]][[ii]]) d[[ii]

Re: [R] rgdal installation with two versions of GDAL

2015-04-22 Thread Roger Bivand
until GDAL 2 stabilises? Neither your motivation nor your affiliation are very convincing. Please never post HTML-formatted mail. Roger > > So, the question is, how to I convince R to use the new library search > path? I'm on xubuntu, and R is installed using apt. &

[R] problem setting default timezone

2015-04-23 Thread Bos, Roger
Period %+% ".txt")$mtime, > tz='America/New_York') [1] "2015-04-21" But my question is why is the as.Date function not using the timezone I have set? Thank you in advance, Roger *** This message and any

Re: [R] Plot by FIPS Code using Shapefiles

2015-05-06 Thread Roger Bivand
Corey Sparks utsa.edu> writes: > > Joining data the way you're doing it is dangerous, Roger Bivand and others > describes a standard way to do this process here: > http://r-sig-geo.2731867.n2.nabble.com/Merging-shapefiles-and-csv-td7586839.html Quite right - the chunks C

Re: [R] R convert pdf/png to html

2015-06-03 Thread Bos, Roger
example, but if you want to produce HTML output in the future, it is good to know: http://rmarkdown.rstudio.com/ Thanks, Roger *** This message and any attachments are for the named person's use only. This message may co

Re: [R] Quantile regression model with nonparametric effect and interaction

2015-06-11 Thread Roger Koenker
The main effect trend seems rather dangerous, why not just estimate the f’s in a loop? url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678

Re: [R] Boxplot using a shapefile

2015-06-16 Thread Roger Bivand
p://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example > - ... and please read the posting guide and don't post in HTML. > Definitely! And note that this is a question that is better suited to the R-sig-geo list. Hope this clarifies, Roger > B. > >

[R] as.numeric looses precision

2015-07-02 Thread Bos, Roger
umeric(p, digits=6) [1] 1087.003 I would be nice if I could keep all the decimal places after the conversion. Thanks, Roger *** This message and any attachments are for the intended recipient's use only. This message may

[R] advice on making HTML tables

2015-07-14 Thread Bos, Roger
This might be a little off topic, but I am starting to produce some HTML reports that contain mostly tables and they look great in Chrome but really bad in IE, so I wanted to see if anyone knows of a better way or an easy fix. One option I have used is to convert to PDF, but sometimes it is nic

Re: [R] advice on making HTML tables

2015-07-14 Thread Bos, Roger
Answering my own question, I was able to make the tables look better in IE using some simple CSS: td { padding: 6px; } -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Bos, Roger Sent: Tuesday, July 14, 2015 11:36 AM To: r-help@r-project.org

[R] using dcast with a function of multiple arguments

2015-07-20 Thread Bos, Roger
y expect that to work, but without an example I don't know what else to try. Any hints would be greatly appreciated. Thanks, Roger *** This message and any attachments are for the intended recipient's use only. This message ma

[R] "se.fit" option to the predict.nls() function

2015-01-05 Thread Roger Coppock
When will "R" implement the "se.fit" option to the predict.nls() function? Is there some schedule? __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://

[R] Reading in an XLS (really XML) file from website

2015-02-27 Thread Bos, Roger
e 14 and Style Error: 1: Opening and ending tag mismatch: Style line 14 and Style Thanks in advance for any help or hints, Roger *** This message and any attachments are for the named person's use only. This message may contain co

Re: [R] R crashes when I run rgeos::gDistance

2015-03-03 Thread Roger Bivand
e you avoid relying on upstream packaging. This will however interfere with other software that may depend on your existing GEOS. > > Please let me know if any additional information would be helpful. > > Thank you for your time, > > Adrian > > [[alternative HTML version

Re: [R] evaluating spatial autocorrelation in a raster

2015-03-14 Thread Roger Bivand
To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. Please do provide a self-contained example. Roger ___

Re: [R] Summarizing R script

2018-09-26 Thread Roger Koenker
I use R CMD BATCH foo which produces a file called foo.Rout and provided the script includes sessionInfo() constitutes a quite sufficient summary for my purposes, it isn’t exactly pretty, but it is informative. > On Sep 26, 2018, at 3:00 PM, Spencer Brackett > wrote: > > R users, > > Is any

Re: [R] Siegel nonparametric regression / mblm package

2019-02-11 Thread Roger Koenker
My first thought was also that this was an artifact of the ties, but dithering the data n <- length(child) child <- child + runif(n,-.5,.5) parent <- parent + runif(n,-.5,.5) and rerunning yields the same discrepancy between the Siegel and other fits. Curiously, both lmsreg and ltsreg from MASS

Re: [R] Siegel nonparametric regression / mblm package

2019-02-11 Thread Roger Koenker
regards. > > P.S.: sorry for my bad english... > > Il giorno lun 11 feb 2019 alle ore 12:54 Roger Koenker <mailto:rkoen...@illinois.edu>> ha scritto: > My first thought was also that this was an artifact of the ties, but > dithering the data > n <- length(child) &

Re: [R] particle count probability

2019-02-20 Thread Roger Koenker
Somewhere, buried in the vast literature on the Wicksell problem, there is probably an answer, or at least a hint. > On Feb 20, 2019, at 11:16 AM, PIKAL Petr wrote: > > Dear all > > Sorry, this is probably the most off-topic mail I have ever sent to this help > list. However maybe somebody

Re: [R] Prediction intervals (i.e. not CI of the fit) for monotonic loess curve using bootstrapping

2014-08-13 Thread Roger Koenker
quite straightforward, although the theory for such bands is rather difficult and still under construction. url:www.econ.uiuc.edu/~roger Roger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax:

Re: [R] Making my own graphics device

2014-08-17 Thread Roger Koenker
In ancient times, ie circa 1981, the S language certainly supported HP pen plotters so there should be code somewhere that could be resuscitated, he said naively. url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217

Re: [R] How to use multi paragraph comment like /* and */ in cpp?

2014-09-09 Thread Bos, Roger
What I do in cases like that is wrap the code in a FALSE statement so it won't get executed: its easy to add and remove as needed. If (FALSE) { ...Code to exclude } *** This message and any attachments are for the named person'

Re: [R] Copying tables from R to Excel

2014-09-23 Thread Bos, Roger
Here is a simple method I saw mentioned on this list a few years ago: toExcel <- function(x, tag=FALSE) {write.table(x, "clipboard-128", sep="\t", row.names=tag)} *** This message and any attachments are for the named person's use on

Re: [R] Workaround for RODBC asymmetric numeric data treatment

2014-10-03 Thread Bos, Roger
packages to get fPortfolio to load. Anyone know why install.packages("fPortfolio",repos="http://R-Forge.R-project.org";) can't install all the dependencies automatically? Thanks, Roger > library(RODBC) > library(fPortfolio) Loading required package: timeSeries

[R] better two column output from rmarkdown

2014-11-05 Thread Bos, Roger
Disclaimer: This question is more about HTML than R, but since I am using rmarkdown, I am not sure of a better forum to ask this. The code example below produces two column HTML output. If you run the example you will see that the "" paragram starts off at the bottom of the first column an

Re: [R] Distance matrix for Spatial Filtering {spdep}

2014-11-09 Thread Roger Bivand
odel. Are you sure that the model you are starting from is well-specified (included variables and their functional forms)? If not, you risk including many eigenvectors that simply mop up other misspecifications. This search will be very time-consuming. Consider using the R-sig-geo list, which

Re: [R] Overlay two shapefiles in same map axes rgdal/maptools

2014-11-18 Thread Roger Bivand
Zilefac Elvis yahoo.com> writes: > > Hello, > > I have two spatial map objects (reproducible example further down) > which I would like to overlay in R. The > ESRI shapefiles were read using: > The reproduction of the data didn't work with dput(). > library(rgdal) > Prairie.Boundaries <- r

Re: [R] Running R in Server

2015-08-18 Thread Bos, Roger
port using RMarkdown and host that on the website. If you really need a web app, I suggest looking into Shiny. You can make a shiny app, host it for free on the cloud. You will have a lot to learn, but it will be time well spent. https://www.rstudio.com/ Than

Re: [R] spatial estimation

2015-09-03 Thread Roger Bivand
s(Q ~ P + D, ~ D + F + A, data=Kmenta)) which explains why I asked for a publication providing something which doesn't spring naturally from the estimation method. Roger Bivand > Thanks for your answer, > > Regards > __ R-help@r-project.

[R] Manova: Data similarly generated are significantly different.

2015-09-04 Thread Roger Xu
Dear R users, Hi. I don't know if my understanding of Manova test is correct. So I test with the following code and got strange results. Any help would be appreciated. y0, y1, and y2 are independently generated by the same method. They are each split into 20 groups by the same method. The summa

Re: [R] Function stslshac {sphet}: heteroskedasticity and autocorrelation consistent (HAC) estimator

2015-09-19 Thread Roger Bivand
using a built-in data set, then your intentions would be clearer. ... > I would be grateful for any reaction. > > Monika > PS. Please post plain text, not HTML Roger Bivand __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,

Re: [R] Quantile Regression without intercept

2015-10-05 Thread Roger Koenker
as for lm() or any other linear model fitting…. rq( y ~ x - 1, … ) url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678

Re: [R] Quantile Regression without intercept

2015-10-06 Thread Roger Koenker
> On Oct 6, 2015, at 7:58 AM, Lorenz, David wrote: > > Did you verify that the correct percentages were above/below the regression > lines? I did a quick check and for example did not consistently get 50% of > the observed response values greater than the tau=.5 line. I did when I > included the

Re: [R] Quantile Regression without intercept

2015-10-06 Thread Roger Koenker
do I interpret the line? As an estimate of the conditional median (quantile) function when constrained to pass through the origin… as with least squares fitting without an intercept, you do this at your peril. > > > On Tue, Oct 6, 2015 at 8:03 AM, Roger Koenker wrote: > > >

Re: [R] Quantile Regression without intercept

2015-10-07 Thread Roger Koenker
me implies 0 expenditure, then all (quantile) Engel curves pass through the origin and one might want to impose this. On the other hand maybe not... > From: Roger Koenker > Sent: ‎06-‎10-‎2015 07:09 PM > To: Lorenz, David > Cc: r-help@r-project.org > Subject: Re: [R] Quantile R

Re: [R] Attaching a pdf file to an email generated with sendmailR?

2015-10-09 Thread Bos, Roger
newwd %+% fnameNotes [1] "//rinnycs0051/research/R_HOME_Research/Markdown/morningNotes/morningNotes_20151009.html" > fnameNotes [1] "morningNotes_20151009.html" Thanks, Roger *** This message and any attachm

Re: [R] quantile regression: warning message

2015-10-13 Thread Roger Koenker
see the output from the quantreg FAQ: FAQ() especially point 2. url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678

Re: [R] algorithmic method quantile regression

2015-10-14 Thread Roger Koenker
Did you read item 1 in the quantreg FAQ()? url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Urbana, IL 61801 > On Oct

Re: [R] GIS Analyses for Economics and Marketing

2016-02-23 Thread Roger Bivand
terials are mostly related to environmental questions, but the techniques used are similar. Roger > Can someone help me with some tutorials, training materials or course > examples in area GIS for Marketing or Economics based on QGIS (with Grass > or any other plugins), and some examples of

Re: [R] rgrass7 problem

2016-02-27 Thread Roger Bivand
h more informative - we can guess that you are using GRASS 7 (but not which release) and rgrass7 (not which version), but have no idea of your OS, or of other output of sessionInfo() in R. Do not post HTML-formatted - plain text is preferred. Roger > > Thank you. > > *___

[R] Functional programming?

2016-03-02 Thread Roger Koenker
t;- function(x) {force(i); return(h[[i]] (x))} } x <- 1:99/10 plot(x, h[[1]](x), type = "l") for(i in 2:5) lines(x, h[[i]](x), col = i) Thanks, Roger __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more

Re: [R] Functional programming?

2016-03-02 Thread Roger Koenker
Thanks, Duncan and Bert, Duncan’s version does replicate my result, Bert’s does something a bit different, now I just need some time to digest what you have done, and try to see how and why. Many thanks!!! Roger url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen

Re: [R] Functional programming?

2016-03-02 Thread Roger Koenker
rying to do. As always the knowledge and generosity of the R-help community is inspiring. url:www.econ.uiuc.edu/~roger Roger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244

Re: [R] Map of Europe at NUTS 2 Level

2016-03-11 Thread Roger Bivand
lim=, ylim= to exclude these, or select based on knowing their NUTS2 codes: plot(nuts2a, axes=TRUE, xlim=c(-11.62699, 19.20174), ylim=c(34.35470, 52.04182)) Roger PS. One gets more, and more rapid replies to this kind of question on R-sig-geo. > > This does produce a plot but its rather &#

[R] weighted regression inside FOREACH loop

2016-10-07 Thread Bos, Roger
I have a foreach loop that runs regressions in parallel and works fine, but when I try to add the weights parameter to the regression the coefficients don’t get stored in the “models” variable like they are supposed to. Below is my reproducible example: library(doParallel) cl <- makeCluster(4)

Re: [R] weighted regression inside FOREACH loop

2016-10-07 Thread Bos, Roger
riginal Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Bos, Roger Sent: Friday, October 07, 2016 9:25 AM To: R-help Subject: [R] weighted regression inside FOREACH loop I have a foreach loop that runs regressions in parallel and works fine, but when I try to add the w

Re: [R] weighted regression inside FOREACH loop

2016-10-07 Thread Bos, Roger
. Thanks again for your help! Roger This message and any attachments are for the intended recipient’s use only. This message may contain confidential, proprietary or legally privileged information. No right to confidential or privileged treatment of this message is waived or lost by an error in

Re: [R] getting data from a webpage

2016-12-19 Thread Bos, Roger
Glenn, R Studio did a webinar on Web Scraping using the rvest package that made it look really easy. I haven't gotten around to using it yet, but the video should be on their website somewhere. The link below is the PDF of the slides. It should be education and will probably give you what yo

[R] Re: OS X solution... was Re: rgdal functions: spTransform and CRS not working

2017-01-16 Thread Roger Bivand
X-specific questions like this, you could try R-sig-mac, for spatial data R-sig-geo. Always show the output of sessionInfo() when reporting issues. -- Roger Bivand Department of Economics, Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; e-mail: roge

[R] selecting colors to be used in a plot

2017-01-25 Thread Bos, Roger
My specific question relates to function charts.PerformanceSummary in package PerformanceAnalytics, but the underlying problem is probably with base plot. In this sample code we see a chart with 10 groupings. Apparently 10 groupings is bigger than the number of default colors (8), so red and bl

Re: [R] selecting colors to be used in a plot

2017-01-25 Thread Bos, Roger
Thanks Richard. That worked beautifully! -Original Message- From: Richard M. Heiberger [mailto:r...@temple.edu] Sent: Wednesday, January 25, 2017 12:58 PM To: Bos, Roger Cc: R Help R Subject: Re: [R] selecting colors to be used in a plot This package uses a nonstandard name colorset

Re: [R] Help with a specific quantile regression problem

2017-02-24 Thread Roger Koenker
. url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Urbana, IL 61801 > On Feb 24, 2017, at 5:10 AM, Helio Santana >

[R] R interprets symbol as q()

2008-07-20 Thread Roger Leenders
ol it interprets it as a quit command! So, rather than continuing, R wants to quit on me. How can I make R see the  symbol as just a particular symbol and NOT as the command to q()? Thanks, Roger ps. I hope the symbol is readable in the various email pro

Re: [R] [R-sig-Geo] Spatial Sample

2008-07-27 Thread Roger Bivand
On Sun, 27 Jul 2008, Paul Hiemstra wrote: Hi, One option is: ?spsample - which implements most of Ripley (1981, reprinted 2004, ch.3), and for even more choice, see functions in the spatstat package for point patterns, and the spsurvey package. Roger cheers, Paul Raphael Saldanha

Re: [R] question about krige code in R

2008-07-31 Thread Roger Bivand
mply will not work without a newdata= argument giving the interpolation coordinates and values of your variables X and Y. Roger PS. The R-sig-geo list is probably better suited to this kind of question. > > when I use this one there is this message (with my data setting), because I >

Re: [R] bug in readRAST6 function in package spgrass6

2008-08-01 Thread Roger Bivand
heuristic should be changed to assign an arbitrary NODATA= (outside the range of the data) using information from r.info or other GRASS commands Roger Thanks Rainer -- Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany) Centre of Excellence fo

Re: [R] variance covariance matrix of parameter estimate using nlrq

2008-08-11 Thread roger koenker
mea culpa: I've not written an extractor for this, so you have to do f <- nlrq(whatever) g <- summary(f) g$cov Note that this is computed by resampling so it varies somewhat depending on the seed. url:www.econ.uiuc.edu/~roger Roger Ko

[R] continuous coloring of a polygon

2008-08-15 Thread Roger Leenders
do the coloring. It is easy to give the polygon only one color (e.g. through polygon(c(x1,x2),c(y1,y2), col="red")), but I need a way in which to color the polygon such that the color moves through the color spectrum from red (left) to green (right). Can anyone help me to achieve t

Re: [R] continuous coloring of a polygon

2008-08-15 Thread Roger Leenders
Hi Aiste, Thank you. It may be slow, but it is exactly what I was looking for! I'd be interested if anyone has a faster solution, but this already does the trick well. Of course there will be additional elements in the figure, but it is the coloring I was struggling with. thanks,

[R] How to remove white spaces

2008-08-15 Thread Bos, Roger
7;t return anything, which I found strange. Any advice? Thanks, Roger ** * This message is for the named person's use only. It may contain confidential, proprietary or legally privileged information. No right to confidential o

Re: [R] continuous coloring of a polygon

2008-08-17 Thread Roger Leenders
ce in which I specifically specify the square size of the plot, but the situation remains in which the figure is wider than it is high. Does anyone know how to solve this final piece of the puzzle? Thanks, Roger example code: #pdf(file = ifelse(onefile, "Rplots.pdf", "Rplot%03d.pd

Re: [R] continuous coloring of a polygon

2008-08-17 Thread Roger Leenders
Thanks Barry, that works beautifully! Roger Barry Rowlingson schreef: 2008/8/17 Roger Leenders <[EMAIL PROTECTED]>: WinXP, R2.7.1 Thanks so much to all who have reponded to my inquery. The solutions are most helpful. There is only one final thing that I can't seem to get ri

Re: [R] Fwd: Parsing XML or KML into CSV /Using R for geo coding , OR problem

2008-08-19 Thread Roger Bivand
ikely to use the rgdal KML driver, now also reading KML files. Install the rgdal package, and use the readOGR() function. Then coerce the imported object to a data frame, and export to csv with write.csv(). Roger > > Which R module with deal with an OR problem (like transportation problem

Re: [R] address (nil), cause 'memory not mapped'

2008-08-20 Thread Roger Bivand
he messages so far do not contain information about what is happening, other than that a system() call from the R side, provoking an issue in running a GRASS module, should not bring R down. I feel that the general workflow is not adequate, and that very pos

Re: [R] spatial probit/logit for prediction

2008-08-20 Thread Roger Bivand
y on contiguities. You could consider asking on the R-sig-geo list, and/or reviewing http://leg.ufpr.br/Rcitrus/, which may prove helpful (especially if you know Portuguese). Roger __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/l

[R] How to send Html using SQL Server db mail

2008-08-20 Thread Bos, Roger
;S", " S") gsub(go, "'S", "S") I understand the problem is that R interprets the single quote as a single quote instead of a character that I want to remove, but I don't know how to tell R that I want to identiy the single quote and remove it. This

Re: [R] How many parameters does my model (gls) have?

2008-08-26 Thread roger koenker
?extractAIC url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Champaign, IL 61820 On Aug 26, 2008, at 2:27 PM, Mark Na wrote

Re: [R] Interpolation Problems

2008-09-02 Thread Roger Bivand
a SpatialPointsDataFrame, with an appropriate CRS("+proj=longlat") coordinate reference system, and define a SpatialPixels object for the new grid, you may be able to use idw() in gstat with a very limited search radius to do a simple interpolation on the sphere. Roger

Re: [R] Non-constant variance and non-Gaussian errors

2008-09-03 Thread roger koenker
You might also consider the additive model fitting with rqss() in the package quantreg. url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678

Re: [R] polygon centroids

2009-08-08 Thread Roger Bivand
, readShapeSpatial() in maptools will return centroids for (multi)-polygon objects. You might have got more response on the R-sig-geo list. Roger Daniel Malter wrote: > > Not knowing the package or the function. Could the round() function help > you? Or do you mean decimals shown in a plot?

Re: [R] Map of UK Counties - to use in R

2009-08-12 Thread Roger Bivand
you can retain only England, and choose the NUTS* boundaries that suit your "counties" - which are not presently well-defined because of boundary and administrative changes. The GISCO shapefile is in geographical coordinates, so you'll be able to overplot points by longitude

Re: [R] Map of UK Counties - to use in R

2009-08-13 Thread Roger Bivand
re large numbers <= 360, subtract 360 from all longitude values over 180 (untried, but does happen, so worth checking). If your csv file has extra data for the county entities, you'll need to identify which are which in order to make thematic maps. Roger Raoul wrote: > > Thanks

Re: [R] Browser and Debug?

2009-08-13 Thread roger koenker
At points of total desperation you can always consider the time-honored, avuncular advice -- RTFM, in this case Section 4.4 of Writing R Extensions. url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558

Re: [R] How to use your own data in gstat and sp?

2009-08-17 Thread Roger Bivand
you are not likely to generate a regular grid in this way. Did you examine the class(), summary() and str() of the input data in the examples? Hope this helps, Roger Bivand PS. consider following this up on the R-sig-geo list. Paul Heinrich Dietrich wrote: > > This seems pretty bas

Re: [R] R2 for SAR and validation

2009-09-01 Thread Roger Bivand
implementations (GeoDa, Matlab spatial econometrics toolbox). It is not obvious what you mean by "validate" in this context. Roger Bivand saniye keser wrote: > > > Hi everyone, > How can i obtain R2 for SAR model? and how can i validate the results, can > i use the coeffici

[R] output layout question

2009-09-08 Thread Roger Leenders
value that is below a given threshold should be printed as empty (thus leaving out the low values from sight) - cells with particular values should be printed in bold face or with a larger font size Is this possible? How would one do such a thing? thanks, Roger

[R] Predicting the Present with R at Google

2009-09-18 Thread Roger Koenker
I thought some R-help readers might enjoy seeing the paper by Hal Varian on predicting the present using R and Google Trends that is linked via the following blog comment: http://googleresearch.blogspot.com/2009/04/predicting-present-with-google-trends.html Apologies in advance if this has alrea

Re: [R] Underlay Mapoutlines with Googlemapspictures

2010-02-19 Thread Roger Bivand
nearest CRAN mirror to load image tiles. Roger Bivand - Roger Bivand Economic Geography Section Department of Economics Norwegian School of Economics and Business Administration Helleveien 30 N-5045 Bergen, Norway -- View this message in context: http://n4.nabble.com/Date-label-lost-while

[R] lty dots pdf issue

2010-02-19 Thread Roger Koenker
t is the setting. I'm viewing all this with Acrobat 9 pro, but my printer produces something quite consistent with what is viewable on the screen. Apologies in advance if this has an obvious resolution, I didn't see anything in the r-help archive. Roger url:www.econ.uiu

Re: [R] map2poly - map lat/long cannot be unconstrained?

2010-03-16 Thread Roger Bivand
() call has set the plotting window and axes. Hope this helps, Roger Euan Reavie wrote: > > Using the maptools function "map2poly" I have created a map/polylist > object > - a map of the Great Lakes. My goal is to passively plot additional data > on > this map. Unfortuna

Re: [R] GeoXp package

2009-06-26 Thread Roger Bivand
ing R sees LD_LIBRARY_PATH. This is a very detailed question, more suited to GRASS or GDAL lists, or R-sig-geo. You'll see when you have the solution because the GDAL utilities will work at the shell, for example gdalinfo --formats - I think that it doesn't at the moment. Roger Bivand

[R] beadarray package

2009-06-30 Thread Vallejo, Roger
t1' not found In addition: Warning message: In readIllumina(arrayNames = NULL, useImages = TRUE, singleChannel = FALSE, : No annotation package was specified. Need to use SetAnnotation later ********** I will apreciate help to troublesh

Re: [R] odd behaviour in quantreg::rq

2009-06-30 Thread roger koenker
something fishy about the data construction that would have allowed us to anticipate this? Regarding the "fn" option, and the non-uniqueness warning, this is covered in the (admittedly obscure) faq on quantile regression available at: http://www.econ.uiuc.edu/~roger/resea

  1   2   3   4   5   6   >