Re: [R] WriteXLS error:Error in get(x, envir = envir) : variable names are limited to 256 bytes

2010-12-10 Thread Marc Schwartz
On Dec 10, 2010, at 4:14 PM, David Winsemius wrote: > > On Dec 10, 2010, at 5:02 PM, David Winsemius wrote: > >> >> On Dec 10, 2010, at 4:39 PM, Patrick McKann wrote: >> >>> Hello all, >>> I don't understand why this won't work. I have entered: >>> >>> WriteXLS(alldata,'test.xls') >> >> I

Re: [R] How to print colorful R output??

2010-12-10 Thread Yihui Xie
Yes (depending on what you mean by "output"), but maybe there is a long way to go. The vignette of the Rd2roxygen package is an example: http://cran.r-project.org/web/packages/Rd2roxygen/vignettes/Rd2roxygen.pdf It makes use of the highlight package and Sweave to output colored code. If the style

[R] randomForest: help with combine() function

2010-12-10 Thread Dennis Duro
I've built two RF objects (RF1 and RF2) and have tried to combine them, but I get the following error: Error in rf$votes + ifelse(is.na(rflist[[i]]$votes), 0, rflist[[i]]$votes) : non-conformable arrays In addition: Warning message: In rf$oob.times + rflist[[i]]$oob.times : longer object lengt

[R] R question: memory usage

2010-12-10 Thread simon lu
Hi I have a large R progress that that is constently running into memory issues, I am trying to rewrite some of the code to make more efficeint However, one thing I have found is that the memory usage shown by gc() is very different from what i see from unix termilnal Garbage collection 2047 = 1

[R] locfit weights not working as expected

2010-12-10 Thread Layla Parast
Hello! I am having a problem understanding what the weights option in the locfit command of the locfit package is doing. I have written a sample program which illustrates the issue (below). The example involves using bootstrap however, that is not my main goal but it illustrates where my problem li

Re: [R] [r] overlap different line in a xyplot (lattice)

2010-12-10 Thread Ben Bolker
Peter Ehlers ucalgary.ca> writes: > > On 2010-12-10 07:04, Francesco Nutini wrote: > > > > dear [R] users, > > is there a way to plot different data (but with the same x-variables) > in the same xyplot window? > > There are already a similar question, but the answer is > not enought explanator

[R] Quantile with discrete types

2010-12-10 Thread Stavros Macrakis
I don't understand why 'quantile' works in this case: > tt <- rep(c('a','b'),c(10,3)) > sapply(0:6/6,function(q) quantile(tt,probs=q,type=1)) 0% 16.7% 33.3% 50% 66.7% 83.3% 100% "a" "a" "a" "a" "a" "b" "b" and also > qua

Re: [R] Stricter read.table?

2010-12-10 Thread Ben Bolker
Stavros Macrakis alum.mit.edu> writes: > > read.table gives idiosyncratic results when the input is formatted > strangely, for example: > > read.table(textConnection( "a'b\nc'd\n"),header=FALSE, fill=TRUE,sep="",quote="'") > => "c'd" "a'b" "c'd" > > > read.table(textConnection( "a'b\nc'd\n

Re: [R] WriteXLS error:Error in get(x, envir = envir) : variable names are limited to 256 bytes

2010-12-10 Thread David Winsemius
On Dec 10, 2010, at 5:02 PM, David Winsemius wrote: On Dec 10, 2010, at 4:39 PM, Patrick McKann wrote: Hello all, I don't understand why this won't work. I have entered: WriteXLS(alldata,'test.xls') I have gotten tripped up by the argument syntax in WriteXLS myself, many times. Please

Re: [R] WriteXLS error:Error in get(x, envir = envir) : variable names are limited to 256 bytes

2010-12-10 Thread David Winsemius
On Dec 10, 2010, at 4:39 PM, Patrick McKann wrote: Hello all, I don't understand why this won't work. I have entered: WriteXLS(alldata,'test.xls') I have gotten tripped up by the argument syntax in WriteXLS myself, many times. Please check the help page for argument names and use them,

[R] How to print colorful R output??

2010-12-10 Thread casperyc
Hi All, I wonder if there is a way to print the R output with COLOR? Not the color plots, but the outputs in the console. Thank. casper -- View this message in context: http://r.789695.n4.nabble.com/How-to-print-colorful-R-output-tp3082750p3082750.html Sent from the R help mailing list archiv

[R] WriteXLS error:Error in get(x, envir = envir) : variable names are limited to 256 bytes

2010-12-10 Thread Patrick McKann
Hello all, I don't understand why this won't work. I have entered: WriteXLS(alldata,'test.xls') and I get this error message: Error in get(x, envir = envir) : variable names are limited to 256 bytes. My variable names are not very long, and are accepted by write.csv. alldata is a list contain

Re: [R] Time Series Row Label

2010-12-10 Thread Gabor Grothendieck
On Fri, Dec 10, 2010 at 3:16 PM, dreadgazebo wrote: > > Simple question... > > I know that when referencing data in a multivariate time series matrix, I > can use the variable name instead of the column number (such as > budget.ts[4,"incometax"]). Is there a way I can use the time unit (say, the >

[R] spatial clusters

2010-12-10 Thread dorina.lazar
Dear all, I am looking for a clustering method usefull to classify the countries in some clusters taking account of: a) the geographical distance (in km) between countries and b) of some macroeconomic indicators (gdp, life expectancy...). Are there some packages in R usefull for this? Thanks

Re: [R] Encoding problem - I fails to read Hebrew text from online

2010-12-10 Thread Tal Galili
Hi Matt and everyone else, Thanks for the help so far. I ended up using the tips provided to create a "dirty hack" based on a translation table between the code and the Hebrew letters. For the future (and for any suggestions), I am attaching this code bellow: Best, Tal # the translation table:

Re: [R] Reorder factor and address embedded escapes

2010-12-10 Thread David Winsemius
On Dec 10, 2010, at 3:14 PM, Rob James wrote: I am trying to reorder a factor variable that has embedded escape characters. The data begins as a csv file with a factor that includes embedded new line characters. By the time read.table has rendered it into a data frame, the variable now h

Re: [R] Reorder factor and address embedded escapes

2010-12-10 Thread Erik Iverson
Does the following help? A = c("A\\nB", "C\\nD") test <-data.frame(A) #access levels directly to change names levels(test$A) <- sub("n", "\n", levels(test$A)) #re-order levels of the factor test$A <- relevel(test$A, "C\nD") Rob James wrote: I am trying to reorder a factor variable that ha

[R] Time Series Row Label

2010-12-10 Thread dreadgazebo
Simple question... I know that when referencing data in a multivariate time series matrix, I can use the variable name instead of the column number (such as budget.ts[4,"incometax"]). Is there a way I can use the time unit (say, the year in an annual time series) instead of the row number? -- Vi

[R] Reorder factor and address embedded escapes

2010-12-10 Thread Rob James
I am trying to reorder a factor variable that has embedded escape characters. The data begins as a csv file with a factor that includes embedded new line characters. By the time read.table has rendered it into a data frame, the variable now has an extra backslash. e.g. "This\nLabel" in the c

Re: [R] survival package - calculating probability to survive a given time

2010-12-10 Thread David Winsemius
On Dec 10, 2010, at 2:07 PM, Andreas Wittmann wrote: Dear R users, i try to calculate the probabilty to survive a given time by using the estimated survival curve by kaplan meier. What is the right way to do that? as far as is see i cannot use the predict-methods from the survival packag

Re: [R] Compare one level of a factor with *all* other non-missing levels

2010-12-10 Thread Peter Ehlers
On 2010-12-10 05:58, deriK2000 wrote: Peter Ehlers wrote: Sounds like you want the Dunnett test procedure which seems to be implemented in a number of packages: multcomp, asd, MCPAN and others. It would probably be a good idea to install package 'sos' and learn how to search with it. Peter

[R] survival package - calculating probability to survive a given time

2010-12-10 Thread Andreas Wittmann
Dear R users, i try to calculate the probabilty to survive a given time by using the estimated survival curve by kaplan meier. What is the right way to do that? as far as is see i cannot use the predict-methods from the survival package? library(survival) set.seed(1) time <- cumsum(rexp(100

Re: [R] Could concurrent R sessions mix up variables?

2010-12-10 Thread Phil Spector
Anthony - I would advise you to use the multicore or snowfall packages to utilize multiple CPUs. As an example using multicore: library(multicore) sim = function(mu)max(replicate(10,max(rnorm(100,mu library(multicore) unlist(mclapply(c(1,5,10,20),sim)) [1] 6.569332 10.268091 15.335

Re: [R] Could concurrent R sessions mix up variables?

2010-12-10 Thread Duncan Murdoch
On 10/12/2010 1:13 PM, Anthony Damico wrote: Hi, I'm working in R 2.11.1 x64 on Windows x86_64-pc-mingw32. I'm experiencing a strange problem in R that I'm not even sure how to begin to fix. I've got a huge (forty-pages printed) simulation written in R that I'd like to run multiple times. When

Re: [R] 45 Degree labels on barplot? Help understanding code previously posted.

2010-12-10 Thread David Winsemius
On Dec 10, 2010, at 10:25 AM, Simon Kiss wrote: Dear colleagues, i found a line or two of code in the help archives from Uwe Ligges about creating slanted x-labels for a barplot and it works well for my purposes (code below). However, I was hoping someone could explain to me precisely wha

Re: [R] help requested

2010-12-10 Thread Petr Savicky
On Fri, Dec 10, 2010 at 07:20:55AM -0800, profaar wrote: > > HI friends, > I have very lengthy graph data in edge list format. I want to convert it > into node list format. > > example: > EDGE LIST FORMAT > 1 2 > 1 3 > 1 4 > 1 5 > 2 3 > 2 4 > 3 2 > 4 1 > 4 3 > 4 5 > 5 2 > 5 4 > > ITS NO

[R] help with RSQLite adding a new column

2010-12-10 Thread Michael D
I'm new to using sql so I'm having difficulties (and worries) in adding a new column of data to a table I have. Its a very large file (around 5 Gb) which is why I'm having to use SQL I have a table with variables ID, IDrec and IDdes and the variables IDrec and IDdes give a mapping of some other va

[R] Could concurrent R sessions mix up variables?

2010-12-10 Thread Anthony Damico
Hi, I'm working in R 2.11.1 x64 on Windows x86_64-pc-mingw32. I'm experiencing a strange problem in R that I'm not even sure how to begin to fix. I've got a huge (forty-pages printed) simulation written in R that I'd like to run multiple times. When I open up R and run it on its own, it works fi

Re: [R] [r] overlap different line in a xyplot (lattice)

2010-12-10 Thread Peter Ehlers
On 2010-12-10 07:04, Francesco Nutini wrote: dear [R] users, is there a way to plot different data (but with the same x-variables) in the same xyplot window? There are already a similar question, but the answer is not enought explanatory... Something like this? x <- rep(1:10, 2) y1 <- rno

Re: [R] Perl "cut" equivalent in R

2010-12-10 Thread Duncan Murdoch
On 10/12/2010 12:05 PM, William Dunlap wrote: > -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Martin Maechler > Sent: Friday, December 10, 2010 8:54 AM > To: Steve Lianoglou > Cc: r-help@r-project.org > Subject: Re: [R]

Re: [R] help requested

2010-12-10 Thread Mike Marchywka
> From: jinyan...@gmail.com > Date: Fri, 10 Dec 2010 17:20:00 +0100 > To: prof...@live.com > CC: r-help@r-project.org > Subject: Re: [R] help requested > > awk '{arr[$1]=arr[$1] " " $2}END{for( i in arr){print i,arr[i]}}' > edgelist.txt | sort -k1

[R] new edition of R Companion to Applied Regression

2010-12-10 Thread John Fox
Dear all, Sandy Weisberg and I would like to announce the publication of the second edition of An R Companion to Applied Regression (Sage, 2011). As is immediately clear, the book now has two authors and S-PLUS is gone from the title (and the book). The R Companion has also been thoroughly rewrit

Re: [R] Perl "cut" equivalent in R

2010-12-10 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Martin Maechler > Sent: Friday, December 10, 2010 8:54 AM > To: Steve Lianoglou > Cc: r-help@r-project.org > Subject: Re: [R] Perl "cut" equivalent in R > > > "SL" == Steve L

Re: [R] Projecting data on a world map using long/lat

2010-12-10 Thread mathijsdevaan
Thanks for the suggestions, but I am not there yet (I'm a real novice). In the code provided by Patrick (see below), I changed the shape input (from sids to world) which I downloaded here: http://thematicmapping.org/downloads/world_borders.php. As a result I also need to change the "CNTY_ID" and "

Re: [R] Perl "cut" equivalent in R

2010-12-10 Thread Martin Maechler
> "SL" == Steve Lianoglou > on Mon, 6 Dec 2010 14:21:59 -0500 writes: >>> if(FALSE) { stuff your don't want executed          } >>> >> > Switching a block of code off/on with editing a single >> character may be done using 0/1 instead of FALSE/TRUE. SL> Or even

Re: [R] Remove 100 years from a date object

2010-12-10 Thread Daniel Brewer
On 10/12/2010 4:17 PM, Barry Rowlingson wrote: > On Fri, Dec 10, 2010 at 3:27 PM, Daniel Brewer > wrote: >> Hello, >> >> I have some data that has dates in the form 27.02.37. I convert them to >> a date object as follows: >> as.Date(data$date,format="%d.%m.%y") >> >> But this gives me years su

[R] Stricter read.table?

2010-12-10 Thread Stavros Macrakis
read.table gives idiosyncratic results when the input is formatted strangely, for example: read.table(textConnection("a'b\nc'd\n"),header=FALSE,fill=TRUE,sep="",quote="'") => "c'd" "a'b" "c'd" read.table(textConnection("a'b\nc'd\nf'\n'\n"),header=FALSE,fill=TRUE,sep="",quote="'") => "f'" "\n

Re: [R] Remove 100 years from a date object

2010-12-10 Thread Gabor Grothendieck
On Fri, Dec 10, 2010 at 10:27 AM, Daniel Brewer wrote: > Hello, > > I have some data that has dates in the form 27.02.37.  I convert them to > a date object as follows: > as.Date(data$date,format="%d.%m.%y") > > But this gives me years such as 2037 when I would like them to be 1937. >  I thought o

Re: [R] Remove 100 years from a date object

2010-12-10 Thread Clint Bowman
There still may be a problem if the dates go back far enough, e.g., 1909. Is '09' 1909 or 2009? No matter what, you have to decide which values need 1900 added and which need 2000. I'd split the date on the delimiter '.', decide whether to add 1900 or 2000, and then paste them together and t

[R] survreg vs. aftreg (eha) - the relationship between fitted coefficients?

2010-12-10 Thread Eleni Rapsomaniki
Dear R-users, I need to use the aftreg function in package 'eha' to estimate failure times for left truncated survival data. Apparently, survreg still cannot fit such models. Both functions should be fitting the accelerated failure time (Weibull) model. However, as Göran Broström points out in

Re: [R] help requested

2010-12-10 Thread Henrique Dallazuanna
Try this: > DF V1 V2 1 1 2 2 1 3 3 1 4 4 1 5 5 2 3 6 2 4 7 3 2 8 4 1 9 4 3 10 4 5 11 5 2 12 5 4 > aggregate(V2 ~ V1, DF, paste, collapse = ' ') V1 V2 1 1 2 3 4 5 2 2 3 4 3 3 2 4 4 1 3 5 5 5 2 4 On Fri, Dec 10, 2010 at 1:20 PM, profa

Re: [R] help requested

2010-12-10 Thread Jinyan Huang
awk '{arr[$1]=arr[$1] " " $2}END{for( i in arr){print i,arr[i]}}' edgelist.txt | sort -k1 On Fri, Dec 10, 2010 at 4:20 PM, profaar wrote: > 1 2 > 1 3 > 1 4 > 1 5 > 2 3 > 2 4 > 3 2 > 4 1 > 4 3 > 4  5 > 5 2 > 5 4 __ R-help@r-project.org mailing list ht

Re: [R] Remove 100 years from a date object

2010-12-10 Thread Barry Rowlingson
On Fri, Dec 10, 2010 at 3:27 PM, Daniel Brewer wrote: > Hello, > > I have some data that has dates in the form 27.02.37.  I convert them to > a date object as follows: > as.Date(data$date,format="%d.%m.%y") > > But this gives me years such as 2037 when I would like them to be 1937. >  I thought of

[R] 45 Degree labels on barplot? Help understanding code previously posted.

2010-12-10 Thread Simon Kiss
Dear colleagues, i found a line or two of code in the help archives from Uwe Ligges about creating slanted x-labels for a barplot and it works well for my purposes (code below). However, I was hoping someone could explain to me precisely what the code is doing. I'm aware it's invoking the text

[R] help requested

2010-12-10 Thread profaar
HI friends, I have very lengthy graph data in edge list format. I want to convert it into node list format. example: EDGE LIST FORMAT 1 2 1 3 1 4 1 5 2 3 2 4 3 2 4 1 4 3 4 5 5 2 5 4 ITS NODE LIST FORMAT SHOULD BE LIKE: 1 2 3 4 5 2 3 4 3 2 4 1 3 5 2 4 Kindly suggest me which package

[R] Remove 100 years from a date object

2010-12-10 Thread Daniel Brewer
Hello, I have some data that has dates in the form 27.02.37. I convert them to a date object as follows: as.Date(data$date,format="%d.%m.%y") But this gives me years such as 2037 when I would like them to be 1937. I thought of trying to take off some time i.e. as.Date(camCD$DoB,format="%d.%m.%y

Re: [R] Textwrangler Languages Folder

2010-12-10 Thread Ben Tupper
Hi, On Dec 10, 2010, at 9:27 AM, Scott Chamberlain wrote: Dear R Community, I recently switched to a Mac (10.6.5), and have installed Textwrangler to run code to R. However, I can't install the syntax highlighting file because I can't find the directory: "~Users/username/Library/Applicati

Re: [R] survival: ridge log-likelihood workaround

2010-12-10 Thread Terry Therneau
-- begin inclusion - Dear all, I need to calculate likelihood ratio test for ridge regression. In February I have reported a bug where coxph returns unpenalized log-likelihood for final beta estimates for ridge coxph regression. In high-dimensional settings ridge regression models usua

Re: [R] (no subject)

2010-12-10 Thread David L Lorenz
Andrija, You should be able to extract the data that you want using a call like this (AD substituted for your c) with(AD, tapply(X2, X1, function(x) sort(x, dec=T)[1:5])) That returns a list like this: $`1` [1] 10 9 8 7 6 $`2` [1] 25 24 23 22 21 Just package it the way that you want. Da

[R] [r] overlap different line in a xyplot (lattice)

2010-12-10 Thread Francesco Nutini
dear [R] users, is there a way to plot different data (but with the same x-variables) in the same xyplot window? There are already a similar question, but the answer is not enought explanatory... Thanks a lot, Francesco [[alternative HTML ve

[R] Textwrangler Languages Folder

2010-12-10 Thread Scott Chamberlain
Dear R Community, I recently switched to a Mac (10.6.5), and have installed Textwrangler to run code to R. However, I can't install the syntax highlighting file because I can't find the directory: "~Users/username/Library/Application Support/TextWrangler/Language Modules/". Is there a different lo

Re: [R] Minimization of the distance

2010-12-10 Thread bluesky
its really help,thanks a lot -- View this message in context: http://r.789695.n4.nabble.com/Minimization-of-the-distance-tp3081345p3081900.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.

Re: [R] ReadWrite.xls problem

2010-12-10 Thread haruo0409
Hans-Peter I have checked the 'library(xlsReadWrite)' startup message. I found that I just failed to 'xls.getshlib()'. Entering 'xls.getshlib()', read.xls() works regularly. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/ReadWrite-xls-problem-tp3078348p3082108.html S

Re: [R] Adding numbers in Outputs

2010-12-10 Thread jim holtman
You should be able to use whatever values you are getting from your script right now. I just did the assignment to match what you were showing on the output. The easiest thing to do is to do 'str(X)' from your data and compare it to the 'x' I created -- str(x). Here is what 'str(x)' gives: > x

[R] melt causes errors when characters and values are used

2010-12-10 Thread Daniel Brewer
Hello, I am finding that the melt function from the reshape library causes errors when applied to a data.frame that contains numeric and character columns. For example, melt(id.vars="ID",data.frame(ID=1:3,date=c("a","b","c"),value=c(1,4,5))) ID variable value 1 1 date a 2 2 date

Re: [R] (no subject)

2010-12-10 Thread jim holtman
try this: > do.call(rbind, lapply(split(x, x$X1), function(.grp){ + .ord <- .grp[order(.grp$X2, decreasing = TRUE),] + .ord[seq(min(5, nrow(.grp))),] + })) X1 X2 1.10 1 10 1.9 1 9 1.8 1 8 1.7 1 7 1.6 1 6 2.25 2 25 2.24 2 24 2.23 2 23 2.22 2 22 2.21 2 21 On Fri, Dec

[R] (no subject)

2010-12-10 Thread andrija djurovic
Hi R-help, I am trying to find a way to select five highest values in data frame according some variable. I will demonstrate: c X1 X2 1 1 1 2 1 2 3 1 3 4 1 4 5 1 5 6 1 6 7 1 7 8 1 8 9 1 9 10 1 10 11 2 11 12 2 12 13 2 13 14 2 14 15 2 15 16 2

Re: [R] Compare one level of a factor with *all* other non-missing levels

2010-12-10 Thread deriK2000
Peter Ehlers wrote: > > > Sounds like you want the Dunnett test procedure which seems > to be implemented in a number of packages: multcomp, asd, MCPAN > and others. > > It would probably be a good idea to install package 'sos' and > learn how to search with it. > > Peter Ehlers > > Thank

Re: [R] Adding numbers in Outputs

2010-12-10 Thread Amelia Vettori
Dear Mr Holtman Sir, Thanks a lot for your great solution. This certainly is helping me achieve what I need to get. However, I shall be hugely thankful to you if you can guide me in one respect. Sir, you have used following commands to assign values to x and y.  > x <- list(40, c(80,160), c(1

Re: [R] Adding numbers in Outputs

2010-12-10 Thread Jinyan Huang
X<-list(40,c(80,160),c(160,80,400)) Y<-list(10,c(10,30),c(5,18,20)) Z<-c(1,2,3) as.data.frame(do.call("rbind",X))->x as.data.frame(do.call("rbind",Y))->y x*y*Z->r r[upper.tri(r)] <- 0 rowSums(r) __ R-help@r-project.org mailing list https://stat.ethz.ch/m

Re: [R] Adding numbers in Outputs

2010-12-10 Thread Peter Ehlers
On 2010-12-10 02:41, Amelia Vettori wrote: two OutputsHello! I am Amelia from Auckland and work for a bank. I am new to R and I have started my venture with R just a couple of weeks back and this is my first mail to R-forum. I need following assistance Suppose my R code generates following ou

Re: [R] Adding numbers in Outputs

2010-12-10 Thread jim holtman
try this: > x <- list(40, c(80,160), c(160,80,400)) > y <- list(10, c(10,30), c(5,18,20)) > z <- c(1,2,3) > mapply(function(a1,a2,a3){ + a3 * sum(a1 * a2) + } + , x + , y + , z + ) [1] 400 11200 30720 On Fri, Dec 10, 2010 at 5:41 AM, Amelia Vettori wrote: > two OutputsHell

Re: [R] Sweave: Setting options with SweaveOpts{} when using driver=RweaveHTML

2010-12-10 Thread Duncan Murdoch
On 10/12/2010 3:40 AM, Søren Højsgaard wrote: When using Sweave in connection with the driver RweaveLatex(), global options can be set with \SweaveOpts{}, e.g. \SweaveOpts{keep.source=T}. Does anybody know if it is possible to set global options in the same way when using Sweave with the driver

Re: [R] Compare one level of a factor with *all* other non-missing levels

2010-12-10 Thread Peter Ehlers
On 2010-12-10 03:24, deriK2000 wrote: Dear list, I try to compare the mean of a variable given a value of a factor with the mean of the same variable for all K-1 other non-missing values of this factor. This procedure I want to repeat for each level of the factor. Having read the recommendatio

[R] Adding numbers in Outputs

2010-12-10 Thread Amelia Vettori
two OutputsHello! I am Amelia from Auckland and work for a bank. I am new to R and I have started my venture with R just a couple of weeks back and this is my first mail to R-forum. I need following assistance Suppose my R code generates following outputs as > X [[1]] [1] 40 [[2]] [1] 80   

Re: [R] Projecting data on a world map using long/lat

2010-12-10 Thread Michael Sumner
Oh, whoops I was looking for the "vote up" button and accidentally hit "Reply All". On Fri, Dec 10, 2010 at 8:25 PM, Barry Rowlingson wrote: > On Fri, Dec 10, 2010 at 2:21 AM, mathijsdevaan > wrote: >> >> Hi, >> >> I have a dataset (CSV) with some counts of firms located around the globe. >>

Re: [R] Using Lagsarlm

2010-12-10 Thread Roger Bivand
This has been answered offlist (the poster also wrote directly to me as package maintainer, but did not post on the R-sig-geo list, as would have seemed natural). The resolution was to read ?formula, and to use either errorsarlm() or lagsarlm() in spdep with formula=y ~ 1. Apparently an insurance

Re: [R] New Installs, Same Trouble Loading doBy and coin Packages

2010-12-10 Thread Peter Ehlers
On 2010-12-10 03:43, Adam Carr wrote: I tried Tal's suggestion of deleting the doBy and coin packages and then reinstalling them from a different mirror. The first install was from the Harvard mirror and the second was from the Case Western Univ. mirror. The new packages generate the same errors

Re: [R] ReadWrite.xls problem

2010-12-10 Thread Hans-Peter Suter
Toby, haruo0409, 2010/12/8 tkdweber : > This is my Error-Message in its German original: > Fehler in .Call("ReadXls", file, colNames, sheet, type, from, rowNames,  : >  Falsche Anzahl von Argumenten (11), erwarte 10 für ReadXls There was a wrong DLL for a short while in the old 1.5.2 version (I f

Re: [R] subset with two factors

2010-12-10 Thread Martin Spindler
Hey Michael, Thank you very much. It works! Best, Martin Original-Nachricht > Datum: Fri, 10 Dec 2010 22:35:56 +1100 > Von: Michael Bedward > An: Martin Spindler > CC: r-help@r-project.org > Betreff: Re: [R] subset with two factors > Hello Martin, > > You were almost ther

[R] Help..Neural Network

2010-12-10 Thread sadanandan
Hi all, I am trying to develop a neural network with single target variable and 5 input variables to predict the importance of input variables using R. I used the packages nnet and RSNNS. But unfortunately I could not interpret the out put properly and the documentation of that packages also not

[R] New Installs, Same Trouble Loading doBy and coin Packages

2010-12-10 Thread Adam Carr
I tried Tal's suggestion of deleting the doBy and coin packages and then reinstalling them from a different mirror. The first install was from the Harvard mirror and the second was from the Case Western Univ. mirror. The new packages generate the same errors when I call them using the library()

Re: [R] [lattice xyplot] Help needed in help in customizing the panel.abline() function

2010-12-10 Thread Girish A.R.
Thanks, Felix! That works. best, -Girish -- View this message in context: http://r.789695.n4.nabble.com/lattice-xyplot-Help-needed-in-help-in-customizing-the-panel-abline-function-tp3079656p3081792.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] subset with two factors

2010-12-10 Thread Michael Bedward
Hello Martin, You were almost there :) T1 <- subset(daten1, Geschlecht=="M" & GG=="A") Hope this helps. Michael On 10 December 2010 22:25, Martin Spindler wrote: > Dear all, > > I have a dataframe of the following strucutre > >  numacc_b coverage_b Geschlecht GG > 1        0          1      

[R] [R-pkgs] pgfSweave 1.1.1 Released

2010-12-10 Thread Cameron Bracken
The next release of pgfSweave is now on CRAN! pgfSweave has seen some significant changes in the past couple of months. The main new features are: - Automatic code highlighting via the highlight package. This can be turned off with the new `highlight` option. - "Tidying" of source code output via

Re: [R] [lattice xyplot] Help needed in help in customizing the panel.abline() function

2010-12-10 Thread Felix Andrews
Hi Girish, Try this: disc <- xyplot(cnt_gt50pct_disc ~ week_num|sku_num, data=DF,type = "h",lwd=2,panel = function(x, y, ...) { panel.abline(v = x[which.max(y)], lty = 2) panel.xyplot(x, y, ...) }) -Felix On 9 December 2010 17:35, Girish A.R. wrote: > > Hi folks, > >

Re: [R] How to enable Arial font for postcript/pdf figure on Windows?

2010-12-10 Thread Camille
Hi Agnes, I converted the Arial font files from ttf to afm using ttf2afm from MikTex complete installation. When used in R with the line recommended by Plos, they seem to give correct Arial font graphics: I checked by opening the ps file with a viewer (gsview), a text editor (notepad++) and Adob

[R] subset with two factors

2010-12-10 Thread Martin Spindler
Dear all, I have a dataframe of the following strucutre numacc_b coverage_b Geschlecht GG 10 1 W A 20 1 M A 30 1 M B 40 1 M B 50 1 W A 60 1

[R] Need help on nnet

2010-12-10 Thread jothy
Hi, Am working on neural network. Below is the coding and the output > library (nnet) > uplift.nn<-nnet (PVU~ConsumerValue+Duration+PromoVolShare,y,size=3) # weights: 16 initial value 4068.052704 final value 3434.194253 converged > summary (uplift.nn) a 3-3-1 network with 16 weights

[R] Compare one level of a factor with *all* other non-missing levels

2010-12-10 Thread deriK2000
Dear list, I try to compare the mean of a variable given a value of a factor with the mean of the same variable for all K-1 other non-missing values of this factor. This procedure I want to repeat for each level of the factor. Having read the recommendations of this list I want to avoid creating

Re: [R] RES: Barplot with "Independent" Lines Y axis

2010-12-10 Thread Jim Lemon
On 12/10/2010 08:48 PM, Rodrigo Aluizio wrote: OK that's it. Working nicely. I sent the final graph with the note of the twoord.plot help attached to it. And some "personal" comments Hi Rodrigo, Although I came rather late to this, I might as well show how I did it: # first I corrected the com

Re: [R] RES: Barplot with "Independent" Lines Y axis

2010-12-10 Thread Rodrigo Aluizio
OK that's it. Working nicely. I sent the final graph with the note of the twoord.plot help attached to it. And some "personal" comments Thank you very much for all the help and remarks. Rodrigo. 2010/12/10 Dennis Murphy > Hi: > > Like Peter Ehlers, I'm not a big fan of multiple response variabl

Re: [R] Minimization of the distance

2010-12-10 Thread Petr Savicky
On Thu, Dec 09, 2010 at 07:43:52PM -0800, bluesky wrote: > > I just contect R,and still learn how to write the code. > I have a problem with argmin sum d(pi,p)/n > for example I have 3 points (a1,b1)(a2,b2)(a3,b3) ,then I want to find > p(x,y) make sure that > (sqrt((a1-x)^2+(b1-y)^2)+sqrt((a2-x

Re: [R] Projecting data on a world map using long/lat

2010-12-10 Thread Barry Rowlingson
On Fri, Dec 10, 2010 at 2:21 AM, mathijsdevaan wrote: > > Hi, > > I have a dataset (CSV) with some counts of firms located around the globe. > Each count is assigned to the longitude and latitude of the specific > location. Now I want to plot these counts on a world map using dots (size of > dots

Re: [R] importing date vector with read.table

2010-12-10 Thread Prof Brian Ripley
On Fri, 10 Dec 2010, Jack Johnson wrote: Hi, I understand this should be an easy task but I am still struggling a bit to read a .txt file with a date vector. My code is as below: data<-read.table(file.choose(),header=TRUE,sep="\t",dec=",",colClasses=c("Date","numeric","numeric","numeric")) B

[R] .por files

2010-12-10 Thread Manderscheid Katharina
hi, i have a spss portable file which i cannot open in R. i tried library(Hmisc) File.spss <- spss.get("File.por", use.value.labels=T) and received: Fehler in read.spss("File.por", : Datei »File.por« ist nicht in einem unterstützten SPSS-Format (in engl.: error, not supported spss format)

[R] importing date vector with read.table

2010-12-10 Thread Jack Johnson
Hi, I understand this should be an easy task but I am still struggling a bit to read a .txt file with a date vector. My code is as below: data<-read.table(file.choose(),header=TRUE,sep="\t",dec=",",colClasses=c("Date","numeric","numeric","numeric")) But I am getting an error: Error in charToDa

[R] Sweave: Setting options with SweaveOpts{} when using driver=RweaveHTML

2010-12-10 Thread Søren Højsgaard
When using Sweave in connection with the driver RweaveLatex(), global options can be set with \SweaveOpts{}, e.g. \SweaveOpts{keep.source=T}. Does anybody know if it is possible to set global options in the same way when using Sweave with the driver RweaveHTML(). Regards Søren [[altern

Re: [R] ReadWrite.xls problem

2010-12-10 Thread peter dalgaard
On Dec 8, 2010, at 15:52 , tkdweber wrote: > > Dear community, > > I have now taken my R-file from lectures and intend to use it at home, but > have a problem > reading the Data from the file. I have installed and loaded the Package > xlsReadWrite so far. xlsReadWrite is a contributed package