Re: [R] scatterplot3d and box3d

2012-11-20 Thread Uwe Ligges
On 19.11.2012 23:34, Ben Bolker wrote: xgao at genome.wustl.edu writes: I created a 3d scatter plot and am trying to change the color of outer box lines with box3d. Anybody can help me to figure out how to do this? My example is: library(scatterplot3d) x=seq(1:6) y=seq(7:12) z=x*2

Re: [R] Is it possible to be sponsored by R?

2012-11-20 Thread Barry Rowlingson
On Mon, Nov 19, 2012 at 9:39 PM, Christophe Genolini cgeno...@u-paris10.fr wrote: Hi the list, I am a member of the organizing comity of the French Statistics Association (SFdS)'s conference. We are looking for sponsors. Some software (SAS, RITME, ...) are represented. Do you know if there

[R] tcl/tk problem with tklistbox,the character and Rcmdr.

2012-11-20 Thread vincent guyader
I everyone, i have a little problem with tklistbox,the character and Rcmdr. Please look at this code require(tcltk) tt-tktoplevel() levels.list2 -tklistbox(tt,selectmode=multiple,exportselection=FALSE, height=4, yscrollcommand=function(...)tkset(levels.list2.scroll,...))

Re: [R] tcl/tk problem with tklistbox,the character and Rcmdr.

2012-11-20 Thread Yves Reecht
Hi Vincent, This is probably because Rcmdr redefine the tclvalue command (check the message frame of the Rcmdr window). You can try: tcltk::tclvalue(tkget(levels.list2,4)) HTH, Yves Le 20/11/2012 11:25, vincent guyader a écrit : I everyone, i have a little problem with tklistbox,the

Re: [R] tcl/tk problem with tklistbox,the character and Rcmdr.

2012-11-20 Thread vincent guyader
very good idea! thanks a lot. Vincent 2012/11/20 Yves Reecht yves.ree...@gmail.com Hi Vincent, This is probably because Rcmdr redefine the tclvalue command (check the message frame of the Rcmdr window). You can try: tcltk::tclvalue(tkget(levels.list2,4)) HTH, Yves Le

Re: [R] help on matrix column removal based on another matrix results

2012-11-20 Thread Rui Barradas
Hello, Try removing the square parenthesis from the left hand side, Vsim - Vsim[, NSErr 0.6] Hope this helps, Rui Barradas Em 20-11-2012 04:07, Irucka Embry escreveu: Hi Rui, how are you? I want to thank you for your assistance again. I'm sorry, but the code that you provided for me did

[R] Classification methods - which one?

2012-11-20 Thread Pedro Silva
Dear Peter, There are several packages that try to address this type of problem (although the remarks made by Max are something that we should always keep in mind), and I also recommend those with perform some form of regularized, penalized or shrunken linear discriminant analysis with a

Re: [R] Help: Meta-analysis with metacor

2012-11-20 Thread Viechtbauer Wolfgang (STAT)
There are multiple issues here. First of all, you should simplify your code to: metacor(cor, n, title=title, complab=comparison, outclab=outcome) The rest of your call to metacor() is just using the defaults, so no need to repeat all of that in your call and apparently one of these arguments is

Re: [R] kinitr

2012-11-20 Thread R. Michael Weylandt
On Tue, Nov 20, 2012 at 1:57 AM, Muhuri, Pradip (SAMHSA/CBHSQ) pradip.muh...@samhsa.hhs.gov wrote: Hello, I am an Intro-level R and ggplot2 user and looking for resources to self teach dynamic report generation in R using knitr. Any advice would be highly appreciated.

Re: [R] github

2012-11-20 Thread R. Michael Weylandt
On Tue, Nov 20, 2012 at 2:07 AM, Muhuri, Pradip (SAMHSA/CBHSQ) pradip.muh...@samhsa.hhs.gov wrote: Hello, I would like to learn how to set up Github/repository and upload/update files and am looking for Github for Dummies. Any help will be appreciated. I believe Hadley has done some

[R] Fit Gumbel Distribution using Method of Moments

2012-11-20 Thread Tonja Krueger
Hi all! I'm sure this is a stupid question but I can't find an answer. How can I fit the Gumbel distribution to my data using The Method of Moments in R? Thank you for helping me, Tonja __ R-help@r-project.org mailing list

Re: [R] lattice: defining grouping variable only for the upper/lower panel with splom

2012-11-20 Thread AnjaM
Thank you very much, this works! -- View this message in context: http://r.789695.n4.nabble.com/lattice-defining-grouping-variable-only-for-the-upper-lower-panel-with-splom-tp4650033p4650144.html Sent from the R help mailing list archive at Nabble.com.

[R] Coefficient of determination for non-linear equations system (nlsystemfit)

2012-11-20 Thread Antti Simola
Hello everyone, I have estimated system of three linear equations with one non-linear restrictions with nlsystemfit. I was wondering how I can calculate the R-squared (or some alternative coefficient of determination) for the whole system. This is automatically given by linear systemfit but

[R] double gaussian with mixdist: what's wrong?

2012-11-20 Thread Sepp Tannhuber
Dear all, I am trying to fit a double gaussian to some data using the mixdist package: --- begin code ---   library(mixdist)   time - seq(673,723)   counts - c(3,12,8,12,18,24,39,48,64,88,101,132,198,253,331,419,563,781,1134,1423,1842,2505,374,6099,9343,13009,                    

[R] Is it possible to draw with function image.asc and show a legend strip in the same window?

2012-11-20 Thread riodementa
I want to represent a Raster map of class asc, and in the same window a list of points with the funcion points(x,y).. I considered two options: 1. With image.plot(map) and then use points(x,y). It works, and the points are in the correct place, but when I maximize the plot the points are

Re: [R] github

2012-11-20 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Thanks, Michael. Pradip Pradip K. Muhuri, PhD Statistician Substance Abuse Mental Health Services Administration The Center for Behavioral Health Statistics and Quality Division of Population Surveys 1 Choke Cherry Road, Room 2-1071 Rockville, MD 20857   Tel: 240-276-1070 Fax: 240-276-1260

Re: [R] kinitr

2012-11-20 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Dear Michael, I really appreciated that you have sent me the link info - Jeromy Anglim's Blog. I was exactly looking for this kind of resources about R Markdown and knitr. All this would be of immense help. Thank you so much. Pradip Pradip K. Muhuri, PhD Statistician Substance Abuse Mental

Re: [R] Fit Gumbel Distribution using Method of Moments

2012-11-20 Thread Bert Gunter
Is this homework? R has a no homework policy. -- Bert On Tue, Nov 20, 2012 at 12:33 AM, Tonja Krueger tonja.krue...@web.de wrote: Hi all! I'm sure this is a stupid question but I can't find an answer. How can I fit the Gumbel distribution to my data using The Method of Moments in

[R] centering gWindow on maximising

2012-11-20 Thread KRISHNA2222
Hi, I am developing a Gui in R for a Statistical Purpose. i am wondering to know whether there is any option in R for centering the objects in a window when a window is maximized. All the gaps, positions and sizes are to be scaled up properly on maximizing, giving an same view as it gives under

[R] lattice density plot: add vertical lines at groupwise medians for all panels

2012-11-20 Thread AnjaM
Suppose you have the following code: ## Start code## data(Chem97, package=mlmRev) densityplot(~gcsescore | factor(score), groups=gender, data=Chem97, auto.key=TRUE, plot.points=FALSE, ref=TRUE, panel=function(x,...){

Re: [R] Performing gage RR study in R w/more than 2 factors

2012-11-20 Thread Matt Jacob
On Mon, Nov 19, 2012, at 16:31, Bert Gunter wrote: I believe that you need to consult a local statistician, as there are likely way too many statistical issues here that you do not fully understand. Alternatively, try posting to a statistical list like stats.stackexchange.com, as I think most

Re: [R] Performing gage RR study in R w/more than 2 factors

2012-11-20 Thread Matt Jacob
On Mon, Nov 19, 2012, at 18:26, David Winsemius wrote: My guess is that you do not understand the meaning of a random factor. I certainly did not when I first encountered it. All my training had been with ordinary regression and analysis of variance. These are methods for what in mixed

Re: [R] help on matrix column removal based on another matrix results

2012-11-20 Thread Irucka Embry
Hi Rui, thank-you very much! That worked perfectly. Irucka -Original Message- From: Rui Barradas [ruipbarra...@sapo.pt] Sent: 11/20/2012 5:42:36 AM To: iruc...@mail2world.com Cc: r-help@r-project.org Subject: Re: [R] help on matrix column removal based on another matrix results

[R] Ordering List Items Chronologically

2012-11-20 Thread Simon Kiss
Dear colleagues, Is there a way to order list items by date? I have a series of surveys in a list where the name of each list item is the date the survey was taken but the list items are out of order. Each data frame has a variable in it with the survey date as well, if that helps. Yours,

Re: [R] Ordering List Items Chronologically

2012-11-20 Thread William Dunlap
o - order(as.Date(names(mylist))) myListInDateOrder - mylist[o] Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Simon Kiss Sent: Tuesday, November 20, 2012 8:41 AM To:

Re: [R] Ordering List Items Chronologically

2012-11-20 Thread Sarah Goslee
As long as all your names of list elements are in the same format, this should do it: mylist[order(names(mylist))] Thanks for the reproducible example, Sarah On Tue, Nov 20, 2012 at 11:40 AM, Simon Kiss sjk...@gmail.com wrote: Dear colleagues, Is there a way to order list items by date? I

[R] Read Write permissions in current directory

2012-11-20 Thread Nick Duncan
Dear All, I was just installing 'linkcomm' and had the following message which I have not seen before and don't really understand: NOTE: To use linkcomm, you require read and write permissions in the current directory (see: help(getwd), help(setwd)) Having looked at these help pages I'm none

[R] correct function formation in R

2012-11-20 Thread Omphalodes Verna
Dear list!   I have question of 'correct function formation'. Which function (fun1 or fun2; see below) is written more correctly? Using ''structure'' as output or creating empty ''data.frame'' and then transform it as output? (fun1 and fun1 is just for illustration).   Thanks a lot, OV   code:

Re: [R] correct function formation in R

2012-11-20 Thread Rui Barradas
Hello, I believe it's a matter of personal taste. I find fun2 more readable, others may not agree. Rui Barradas Em 20-11-2012 17:39, Omphalodes Verna escreveu: Dear list! I have question of 'correct function formation'. Which function (fun1 or fun2; see below) is written more correctly?

[R] Puzzling Execution halted

2012-11-20 Thread Gang Chen
I wrote an R program that does heavy computations with hundreds of lines of code. It's running fine both interactively and in batch mode on a Mac OS X computer. The program also has no problem running on a Linux system (Fedora 14) interactively. However, when I try it on the terminal in batch mode

[R] Multilevel analysis using nlme (lme) . Error using z-scores

2012-11-20 Thread Mikołaj Hnatiuk
Hi, i am trying to learn something about multilevel analysis using a great Discovering statistics using R. I constructed some sample data and then tried to fit a model. Generally model fits well, however when trying to fit the same model using z-score (standarizded) variables i got an error:

[R] Deviation standard variance

2012-11-20 Thread Elli
How I can plot the standard deviation and variance separately or is there a graph where I can show the two? -- View this message in context: http://r.789695.n4.nabble.com/Deviation-standard-variance-tp4650181.html Sent from the R help mailing list archive at Nabble.com.

[R] Controlling the number of interactions of a lme

2012-11-20 Thread cleberchaves
Hello everybody, i'm trying to run a lme but am having a problem. My model have many response variables and when i run the anova, the number of interactions (up to six) is great and the p-values of all variables not appear. I wanted to know if i could to control the number of interactions of the

Re: [R] Puzzling Execution halted

2012-11-20 Thread Milan Bouchet-Valat
Le mardi 20 novembre 2012 à 12:54 -0500, Gang Chen a écrit : I wrote an R program that does heavy computations with hundreds of lines of code. It's running fine both interactively and in batch mode on a Mac OS X computer. The program also has no problem running on a Linux system (Fedora 14)

Re: [R] Read Write permissions in current directory

2012-11-20 Thread Jeff Newmiller
This question is at best about how R is configured on your computer, but far more likely the problem lies in your operating system, which you have not identified and details of which are off-topic here anyway. If you use getwd() to identify which directory you are in, or pay attention to which

[R] Updating Tom Short's R Reference Card, contacting him?

2012-11-20 Thread Matthew Baggott
I made an update/reboot of Tom Short's classic and public domain R Reference Card. His is from late 2004 and I've found myself giving it to new R users with additional notes about packages. If anyone knows how to reach Tom, that would be great. I am titling this reboot Short R Reference, in a

Re: [R] Can't install rJava

2012-11-20 Thread josueao
For me this worked: install.packages('rJava',,'http://www.rforge.net/') -- View this message in context: http://r.789695.n4.nabble.com/Can-t-install-rJava-tp3948460p4650191.html Sent from the R help mailing list archive at Nabble.com. __

[R] [solved!] Re: [lattice] format and rotation of strip text

2012-11-20 Thread Tom Roche
https://stat.ethz.ch/pipermail/r-help/2012-November/329509.html - without the line above commented out, strip values are (correctly) typo, should be incorrectly ^ all 1.123 + with the line above commented out, strip values are (correctly) 1.123

Re: [R] Puzzling Execution halted

2012-11-20 Thread Gang Chen
Thanks for the pointer! Unfortunately it doesn't look like a memory issue: /var/log/messages contains nothing about memory problems. And the Linux system has enough RAM for this computation. Gang On Tue, Nov 20, 2012 at 1:03 PM, Milan Bouchet-Valat nalimi...@club.fr wrote: Le mardi 20 novembre

Re: [R] Deviation standard variance

2012-11-20 Thread Jessica Streicher
At least for me this is not enough information to go on Elli. What exactly is the problem? Can you post example data? On 20.11.2012, at 18:39, Elli wrote: How I can plot the standard deviation and variance separately or is there a graph where I can show the two? -- View this message

Re: [R] Puzzling Execution halted

2012-11-20 Thread Duncan Murdoch
On 20/11/2012 12:54 PM, Gang Chen wrote: I wrote an R program that does heavy computations with hundreds of lines of code. It's running fine both interactively and in batch mode on a Mac OS X computer. The program also has no problem running on a Linux system (Fedora 14) interactively. However,

Re: [R] correct function formation in R

2012-11-20 Thread Duncan Murdoch
On 20/11/2012 12:39 PM, Omphalodes Verna wrote: Dear list! I have question of 'correct function formation'. Which function (fun1 or fun2; see below) is written more correctly? Using ''structure'' as output or creating empty ''data.frame'' and then transform it as output? (fun1 and fun1 is

Re: [R] Updating Tom Short's R Reference Card, contacting him?

2012-11-20 Thread Jeff Newmiller
Some issues: a) I would advocate marking - as not recommended. Anyone who needs this reference probably should avoid this operator. b) The $ operator does an approximate match with the name. [[name]] does an exact match. c) Matrices are not data frames. Their indexing by the [ operator is

Re: [R] Is it possible to draw with function image.asc and show a legend strip in the same window?

2012-11-20 Thread David Winsemius
On Nov 20, 2012, at 3:15 AM, riodementa wrote: I want to represent a Raster map of class asc, and in the same window a list of points with the funcion points(x,y).. I considered two options: 1. With image.plot(map) and then use points(x,y). It works, and the points are in the correct

Re: [R] Deviation standard variance

2012-11-20 Thread Elli
I have a data set consisting of 364 variables and I need to graph the standard deviation and variance of each variable as only the function sd () or var () is not enough. And I want to know how I can graph this? -- View this message in context:

[R] parApply computing

2012-11-20 Thread winston87
Hi All. I'm using iparApply/i() function in snow package for parallel computing (boostrapping repeating calculation), as the follows: MyBoostrapping - function( i ){ . } cl - makeCluster( ncluster, type = SOCK) i.circle - as.matrix( 1:128, , 1) parApply( cl, i.circle, 1, FUN =

[R] Using if

2012-11-20 Thread Virgile Capo-Chichi
Dear R users, As a new comer to R, I would like to create a new variable using if statements but don't know how to do it. Basically, I have two variables (EvHint and MinTex). I want to create a third variable called RiskTest. In SPSS, my syntax would look like Compute RiskTest=0. if (EvHint=1

Re: [R] Using if

2012-11-20 Thread Jeff Newmiller
?ifelse Quite distinct from if () {} else {}. --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

Re: [R] Using if

2012-11-20 Thread Sarah Goslee
Hi, Thanks for providing a reproducible example. On Tue, Nov 20, 2012 at 2:08 PM, Virgile Capo-Chichi vcapochi...@gmail.com wrote: Dear R users, As a new comer to R, I would like to create a new variable using if statements but don't know how to do it. Basically, I have two variables

Re: [R] Using if

2012-11-20 Thread William Dunlap
In SPSS, my syntax would look like Compute RiskTest=0. if (EvHint=1 MinTex=1) RiskTest=1. RiskTest - EvHint == 1 MinTex == 1 will create a logical vector, one with values TRUE and FALSE. If you really want 1 and 0 then add RiskTest - as.integer(RiskTest), but using logicals instead of

Re: [R] Using if

2012-11-20 Thread jim holtman
Try this (convert the logical TRUE to the value 1): EvHint-c(0, 0, 0, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1, 0, 0) MinTex-c(0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0) RiskTest - (EvHint == 1 MinTex == 1) + 0L rbind(EvHint, MinTex, RiskTest) [,1] [,2] [,3]

Re: [R] Deviation standard variance

2012-11-20 Thread Jessica Streicher
Replying to post on Nabble: so, basically first you need these values for each variable, yes? m-data.frame(a=1:3,b=4:6,c=7:9) variances-sapply(m,var) sdevs-sapply(m,sd) but i still don't know what you want to show. Just plot the values per variable? That won't be properly visible in just one

[R] data after write() is off by 1 ?

2012-11-20 Thread Brian Feeny
I am new to R, so I am sure I am making a simple mistake. I am including complete information in hopes someone can help me. Basically my data in R looks good, I write it to a file, and every value is off by 1. Here is my flow: str(prediction) Factor w/ 10 levels 0,1,2,3,..: 3 1 10 10 4 8 1

Re: [R] Deviation standard variance

2012-11-20 Thread Jessica Streicher
Anyway, need to go home. If you really only want to plot the values of all variables, and its for some reason not important what value belongs to what exact variable: plot(variances,ylim=range(variances,sdevs)) # black circles points(sdevs,col=2,pch=3) #red crosses On 20.11.2012, at 19:32,

Re: [R] [solved!] Re: [lattice] format and rotation of strip text

2012-11-20 Thread David Winsemius
On Nov 20, 2012, at 10:24 AM, Tom Roche wrote: https://stat.ethz.ch/pipermail/r-help/2012-November/329509.html - without the line above commented out, strip values are (correctly) typo, should be incorrectly ^ all 1.123 + with the line above

Re: [R] data after write() is off by 1 ?

2012-11-20 Thread Duncan Murdoch
On 20/11/2012 2:30 PM, Brian Feeny wrote: I am new to R, so I am sure I am making a simple mistake. I am including complete information in hopes someone can help me. Basically my data in R looks good, I write it to a file, and every value is off by 1. Here is my flow: str(prediction)

Re: [R] data after write() is off by 1 ?

2012-11-20 Thread Brian Feeny
A followup to my own post, I believe I figured this out, but if I should be doing something different please correct: prediction.out - levels(prediction)[prediction] write(prediction.out, file=prediction.csv) This gives me my correctly adjusted values Brian On Nov 20, 2012, at 2:30 PM,

Re: [R] Using if

2012-11-20 Thread Peter Alspach
Tena koe If the vectors are long and time is an issue, it might be better to avoid the use of ifelse(). Try: set.seed(0) EvHint - rbinom(10^7, 1, 0.5) MinTex - rbinom(10^7, 1, 0.5) system.time(RiskTest0 - ifelse(EvHint==1 MinTex==1, 1, 0)) system.time({ RiskTest1 - rep(0, length(EvHint))

Re: [R] data after write() is off by 1 ?

2012-11-20 Thread Rui Barradas
Hello, You are seeing the levels of a factor but saving its values. Internally, factors are coded as consecutive integers starting at 1, and that's what is saved to file using write.table. To have the levels 0, 1, etc and not the corresponding values 1, 2, etc, try

Re: [R] Using if

2012-11-20 Thread arun
HI, May be this: RiskTest-ifelse(EvHint==1 MinTex==1,1,0) RiskTest # [1] 0 0 0 1 0 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 A.K. - Original Message - From: Virgile Capo-Chichi vcapochi...@gmail.com To: r-help@r-project.org Cc: Sent: Tuesday, November 20, 2012 2:08 PM Subject: [R] Using if Dear

Re: [R] Deviation standard variance

2012-11-20 Thread Mark Lamias
What about something like this: #Create dummy data mydf=data.frame(matrix(rnorm(1000), ncol=100)) #Obtain summary stats my.vars=apply(mydf, MARGIN=2, FUN=var) my.std=sqrt(my.vars) #Plot summary stats plot(my.vars, col=red) points(my.std, col=blue) From:

[R] parApply computing

2012-11-20 Thread George
I'm using /parApply/() function in snow package for parallel computing (boostrapping repeating calculation), as the follows: MyBoostrapping - function( i ){ . } cl - makeCluster( ncluster, type = SOCK) i.circle - as.matrix( 1:128, , 1) parApply( cl, i.circle, 1, FUN = MyBoostrapping)

Re: [R] data after write() is off by 1 ?

2012-11-20 Thread Prof Brian Ripley
On 20/11/2012 19:46, Duncan Murdoch wrote: On 20/11/2012 2:30 PM, Brian Feeny wrote: I am new to R, so I am sure I am making a simple mistake. I am including complete information in hopes someone can help me. Basically my data in R looks good, I write it to a file, and every value is off by

Re: [R] data after write() is off by 1 ?

2012-11-20 Thread Duncan Murdoch
On 12-11-20 4:35 PM, Prof Brian Ripley wrote: On 20/11/2012 19:46, Duncan Murdoch wrote: On 20/11/2012 2:30 PM, Brian Feeny wrote: I am new to R, so I am sure I am making a simple mistake. I am including complete information in hopes someone can help me. Basically my data in R looks good, I

[R] about MCMC

2012-11-20 Thread Tania Patiño
Hello all, could you tell wehere I can find information related to MCMC (Monte Carlo). And some examples about this topic. Thanks, Tania [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] about MCMC

2012-11-20 Thread Bert Gunter
Surely you jest! -- Bert On Tue, Nov 20, 2012 at 1:52 PM, Tania Patiño taniu...@gmail.com wrote: Hello all, could you tell wehere I can find information related to MCMC (Monte Carlo). And some examples about this topic. Thanks, Tania [[alternative HTML version deleted]]

[R] Greek letters on title

2012-11-20 Thread bgnumis
Hi all, I want to plot a series x and I want to put on title Time evolution with alpha=0.2, e.g., plot(x, main=Time evotion with alpha=0.2,col=orange) Is there any posibility to put the greek letter and avoid written the word? Many thanks in advance. [[alternative HTML version

[R] Help with glht function and binary data!

2012-11-20 Thread Jessica McLachlan
Hi, I carried out an experiment, using a repeated measures design, in which I broadcast 6 playbacks to 15 nests. I am trying to run a posthoc pairwise comparison on a binomial GLMM with missing data points to determine which playback treatments differ. The response variable (Response) is

Re: [R] kinitr

2012-11-20 Thread Mark Lamias
This is how I learned everything about knitr: http://yihui.name/knitr/ Yihui is great and his site gives  you pretty much all the information you need to get you started. From: Muhuri, Pradip (SAMHSA/CBHSQ) pradip.muh...@samhsa.hhs.gov To: 'R. Michael

[R] Removing columns that are na or constant

2012-11-20 Thread Brian Feeny
I have a dataset that has many columns which are NA or constant, and so I remove them like so: same - sapply(dataset, function(.col){ all(is.na(.col)) || all(.col[1L] == .col) }) dataset - dataset[!same] This works GREAT (thanks to the r-users list archive I found this) however, then

Re: [R] Greek letters on title

2012-11-20 Thread Uwe Ligges
On 20.11.2012 21:47, bgnumis wrote: Hi all, I want to plot a series x and I want to put on title Time evolution with alpha=0.2, e.g., plot(x, main=Time evotion with alpha=0.2,col=orange) Is there any posibility to put the greek letter and avoid written the word? Yes, see ?plotmath

Re: [R] Greek letters on title

2012-11-20 Thread arun
HI, Try this: plot(1,main= expression(paste(Time evolution with , alpha,= 0.2))) A.K. - Original Message - From: bgnumis bgnu...@gmail.com To: r-help@r-project.org Cc: Sent: Tuesday, November 20, 2012 3:47 PM Subject: [R] Greek letters on title Hi all, I want to plot a series x and

Re: [R] Removing columns that are na or constant

2012-11-20 Thread Rui Barradas
Hello, Inline. Em 20-11-2012 22:03, Brian Feeny escreveu: I have a dataset that has many columns which are NA or constant, and so I remove them like so: same - sapply(dataset, function(.col){ all(is.na(.col)) || all(.col[1L] == .col) }) dataset - dataset[!same] This works GREAT (thanks

Re: [R] about MCMC

2012-11-20 Thread Suzen, Mehmet
Hello Tania, You may want to read MCMCpack's paper; http://www.jstatsoft.org/v42/i09/paper and one of the best book in the field: http://www.amazon.com/Markov-Practice-Chapman-Interdisciplinary-Statistics/dp/0412055511 Also look on CRAN view on Bayesian Inference Best, -m On Tue, Nov 20,

Re: [R] Greek letters on title

2012-11-20 Thread Mark Lamias
Try this: titleline=expression(paste(Time evolution with , alpha, =0.2)) plot(x, main=titleline,col=orange) --Mark J. Lamias From: bgnumis bgnu...@gmail.com To: r-help@r-project.org Sent: Tuesday, November 20, 2012 3:47 PM Subject: [R] Greek letters on

Re: [R] kinitr

2012-11-20 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Mark- Thank you for your help. Pradip Pradip K. Muhuri, PhD Statistician Substance Abuse Mental Health Services Administration The Center for Behavioral Health Statistics and Quality Division of Population Surveys 1 Choke Cherry Road, Room 2-1071 Rockville, MD 20857 Tel: 240-276-1070 Fax:

Re: [R] Greek letters on title

2012-11-20 Thread Mark Lamias
May may also find this helpful, especially if you know LaTex: http://yihui.name/en/2011/04/produce-authentic-math-formulas-in-r-graphics/ To: bgnumis bgnu...@gmail.com; r-help@r-project.org r-help@r-project.org Sent: Tuesday, November 20, 2012 6:22 PM

[R] [lattice] how to overlay a geographical map on a levelplot?

2012-11-20 Thread Tom Roche
r-help lattice adepts: I have a question which is somewhat geospatial, so I posted to r-sig-geo rather than here: https://stat.ethz.ch/pipermail/r-sig-geo/2012-November/016757.html summary: How to overlay a geographical map on each panel in a lattice (or Trellis), e.g., of levelplot's? Note I

Re: [R] [lattice] how to overlay a geographical map on a levelplot?

2012-11-20 Thread Pascal Oettli
Hello, Your image is projected, not your panels. Does it match what you are looking for? wld - map('world', xlim=c(west.lon.deg,east.lon.deg), ylim=c(south.lat.deg,north.lat.deg),plot=FALSE) wld - data.frame(lon=wld$x, lat=wld$y) state - map('state', xlim=c(west.lon.deg,east.lon.deg),

Re: [R] [lattice] how to overlay a geographical map on a levelplot?

2012-11-20 Thread Tom Roche
https://stat.ethz.ch/pipermail/r-sig-geo/2012-November/016757.html summary: How to overlay a geographical map on each panel in a lattice (or Trellis), e.g., of levelplot's? https://stat.ethz.ch/pipermail/r-help/2012-November/329714.html Does [this] match what you are looking for? Alas, no:

Re: [R] [lattice] how to overlay a geographical map on a levelplot?

2012-11-20 Thread Pascal Oettli
Hello, I think you need to launch the library 'latticeExtra'. I forgot to mention it. Sorry. library(latticeExtra) Regards, Pascal Le 21/11/2012 14:29, Tom Roche a écrit : https://stat.ethz.ch/pipermail/r-sig-geo/2012-November/016757.html summary: How to overlay a geographical map on

[R] [solved!] [lattice] how to overlay a geographical map on a levelplot?

2012-11-20 Thread Tom Roche
https://stat.ethz.ch/pipermail/r-sig-geo/2012-November/016757.html summary: How to overlay a geographical map on each panel in a lattice (or Trellis), e.g., of levelplot's? https://stat.ethz.ch/pipermail/r-help/2012-November/329716.html you need to launch the library 'latticeExtra' That's

[R] I don't know the difference between rank and order

2012-11-20 Thread li1127217ye
I don't know the difference between rank and order.For example: x=c(10,30,30,20,10,20) x[rank(x,ties.method=first)] [1] 10 10 20 30 30 20 x[order(x)] [1] 10 10 20 20 30 30 the result is quite different, x[rank(x,ties.method=first)] [1] 10 10 20 30 30 20 It is not sorted,why? -- View this

Re: [R] Discrete trait Ornstein–Uhlenbeck in R?

2012-11-20 Thread KRAmazon
Thanks for the feedback. By discrete, I meant a categorical trait (diel activity). I've looked at the discrete models in geiger, but I had to ask about Ornstein-Uhlenbeck because one of my advisors suggested it. -- View this message in context:

[R] quantifying the latent variable of an ordered probit

2012-11-20 Thread ak13
hi, i have an ordered probit model where i estimate the probability of survey participants to have selected a certain category of a censored latent variable from a survey (the latent variable is the evolution of consumer prices). the latent variable is of course continuous, the survey however

Re: [R] e1071 SVM: Cross-validation error confusion matrix

2012-11-20 Thread signal
Did you ever receive a response to this? I did not see one public. I would think that if your dataset was of a large enough size, that 10-fold validation would show an improvement over N:N. Also, any ideas if there is any difference really in using fitted() vs. predict() in your second step? I

[R] McNemar's Test Question

2012-11-20 Thread calicyc
Hi all, I'm wondering if there's a way that I can perform McNemar's test for a small dataset easier. I know this may be different from the usual way, but I'm hoping to throw it out there for help. For now, my orginnal dataset looks something like this... Chemical shifts Pre_yes Pre_no

[R] Function storing error messages in 32 bit R-2.15.2 version

2012-11-20 Thread Maulik Shah
I am using 32 bit R - 2.15.2 version and working on package ltm. The program does not give me any warning message when I run tpm command for a dataset. If run tpm command for the same dataset in 64-bit R 2.15.1 version, I get the warning message shown below. Can someone suggest why is no warning

Re: [R] Updating Tom Short's R Reference Card, contacting him?

2012-11-20 Thread peter dalgaard
On Nov 20, 2012, at 19:45 , Jeff Newmiller wrote: a) I would advocate marking - as not recommended. Anyone who needs this reference probably should avoid this operator. Could we say not for beginners? I mean, it has legitimate uses. If we say not recommended, people will remember it as

Re: [R] e1071 SVM: Cross-validation error confusion matrix

2012-11-20 Thread Brian Feeny
responding to my own question, I see in ?svm man it states fitted() and predict() can do the same thing: # test with train data pred - predict(model, x) # (same as:) pred - fitted(model) On Nov 21, 2012, at 1:08 AM, signal bfe...@mac.com wrote: Did you ever receive a response to this? I did

[R] BiodiversityR and box cox transformation

2012-11-20 Thread Katherine Lizama Allende
Hi all: I am trying to use the BiodiversityR package in Rcommander since I need to box-cox transform some data. However, even when trying with the example files, I get this error just after selecting the box cox transformation in the Menu, and Rcommander and Rstudio just shut down. What is it?

[R] Two way manova

2012-11-20 Thread DusanP
Hello everyone, I would like to perform a 2-way manova test, but I'm having some issues. I implemented like this Y-cbind(Resp1,Resp2,Resp3,,Respn) model-manova(Y tilda FactorA*FactorB) summary.aov(model) 1. I don't know at what level I have to do the Type I error correction. Is it on