Re: [R] map and shapefile help

2012-03-31 Thread Michael Sumner
Try the R-Sig-Geo mailing list for a better target community, but if these are shapefiles there is read support in rgdal (or maptools for a simpler alternative) and if these are in long-lat like the maps package then plot(x, add=TRUE) will be most of the way there. If you need to transform either

Re: [R] How to adjust the stack size of R

2012-03-31 Thread Xiaobo Gu
2012/3/31 Uwe Ligges > > > On 30.03.2012 03:16, Xiaobo Gu wrote: > >> Hi, >> >> I got a stack overflow error when training a glm model with a very long >> formula. >> > > I just tried with a formula of length 1000. How long was yours? > Which version of R? Where is the repdroducible example? > >

[R] map and shapefile help

2012-03-31 Thread chuck.01
Hi, I want to use map("state") and have the ecoregion shape (please see link) file projected onto this. ftp://ftp.epa.gov/wed/ecoregions/cec_na/NA_CEC_Eco_Level3.zip ftp://ftp.epa.gov/wed/ecoregions/cec_na/NA_CEC_Eco_Level3.zip Could someone please show me how; I have never messed with this so

[R] Degrees of Freedom for lme.

2012-03-31 Thread SHouston
Hi, I am trying to run a linear mixed effect model on data. I have 17 longitudinal subjects and 36 single subjects, and this is the code I'm using (below). So, INDEX1 is the column with brain volumns, and the predictors are gort and age, by time ID (time they were seen). I believe my data is

[R] Filling empty List in a FOR LOOP

2012-03-31 Thread michaelyb
Hello all, I am trying to automate a list to be filled in a FOR LOOP. Basically I need to load info from a third source and store it into a list. This list goes to object "A"; The next list brought in by the LOOP goe to object "B".. os on ans so forth. When the iteration on the LOOP goes well (wh

Re: [R] simulate correlated binary, categorical and continuous variable

2012-03-31 Thread Greg Snow
Your explanation below has me more confused than before. Now it is possible that it is just me, but it seems that if others understood it then someone else would have given a better answer by now. Are you restricting your categorical and binary variables to be binned versions of underlying normal

[R] How do I get a rough quick utility plot of a time series?

2012-03-31 Thread Hurr
The following is an R run showing that it has read a 2000 point time series having two values at each time. What is the simplest way to plot the data? It seems like the whatever plot routine should be able to make a rough utility plot from just the TS matrix. Note that I am already using spec.p

[R] trouble with small multiples on a date variable

2012-03-31 Thread davideps
Hello, I am loading a DBF file into R via JGR and am having trouble creating small multiple histograms on a date variable. Hist() handles the variable correctly. But I've been unable to work with lattice or ggplot2 despite trying format() and as.POSIXct(). Dates are in the format "2010-05-15" and

Re: [R] Adding text for written comments to bottom of graphs

2012-03-31 Thread baptiste auguie
Hi, On 1 April 2012 03:41, Paul Miller wrote: > Hello Baptiste, > > What you've done is very interesting. Went through and tried to understand > all the steps. Reminded me of studying languages in years gone by. Always > found it easier to read and understand a sentence than to construct a > s

Re: [R] simulate correlated binary, categorical and continuous variable

2012-03-31 Thread David Duffy
Burak Aydin asked: Lets say I know Pearson covariance matrix. When I use rmvnorm to simulate 9 variables and then dichotomize/categorize them, I cant retrieve the population covariance matrix. library(polycor) sim1 <- function(thresh=0.5, r=0.3) { x <- rmvnorm(1000,c(0,0),matrix(c(1,r,r,1)

Re: [R] trouble compiling on Windows

2012-03-31 Thread Prof Brian Ripley
On 31/03/2012 10:38, Joshua Wiley wrote: I may be putting my foot in my mouth, but I want to say that Rtools changed at 2.15 from mingw to gcc. I would check that you have the most recent version of Rtools installed, make sure you have the right directories in your path environment variable (wh

Re: [R] help interpreting aov results

2012-03-31 Thread Rolf Turner
On 01/04/12 06:22, Peter Ehlers wrote: On 2012-03-31 10:09, peter dalgaard wrote: On Mar 31, 2012, at 18:57 , David Winsemius wrote: I some how got idea (after studying) about QR, but here how should I interpret these values? Trying to teach yourself statistics this way is a bit like smashi

Re: [R] subtract a list of vectors from a list of data.frames in an elegant way

2012-03-31 Thread Gabor Grothendieck
On Thu, Mar 29, 2012 at 4:28 AM, Eric Fail wrote: > Dear R experts, > > I've realized that it might not be possible to define a negative SELCET > statement in a SQL call so now I'm looking for the smoothest way to generate > a list of what I would like from my large database by first pulling all

Re: [R] gradient

2012-03-31 Thread David Winsemius
On Mar 31, 2012, at 3:58 PM, mariam behboudi wrote: Hello In matlab we have " gradient(F,h) " where h is a scalar uses h as the spacing between points in each direction. Now I need to use this function in R. and I dont know how can I should define my function that I haave "h" in R?

Re: [R] subtract a list of vectors from a list of data.frames in an elegant way

2012-03-31 Thread Henrique Dallazuanna
Try this: mapply('[', DBquery, mapply(setdiff, lapply(DBquery, names), lookup)) On Sat, Mar 31, 2012 at 5:50 PM, Eric Fail wrote: > a working solution to the problem, > >  a <- DBquery[names(lookup)] > > mother.of.lookup <- list() > for(string in names(a)) { >   a[[string]] <- names(a[[string]])

Re: [R] gradient

2012-03-31 Thread Joshua Wiley
Hi, Have you looked at: ?grad in the numDeriv package? It works very similarly and there are examples. Hope this helps, Josh On Sat, Mar 31, 2012 at 12:58 PM, mariam behboudi wrote: > Hello > In matlab we have "   gradient(F,h)   "  where h is a scalar uses h as the > spacing between points

[R] gradient

2012-03-31 Thread mariam behboudi
Hello In matlab we have " gradient(F,h) " where h is a scalar uses h as the spacing between points in each direction. Now I need to use this function in R. and I dont know how can I should define my function that I haave "h" in R? thanks [[alternative HTML version deleted]]

Re: [R] subtract a list of vectors from a list of data.frames in an elegant way

2012-03-31 Thread Eric Fail
a working solution to the problem, a <- DBquery[names(lookup)] mother.of.lookup <- list() for(string in names(a)) { a[[string]] <- names(a[[string]]) mother.of.lookup[[string]] <- setdiff(a[[string]], lookup[[string]]) } identical(mother.of.lookup, result) It might not be the most elegan

Re: [R] basic subset question of matrix

2012-03-31 Thread jim holtman
I think you want: index <- row.names(myMatrix) %in% list_to_keep[,1] This will give a boolean expression as long as the number of rows in myMatrix; your original only had a list as long as list_to_keep On Sat, Mar 31, 2012 at 12:26 PM, Peter Davidsen wrote: > Dear list, > > I would like to sub

[R] basic subset question of matrix

2012-03-31 Thread Peter Davidsen
Dear list, I would like to subset a large expression matrix based on rownames. That is, I have a list (as a txt-file) with gene names that matches some of the rows in my matrix. I've loaded my matrix as well as gene list using the read.table() command. myMatrix <- read.table("name_of_file.txt",

Re: [R] A introductory question about Zips law (Newbie to statistics)

2012-03-31 Thread Daniel Nordlund
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of ali_protocol > Sent: Friday, March 30, 2012 11:04 PM > To: r-help@r-project.org > Subject: [R] A introductory question about Zips law (Newbie to statistics) > > Hi everyone. > >

Re: [R] help interpreting aov results

2012-03-31 Thread Peter Ehlers
On 2012-03-31 10:09, peter dalgaard wrote: On Mar 31, 2012, at 18:57 , David Winsemius wrote: I some how got idea (after studying) about QR, but here how should I interpret these values? Trying to teach yourself statistics this way is a bit like smashing a watch with a hammer and then tryin

Re: [R] help interpreting aov results

2012-03-31 Thread peter dalgaard
On Mar 31, 2012, at 18:57 , David Winsemius wrote: >> I some how got idea (after studying) about QR, but here how should I >> interpret these values? > > Trying to teach yourself statistics this way is a bit like smashing a watch > with a hammer and then trying to understand how it works by loo

Re: [R] POSIXlt vs POSIXct

2012-03-31 Thread Uwe Ligges
On 29.03.2012 23:44, Steven R Corsi wrote: Thanks for this feedback. I was able to get that example to work, but I have gotten conflicting results with factors converted via POSIXct using the format= and tz= options to bring in excel-type dates and times. Here is what I have tried: dt.factor<-

Re: [R] Error in use of "gwindow" and ".First" function

2012-03-31 Thread Uwe Ligges
On 30.03.2012 05:48, mrzung wrote: Hi, I saved an image and planed to open it whenever I want to load it. However, when I open it, It doesn't work. The code is following, and saved it as an image. .First<-function(){ require(tcltk) require(TeachingDemos) library(gWidge

Re: [R] help interpreting aov results

2012-03-31 Thread David Winsemius
On Mar 31, 2012, at 4:52 AM, Kamakshaiah wrote: Dear Friends, I had performed anova test on certain data frame (Health Care Management) and got results [summary(aov)]. I am new to R and also some extent to Statistics. Can somebody help me how should I interpret these figures. I feel difficulty

Re: [R] Adding text for written comments to bottom of graphs

2012-03-31 Thread Paul Miller
Hello Baptiste, What you've done is very interesting. Went through and tried to understand all the steps. Reminded me of studying languages in years gone by. Always found it easier to read and understand a sentence than to construct a sentence of one's own. This is particularly true when you're

Re: [R] list assignment syntax?

2012-03-31 Thread ivo welch
"what is the problem you are trying to solve?" elegance, ease, and readability in my programs. R has morphed from a data manipulation, graphics, and stats program into my mainstay programming language. most of this has been a huge gain. the addition of the parallel package was another recent bi

Re: [R] unable to move temporary installation

2012-03-31 Thread Uwe Ligges
On 31.03.2012 16:15, Drew Tyre wrote: Hi all, I'm having a strange error that prevents me from installing new packages, or updating packages after reinstalling. The error message is Warning: unable to move temporary installation ‘C:\Program Files\R\R-2.14.2\library\file15045004ac2\sandwich’ to ‘

Re: [R] R Help

2012-03-31 Thread R. Michael Weylandt
This sounds a little bit like homework (since you know the actual model), so I'll just point you in the right direction: I'd do something like: n <- 2^(1:10) y <- 2*n^2 - 5*n + 3 dtf <- data.frame(y = y, n=n) lm( **, data = dtf) The * should be replaced by a "formula" object: to fin

Re: [R] How to adjust the stack size of R

2012-03-31 Thread Uwe Ligges
On 30.03.2012 03:16, Xiaobo Gu wrote: Hi, I got a stack overflow error when training a glm model with a very long formula. I just tried with a formula of length 1000. How long was yours? Which version of R? Where is the repdroducible example? Uwe Ligges Regards, Xiaobo Gu [[al

Re: [R] clear console

2012-03-31 Thread Gabor Grothendieck
On Sat, Mar 31, 2012 at 9:53 AM, Rui Barradas wrote: > Hello, > > Try > > cls2 <- function(lines=25) cat(rep("\n", lines)) > cls2() > > It's simpler, and doesn't need any special package. > (In my system, R 2.14.1/Windows 7, i386 or x64,  '\f' didn't work.) > This may or may not be important but

Re: [R] clear console

2012-03-31 Thread Bert Gunter
Will cat(rep("\n",30)) do? This does not remove any objects from the global environment. That requires ?rm -- Bert On Sat, Mar 31, 2012 at 5:09 AM, mrzung wrote: > hi, > > I use  R - 2.15(32bit), and want to make a code to clear a console. > > Actually, I used to run following code to do that

[R] R Help

2012-03-31 Thread Jaymin Shah
Hi, I have a polynomial of 2n^2-5n+3 and I have my n values going up in powers of 2 i.e. n=2,4,8,16….. I wanted to fit this curve to the function A*n*log2(n) +B*n where A and B are to be found. How would i do this? Thank you Jaymin __ R-help@r-

[R] unable to move temporary installation

2012-03-31 Thread Drew Tyre
Hi all, I'm having a strange error that prevents me from installing new packages, or updating packages after reinstalling. The error message is Warning: unable to move temporary installation ‘C:\Program Files\R\R-2.14.2\library\file15045004ac2\sandwich’ to ‘C:\Program Files\R\R-2.14.2\library\sand

Re: [R] clear console

2012-03-31 Thread Rui Barradas
Hello, Try cls2 <- function(lines=25) cat(rep("\n", lines)) cls2() It's simpler, and doesn't need any special package. (In my system, R 2.14.1/Windows 7, i386 or x64, '\f' didn't work.) Hope this helps, Rui Barradas -- View this message in context: http://r.789695.n4.nabble.com/clear-con

Re: [R] clear console

2012-03-31 Thread Gabor Grothendieck
On Sat, Mar 31, 2012 at 8:09 AM, mrzung wrote: > hi, > > I use  R - 2.15(32bit), and want to make a code to clear a console. > > Actually, I used to run following code to do that but after update the > version of R from 2.14 to 2.15, it doesn't work. > > cls <- function (t) { >        require(RDCO

Re: [R] R warning and error messages taking a long time and generating pop-up windows

2012-03-31 Thread Uwe Ligges
On 30.03.2012 07:18, achernia wrote: Balko, Justin wrote Hi all, thanks for your help in advance. I was running some scripts and suddenly R began taking a very long time to process mundane errors, which I often generate since I am new to R...such as forgetting to add in my quote getOption(w

Re: [R] how to avoid farword slashes in RJSONIO output

2012-03-31 Thread jim holtman
The forward slashes are to escape quotes and the new lines. Here is an example showing that when you output the file with 'cat' you probably get the results you are expecting: > x <- "({\"total\":5,\"results\":[\n {\n \"id\": 1,\n\"name\": \"Monkey D + Luffy\",\n\"position\": \"Captain\",\n\"amb

[R] How can I extract the weights (for feature selection) from linear KSVM?

2012-03-31 Thread Reza Salimi-Khorshidi
Hello all, I need some help with extracting the feature-wise weights from the SVM object that is the output of *ksvm* function in R (from *kernlab* package). For instance, if I use the same data/example as the help page does, i.e., data(spam) index <- sample(1:dim(spam)[1]) spam.train <- spam[in

Re: [R] Not getting correct graphs

2012-03-31 Thread Uwe Ligges
On 31.03.2012 15:03, Bryan Hanson wrote: Attachments don't come through on this list. It would be helpful to know what you mean by "junk graphs". Your coding style is a little hard to follow as it is a mix of -> and<- but I think you definitely have a problem in for (i in refid) which s

Re: [R] Not getting correct graphs

2012-03-31 Thread Bryan Hanson
Attachments don't come through on this list. It would be helpful to know what you mean by "junk graphs". Your coding style is a little hard to follow as it is a mix of -> and <- but I think you definitely have a problem in > for (i in refid) which should probably be for (i in 1:length(refid)

[R] clear console

2012-03-31 Thread mrzung
hi, I use R - 2.15(32bit), and want to make a code to clear a console. Actually, I used to run following code to do that but after update the version of R from 2.14 to 2.15, it doesn't work. cls <- function (t) { require(RDCOMClient) wsh <- COMCreate("Wscript.Shell")

[R] Is this is correct way to calculate AIC value for time series ?

2012-03-31 Thread sagarnikam123
i have attached time series file http://r.789695.n4.nabble.com/file/n4521545/1A1I_A_phi_psi_pot_r_k.txt 1A1I_A_phi_psi_pot_r_k.txt i am not understanding formula() function,but what i read,from that 1 st method:- > t<-read.table(file.choose()) #choosing attached file > lm<-lm(t$V1~ 1) # i used

[R] Not getting correct graphs

2012-03-31 Thread Sri krishna Devarayalu Balanagu
Hi all, Can anybody debug the following programme, as I am getting some Junk graphs in the pdf. Please find the attached raw data file. Thank you Regards Rayalu library(ggplot2) setwd("D:\\General Check list") library(RODBC) conn <- odbcConnectExcel ("Book1.xls") Orange1 <- sqlFetch (conn, "Shee

Re: [R] avoiding expression evaluation when calling a function

2012-03-31 Thread Rui Barradas
Hello, btc1 wrote > > Hey thanks, that worked. So you're right, I'd like to run it for multiple > values of i. As it's written, I'm doing it in a for loop, as > > plotter<-function(i,fram,framvec,obj,form1,form2){ > temp.i<-fram[framvec <=(i*.10),] > plot(form1, data=temp.i, xlim=c(0,1500), yl

Re: [R] Is it there any std pattern in R which show same representation style of ppt found on net

2012-03-31 Thread Michael Bibo
sagarnikam123 gmail.com> writes: > > i m searching on net,material related to R, > whenever i get ppt presentation of some topic ,i found it in speical format > white backgroud & blue bars,having links for further topics > see link > > 1) https://www.rmetrics.org/files/Meielisalp2007/Presentati

Re: [R] Is this is correct way to calculate AIC value for time series ?

2012-03-31 Thread sagarnikam123
i want aic value to find best model from AR & ARMA model i calculate ARMA model via autofit() function from "itsmr" package when i get aic of AR model,then i will compare with ARMA's aicc value is my way logically correct? -- View this message in context: http://r.789695.n4.nabble.com/Is-this

[R] help interpreting aov results

2012-03-31 Thread Kamakshaiah
Dear Friends, I had performed anova test on certain data frame (Health Care Management) and got results [summary(aov)]. I am new to R and also some extent to Statistics. Can somebody help me how should I interpret these figures. I feel difficulty in interpreting values and respective rows and colu

Re: [R] Is it there any std pattern in R which show same representation style of ppt found on net

2012-03-31 Thread Sarah Goslee
Hi, It's not a standard format for R presentations, but you can certainly use it. The blue and white style you referenced is one of the standard Beamer themes. Those people are using LaTeX and Sweave to write their presentations. You'd need to install LaTeX for your OS and do some reading, but

Re: [R] lubridate:ymd_hm and coercion of class POSIXct. Smooth way to restore the date format.

2012-03-31 Thread Henrik Pärn
Good old base R... Thanks Michael for your rapid response. -- Henrik Pärn Centre for Conservation Biology Department of Biology Norwegian University of Science and Technology NO-7491 Trondheim NORWAY Office: +47 735 96084 Mobile: +47 909 89 255 Fax: +47 735 96100

Re: [R] plot points using circles filled half in red and half in blue.

2012-03-31 Thread Jim Lemon
On 03/28/2012 01:49 PM, alan wrote: I want to plot many points and want to use circles. The filling color depends on variable a. if a=1, then not fill if a=2 then fill with red, if a=3 then fill with blue, if a=4, fill half with red and half with blue. Can anyone tell me how to plot the case "a=4

Re: [R] trouble compiling on Windows

2012-03-31 Thread Joshua Wiley
I may be putting my foot in my mouth, but I want to say that Rtools changed at 2.15 from mingw to gcc. I would check that you have the most recent version of Rtools installed, make sure you have the right directories in your path environment variable (which also changed, I believe), and make su

Re: [R] trouble compiling on Windows

2012-03-31 Thread Erin Hodgess
It is the OPENMP switch. I changed it from -fopenmp to blank and all was well. thanks, Erin On Sat, Mar 31, 2012 at 4:38 AM, Joshua Wiley wrote: > I may be putting my foot in my mouth, but I want to say that Rtools changed > at 2.15 from mingw to gcc.  I would check that you have the most rec

[R] trouble compiling on Windows

2012-03-31 Thread Erin Hodgess
Dear R People: I'm trying to compile R-2.15.0 on a 64 bit Windows machine. However, I'm getting the following error: c:\R64\R-2.15.0\src\gnuwin32>set TMPDIR=c:\temp set TMPDIR=c:\temp c:\R64\R-2.15.0\src\gnuwin32>make all recommended make all recommended make[1]: `MkRules' is up to date. make[

[R] A introductory question about Zips law (Newbie to statistics)

2012-03-31 Thread ali_protocol
Hi everyone. Newbie to statistics. I have 40 matrices of ~400 values. how may I determine whether the distribution follows zips law? response <-sample (1:20,400*4, replace= TRUE) Thank you vry much. -- View this message in context: http://r.789695.n4.nabble.com/A-introductory-question-about-

[R] Is it there any std pattern in R which show same representation style of ppt found on net

2012-03-31 Thread sagarnikam123
i m searching on net,material related to R, whenever i get ppt presentation of some topic ,i found it in speical format white backgroud & blue bars,having links for further topics see link 1) https://www.rmetrics.org/files/Meielisalp2007/Presentations/Pfaff.pdf 2)http://grapevine.com.au/~yanchang/