Re: [R] how to combine presence only data sets to one presence/absence table

What about the opposite of this, This has been very helpful for me, but at the same time, I needed the opposite of this to.. ie, this to spl_A spl_B spl_C spcs1 1 1 0 spcs2 1 0 1 spcs3 0 1 1 this spl_A spl_B spl_C spcs1 spcs1spc

Re: [R] Help needed!

On 2011-04-20 14:33, Shuangyan Xiong wrote: Hi everyone, I have a question. Now I am reading the resource code of the package "ssfcov". The resource code is as following. I cannot find the resource code of the function "myss2d" anywhere in the package. Can anyone give me a hint how to find it in

[R] automatic font selection, please help

Hi there, With the helps from this list, I can set specific CJK fonts for character string using text() function. for example: song <- CIDFont("SimSun", "GBK-EUC-H", "GBK", "") postscriptFonts(song = song) postscript("test.ps", height = 7, width =7, family = "Times", fonts = c("song"), hori

[R] Stymied by plyr

Hello, This is my first time trying to use plyr, and I'm getting nowhere. I have teacher ratings data (1:4), on 10 components, by external observers and internal observers, in schools in areas. I want to calculate the percentage of each rating given on each component, by each type of observer, with

Re: [R] GLM output for deviance and loglikelihood

As you mentioned, the deviance does not always reduce to: D = -2(loglikelihood(model)) It does for ungrouped data, such as for binary logistic regression. So let's stick with the original definition. In this case, we need the log-likelihood for the saturated model. x = rnorm(10) y = rpois(10,l

[R] Help with matching rows

Dear Sir, Please excuse my akwardness as I a new to R and computers, but would kindly appreciate help. { a <- sample (1:10,100,replace=T ) b <-sample(10:20,100,replace=T) c <- sample(20:30,100,replace=T) d <- sample(30:40,100,replace=T) e <- sample(40:50,100,replace=T) } d1 <- a d2 <- b d3 <-c d4

[R] override default arguments in nested function

Hi All, I think I already know the answer, but I am hoping I am missing something. I am using function omega from the psych_1.0-96 in R version 2.13.0. I would like to override one of the default arguments of a function that is eventually called (to fix convergence issues), the problem is there

[R] Error running pvals.fnc in R version 2.13.0

Dear R-help: I've been trying to run pvals.fnc in the newest version of R (2.13.0). The function lmer worked fine, but when I tried to use pvals.fnc on the lmer object, I got the following error: "Error in pvals.fnc(elogr.subj.dys.sum.3x3.p, nsim = 1) : trying to get slot "coefs" from an ob

[R] question on automated content analysis program

Hello. I am trying to replicate “Automated Content Analysis” program, which is developped by Gary King. As was instructed in “readme.pdf”, I put the command “undergrad.results <- undergrad(sep = ",")” in R. It gave me the following message: “[1] "python C:/PROGRA~1/R/R-212~1.2/library/ReadMe/m

Re: [R] invoking python from R

On 21/04/11 10:43, 엄기홍 wrote: Dear colleagues, I was wondering whether anyone may give me some help on a python issue in R. I tried to invoke Python from R, but failed. I have set up a system path to call python: "C:\python27" in path file. In addition, I have confirmed python interpreter proper

Re: [R] grid.table + splom: how to nicely align panel entries

Dear Baptiste, *fantastic*, thank you very much, *precisely* what I was looking for! Cheers, Marius On 2011-04-21, at 01:31 , baptiste auguie wrote: > On 21 April 2011 09:54, Marius Hofert wrote: >> Dear Baptiste, >> >> great, many thanks! >> One last thing: Do you know why the gpar(cex=0.1)

Re: [R] random value generation in a row

On Apr 20, 2011, at 7:30 PM, Rujealous wrote: the plan is i have to create a column with random numbers from 1 to 5 forming a vector with 200 numbers. these numbers should have the same proportion, like 40 "1", 40 ''2'', etc. So you really want a random permutation of a fixed set of 200

[R] random value generation in a row

the plan is i have to create a column with random numbers from 1 to 5 forming a vector with 200 numbers. these numbers should have the same proportion, like 40 "1", 40 ''2'', etc. thx, pvm -- View this message in context: http://r.789695.n4.nabble.com/random-value-generation-in-a-row-tp3464503

Re: [R] Sqldf INSERT INTO

On Wed, Apr 20, 2011 at 12:39 PM, new2R wrote: >  Hi, > > I am new to R and trying to migrate from SAS. I am trying to copy data from > one table to another table which have same columns using sqldf. but not > working and showing "NULL" > > I wrote statement as sqldf("INSERT INTO new select * from

Re: [R] taking rows from data.frames in list to form new data.frame?

Dennis, Thanks, the first example works perfectly. > do.call(rbind, lapply(database, function(df) subset(df, Symbol == 'IBM'))) I haven't tried the second, but will look into plyr. Thanks, again, James On Wed, Apr 20, 2011 at 6:36 PM, Dennis Murphy wrote: > Hi: > > Perhaps you're looking f

Re: [R] taking rows from data.frames in list to form new data.frame?

Hi: Perhaps you're looking for subset()? I'm not sure I understand the problem completely, but is do.call(rbind, lapply(database, function(df) subset(df, Symbol == 'IBM'))) or library(plyr) ldply(lapply(database, function(df) subset(df, Symbol == 'IBM'), rbind) in the vicinity of what you're l

Re: [R] grid.table + splom: how to nicely align panel entries

On 21 April 2011 09:54, Marius Hofert wrote: > Dear Baptiste, > > great, many thanks! > One last thing: Do you know why the gpar(cex=0.1) argument is ignored? > Yes – the theme overrides it, you need to include it in the theme.list(). baptiste > Cheers, > > Marius > > library(lattice) > librar

[R] R with X_Trader XTAPI

Is there a way or an existing packages to connect R to X_Trader using XTAPI ? I know there is a way of connecting R to TWS (Interactive broker).. but I was wondering if anyone has done connecting R to X_Trader ?? Thank you. -- View this message in context: http://r.789695.n4.nabble.com/R-with-X-

[R] Help needed!

Hi everyone, I have a question. Now I am reading the resource code of the package "ssfcov". The resource code is as following. I cannot find the resource code of the function "myss2d" anywhere in the package. Can anyone give me a hint how to find it in the package. Thanks a lot!!bv > ssfcov

Re: [R] Can R replicate this data manipulation in SAS?

Oops, I missed the HAART part. Fortunately that translates straightforwardly: n.dat$HAART <- with(n.dat, ifelse((NRTI >= 3 & NNRTI==0 & PI==0) | (NRTI >= 2 & (NNRTI >= 1 | PI >= 1)) | (NRTI == 1 & NNRTI >= 1 & PI >= 1),

[R] invoking python from R

Dear colleagues, I was wondering whether anyone may give me some help on a python issue in R. I tried to invoke Python from R, but failed. I have set up a system path to call python: "C:\python27" in path file. In addition, I have confirmed python interpreter properly. The code I used is typical

[R] multiple imputation

I have missing values from a few subjects due to instrumentation not working. My data set is N=283 data points. For some subjects i have 60 data points missing max. I tried to use Amelia 2 to impute the missing values but i am getting a negative number and i am sure this is wrong because its bi

[R] taking rows from data.frames in list to form new data.frame?

Hi, I am having a problem figuring out how to extract a subset of rows. I have a list with 68 similar data.frames. Each data.frame is 500 rows by 5 columns. I want to take one row from each data.frame based upon the data in a particular column (i.e. it matches a symbol). For example: > str(da

Re: [R] BMA, logistic regression, odds ratio, model reduction etc

I think it's OK. You can also use the Hmisc package's varclus function. Frank 細田弘吉 wrote: > > Dear Prof. Harrel, > > Thank you very much for your quick advice. > I will try rms package. > > Regarding model reduction, is my model 2 method (clustering and recoding > that are blinded to the out

Re: [R] question regarding qmvnorm

On 21/04/11 09:59, li li wrote: Dear all, I wrote the following function previously. It worked fine with the old mvtnorm package. Somehow with the updated package, I got a error message when trying to use the function. I need some help. It is sort of urgent. Can anyone please take a look. Th

Re: [R] succession time series graph

On Apr 20, 2011, at 6:23 PM, David Bird wrote: Dear gracious R community, I would like to produce charts of phytoplankton biomass changes through time. Each species has a line, and the biomass varies in mirror form along the line for each species along the X time axis. Here is an example of wh

Re: [R] question regarding qmvnorm

If you had told us what the error message was, my job would have been easier. But, at least you provied the code, so it was not hard for me to see where the problem was. There is a problem with the strategy used by `qmvnorm' to locate the initial interval in which the quantile is supposed to l

Re: [R] succession time series graph

Hi David, On Wed, Apr 20, 2011 at 6:23 PM, David Bird wrote: > Dear gracious R community, > > I would like to produce charts of phytoplankton biomass changes > through time. Each species has a line, and the biomass varies in > mirror form along the line for each species along the X time axis. > H

[R] succession time series graph

Dear gracious R community, I would like to produce charts of phytoplankton biomass changes through time. Each species has a line, and the biomass varies in mirror form along the line for each species along the X time axis. Here is an example of what I'd like to do: http://www.er.uqam.ca/nobel/r30

Re: [R] Integrate na.rm in own defined functions

On 20/04/11 22:25, vioravis wrote: This should work!! rmse<-function (x){ dquared<-x^2 sum1<-sum(x^2,na.rm=TRUE) rmse<-sqrt((1/length(x))*sum1) rmse} Shouldn't the divisor be the number of non-missing values in x? Rather than the length of x? (Like, e.g. sum(!is.na(x)) ?) cheers,

Re: [R] Rcmdr vs SPSS

Hi, maybe my eyes are not very well, but I can't see any difference between both results up to different rounding policies... Am 20.04.2011 23:36, schrieb Tamas Barjak: > Hy all! > > Excuse me for the inaccurate composition, but I do not speak well in > English. > > I noticed a mistake in Rcmd

Re: [R] Rcmdr vs SPSS

What's the mistake? They look like the same numbers to me. (Although I didn't check them all). Oh, hang on, are you saying that they're different kinds of residuals, but they are the same? This is because SPSS names its residuals wrongly. SPSS has standardized residuals, these are residuals di

[R] question regarding qmvnorm

Dear all, I wrote the following function previously. It worked fine with the old mvtnorm package. Somehow with the updated package, I got a error message when trying to use the function. I need some help. It is sort of urgent. Can anyone please take a look. The function is the following. Thank

Re: [R] grid.table + splom: how to nicely align panel entries

Dear Baptiste, great, many thanks! One last thing: Do you know why the gpar(cex=0.1) argument is ignored? Cheers, Marius library(lattice) library(grid) library(gridExtra) ## function for correct digit alignment align.digits <- function(l){ sp <- strsplit(as.character(l), "\\.") chars

[R] Rcmdr vs SPSS

Hy all! Excuse me for the inaccurate composition, but I do not speak well in English. I noticed a mistake in Rcmdr (?) -- Models menu --- Add observation statistics to data --- Studentized residuals. My output : (Rcmdr !!!) rstudent.RegModel.1 (= *Studentized residuals*) -1.5690952 -0.0697492

Re: [R] Can R replicate this data manipulation in SAS?

I think this is kind of like asking "will your Land Rover make it up my driveway?", but I'll assume the question was asked in all seriousness. Here is one solution: ## Read in test data; dat <- read.table(textConnection("iddrug start stop 1004NRTI 07/24/9501/05/99

Re: [R] Pattern match

Hi: This is a bit of a roundabout approach; I'm sure that folks with regex expertise will trump this in a heartbeat. I modified the last piece of the string a bit to accommodate the approach below. Depending on where the strings have line breaks, you may have some odd '\n' characters inserted. #

Re: [R] Include C++ DLL, error in ...: C symbol name not in load table

On 20/04/2011 4:38 PM, Dirk Eddelbuettel wrote: On 20 April 2011 at 16:24, Duncan Murdoch wrote: | On 20/04/2011 4:06 PM, Sascha Vieweg wrote: |> Hello R experts |> |> I am googling and reading around, however, I can't get it working |> (perhaps because I do not understand much C, however, I'

Re: [R] get cells by the combination of their column and row names

On Apr 20, 2011, at 3:52 PM, yoav baranan wrote: Hi, I have a (correlation) matrix and I want to select a subset of its cells depending on the combination of their column and row names. This illustrates my problem: mtrx <- matrix(c(1,2,3,4,5,6,7,8,9), nrow=3, ncol=3, dimnames = list(c('c13

Re: [R] grid.table + splom: how to nicely align panel entries

Try this, align.digits = function(l) { sp <- strsplit(as.character(l), "\\.") chars <- sapply(sp, function(x) nchar(x)[1]) n = max(chars) - chars l0 = sapply(n, function(x) paste(rep("0", x), collapse="")) labels = sapply(seq_along(sp), function(i) { point <- if(is.na(sp[[i]][2])) NULL else quo

Re: [R] Include C++ DLL, error in ...: C symbol name not in load table

On 20 April 2011 at 16:24, Duncan Murdoch wrote: | On 20/04/2011 4:06 PM, Sascha Vieweg wrote: | > Hello R experts | > | > I am googling and reading around, however, I can't get it working | > (perhaps because I do not understand much C, however, I'll give it | > a try). I am trying to include C++

Re: [R] get cells by the combination of their column and row names

Yoav - Here's one possibility: wh = outer(rownames(mtrx),colnames(mtrx), + function(x,y)substr(x,nchar(x)-2,nchar(x)) == substr(y,nchar(y)-2,nchar(y))) mtrx[wh] [1] 1 6 8 If you knew that all of the row and column names were 4 characters long, it would simplify to mtrx[out

Re: [R] Include C++ DLL, error in ...: C symbol name not in load table

On 20/04/2011 4:06 PM, Sascha Vieweg wrote: Hello R experts I am googling and reading around, however, I can't get it working (perhaps because I do not understand much C, however, I'll give it a try). I am trying to include C++ code into an R routine, where the C++ code looks: #include using na

[R] get cells by the combination of their column and row names

Hi, I have a (correlation) matrix and I want to select a subset of its cells depending on the combination of their column and row names. This illustrates my problem: mtrx <- matrix(c(1,2,3,4,5,6,7,8,9), nrow=3, ncol=3, dimnames = list(c('c132','c432', 'c233'), c('r132','r233', 'r432')))> mt

[R] Include C++ DLL, error in ...: C symbol name not in load table

Hello R experts I am googling and reading around, however, I can't get it working (perhaps because I do not understand much C, however, I'll give it a try). I am trying to include C++ code into an R routine, where the C++ code looks: #include using namespace std; void foo (double* x, double

Re: [R] Error in dimnames(x) for Poisson EWMA model

On Apr 20, 2011, at 3:05 PM, Ashley E. Jochim wrote: I am attempting to run a Poisson EWMA model using Patrick Brandt's source code. I get the following error when I run the code: Error in dimnames(x) <- dn : length of 'dimnames' [1] not equal to array extent Dimnames(x) looks like this:

[R] Error in dimnames(x) for Poisson EWMA model

I am attempting to run a Poisson EWMA model using Patrick Brandt's source code. I get the following error when I run the code: Error in dimnames(x) <- dn : length of 'dimnames' [1] not equal to array extent Dimnames(x) looks like this: [[1]] NULL [[2]] [1] "mip" "div" "nom"

[R] output usage info to stdout

Hi all, I want to run a R code in the batch mode under UNIX system. Inside that code, I have a usage() function to give the hints regarding parameters. For example usage = function() { msg = "R CMD BATCH --save \"--args path=\\\"input_path\\\" file=\\\"input_file\\\" [seed=

[R] avoiding if-then statements for looped chi-square tests

Hi, I am trying to test for pairwise associations between genotypes ( Rows=individuals, Columns =genes, data are up to 4 genotypes per gene, some with 2,3 or 4) where each chisquare comparison is different depending on the genes tested. The test is the observed multilocus (across columns for each

Re: [R] Random Relabelling

Hi: How about y <- rnorm(4000) ymat <- rowMeans(replicate(1000, y[sample(4000)])) hist(ymeans) system.time({y <- rnorm(4000); yy <- rowMeans(replicate(1000, y[sample(4000)]))}) user system elapsed 0.190.030.22 HTH, Dennis On Wed, Apr 20, 2011 at 7:04 AM, kmatthews wrote: > I ha

Re: [R] Random Relabelling

I have a map of Iowa of with 4000 locations. At each location, I have a cancer mortality rate. I need to test my null hypothesis; that the spatial distribution of the mortality rates is random. For this test, I need to establish a spatial reference distribution. My reference distribution will

[R] Pattern match

Hi ALL, I have very simple question regarding pattern matching. Could anyone tell me how to I can use R to retrieve string pattern from text file. for example my file contain following information SpeciesCommon=(Human);SpeciesScientific=(Homo sapiens);ReactiveCentres=(N,C,C,C,+ H,O,C,C,C,C,O,H);

[R] Can R replicate this data manipulation in SAS?

[*** PLEASE NOTE: I am sending this message on behalf of Paul Miller: Paul Miller (to whom this message has also been copied). He has been trying to send it, but it has never got through. Please do not reply to me, but either to the list and/or to Paul at that address ***] =

[R] Sqldf INSERT INTO

Hi, I am new to R and trying to migrate from SAS. I am trying to copy data from one table to another table which have same columns using sqldf. but not working and showing "NULL" I wrote statement as sqldf("INSERT INTO new select * from data") but showing NULL Please help me in this regard. Th

[R] Polynomial question solved.

Dear R People: Here is the solution: > polynomial(c(-2,1))*polynomial(c(-1,1)) 2 - 3*x + x^2 > Sorry for the trouble. Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodg...@gmail.com

[R] polynomial question

Dear R People: Suppose I have (x-1)*(x-2) and would like to produce x^2 - 3x + 2. I have tried > polynomial(c(1,1)*c(2,1)) 2 + x from the polynom library, but as you can see, that doesn't do it. Is there a function to do this please? I'm sure that I've done this before, but I can't remember th

[R] Crouts algorithm

R-listers Quick question for the group. Is there any LU decomposition that makes use of Crout's algorithm in R. I've been looking for it and I really haven't seen it among the R packages. A -- Adrienne Wootten Graduate Research Assistant State Climate Office of North Carolina Department of Ma

Re: [R] Random Relabelling

Kevin, The following follows John's suggestion, but without the loop. It's quick for me. Jeremy Jeremy T. Hetzel Boston University ## Generate sample data n <- 4000 rep <- 1000 rate <- rnorm(n, mean = 15, sd = 2) / 10 # Mortality rates around 15/100k ## Create an empty matrix with ap

Re: [R] Matrix package transpose

On Wed, Apr 20, 2011 at 8:37 AM, Tobias Abenius wrote: > Since I installed R 2.13 I cannot use the transpose method "t" on sparse > matrices inside my package. Outside the package works. Is there something > new that I have to import methods? Can I then import everything from the > Matrix package

Re: [R] Random Relabelling

There is probably a better way to do this but a for loop like this should work. You would just need to change the numbers to yours and then add on the locations = scores  <- 1:5 mydata <- matrix(data=NA, nrow=5, ncol=10) for(i in 1:10) {

Re: [R] 'Record' row values every time the binary value in a collumn changes

Here's one way to do part 1: rr = rle(Table[,'binary']) cc = cumsum(rr$lengths)+1 thestarts = c(1,cc[cc<=nrow(Table)]) theends = cc-1 answer = cbind(Table[thestarts,'Chromosome'],Table[thestarts,'start'],Table[theends,'start'],rr$values) answer [,1] [,2] [,3] [,4] [1,]1 12 18

Re: [R] Two Questions

When running a large number of commands from a script that produces multiple plots it is often best to send the plots to the pdf device (or other system) that you can then page through after it is finished. You could also specify par(ask=TRUE) then you would be prompted before changing the plot

Re: [R] 'Record' row values every time the binary value in a collumn changes

Hi: Here are a couple more options using packages plyr and data.table. The labels in the second part are changed because they didn't make sense in a 2M line file (well, mine may not either, but it's a start). You can always change them to something more pertinent. # Question 1: Table <- data.fram

Re: [R] How can I 'predict' from an nls model with a fit specified for separate groups?

On Apr 20, 2011, at 17:04 , Stuart Rosen wrote: > Following an example on p 111 in 'Nonlinear Regression with R' by Ritz & > Streibig, I have been fitting nls models using square brackets with the > grouping variable inside. In their book is this example, in which 'state' is > a factor indicat

Re: [R] 'Record' row values every time the binary value in acollumn changes

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of jim holtman > Sent: Wednesday, April 20, 2011 9:59 AM > To: baboon2010 > Cc: r-help@r-project.org > Subject: Re: [R] 'Record' row values every time the binary > value in acollumn

[R] Expanding a VCV Matrix

Dear all, I have special task to expand a given VCV matrix, however could not accomplice yet. Let say I have following VCV matrix > mat <- matrix(c(1,2,0,2,5,0.5,0,0.5,3), 3, 3) > colnames(mat) <- rownames(mat) <- paste("variable", 1:3) > mat variable 1 variable 2 variable 3 variable 1

Re: [R] 'Record' row values every time the binary value in a collumn changes

Here is an answer to part 1: > binary<-c(1,1,1,0,0,0,1,1,1,0,0) > Chromosome<-c(1,1,1,1,1,1,2,2,2,2,2) > start<-c(12,17,18,20,25,36,12,15,16,17,19) > Table<-cbind(Chromosome,start,binary) > # determine where the start/end of each group is > # use indices since the size is large > startEnd <- lappl

[R] Fibonacci

Thank-you all!! Very helpful. -- Forwarded message -- From: Bart Joosen Date: 20 April 2011 15:46 Subject: Re: [R] Fibonacci To: r-help@r-project.org Another solution: while (Fibonacci[1] < 500) Fibonacci <- c(sum(Fibonacci[c(1,2)]), Fibonacci) While this adds the sum before

[R] Partial Least Squares Regression independent variables

Hi, I'm running PLSR on the PLS package. I have 507 independent and one dependent variable. model<-plsr(y~x1+x2+x3.., data = mydata, validation = "CV") the problem with this is writing in 507 variable names is not realistic as I run out of line space in the command window. I cannot run the

[R] survexp with weights

Hello, I probably have a syntax error in trying to generate an expected survival curve from a weighted cox model, but I can't see it. I used the help sample code to generate a weighted model, with the addition of a "weights=albumin" argument (I only chose albumin because it had no missing values,

Re: [R] Two Questions

> From: Stephen P Molnar > Subject: [R] Two Questions > To: "R-help" > Received: Wednesday, April 20, 2011, 9:23 AM > 1.        What is a really good book on > R for a nonprogrammer? Have a look at the books listed on the R website. Books by Peter Dalgaard, Phil Spector, Michael Crawley & Joh

[R] Class htest with non-numeric p-values

Hi everyone, For some tests, tables of critical values are readily available while the algorithm to compute those critical values is not. In such cases, only a range for the p-value can be evaluated (e.g. 0.05 < p < 0.1) from the table of critical values and the statistic. Is there anyway to incl

Re: [R] Sweave

Hi, On Wed, Apr 20, 2011 at 12:11 PM, R Heberto Ghezzo, Dr wrote: > Hello, > I never used Sweave before but now I try and got: Have you installed LaTeX on your computer? It doesn't come installed on Windows machines by default. I would suggest MiKTeX http://www.miktex.org/ but people who use Win

Re: [R] Sweave

Hello, I never used Sweave before but now I try and got: > > rnwfile <- system.file("Sweave", "example-1.Rnw", package = "utils") > Sweave(rnwfile) Writing to file example-1.tex Processing code chunks with options ... 1 : echo term verbatim 2 : term verbatim pdf You can now run (pdf)latex on '

Re: [R] Random Relabelling

Can you explain this a bit more. At the moment I don't see what you are trying to achieve. "calculate the mean of the 1000 values at each of the 4000 points" does not seem to make sense. --- On Wed, 4/20/11, kmatthews wrote: > From: kmatthews > Subject: [R] Random Relabelling > To: r-help@r

[R] Test email #2: Please ignore

Having trouble posting. This is a second test email to help determine if the problem has been resolved. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-g

[R] How can I 'predict' from an nls model with a fit specified for separate groups?

Following an example on p 111 in 'Nonlinear Regression with R' by Ritz & Streibig, I have been fitting nls models using square brackets with the grouping variable inside. In their book is this example, in which 'state' is a factor indicating whether a treatment has been used or not: > Puromyci

Re: [R] Simple question about symbols()

Thank you for the answer and sorry about the bad post i'll remember that in the future. By the way, the line code i used to read the data was inv <- read.csv("data.csv", header=TRUE, sep=";") I tried before to use the bg, but for some reason it wasn't working out for me. But now i got it. Th

Re: [R] Unlist command drops all my column names in the first row and adopts NAs

Why not just data.matrix(UN2010) ? --- On Wed, 4/20/11, Haillie wrote: > From: Haillie > Subject: [R] Unlist command drops all my column names in the first row and > adopts NAs > To: r-help@r-project.org > Received: Wednesday, April 20, 2011, 11:05 AM > Hi Everyone, > > I am having trouble t

Re: [R] How to get R plots with FastRweb

Devi, FastRWeb doesn't use Java - it is entirely R based so all you need is a web server with either CGI or PHP. The client it uses is either a C++ client (part of FastRWeb in the cgi-bin directory of the installed package - just copy to you server's cgi-bin) or a PHP client (in Rserve/src/clie

Re: [R] grid.table + splom: how to nicely align panel entries

Dear Baptiste, very nice, indeed! Two minor issues that remain, are: (1) I tried to omit the decimal dot for those numbers that do not have digits after the decimal dot. But somehow it does not work... (2) Do you know how one can decrease the text size for the text appearing in the lo

[R] Unlist command drops all my column names in the first row and adopts NAs

Hi Everyone, I am having trouble turning my data.frame into a matrix format. Because I wanted to change my data.frame with mostly factor variables into a numeric matrix, I used the following code --> UN2010frame<-data.matrix(lapply(UN2010,as.numeric)) However when i checked the mode of the UN201

[R] Matrix package transpose

Hi, Since I installed R 2.13 I cannot use the transpose method "t" on sparse matrices inside my package. Outside the package works. Is there something new that I have to import methods? Can I then import everything from the Matrix package? The problem is that R tries to use t.default which of

Re: [R] Fibonacci

Another solution: while (Fibonacci[1] < 500) Fibonacci <- c(sum(Fibonacci[c(1,2)]), Fibonacci) While this adds the sum before the existing values, the length or tail function or avoided, but even with reordering, its faster (Fibonacci[length(Fibonacci):1]) Best regards Bart -- View this mess

[R] Random Relabelling

I have 4000 observations that I need to randomly relabel 1000 times and then calculate the mean of the 1000 values at each of the 4000 points. Any ideas for where to begin? Thanks Kevin -- View this message in context: http://r.789695.n4.nabble.com/Random-Relabelling-tp3463100p3463100.html Se

Re: [R] BMA, logistic regression, odds ratio, model reduction etc

Dear Prof. Harrel, Thank you very much for your quick advice. I will try rms package. Regarding model reduction, is my model 2 method (clustering and recoding that are blinded to the outcome) permissible? Sincerely, -- KH (11/04/20 22:01), Frank Harrell wrote: Deleting variables is a bad i

Re: [R] Extrapolating data points for individuals who lack them

> From: de...@exeter.ac.uk > To: r-help@r-project.org > Date: Wed, 20 Apr 2011 12:41:29 +0100 > Subject: [R] Extrapolating data points for individuals who lack them > > Hi, > > We have an experiment where individuals responses were measured over 5 days.

Re: [R] Simple question about symbols()

On Apr 20, 2011, at 8:45 AM, murilofm wrote: The link to the csv file is http://www.filedropper.com/data_5 I use the "d" variable to create the radius: radius <- sqrt( inv$d/ pi ) and i tried symbols(inv$a, inv$b, circles=radius, inches=0.35, fg="white", bg="red", xlab="aa", ylab="bb",

Re: [R] Yearly aggregates and matrices

Thanks for clarifying that. Best Gabor Grothendieck wrote: > > On Wed, Apr 20, 2011 at 5:49 AM, mathijsdevaan > wrote: >> As a follow up on this post, I am trying to slightly adjust the solution >> kindly provided by Gabor. However, I am getting some results that I do

[R] Test email: Please disregard

Having trouble posting. Thought it might have something to do with the particular message I'm sending. Sending this message to test that possibility. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the po

Re: [R] Two Questions

On Apr 20, 2011, at 9:23 AM, Stephen P Molnar wrote: Sorry for the somewhat nondescript subject line, but I have two questions: 1.What is a really good book on R for a nonprogrammer? 2. How do I open more than one R Graphics: Device 2(ACTIVE). That what is the R command

Re: [R] Two Questions

On 20/04/2011 9:23 AM, Stephen P Molnar wrote: Sorry for the somewhat nondescript subject line, but I have two questions: 1.What is a really good book on R for a nonprogrammer? Any book that teaches you the basics of programming would be good, it doesn't need to be about R. If you

Re: [R] SPDEP Package, neighbours list for Moran's I on large grid dataset

Laurent Jégou wrote: > Hello list members, i'd like to calculate the Moran I on a large dataset, > a 8640x3432 grid of values. > > When i try to create the neighbours list with the cell2nb function, on > such a scale, R works for several hours and seems to crash. I'm using the last > version (2.1

Re: [R] Simple question about symbols()

The link to the csv file is http://www.filedropper.com/data_5 I use the "d" variable to create the radius: radius <- sqrt( inv$d/ pi ) and i tried symbols(inv$a, inv$b, circles=radius, inches=0.35, fg="white", bg="red", xlab="aa", ylab="bb", col=c("blue","red")[inv$c+1]) Thanks for the

Re: [R] Simple question about symbols()

David Winsemius comcast.net> writes: > > > On Apr 19, 2011, at 10:51 PM, murilofm wrote: > > > Thanks for the answer; I see that col=c("blue","red")[inv$c+1] > > creates a > > vector of "red" and "blue" associated with the binnary c. > > But still I got everything red. > > If you want teste

[R] Two Questions

Sorry for the somewhat nondescript subject line, but I have two questions: 1.What is a really good book on R for a nonprogrammer? 2. How do I open more than one R Graphics: Device 2(ACTIVE). That what is the R command that I can use to keep more than one plot open. I am running

Re: [R] lapply sequence

On 20/04/2011 7:26 AM, Dean Marks wrote: Good day, My question is: Does the lapply function guarantee a particular sequence in which elements are mapped? And, are we guaranteed that lapply will always be sequential (i.e. never map elements in parallel) ? No. The reason I ask is if I use lapp

Re: [R] BMA, logistic regression, odds ratio, model reduction etc

Deleting variables is a bad idea unless you make that a formal part of the BMA so that the attempt to delete variables is penalized for. Instead of BMA I recommend simple penalized maximum likelihood estimation (see the lrm function in the rms package) or pre-modeling data reduction that is blinde

Re: [R] problem reading csv file

Thank you. r-help-boun...@r-project.org napsal dne 19.04.2011 22:00:56: > Not attached. You might succeed if you rename the file with a .txt > extension and re-post. > > Almost surely an encoding issue. We may need your session Info to get > "locale". Maybe you are right. After some further

Re: [R] Yearly aggregates and matrices

On Wed, Apr 20, 2011 at 5:49 AM, mathijsdevaan wrote: > As a follow up on this post, I am trying to slightly adjust the solution > kindly provided by Gabor. However, I am getting some results that I do not > understand. Example: > > # devel version of zoo > install.packages("zoo", repos = "http://

  1   2   >