Re: [R] Help on predict.lm

2012-03-28 Thread Nederjaard
Hello all, Thanks for all your replies. I have studied on it some more in the meantime, and found indeed out that what I was trying to do was not correct to begin with. Sorry to have wasted your time, but thanks for the comments. -- View this message in context:

Re: [R] Get data from Access 2010 database

2012-03-28 Thread Prof Brian Ripley
On 28/03/2012 06:28, Jeff Newmiller wrote: AFAIK you need to use the 32bit version of R for getting the data. Then you can save it into a more versatile format and re-read it if you really need to run R in 64 bit mode. Not quite: you need to run R of the same architecture as your ODBC

[R] Prediction using dynlm

2012-03-28 Thread jpm miao
Hello, I intend to perform FORECASTING with the dynlm package since I have a bunch of time series data, and my model is a little deviating from the standard AR. Therefore my first attempt is to perform AR forecasting based on dynlm package. First of all, I generate an ARIMA model and then

[R] reading files from two folders

2012-03-28 Thread uday
I would like to read data from two different folder and then combine this together the code which I have tried are as follows setwd(/Groups/data_first/) file_was - list.files(path = ., pattern = v2.0.2.was, all.files = FALSE, full.names = FALSE, recursive = FALSE, ignore.case = FALSE)

[R] Getting a complete vector of Theta estimates from Package LTM

2012-03-28 Thread Jeff
I'm new to R and this list. I'm using the LTM package for the first time to estimate Rasch and grm models. I've searched the web and the list archives and have found an almost identical question to my current one, but I haven't found an answer. I'm going to repeat the question

Re: [R] Exporting a data.frame to excel using sqlSave - adds a character ' to values

2012-03-28 Thread Juliette Fabre
Hello, Sorry for not providing my data: dput(data_col) structure(list(Parameter_Name = c(Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium, Cd Cadmium), Unit = c(ppb , ppb ,

[R] FlexMix for clustering of gene expression data

2012-03-28 Thread Tyler Funnell
Hello, I'm trying to cluster gene expression data as in http://bioinformatics.oxfordjournals.org/content/28/2/222.full It looks like they are using some combination of FlexMix and mgcv but I haven't been able to figure out the exact steps necessary to perform the clustering. Has anyone else

[R] Problems with R Commander version 1.8-3

2012-03-28 Thread petretta
Dear all, I use R 2.14.2 for Windows XP I have no problem with R Commander version 1.7-2, but now I have a problem with R Commander (Rcmdr) version 1.8-3. After I lunch R Commander (sdi as well standard) if I try to import data (i. from Stata datafile) I have this Warning message: Error

Re: [R] Exporting a data.frame to excel using sqlSave - adds a character ' to values

2012-03-28 Thread Ivan Calandra
Hi, I cannot test with RODBC because it doesn't work on my Mac, but if RODBC is not necessary to you, you could do the same using xlsx::write.xlsx(), or xlsReadWrite::wirte.xls(). For example: library(xlsx) write.xlsx(your_df, file=file.xlsx) HTH, Ivan -- Ivan CALANDRA Université de

Re: [R] reading files from two folders

2012-03-28 Thread Gerrit Eichner
Hello, uday, there's presumably a typo in your code because you use path1t in data2 - read.table(paste(path1t, file_wasaux2[i],header=TRUE)) and not path which you defined above. Hth -- Gerrit - Dr. Gerrit Eichner

Re: [R] Test Normality

2012-03-28 Thread Eik Vettorazzi
Hi Sindy, you might try Snows penultimate normality test from the TeachingDemos package. But read the help file carefully. http://www.inside-r.org/packages/cran/TeachingDemos/docs/SnowsPenultimateNormalityTest cheers. Am 28.03.2012 02:32, schrieb Sindy Carolina Lizarazo: Good Night I made

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

2012-03-28 Thread Rolf Turner
On 28/03/12 15:49, 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] SVM. How to use categorical attributes?

2012-03-28 Thread Alekseiy Beloshitskiy
Thank you, Steve, I was thinking about smth like this. Just not sure about the efficiency of using several thousands of additional variables. And the second problem will be time-consumption for managing all these data in memory. Here I posted more brief description:

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

2012-03-28 Thread Eik Vettorazzi
Hi Alan, on an UTF-8 locale you can use \u25D6 and \u25D7, but you have to plot group a=4 twice. On Windows in a non-utf locale you can use a special Windows font: windowsFonts(wdg2=windowsFont(Wingdings 2)) plot(0:1,0:1) strh-strheight(x)/2.02 points(x = 0.5, y =

[R] Nested brew call yields Error in .brew.cat(26, 28) : unused argument(s) (26, 28)

2012-03-28 Thread Chris Beeley
I am writing several webpages using the brew package and R2HTML. I would like to work off one script so I am using nested brew calls. The documentation for brew states that: NOTE: brew calls can be nested and rely on placing a function named ’.brew.cat’ in the environment in which it is

[R] Help on reading data.frame in RJSONIO

2012-03-28 Thread MLSC
Hi friends, I was trying to read a MySQL table and then converting to data.frame to JSON format. Following is the code library(RJSONIO); con - dbConnect(MySQL, host=localhost, user =root, dbname = myDB) rs1 - dbGetQuery(con, select * from strawhats;) print(rs1) json_rs - toJSON(rs1);

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

2012-03-28 Thread Barry Rowlingson
On Wed, Mar 28, 2012 at 11:02 AM, Eik Vettorazzi e.vettora...@uke.de wrote: Hi Alan, on an UTF-8 locale you can use \u25D6 and \u25D7, but you have to plot group a=4 twice. The problem with that is that by default the two half-circles will be centred on the point location, so won't join look

Re: [R] Urgent - I really need some help lme4 model avg Estimates

2012-03-28 Thread Mitchell Maltenfort
You said it right in your first letter. You want to keep it at a level that is comprehensible. Not just to you, but your colleagues, reviewers, readers... Remember the sage wisdom that all models are wrong, but some of them are useful. Rather than try to forge the Excalibur of statistical

[R] Transferring a dataset from WINXP to Mac - problems with umlauts

2012-03-28 Thread Frank Bloos
Dear R-Help, is it possible to correctly display umlauts in a dataset, which has been created under WINXP, on a Mac-system? This is the problem: I have imported a German MS-Access database into R 2.14.1 by using RODBC on a WINXP-System. It runs with a Latin-1 encoding (codepage 1252). However,

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

2012-03-28 Thread Gabor Grothendieck
On Tue, Mar 27, 2012 at 10:49 PM, alan alan.wu2...@gmail.com 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

[R] Spatstat users.

2012-03-28 Thread Rolf Turner
This is a cry for help. My apologies for taking up bandwidth with an issue that is not really on topic. But I really do want to acquire the requested information. In the course of preparing my PBRF (performance based research funds) portfolio for the upcoming PBRF exercise here in New

Re: [R] SVM. How to use categorical attributes?

2012-03-28 Thread Ulrich Bodenhofer
Alex, To avoid the memory issue, you can directly use a bag of words kernel (which corresponds to using the linear kernel on the sparse bag of words matrix Steve suggested). Just a little toy example how this is done for two : x1 - c(how, to, grow, tree) x2 - c(where, to, go, weekend, cinema)

Re: [R] SVM. How to use categorical attributes?

2012-03-28 Thread Ulrich Bodenhofer
Sorry, I forgot to mention the following: all I wrote is only valid as long as your number of samples is smaller than the number of different words. If the number of samples exceeds the total number of different words, you should better use the explicit matrix representation and use some kernel

Re: [R] R extract parts

2012-03-28 Thread MSousa
Thank you, Works perfectly. When the column sector, is the value of lines in the example has six sectors -- View this message in context: http://r.789695.n4.nabble.com/R-extract-parts-tp4509042p4511689.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] matrix(unlist(strsplit())) 'missing value' issue

2012-03-28 Thread MaartenJacobs
Nobody any solution for my problem?? -- View this message in context: http://r.789695.n4.nabble.com/matrix-unlist-strsplit-missing-value-issue-tp4509065p4511668.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

Re: [R] reading files from two folders

2012-03-28 Thread Roland Rau
Hi, On 03/28/2012 08:48 AM, uday wrote: I would like to read data from two different folder and then combine this together the code which I have tried are as follows setwd(/Groups/data_first/) file_was - list.files(path = ., pattern = v2.0.2.was, all.files = FALSE, full.names = FALSE,

[R] Generate distribution from asymmetrical confidence interval and mean

2012-03-28 Thread Andreas Emanuelsson
Hi everyone, does anyone know a god way or package to generate a distribution satisfying a given Mean + asymmetrical Upper and Lower 95% confidence interval? Approximately lognormal distributions but not necessary. In my case I have a set of 42 means with UCL and LCL limits, and the data is

[R] Lacking a dataset called nztrees?

2012-03-28 Thread mandova
I'm looking into the help page on some functions of the {spatstat} packages and they're exemplified by a dataset called nztrees. But it does not exist, and it's not listed on R Datasets Package (http://stat.ethz.ch.sixxs.org/R-manual/R-devel/library/datasets/html/00Index.html). Why is this

Re: [R] Is it possible to de-select with sqlQuery from the RODBC library?

2012-03-28 Thread Bart Joosen
What you can do: SELECT top 1 * FROM your_table; Use this selection to find all your column names in R then paste everything together without the names you don't want and then run your query. Bart -- View this message in context:

Re: [R] reading files from two folders

2012-03-28 Thread uday
Hi Hth -- Gerrit , I found the problem I added sep and one bracket was missing the new codes are as follows setwd(/Groups/data_first/) file_was - list.files(path = ., pattern = v2.0.2.was, all.files = FALSE, full.names = FALSE, recursive = FALSE, ignore.case = FALSE)

[R] customizing plot, Vegan CCA

2012-03-28 Thread millerjeremya
I’m working on a canonical correspondence analysis using Vegan. I’m new to this. I would like to control the shapes and colors for each of my individual display points and can’t figure out the syntax. Here is the script I am using. It comes from Numerical Ecology in R by Bocard et al.

Re: [R] reading files from two folders

2012-03-28 Thread uday
Hi Hth -- Gerrit , thanks for reply sorry it was typo mistake , but I corrected that problem before . But the error is still same. setwd(/Groups/data_first/) file_was - list.files(path = ., pattern = v2.0.2.was, all.files = FALSE, full.names = FALSE, recursive = FALSE, ignore.case =

[R] xyplot lattice fine control of axes limits and thick marks (with log scale)

2012-03-28 Thread maxbre
After a long and winding road (sorry but I'm a novice) I get to a final result which is quite close to what I need; nevertheless I would like to tweak a little further the xyplot so that I can get ***for each single panel defined by variable z*** a finer control over: -the x and y the limits: I

Re: [R] Problems with R Commander version 1.8-3

2012-03-28 Thread John Fox
Dear Mario, This problem is due to an error in the Italian translation file for the Rcmdr. It's fixed in the development version of the package on R-Forge, which eventually will be moved to CRAN. Your options at present are (1) to install the Rcmdr from R-Forge, via the command

Re: [R] What error distribution should I use?

2012-03-28 Thread Lívia Dorneles Audino
Thank's a lot for your help Ben. I already posted my question on r-sig-mixed-models! However you already answered what I needed to know. 2012/3/27 Ben Bolker bbol...@gmail.com Lívia Dorneles Audino livia.audino at gmail.com writes: I'm trying to make a glmm to identify the relationship

Re: [R] SVM. How to use categorical attributes?

2012-03-28 Thread Alekseiy Beloshitskiy
Thank you so much, Ulrich, Will play with this. Best, -Alex From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] on behalf of Ulrich Bodenhofer [bodenho...@bioinf.jku.at] Sent: 28 March 2012 14:40 To: r-help@r-project.org Subject: Re: [R]

Re: [R] What error distribution should I use?

2012-03-28 Thread Lívia Dorneles Audino
Thank's for your help Peter! This post is really very informative. It helped a lot! 2012/3/27 Peter Ehlers ehl...@ucalgary.ca On 2012-03-27 15:11, Ben Bolker wrote: Lívia Dorneles Audinolivia.audinoat gmail.com writes: I'm trying to make a glmm to identify the relationship between insect

Re: [R] What error distribution should I use?

2012-03-28 Thread Lívia Dorneles Audino
Hi Chuck! Actually my data has a lot of zeros because I counted the number of species that is exclusive to each studied fragment. So, my response variable is the number of exclusive species. That's why it has a lot of zeros. I have fragments that do not present exclusive species. Thank's for

[R] getopt does not work as expected!

2012-03-28 Thread Juergen Rose
I have the following script (also attached): #!/usr/bin/Rscript spec=matrix(c( 'verbose', 'v', 1, integer, 'help' , 'h', 0, logical ),ncol=4, byrow=TRUE) spec.dim=dim(spec) spec.opt.long=spec[,1] spec.opt.short=spec[,2] spec.opt.l - spec.dim[1] infile - test.dat args=commandArgs(TRUE);

Re: [R] Transferring a dataset from WINXP to Mac - problems with umlauts

2012-03-28 Thread David Winsemius
On a Mac running with a US locale, when I first type option-u (in Mail.app that shows a highlighted naked umlaut) and then type 'u' I get: ü Similarly with opt-u, 'a' ä -- David. On Mar 28, 2012, at 7:15 AM, Frank Bloos wrote: Dear R-Help, is it possible to correctly display umlauts in

[R] What are the color's name in heat.color(5)

2012-03-28 Thread YN Kim
Hi all, I'm trying to using heat.color function in my data visualization. So, when I set col = heat.color(5), five different colors are displayed from red to weak yellow. But, now I want to know the color's keys or the names of the displayed heat.color. For instance, I think one of the five

Re: [R] Spatstat users.

2012-03-28 Thread David Winsemius
No I am not an active spatstat user. David Winsemius, MD West Hartford, CT __ 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-guide.html and provide

Re: [R] matrix(unlist(strsplit())) 'missing value' issue

2012-03-28 Thread Petr PIKAL
What problem? Nabble is not available to all and here is not much to cook from. Nobody any solution for my problem?? -- View this message in context: http://r.789695.n4.nabble.com/matrix-unlist- strsplit-missing-value-issue-tp4509065p4511668.html Sent from the R help mailing list

Re: [R] Transferring a dataset from WINXP to Mac - problems with umlauts

2012-03-28 Thread Prof Brian Ripley
On 28/03/2012 13:27, David Winsemius wrote: On a Mac running with a US locale, when I first type option-u (in Mail.app that shows a highlighted naked umlaut) and then type 'u' I get: ü Similarly with opt-u, 'a' ä I don't think that was the problem: he very likely has a German keyboard in

Re: [R] Test Normality

2012-03-28 Thread Michael Friendly
On 3/27/2012 8:32 PM, Sindy Carolina Lizarazo wrote: Good Night I made different test to check normality and multinormality in my dataset, but I don´t know which test is better. To verify univariate normality I checked: shapiro.test, cvm.test, ad.test, lillie.test, sf.test or jaque.bera.test

Re: [R] What are the color's name in heat.color(5)

2012-03-28 Thread Kevin Wright
They have hex RGB values instead of names: R heat.colors(5) [1] #FFFF #FF5500FF #FFAA00FF #00FF #80FF Kevin On Wed, Mar 28, 2012 at 7:24 AM, YN Kim y2sile...@gmail.com wrote: Hi all, I'm trying to using heat.color function in my data visualization. So, when I set col =

[R] how to avoid farword slashes in RJSONIO output

2012-03-28 Thread MLSC
Hi All, I wanted know how to get RJSONIO output with out forward slashes? for example: ({\total\:5,\results\:[\n {\n \id\: 1,\n\name\: \Monkey D Luffy\,\n\position\: \Captain\,\n\ambition\: \I Will become the pirate king\ \n},\n{\n \id\: 2,\n\name\: \Roronoa zoro\,\n\position\:

Re: [R] fitted values with locfit

2012-03-28 Thread Liaw, Andy
I believe you are expecting the software to do what it did not claim being able to do. predict.locfit() does not have a type argument, nor can that take on terms. When you specify two variables in the smooth, a bivariate smooth is done, so you get one bivariate smooth function, not the sum of

Re: [R] Lacking a dataset called nztrees?

2012-03-28 Thread R. Michael Weylandt
Oh ye of little faith... That's because it's in the spatstat package: type data(nztrees) to load it. Michael On Wed, Mar 28, 2012 at 4:36 AM, mandova mandova1...@yahoo.cn wrote: I'm looking into the help page on some functions of the {spatstat} packages and they're exemplified by a dataset

Re: [R] What are the color's name in heat.color(5)

2012-03-28 Thread David Winsemius
On Mar 28, 2012, at 9:08 AM, Kevin Wright wrote: They have hex RGB values instead of names: R heat.colors(5) [1] #FFFF #FF5500FF #FFAA00FF #00FF #80FF Kevin The first one is the same as 'red' as can be seen by parsing the RGB values above: col2rgb(red) [,1] red

Re: [R] SVM. How to use categorical attributes?

2012-03-28 Thread Steve Lianoglou
Hi, These suggestions still require you to explicitly compute your feature space or kernel matrix first, which might kill you memory wise. You might consider taking a look at the shogun toolbox: http://www.shogun-toolbox.org/ With some digging, I'm pretty sure you'll find a bag-of-words type

Re: [R] SVM. How to use categorical attributes?

2012-03-28 Thread Steve Lianoglou
Sorry -- I should add that I'm pointing out the potential shogun implementation because I suspect their implementation of a bag-of-words -like kernel would use the kernel trick, so you won't have to map all of your data explicitly into some huge feature space that will blow your memory away. I'm

Re: [R] rgl package broke with R 2.14.2

2012-03-28 Thread Duncan Murdoch
On 12-03-27 6:31 PM, Grimes Mark wrote: Dear People I can't figure out how to fix this problem: rgl won't run under R 2.14.2 (it was working for me before under 2.14.0). The error message is: rgl is currently changing fairly rapidly. I'd suggest trying to install again (the current

Re: [R] R_unload_[PACKAGENAME] not called

2012-03-28 Thread Henrik Alsing Friberg
I experimented with this earlier, but it is indeed tricky business. Last I tried my registered finalizer addressed unmapped memory, because garbage collection happened after my package (and thus the clean-up stuff) was unloaded. This happened using the R_RegisterFinalizerEx function from C++ on

Re: [R] R_unload_[PACKAGENAME] not called

2012-03-28 Thread Duncan Murdoch
On 28/03/2012 10:25 AM, Henrik Alsing Friberg wrote: I experimented with this earlier, but it is indeed tricky business. Last I tried my registered finalizer addressed unmapped memory, because garbage collection happened after my package (and thus the clean-up stuff) was unloaded. This happened

Re: [R] xyplot lattice fine control of axes limits and thick marks (with log scale)

2012-03-28 Thread ilai
On Wed, Mar 28, 2012 at 4:21 AM, maxbre mbres...@arpa.veneto.it wrote: After a long and winding road (sorry but I'm a novice) I get to a final result which is quite close to what I need; nevertheless I would like to tweak a little further the xyplot Without dput(mydata) you are the only one

Re: [R] inefficient ifelse() ?

2012-03-28 Thread Manta
I have a similar problem. I have a dataset and an element. If the element is equal to YY, I want to take the first column of the dataset, otherwise I want to take the second column. The following does not work, as it only evaluates the first element. Any idea? a=c(AAAXXX,BBBXXX)

Re: [R] What are the color's name in heat.color(5)

2012-03-28 Thread Eik Vettorazzi
Hi, just one more note, col2rgb works for color codes as well: (HC.m - col2rgb(heat.colors(5))) and to convert named colors to color codes, something like tabl-apply(col2rgb(colors()),2,function(x)do.call(rgb,c(as.list(x),alpha=255,maxColorValue=255))) names(tabl)-colors() tabl should work.

Re: [R] inefficient ifelse() ?

2012-03-28 Thread R. Michael Weylandt
Just use a regular if statement: X - if(b == YY) a else substr(a, 1,3) Michael On Wed, Mar 28, 2012 at 11:19 AM, Manta mantin...@libero.it wrote: I have a similar problem. I have a dataset and an element. If the element is equal to YY, I want to take the first column of the dataset, otherwise

Re: [R] matrix(unlist(strsplit())) 'missing value' issue

2012-03-28 Thread ilai
On Wed, Mar 28, 2012 at 6:49 AM, Petr PIKAL petr.pi...@precheza.cz wrote: What problem? Nabble is not available to all and here is not much to cook from. Indeed. Also the OP actually provided their own solution, just 5 more minutes of googling to find ?sub. bankoffer.3 -

Re: [R] inefficient ifelse() ?

2012-03-28 Thread Berend Hasselman
On 28-03-2012, at 17:19, Manta wrote: I have a similar problem. I have a dataset and an element. If the element is equal to YY, I want to take the first column of the dataset, otherwise I want to take the second column. The following does not work, as it only evaluates the first element. Any

Re: [R] Urgent - I really need some help lme4 model avg Estimates

2012-03-28 Thread Dragonwalker
Thank you Mitchell, I will try that. So I presume that the initial paper where they showed the estimates AND the intercept from a model averaging procedure may have been done using a different method? Would it still be prudent to use a global model and then perhaps show the top so many, perhaps

[R] Local linear regression

2012-03-28 Thread Dimitris.Kapetanakis
Dear all, I wrote a code which supposedly should create a function which performs a local linear regression at point x (lower case). The issue is that I have no reference to check if indeed my code performs correctly a local linear regression. Is there any function in R that performs a local

[R] fun.zero.omit {GLDEX}

2012-03-28 Thread Nelson Nomura
Hello. Sorry for the newbie question. When I run the function fun.zero.omit (from the GLDEX package), I get a matrix and some descriptive information. Can someone tell me how to strip out the descriptive information so that I'm only left with the matrix? Here's an example:

[R] How to move files from one folder to another

2012-03-28 Thread sagarnikam123
How to move files from one folder to another i.e. cutting/copying from one location pasting to another -- View this message in context: http://r.789695.n4.nabble.com/How-to-move-files-from-one-folder-to-another-tp4512087p4512087.html Sent from the R help mailing list archive at Nabble.com.

[R] blocked MRPP and ISA

2012-03-28 Thread Whitby, Michael D
I have been trying to find a way to do a blocked MRPP (MRBP) and blocked ISA in R but have not found anything (in fact I found some things that say it is not possible). Does anyone know of a package, or written a script, that can do either of these? Thank you for any help.

Re: [R] Lacking a dataset called nztrees?

2012-03-28 Thread mandova
yeah I forgot to library(XX) first... Didn't know packages have their own datasets. Thx! -- View this message in context: http://r.789695.n4.nabble.com/Lacking-a-dataset-called-nztrees-tp4511692p4512375.html Sent from the R help mailing list archive at Nabble.com.

[R] POSIXlt vs POSIXct

2012-03-28 Thread Steven R Corsi
Hello R users I am searching for a descriptive summary of the use of POSIXlt as compared to POSIXct date/time formats. I have been using them extensively for different purposes, but still can't quite understand when to use which one for the most efficient coding and use. I typically use them

[R] Logistic regression

2012-03-28 Thread carlb1
Hi does anyone know how to do a logistic regression in R? any help appreciated Carl -- View this message in context: http://r.789695.n4.nabble.com/Logistic-regression-tp4512658p4512658.html Sent from the R help mailing list archive at Nabble.com.

[R] Extract index column the max

2012-03-28 Thread MSousa
Hello, I have this data.frame, and i need the value the column where the value in matrix is maximum res3 77 86 93 106 161 77 20 0 0 1 0 86 1 12 1 2 1 93 0 2 3 3 0 106 0 4 3 18 0 161 0 0 0 1 55 I use this d-which(res3==max(res3),

[R] Need Help

2012-03-28 Thread Partha Sinha
Sir, Which modules/package should I use to solve the following statistical models: 1. Multiple Objective Programming (MOP) 2. Data Envelopment Analysis (DEA) 3. Lab Experiments 4. Multiple Criteria Decision Making (MCDM) Regards, PS __

[R] Making Knitr work

2012-03-28 Thread Jessica Streicher
Might not be the best place to ask, but i could get lucky.. I have setup an eclipse environment to write sweave files lately and wanted to switch to knitr. I could get it to work on easy files, but my earlier written sweave file fails to be knit properly. Here is the error message: Quitting

[R] p-value using survival ROC

2012-03-28 Thread Anamaria Crisan
It is possible to calculate the c-index for time dependent outcomes (such as disease) using the survivalROChttp://cran.r-project.org/web/packages/survivalROC/survivalROC.pdf package in R. My question is : is it possible to produce a p-value for the c-index that is calculated (at a specific point

Re: [R] matrix(unlist(strsplit())) 'missing value' issue

2012-03-28 Thread MaartenJacobs
Thank you!! Been googling for hours, but kind of hard to find something if you don't know how to look for it. So thanks again!! Greetings Maarten ilai-2 wrote On Wed, Mar 28, 2012 at 6:49 AM, Petr PIKAL lt;petr.pikal@gt; wrote: What problem? Nabble is not available to all and here is not

Re: [R] Logistic regression

2012-03-28 Thread Sarah Goslee
On Wed, Mar 28, 2012 at 11:49 AM, carlb1 carl19...@hotmail.com wrote: Hi does anyone know how to do a logistic regression in R? I do. And doubtless many other folks on R-help do too. any help appreciated You should probably start here:

Re: [R] How to move files from one folder to another

2012-03-28 Thread Sarah Goslee
Within R? On what OS? You could presumably use system() with the OS-specific command, which might be cp or mv or copy or move depending on your OS. But I'm not sure why you'd want to do this from R. On Wed, Mar 28, 2012 at 8:07 AM, sagarnikam123 sagarnikam...@gmail.com wrote: How to move files

Re: [R] Need Help

2012-03-28 Thread Sarah Goslee
On Wed, Mar 28, 2012 at 11:49 AM, Partha Sinha pnsinh...@gmail.com wrote: Sir, I am not a Sir, but nonetheless I'd suggest you start with http://www.rseek.org and do some reading on the packages and functions available for those methods. The Optimization task view might also be of interest. I'm

Re: [R] POSIXlt vs POSIXct

2012-03-28 Thread Duncan Murdoch
On 28/03/2012 10:06 AM, Steven R Corsi wrote: Hello R users I am searching for a descriptive summary of the use of POSIXlt as compared to POSIXct date/time formats. I have been using them extensively for different purposes, but still can't quite understand when to use which one for the most

Re: [R] Extract index column the max

2012-03-28 Thread David Winsemius
On Mar 28, 2012, at 11:18 AM, MSousa wrote: Hello, I have this data.frame, and i need the value the column where the value in matrix is maximum res3 77 86 93 106 161 77 20 0 0 1 0 86 1 12 1 2 1 93 0 2 3 3 0 106 0 4 3 18 0 161 0 0 0 1 55 I use

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

2012-03-28 Thread Paul Miller
Hello All, Figured out how to get the gaps in the bars when a drug stops and the starts again (see below). Made the graph overlaid/superimposed/stacked and got the desired result. Not sure this is how an expert would do it. But it's simple and it works. Paul setwd(N:/Regimen

Re: [R] Making Knitr work

2012-03-28 Thread Peter Langfelder
On Wed, Mar 28, 2012 at 9:09 AM, Jessica Streicher j.streic...@micromata.de wrote: Might not be the best place to ask, but i could get lucky.. I have setup an eclipse environment to write sweave files lately and wanted to switch to knitr. I could get it to work on easy files, but my earlier

Re: [R] blocked MRPP and ISA

2012-03-28 Thread Rachel W
I am not sure why it wouldn't let me see your thread, so I am replying directly. Look at the BLOSSOM software. It comes with a pretty clear manual too. http://www.fort.usgs.gov/products/software/blossom/ Hope this helps. From: mdwhi...@bsu.edu To: r-help@r-project.org Date: Wed, 28

Re: [R] Making Knitr work

2012-03-28 Thread Yihui Xie
You must have set the package option 'base.dir' (opts_knit$set(base.dir = )), and this directory does not exist. I do not recommend you to use this option unless you really understand it. Regards, Yihui -- Yihui Xie xieyi...@gmail.com Phone: 515-294-2465 Web: http://yihui.name Department of

[R] discrepancy between paired t test and glht on lme models

2012-03-28 Thread Rajasimhan Rajagovindan
Hi folks, I am working with repeated measures data and I ran into issues where the paired t-test results did not match those obtained by employing glht() contrasts on a lme model. While the lme model itself appears to be fine, there seems to be some discrepancy with using glht() on the lme

[R] lapply and paste

2012-03-28 Thread Ed Siefker
I have a list of suffixes I want to turn into file names with extensions. suff- c(C1, C2, C3) paste(filename_, suff[[1]], .ext, sep=) [1] filename_C1.ext How do I use lapply() on that call to paste()? What's the right way to do this: filenames - lapply(suff, paste, ...) ? Can I have lapply()

Re: [R] Is it possible to de-select with sqlQuery from the RODBC library?

2012-03-28 Thread Eric Fail
Thank you Bart for your idea, the thing is that I have a large number of tables and I would like to avoid having to pull them at all. I currently have a list that I use as a lookup table in a loop with an if else statement to sort between tables I want to sqlFetch (take everything) and tables

Re: [R] lapply and paste

2012-03-28 Thread R. Michael Weylandt
I think you're confused about the need for lapply -- paste is vectorized so this paste(filename_, suff, .ext, sep = ) will work. But if you want to use lapply (for whatever reason) try this: lapply(suff, function(x) paste(filename_, x, .ext, sep = ) Michael On Wed, Mar 28, 2012 at 2:31 PM, Ed

Re: [R] lapply and paste

2012-03-28 Thread Ed Siefker
Thank you, I was confused about that. What exactly is lapply for then, if R handles this kind of thing automatically? Are there functions that are not vectorized? On Wed, Mar 28, 2012 at 1:37 PM, R. Michael Weylandt michael.weyla...@gmail.com wrote: I think you're confused about the need for

Re: [R] lapply and paste

2012-03-28 Thread Sarah Goslee
suff isn't a list, so lapply() isn't the right choice. How about instead: suff- c(C1, C2, C3) sapply(suff, function(x)paste(filename_, x, .ext, sep=)) C1C2C3 filename_C1.ext filename_C2.ext filename_C3.ext On Wed, Mar 28, 2012 at 2:31 PM, Ed

Re: [R] lapply and paste

2012-03-28 Thread R. Michael Weylandt michael.weyla...@gmail.com
Yes, there are non-vectorized functions e.g, integrate(), or you can use lapply() to apply a vectorized function to each element of a list (which is not what suff was) individually: x - list(1:3, 1:4, 1:5, 2:7) mean(x) # bad lapply(x, mean) #good Michael On Mar 28, 2012, at 2:44 PM, Ed

Re: [R] lapply and paste

2012-03-28 Thread Thomas Lumley
On Thu, Mar 29, 2012 at 7:44 AM, Ed Siefker ebs15...@gmail.com wrote: Thank you, I was confused about that.  What exactly is lapply for then, if R handles this kind of thing automatically?  Are there functions that are not vectorized? There are, especially ones you write yourself that don't

[R] CI with confint

2012-03-28 Thread Benjamin Caldwell
Hello all, I'm trying to use confint from the MASS package to compute confidence intervals for an nls object. When I plot the results, however, they don't make sense - lines cross over the fitted model or just don't match the data. Code is : Thanks for help dat-data.frame(a,b) with(dat,

[R] length(unique(v))

2012-03-28 Thread Sam Steingold
is there a built-in synonym for function(v) length(unique(v)) ?? thanks -- Sam Steingold (http://sds.podval.org/) on Ubuntu 11.10 (oneiric) X 11.0.11004000 http://www.childpsy.net/ http://pmw.org.il http://iris.org.il http://ffii.org http://honestreporting.com http://openvotingconsortium.org Two

[R] scanning data into r

2012-03-28 Thread joel.green
Hey I am having trouble importing data into R, my data field looks like this 21 TEST DATA 32 year:2012 33 34 5 36 I require the the number at the start of each line however the text is not needed, i am struggling to get R to import the data with out changing the file itself? how do i

[R] How to get all possible combinations?

2012-03-28 Thread Somnath Dhar
Dear all, suppose I have a vector with elements as: Vec - c(2,3,4,5,6) Now I want to have all possible combination of length 3 using those elements and without any repetition. Like, I want to have 1 possibility like 2-3-4 but not 3-2-4. Can somebody guide me how to achieve that in R? Thanks

[R] Why does this work? plyr within-subset normalization

2012-03-28 Thread z2.0
Working code that normalize each row's value against the subset's maximum. Does the invocation of max() somehow instruct R to 'step back' and evaluate the subset? Thanks, Zack -- View this message in context:

Re: [R] scanning data into r

2012-03-28 Thread z2.0
Just read into a data.frame with read.table and then subset to use the first column. e.g., your_desired_data - data.frame(read.table(path_to_file, sep = , fill = T)) your_desired_data - your_desired_data[,1] -- View this message in context:

Re: [R] How to get all possible combinations?

2012-03-28 Thread Jorge I Velez
Perhaps this? t(combn(c(2,3,4,5,6), 3)) HTH, Jorge.- On Wed, Mar 28, 2012 at 3:17 PM, Somnath Dhar wrote: Dear all, suppose I have a vector with elements as: Vec - c(2,3,4,5,6) Now I want to have all possible combination of length 3 using those elements and without any repetition.

[R] how to match exact phrase using gsub (or similar function)

2012-03-28 Thread Markus Weisner
trying to switch out addresses that have double directions, such as the following example: a = S S Main St Interstate 95 a = gsub(pattern=S S , replacement=S , a) … the problem is that I don't want to affect instances where this might be a correct address such as the following: 3421 BIGS St

Re: [R] scanning data into r

2012-03-28 Thread David Winsemius
On Mar 28, 2012, at 2:23 PM, joel.green wrote: Hey I am having trouble importing data into R, my data field looks like this 21 TEST DATA 32 year:2012 33 34 5 36 read.table(text=21 TEST DATA + 32 year:2012 + 33 + 34 + 5 + 36, fill=TRUE)[1] V1 1 21 2 32 3 33 4 34 5 5 6 36 I

  1   2   >