Re: [R] Linux editor for R+LaTeX, but not Emacs

Christian Salas wrote: > Using Tinn-R (in windows) is possible to run latex and R from the same > editor, which was great. Now, I am using Ubuntu-linux, which has been > better than WinXP. Currently, I am using Emacs (and then install ESS) > for running LaTeX and R from a same editor-program (li

Re: [R] special characters in linux using dev.print

Thanks Brian, I got lost with the coding between windows and linux. I used the ± (Windows: Alt-0177, instead of "\u00B1") in a code writen under windows and than used on linux. eg: plot(1:10) text(3,8,paste("±", "alt-numeric 0177 works only in windows",sep="")) text(8,3,paste("\u00B1"," \\u00B

Re: [R] Errors in R-2-6-0?

francogrex wrote: > I re-installed version 2-5-0 and it works well. I believe this problem is due > to an incompatibility of BRugs with R-2-6-0 > > > francogrex wrote: >> I was happy today to install the new version of R 2-6-0 But I ran into >> problems I did not have before: >> >> >> setwd("C

Re: [R] tcltk scrollbar

Maybe I do... I think the OP means the first scrollbar isn't tight to the first listbox. This is due to the length of the "Basic Manipulation" label. For example, try to pack this label as follows: tkgrid(lbl.MainT,columnspan=2,sticky="nw") That helps? Jonne. On Mon, 2007-10-08 at 08:37 +0200

Re: [R] Errors in R-2-6-0?

Uwe Ligges wrote: > > The recent version of BRugs works for me. > The file \buffer.txt is required for some inetraction between > OpenBUGS and R and is in R's temporary directory. > No idea why it does not work for you. > That's because I was using the older version of BRugs which was ok

Re: [R] R-2.6.0 and RWinEdt

Installed and tested right now. Works fine, no problem. Thanks, Patrick Uwe Ligges a écrit : > > > Patrick Giraudoux wrote: >> Dear Listers, >> >> I have just installed R-2.6.0 and the RWinEdt package 1.7-6 under >> Windows XP. > > wait for 1.7-7 which should appear on CRAN real soon now. > > Uw

Re: [R] Trust p-values or not ?

Hi John,Peter Thanks for your quick response. 1) I plotted the residual vs fitted values graph and it was fairly random, there was no observable pattern in it - So I guess we are OK on that front. The Q-Q plot was stretched 'S' shaped with the middle part being mostly on the straight line. I trie

Re: [R] tcltk scrollbar

Jonne Zutt wrote: > Maybe I do... > I think the OP means the first scrollbar isn't tight to the first > listbox. This is due to the length of the "Basic Manipulation" label. > > For example, try to pack this label as follows: > tkgrid(lbl.MainT,columnspan=2,sticky="nw") > > Aha. Now I see it t

[R] Incompatible methods ("-.POSIXt", "Ops.difftime") for "-"

Dear all, according to the Help-page of DateTimeClasses {base} I should be able to do time - z with timedate-time objects z a numeric vector (in seconds) or an object of class "difftime". However, on R version 2.6.0 (Windows XP) I get > Sys.time() - as.difftime(c("0:3:20", "11:23:15

[R] Residuals for binomial lmer fits

Dear all, I would like to use the residuals in a general linear mixed effect model to diagnose model fit. I know that the resid function has been implemented for linear mixed models but not yet for general linear mixed effects. Is there a way to get them out of lmer fit objects? I tried sear

[R] How transform log(x+1) all data matrix?

I just have started using R. I am working with abundance data (57 parasite species found in cod) and I would like to know if there is possible to transform all the data matrix. I read that in LDA it is possible to do it with log transformation but I need log (x+1) transformation because I have man

[R] heatmap

Hi, I am having troubles with heatmap(). I have a matrix containing pairwise distance values and I want to plot this matrix with heatmap. I wonder now what distfun is for. Is a distance matrix computed from my initial matrix again and then plotted? cheers, andreas

[R] How to remove legend?

Hello, I want to do some simple bar plots on email arrival data, but I do not quite get the point how to remove the legend. I use the following code: > dow <- read.csv(file="j:/uni_lernen/da_stockspam/svn/data/c411_weekday.txt", header=TRUE) > dow$dayofweek <- factor(dow$dayofweek) > days <- c("S

Re: [R] heatmap

If you just want a coloured representation of your distance matrix, use image() ?image() If you want to cluster your original data and show the original data (not the distances) as a heatmap then use heatmap(), but you should use something like hclust() to cluster the data first. ___

Re: [R] How to remove legend?

Try somelike this library(ggplot2) dow <- data.frame(dayofweek = factor(1:7), p = runif(7), w = c(0, 1, 1, 1, 1, 1, 0)) sc <- scale_fill_continuous() sc$legend <- FALSE qplot(dayofweek, p, data=dow, geom="bar", xlab="Day", ylab="Arrival Rate", main="Spam by Weekday", fill=w) + sc HTH, Thierry -

Re: [R] comparing matched proportions using glm

Thanks very much for your reply Chuck, I have a quick follow up question. You mention putting the data into a 2x2x3 for log-linear model, however my lists have many more than 3 strata, actually thousands. I am trying to work out whether the proportions in list 1 tend to be equal to the proportions

Re: [R] tcltk scrollbar

This is approximately what I get: -- --- ! !!S! ! !S! ! !!C! ! !C! ! !!R! !

Re: [R] a function to compute the cumulative distribution function (cdf) of the gamma

elodie gillain wrote: > > Dear Forum, > > I would need to write a function to calculate the cumulative distribution > function (cdf) of the gamma random variable. I am allowed to use the > following precoded functions: rgamma, pgamma, and dgamma. I am a little > overwhelmed by this project, I

[R] How to change row value based on previous row value

Hi, I would like to fix the data in the following data.frame, and I am having trouble coding this in R - help please! > x AB xy z 1 1 10.0 100 1000 1 2 2 19.8 200 2000 2 3 3 20.1 300 3000 3 4 4 20.3 400 4000 4 5 5 30.0 500 5000 5 Column B is the problem. The

Re: [R] How to remove legend?

Hello Thierry, thanks for your fast response. This does what I want, thank you. Where can I get more documentation for this scale_fill_continuous() method and similar ones? I figured out how to choose the colors but I don't quite get how to add a solid border to the bars. Regards, Christoph -

[R] Mouse Location on Plots

Hi, Does anyone know of a function that dynamically updates a variable with the (x,y) co-ordinates of the mouse pointer on a plot? locator and identify only return the values when you click 'stop' on the graphics device. Any help would be gratefully appeciated. Kind regards, Sam. [[al

[R] Odp: How transform log(x+1) all data matrix?

Hi [EMAIL PROTECTED] napsal dne 08.10.2007 13:20:56: > > I just have started using R. I am working with abundance data (57 parasite > species found in cod) and I would like to know if there is possible to > transform all the data matrix. I read that in LDA it is possible to do it By data matr

Re: [R] R and FDA trials

On Sun, 2007-10-07 at 21:25 -0500, Frank E Harrell Jr wrote: > Ricardo Pietrobon wrote: > > Yesterday I just noticed the new document on R and regulatory aspects > > for biomedical research posted at > > http://www.r-project.org/doc/R-FDA.pdf > > > > Coming from an institution that performs a larg

Re: [R] How to remove legend?

Hi Christophe, Although I'm a fan of ggplot2, it's lack of documentation is a drawback. I was having a similar problem and Hadley suggested the "sc$legend <- FALSE" solution. AFAIK it isn't documented (yet). Have a look at Hadley's website and the archive of this list for more documentation. This

[R] tcltk scrollbar

For example, try to pack this label as follows: tkgrid(lbl.MainT,columnspan=2,sticky="nw") This worked once I changed columnspan=5, can you explain why. It seems i need to do the same for all widgets located above the ListBoxes. Thanks [[alternative HTML version deleted]] _

[R] Odp: How to change row value based on previous row value

Hi [EMAIL PROTECTED] napsal dne 08.10.2007 14:24:13: > Hi, > > I would like to fix the data in the following data.frame, and I am > having trouble coding this in R - help please! > > > x > AB xy z > 1 1 10.0 100 1000 1 > 2 2 19.8 200 2000 2 > 3 3 20.1 300 3000 3 > 4

[R] Specify plot size

Hello, I have another (possibly easy) question: How to specify the size of a plot? When I draw a plot, I can freely change the size of the window, which is nice for single plots to find the best height/width ratio, but as I need a lot of plots in my work, I want to look them all the same, so I ne

Re: [R] Specify plot size

Have you tried using pdf(), postscript(), or jpeg()? All of these have arguments to specify the height and width of the device. On 10/8/07, Christoph Krammer <[EMAIL PROTECTED]> wrote: > Hello, > > I have another (possibly easy) question: > > How to specify the size of a plot? When I draw a plot

[R] Odp: Specify plot size

Hi [EMAIL PROTECTED] napsal dne 08.10.2007 15:15:58: > Hello, > > I have another (possibly easy) question: > > How to specify the size of a plot? When I draw a plot, I can freely change > the size of the window, which is nice for single plots to find the best > height/width ratio, but as I nee

Re: [R] Ambiguities in vector

Hello James, all of your suggestions work very well except of this: FemMal <- cbind(FemV1gezählt[2,], MalV1gezählt[2,]) colnames(FemMal) <- ("Females", "Males") Fehler: syntax error FeMMal [,1] [ ,2] 1 133 79 2 203 237 3 51 76 But it works if I do that: Namen<-c("Female","Mal

Re: [R] Ambiguities in vector

Hi [EMAIL PROTECTED] napsal dne 08.10.2007 15:35:01: > Hello James, > > all of your suggestions work very well except of this: > > FemMal <- cbind(FemV1gezählt[2,], MalV1gezählt[2,]) > > colnames(FemMal) <- ("Females", "Males") > Fehler: syntax error Syntax error means something is missing an

Re: [R] Ambiguities in vector

On Mon, 2007-10-08 at 15:35 +0200, Birgit Lemcke wrote: > Hello James, > > all of your suggestions work very well except of this: > > FemMal <- cbind(FemV1gezählt[2,], MalV1gezählt[2,]) > > colnames(FemMal) <- ("Females", "Males") > Fehler: syntax error The OP missed out c() above, hence the sy

Re: [R] How to remove legend?

On 10/8/07, ONKELINX, Thierry <[EMAIL PROTECTED]> wrote: > Hi Christophe, > > Although I'm a fan of ggplot2, it's lack of documentation is a drawback. > I was having a similar problem and Hadley suggested the "sc$legend <- > FALSE" solution. AFAIK it isn't documented (yet). Have a look at > Hadley'

[R] Applying function to data.frame

Hi, If I have the following data.frame >y time val 1 08:00:05.834 1 2 08:03:13.345 2 3 08:10:12.443 3 > and the following function which converts the time string to the number of milliseconds since midnight > str_to_millis function( s ) { a <- as.numeric( unlist( strsplit(s,":

Re: [R] tcltk scrollbar

kapo coulibaly wrote: > This is approximately what I get: > -- --- > ! !!S! ! !S! > ! !!C! ! !C! > ! !

Re: [R] .C("myHadamaProduct", xx, yy, nr, nc, out=as.double(rep(0.0, n)))

> > Konw anybody as run the function psi on > > > http://www.alglib.net/translator/dl/specialfunctions.psi.csharp.zip > > or > > http://www.alglib.net/translator/dl/specialfunctions.psi.cpp.zip > > using .C with R. > ? > > Thank for your attention. > > Bernardo. > > >> If anybody has an idea what i

Re: [R] heatmap

From: michael watson > > If you just want a coloured representation of your distance > matrix, use image() > > ?image() W/o the parens... > If you want to cluster your original data and show the > original data (not the distances) as a heatmap then use > heatmap(), but you should use so

Re: [R] Applying function to data.frame

try this: > x <- read.table(textConnection(" time val + 1 08:00:05.834 1 + 2 08:03:13.345 2 + 3 08:10:12.443 3"), header=TRUE) > z <- strsplit(as.character(x$time), ":") > newtime <- lapply(z, function(a) as.numeric(a) %*% c(360, 6, 1000)) > x$time <- unlist(newtime) > x

Re: [R] Odp: How to change row value based on previous row value

Petr PIKAL wrote: > Hi > > [EMAIL PROTECTED] napsal dne 08.10.2007 14:24:13: > > >> Hi, >> >> I would like to fix the data in the following data.frame, and I am >> having trouble coding this in R - help please! >> >> >>> x >>> >> AB xy z >> 1 1 10.0 100 1000 1 >>

[R] RTable width sum colum

Hello I need to print table width "colum sum" like XV1V2TOT A235 B167 C46 10 how i can do it by table function? Thanks Alessandra __ R-help@r-project.org mailing list https://stat.e

[R] "R graphical manuals" updated.

Dear all I have updated "R graphical manuals" homepage. http://cged.genes.nig.ac.jp/RGM2/index.php Since obsolute homepage http://bg9.imslab.co.jp/Rhelp/ is no longer maintained, please change your bookmark. Osamu Ogasawara [EMAIL PROTECTED] [EMAIL PROTECTED] ___

[R] semivariogram

I need the most straightforward way to build semivariograms within R i am currently using s-gems software but i would like to even test R Thanks D [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

[R] data.frame's and [2, c(T, F)] in R 2.6.0

In R prior to 2.6.0 with matrizes as well as with data.frames it is possible to do: a[2, c(TRUE, FALSE)] In R 2.6.0 the mixed indexing works only for matrizes. Is this the intention, a bug, or did I get something wrong? Best regards Meinhard Meinhard Ploner Althingstrasse 15 39031 Brun

Re: [R] (no subject)

This is the same issue as FAQ 7.19 How do I produce PNG graphics in batch mode? The advice there may be helpful. -thomas On Fri, 5 Oct 2007, Zhang, Yongqing (NIH/NIA/IRP) [C] wrote: > Hi, All, > > I try to accessing the x11() device over cgi - namely, when trying to > create a graphi

[R] sort time

Hello, I have got the following problem: > times <- c("02.07.2007", "03.07.2007","03.09.2007", "04.07.2007","05.07.2007") > mode(times) [1] "numeric" > tim <- as.character(times) > mode(tim) [1] "character" > sort(times) [1] "02.07.2007" "03.07.2007" "03.09.2007" "04.07.2007" "05.07.2007" I

[R] .C("myHadamaProduct", xx, yy, nr, nc, out=as.double(rep(0.0, n)))

Konw anybody as run the function psi on >> >> >> http://www.alglib.net/translator/dl/specialfunctions.psi.csharp.zip >> >> or >> >> http://www.alglib.net/translator/dl/specialfunctions.psi.cpp.zip >> >> using .C with R. >> ? >> >> Thank for your attention. >> >> Bernardo. >> >> >>> If anybody has a

Re: [R] semivariogram

Have a look at the gstat package and it's documentation. Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie, methodologie en kwaliteitszorg / Sec

Re: [R] Mouse Location on Plots

On 08/10/2007 8:42 AM, Samuel Kemp wrote: > Hi, > > Does anyone know of a function that dynamically updates a variable with the > (x,y) co-ordinates of the mouse pointer on a plot? > > locator and identify only return the values when you click 'stop' on the > graphics device. > > Any help would

Re: [R] data.frame's and [2, c(T, F)] in R 2.6.0

data-ploner.com wrote: > In R prior to 2.6.0 with matrizes as well as with data.frames it is > possible to do: > > a[2, c(TRUE, FALSE)] > > In R 2.6.0 the mixed indexing works only for matrizes. Is this the > intention, a bug, or did I get something wrong? > Best regards > > Meinhard > > It

Re: [R] Odp: How to change row value based on previous row value

Peter Dalgaard <[EMAIL PROTECTED]> napsal dne 08.10.2007 16:31:28: > Petr PIKAL wrote: > > Hi > > > > [EMAIL PROTECTED] napsal dne 08.10.2007 14:24:13: > > > > > >> Hi, > >> > >> I would like to fix the data in the following data.frame, and I am > >> having trouble coding this in R - help please!

Re: [R] Bug or not

Darius Kasiulevičius wrote: > > I have this sample from help. I try use command erase.screen() but in > split mode sreen dont erase. Why? > > Not a bug. From the documentation: 'erase.screen' will appear not to work if the background colour is transparent (as it is by default on m

Re: [R] linewidth in Rgraphviz

Robert Gentleman wrote: > > Hi Ben, >Sorry for the slow reply, but something like this works, > > [...] > Sorry not to have responded sooner. I had actually worked it out for myself (after finding another e-mail exchange on the Bioconductor list that indicated the problem had been

Re: [R] Odp: How to change row value based on previous row value

Petr PIKAL wrote: > Peter Dalgaard <[EMAIL PROTECTED]> napsal dne 08.10.2007 16:31:28: > > >> Petr PIKAL wrote: >> >>> Hi >>> >>> [EMAIL PROTECTED] napsal dne 08.10.2007 14:24:13: >>> >>> >>> Hi, I would like to fix the data in the following data.frame, and I am h

Re: [R] data.frame's and [2, c(T, F)] in R 2.6.0

data-ploner.com wrote: > On Oct 8, 2007, at 4:58 PM, Peter Dalgaard wrote: >> >> >> As far as I can tell, It happens only when drop=TRUE and exactly one >> column is selected using a logical index. > > Thank you. > Maybe it is better letting [] work as in previous R ... ? Well, presumably the code

Re: [R] sort time

Samuel Okoye wrote: > > Hello, I have got the following problem: > [snip] > > Assuming that these are in month/day/year format (which I conclude from your sorted values): times <- c("02.07.2007", "03.07.2007","03.09.2007", "04.07.2007","05.07.2007") times2 <- as.Date(times,format="%m.%d.%Y

Re: [R] sort time

On Mon, 2007-10-08 at 07:26 -0700, Samuel Okoye wrote: > Hello, I have got the following problem: > > times <- c("02.07.2007", "03.07.2007","03.09.2007", > "04.07.2007","05.07.2007") > > mode(times) > [1] "numeric" > > tim <- as.character(times) > > mode(tim) > [1] "character" > > sort(times)

Re: [R] RTable width sum colum

On Mon, 2007-10-08 at 16:18 +0200, Alessandra Marmo wrote: > Hello > I need to print table width "colum sum" like > XV1V2TOT > A235 > B167 > C46 10 > > how i can do it by table function? > > Thanks > > Alessandra Presuming th

[R] pruning of Additive Regression Trees

Hi, Does R provide functions for pruning of Additive Regression Tree Models? Suppose I have grown a large set of trees, say 3000 trees, and now want to find a 'computationally more economic' approximation of the function defined by the trees. In the case of a single tree this is pruning. Is t

[R] Error message on script execution

Dear R_users, I have some troubles with a visual basic application I have recently created. This application automatically creates and executes an R script based on a dataset and settings defined by the user. The .r file (Create_Diagnostic_plots.r) is saved first in a given folder and then cal

Re: [R] Incompatible methods ("-.POSIXt", "Ops.difftime") for "-"

On Mon, 8 Oct 2007, Albrecht, Dr. Stefan (APEP) wrote: Dear all, according to the Help-page of DateTimeClasses {base} I should be able to do time - z with timedate-time objects z a numeric vector (in seconds) or an object of class "difftime". However, on R version 2.6.0 (Windows X

Re: [R] Specify plot size

On Mon, 8 Oct 2007, Eric Thompson wrote: > Have you tried using pdf(), postscript(), or jpeg()? All of these have > arguments to specify the height and width of the device. And so does every other graphics device. Note though that the height/width of the device is not that of a plot (nor necess

Re: [R] R and FDA trials

Ricardo, Tony Rossini (Novartis) provided a very good post on point (3) below (which you may have already read): http://tolstoy.newcastle.edu.au/R/e2/help/07/08/23820.html. Regards, -Cody -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ricardo Pietr

Re: [R] A rebel boxplot question

Thanks to all for the answers --- Thierry : a very interesting answer... but I'm still trying to figure how to install ggplot2 :-) install.packages("ggplot") Warning in install.packages("ggplot2") : argument 'lib' is missing: using '...R/i486-pc-linux-gnu-library/2.5' ?? also installing

[R] as.ordered

Hello Members, I try to convert variables in a data.frame (bract.awn) in the class ordered. str(bract.awn) 'data.frame': 348 obs. of 2 variables: $ bracts.length.relative.to.flower...Min: Factor w/ 4 levels "1","2","3","4": 2 3 3 3 3 2 1 4 3 2 ... $ bract.awn.relative.to.body..

Re: [R] comparing matched proportions using glm

On Mon, 8 Oct 2007, Corry Gellatly wrote: > > Thanks very much for your reply Chuck, I have a quick follow up > question. You mention putting the data into a 2x2x3 for log-linear > model, however my lists have many more than 3 strata, actually > thousands. I am trying to work out whether the propo

Re: [R] RTable width sum colum

If it really is of class "table" (or is an array) or you convert it to such you can use addmargins. You may need to read ?addmargins carefully. On 10/8/07, Alessandra Marmo <[EMAIL PROTECTED]> wrote: > Hello > I need to print table width "colum sum" like > XV1V2TOT > A23

Re: [R] Applying function to data.frame

Your function needs to be able to work on vector input, e.g. transform( y, time=Vectorize(str_to_millis)( as.character(time) ) ) or just do this: library(chron) y$time <- 24 * 60 * 60 * 1000 * as.numeric(times(y$time)) On 10/8/07, Rees, David <[EMAIL PROTECTED]> wrote: > Hi, > > If I have the

[R] Reading distance matrices

Hi all, I've been having trouble reading in distance matrices (values in lower triangle only) that have been created by other programs. They load with NA values in the upper triangle. Is there an option to read the matrices properly? Thanks, Muri _

[R] estfun & df

Hello EVERYONE, I need an URGENT help from you please! How can I see the "estfun" (empirical estimating function) and "df" (degree of freedom) from the following mixed-model please? (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) Many thanks in advance for your kind help.

Re: [R] data.frame's and [2, c(T, F)] in R 2.6.0

On Oct 8, 2007, at 4:58 PM, Peter Dalgaard wrote: > data-ploner.com wrote: >> In R prior to 2.6.0 with matrizes as well as with data.frames it is >> possible to do: >> >> a[2, c(TRUE, FALSE)] >> >> In R 2.6.0 the mixed indexing works only for matrizes. Is this the >> intention, a bug, or did I get

[R] problem with times

Hello, I have got the following problem: > times <- c("02.07.2007", "03.07.2007","03.09.2007", "04.07.2007","05.07.2007") > mode(times) [1] "numeric" > tim <- as.character(times) > mode(tim) [1] "character" > sort(times) [1] "02.07.2007" "03.07.2007" "03.09.2007" "04.07.2007" "05.07.2007" Is it p

[R] CompetingRiskFrailty modeling

Hi all can anybody give info what I do wrong when modeling a survival function with the CompetingRiskFrailty package in the following way? "accumulate", "reduce" and "closed" are dummy vectors with 0s indicating no event and 1s indicating the respective event in the appropriate column. ac

Re: [R] Reading distance matrices

Muri If it is lower triangular, then the upper triangle is empty. What do you expect it to have? We don't have any information on how you're reading in your data. But, read.table() has an argument for na.strings that you should look at. > -Original Message- > From: [EMAIL PROTECTED] > [m

Re: [R] estfun & df

On Mon, 8 Oct 2007, Abdus Sattar wrote: > Hello EVERYONE, > > I need an URGENT help from you please! This type of requests is not considered to be very polite, please have a look at the posting guide. > How can I see the "estfun" (empirical estimating function) I guess (because you are not tel

Re: [R] estfun & d

On 10/8/07, Abdus Sattar <[EMAIL PROTECTED]> wrote: > Hello EVERYONE, > I need an URGENT help from you please! > How can I see the "estfun" (empirical estimating function) and "df" (degree > of freedom) from the following mixed-model please? > (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sle

[R] Dice simulation: Getting rep to re-evaluate sample()?

I'm trying to get R to simulate the sum of the values on 10 fair dice (yes, it's related to a homework problem, but is not the problem itself). I tried to do this: > rep(sum(sample(1:6,100,replace=T)), times=10) [1] 341 341 341 341 341 341 341 341 341 341 and noticed that sum(sample()) seems to

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

See ?replicate which I think is what you are after. Chuck On Mon, 8 Oct 2007, Zembower, Kevin wrote: > I'm trying to get R to simulate the sum of the values on 10 fair dice > (yes, it's related to a homework problem, but is not the problem > itself). I tried to do this: >> rep(sum(s

[R] another try with Rmpi on Linux

Dear R People: Here are my latest attempts for Rmpi: This is the list of locations for mpi.h: /home/faculty/hodgess/mpich2-1.0.5p4/src/include/mpi.h /home/faculty/hodgess/include/mpi.h /home/faculty/hodgess/lam-7.1.4/romio/adio/sgi/mpi3.1/mpi.h /home/faculty/hodgess/lam-7.1.4/share/incl

Re: [R] problem with times

On Mon, 2007-10-08 at 07:38 -0700, Samuel Okoye wrote: > Hello, I have got the following problem: > > times <- c("02.07.2007", "03.07.2007","03.09.2007", > > "04.07.2007","05.07.2007") > > mode(times) > [1] "numeric" > > tim <- as.character(times) > > mode(tim) > [1] "character" > > sort(times) >

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

See ?replicate, e.g. replicate(sum(sample(1:6,100,replace=TRUE)), n=10) Function arguments are (in gneral) evaluated as if they were evaluated when passed to the function. On Mon, 8 Oct 2007, Zembower, Kevin wrote: > I'm trying to get R to simulate the sum of the values on 10 fair dice > (yes,

Re: [R] estfun & df

On 10/8/07, Achim Zeileis <[EMAIL PROTECTED]> wrote: > On Mon, 8 Oct 2007, Abdus Sattar wrote: > > > Hello EVERYONE, > > > > I need an URGENT help from you please! > > This type of requests is not considered to be very polite, please have a > look at the posting guide. > > > How can I see the "estf

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

Thanks so much, Chuck and Mark. Here's my script to simulate 10,000 rolls of 100 fair dice to demonstrate their conformity to a normal curve: > x<-replicate(1, sum(sample(1:6,100,replace=T))) > sdx<-sd(x) > sdx [1] 17.13966 > meanx<-mean(x) > meanx [1] 350.0451 > hist(x, freq=FALSE) > curve(dno

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

rep( sum(sample(1:6,100,replace=T)), times=10) Read carefully what you wrote. You asked R to sum a single sample, then make ten copies of the sum. You need to take ten samples, then sum each. apply() would be helpful. -Original Message- I'm trying to get R to simulate the sum

Re: [R] print Text on device

Look at the textplot function in the gplots package. Does that do what you want? -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of

[R] coxph() command design and data setup

Hello all: I'm attempting to run a Cox proportional hazards function on survival data from insects and I have a few questions. My current command that I'm using to call the model is as follows (using coxph() from the survival library): coxph(Surv(day, censor) ~ treatment + room + chamber %in% tr

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

Kevin Zembower wrote: > > I'm trying to get R to simulate the sum of the values on 10 fair dice > (yes, it's related to a homework problem, but is not the problem > itself). I tried to do this: > > rep(sum(sample(1:6,100,replace=T)), times=10) > [1] 341 341 341 341 341 341 341 341 341 341 > rep(s

Re: [R] as.ordered

Hello, (Warning. This might not be the most complete or elegant solution ...) If you want a sorted dataframe: look here for example http://tolstoy.newcastle.edu.au/R/help/05/02/12391.html To convert the factors from a data frame, you have to call as.ordered for each factor separately (not for

Re: [R] Simulate data based on correlation coefficient

Similar questions have been asked here before. See http://finzi.psych.upenn.edu/R/Rhelp02a/archive/110042.html for one possible answer (your case will be simpler than this). -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 >

Re: [R] as.ordered

Hello, (Warning. This might not be the most complete or elegant solution ...) If you want a sorted dataframe: look here for example http://tolstoy.newcastle.edu.au/R/help/05/02/12391.html To convert the factors from a data frame, you have to call as.ordered for each factor separately (not for t

Re: [R] Tart charts

Michael, > Try this alternative: > > # from http://research.microsoft.com/users/lamport/pubs/hair.pdf > hairsex <- matrix( >c(46, 45, 13, 12, > 1, 101, 0, 20), 2, 4, byrow=TRUE) > dimnames(hairsex) <- list("Gender"=c("Female", "Male"), > "Hair color"=c("Blond", "Brown", "Re

Re: [R] Mouse Location on Plots

Another option (which should work on all platforms) is to use tcltk and bind a function to the "Motion" event. The playSudoku function in the sudoku package has an example of this (as well as using getGraphicsEvent). (wow, could that package actually be useful?) -- Gregory (Greg) L. Snow Ph.D.

[R] get more leads with these healthcare lists

Certified MDs in the United States 788,437 in total – 17,880 emails Many popular specialties like Emergency Medicine, Plastic Surgery, OBGYN, Oncology, Pediatrics and more Many unique fields like 'medical school attended' and 'location of residency training' Special Price: $321 *** FR

[R] variance explained by each term in a GAM

Hello fellow R's, I do apologize if this is a basic question. I'm doing some GAMs using the mgcv package, and I am wondering what is the most appropriate way to determine how much of the variability in the dependent variable is explained by each term in the model. The information provided by

Re: [R] R and FDA trials

Hi All, I'm excited to see that R-Core has released the document "R: Regulatory Compliance and Validation Issues A Guidance Document for the Use of R in Regulated Clinical Trial Environments" (http://www.r-project.org/doc/R-FDA.pdf). I know it represents a great deal of effort on the part of

Re: [R] Mouse Location on Plots

Greg Snow wrote: > > Another option (which should work on all platforms) is to use tcltk and > bind a function to the "Motion" event. The playSudoku function in > the sudoku package has an example of this (as well as using > getGraphicsEvent). > > (wow, could that package actually be useful?)

[R] embedFonts rotates figure

Dear All, Consider the following code: pdf(file="figure.pdf",family="URWPalladio") curve(dlnorm(x,0,1.5),0,10,xlim=c(0,10),ylim=c(0,0.85),axes=F,xlab="",ylab="f") segments(exp(-1.5^2),0,exp(-1.5^2),dlnorm(exp(-1.5^2),0,1.5),lty="dashed") segments(1,0,1,dlnorm(1,0,1.5),lty="dashed") segments(exp(

[R] lmList function in lme4

I don¹t understand the following error message in lme4. > lmList(formula = oaas.bin ~ oaas.prob | subject, data = inductionPropRemiAll.df, family = binomial) Error in FUN(X[[1L]], ...) : unused argument(s) (family = function (link = "logit") > lmList(formula = oaas.bin ~ oaas.prob | subject, dat

Re: [R] embedFonts rotates figure

Paul Smith wrote: > Dear All, > > Consider the following code: > > pdf(file="figure.pdf",family="URWPalladio") > > curve(dlnorm(x,0,1.5),0,10,xlim=c(0,10),ylim=c(0,0.85),axes=F,xlab="",ylab="f") > segments(exp(-1.5^2),0,exp(-1.5^2),dlnorm(exp(-1.5^2),0,1.5),lty="dashed") > segments(1,0,1,dlnorm(1,0

Re: [R] Getting intervals for within-group standard errors for each group using nlme and varIdent

Just to follow-up on this for the archive and anyone else who is interested in this problem. I ended up solving this by running lme seperately on subsets of the data (one for group=='A' and one for group=='B'), the intervals for the statistics computed both ways were within rounding of each other,

[R] "xlab" and "ylab" in clusplot

Hi there, How can I change the "xlab" and "ylab" in clusplot instead of "Component 1" and "Component 2"? It always gives me the following error message. Error in plot.default(x1[, 1], x1[, 2], xlim = xlim, ylim = ylim, xlab = "Component 1", : formal argument "ylab" matched by multi

  1   2   >