[R] Words appear to be bolded in the PDF output

Hi all, I have written a note near each of my graphs using mtext. mtext(text,side=1,line=4,cex=0.5,adj=0) Then I have exported the graphs as a PDF file. pdf(file=name,paper='a4',width=7.27,height=10.69) The mtext appears OK in R. But it looks like it is bolded in the PDF file. http://n4.nabble.

[R] : Followup Error in library(gplots) : there is no package called 'gplots' : 64bits

Dear James and everybody interested, James'idea was correct : the problem was related to the 64bits version of R. Trying on the same computer a 32bits version of R did not produce this 'gplots' error (but many other conflicts linked to 64 / 32 bits) when installing RODBC package. I will ins

Re: [R] How to embed italic Greek letters in a eps file?

Hi, Another option might be the tikzDevice package, which uses LaTeX to process the fonts, library(tikzDevice) tikz(standAlone=T) plot(1,1, type = 'n') mtext(side = 3, line = 2, "$\\mu$") dev.off() ## system("/usr/texbin/pdflatex Rplots.tex") HTH, baptiste On 20 April 2010 07:30, Prof Brian R

Re: [R] How to embed italic Greek letters in a eps file?

This is discussed in detail on the help page for postcript! Default PostScript with the 14 standard fonts does not cover Greek at all (but it does cover mathematical symbols such as mu, in a different typeface). See the section 'Encodings': to display Greek you need to have an encoding which

Re: [R] Help: coxph() in {survival} package

Hi Xin, to answer your question: say you have your regression reg = coxph(... Then you can assess the log likelihood by reg$loglik This vector contains typically two values, the null log likelihood of the restricted model that excludes the fixed effects and the log likelihood of the unrestrict

Re: [R] Help: coxph() in {survival} package

Xin Ge wrote: Hi All, I'm runnning coxph() on 90 different datasets (in a loop). 1. I'm wondering how can I get "log-likelihood" value from coxph() output. Currently I can only see following: Likelihood ratio test = Wald test = Score (logrank) test = 2. Once I have likelihood value, I would

[R] Help: coxph() in {survival} package

Hi All, I'm runnning coxph() on 90 different datasets (in a loop). 1. I'm wondering how can I get "log-likelihood" value from coxph() output. Currently I can only see following: Likelihood ratio test = Wald test = Score (logrank) test = 2. Once I have likelihood value, I would like to extract

Re: [R] converting a zoo or an xts to a data frame

If z is a zoo series as.data.frame(z) converts it to a data frame. If that is not what you are asking please clarify. Also read the last line to every message to r-help asking for complete self contained code. Yours has undefined variables. Also what is the purpose of the code? Normally its unn

Re: [R] bwplot puts the bars in the wrong place

James, It's actually the bars for hour 3 (which don't exist) that are missing. You still need the 'drop.unused.levels=FALSE' and if you make 'OnHour' into a factor then you won't need the 'horizontal=FALSE'. -Peter Ehlers On 2010-04-19 8:27, James Rome wrote: You were right about the gdf$. So

[R] converting a zoo or an xts to a data frame

Dear R People: I have the following code that I use to convert a monthly zoo object to a data.frame, and it works perfectly: library(tseries) z <- get.hist.quote(instrument=inst1, start=start1,end=end1, quote=quot1,comp = "m") y <- as.ts(aggregate(z, as.yearmon, tail, 1)

Re: [R] Problem installing RODBC

Marc Schwartz wrote: On Apr 19, 2010, at 7:02 PM, David Scott wrote: I have a problem with some missing headers when trying to install RODBC on my linux box. install.packages("RODBC", + lib="/usr/local/david/R") trying URL 'http://cran.stat.auckland.ac.nz/src/contrib/RODBC

Re: [R] How to set proxy settings for R

Hi: One alternative to sneaking around the proxy restrictions is to load R onto a jump drive and run it off of that. I know of several instructors who are using that approach with their students for similar reasons. Unless you have a compelling reason to install every single package, you should be

[R] Solved: DCA and centroid analysis in Vegan

 Sorry guys,  I just have it figured out. But thanks anyway. Fan [[alternative HTML version deleted]] __ 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/po

[R] DCA and centroid analysis in Vegan

 Hi everyone, I am performing some correspondence analysis. In the Vegan package, all centroid analysis was performed based on CCA. I tried to use the sample data set (dune & dune.env) to run a centroid analysis on its DCA output but failed.  So my question is: Is it possible to perform centr

Re: [R] Problem installing RODBC

On Apr 19, 2010, at 7:02 PM, David Scott wrote: > I have a problem with some missing headers when trying to install RODBC on my > linux box. > > > install.packages("RODBC", > + lib="/usr/local/david/R") > trying URL 'http://cran.stat.auckland.ac.nz/src/contrib/RODBC_1.3-1.tar.g

Re: [R] Type-I v/s Type-III Sum-Of-Squares in ANOVA

Dear Kevin and Daniel, This question comes up, in one form or another, with some frequency. Here's a closely related reply from last month: . Regards, John John Fox Senator William McMaster Profess

[R] Problem installing RODBC

I have a problem with some missing headers when trying to install RODBC on my linux box. > install.packages("RODBC", + lib="/usr/local/david/R") trying URL 'http://cran.stat.auckland.ac.nz/src/contrib/RODBC_1.3-1.tar.gz' Content type 'application/x-gzip' length 990220 bytes (96

Re: [R] Week number calculation

If dd is of Date class then this code finds the number of days since May 1st of the same year and then does an integer divide by 7 plus 1 so that the if dd lies in May 1 - May 7 it returns 1 as the week number, etc. dd <- as.Date("2000-06-01") # sample input as.numeric(dd - as.Date(sub("-.*", "-05

[R] label the bars by the percentage values in the conditional histogram?

HI, Dear R community, HOW to LABEL the bars by the percentage values in the conditional histogram: Thanks so much! h<- sample(1:14, 319, rep=T) c<- sample(1:14, 608, rep=T) n<- sample(1:14, 1140, rep=T) vt<-c(h, c, n) ta<-rep(c("h", "c", "n"), c(319, 608, 1140)) to<-data.frame(vt,ta) library(

[R] How to embed italic Greek letters in a eps file?

Hi, I need to add on a plot text containing italic Greek characters using the function mtext (i.e. I cannot use Hershey vectors). The characters are nicely displayed when the file is saved as png but not when saved as eps. See code below as example: #postscript('test.eps') png('test.png') p

Re: [R] Efficiency of C Compiler in "R CMD SHLIB"

do you know how to increase the optimization level of gcc in R? -- View this message in context: http://n4.nabble.com/Efficiency-of-C-Compiler-in-R-CMD-SHLIB-tp1934429p2016646.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-p

[R] Week number calculation

R community, How could I extract the week number from a date vector (in class Date) such that week numbering (week 1...2...) begins (May 01) and ends (October 31) on the same specific dates each year? Thank you, Mike __ R-help@r-project.org mailing

Re: [R] nls minimum factor error

Hi Dennis, Thank you for your suggestion. Firstly, the below minimum step size is only a warning because I told it to be (warnOnly=T), otherwise it's thrown as an error. I'm not sure if this is completely arbitrary. About segmented regression, I have tried to use it, but it doesn't seem to work ri

Re: [R] Overlay of barchart and xyplot

Huapeng, It sounds like you have different y scales on the two plots you want to overlay? If so you could try doubleYScale() from the latticeExtra package. Please post a minimal, reproducible example (with data-generating code that can be pasted in to R) if you want more help. Regards -Felix O

[R] What is mclust up to? Different clusters found if x and y interchanged

Hello All... I gave a task to my students that involved using mclust to look for clusters in some bivariate data of isotopes vs various mining locations. They discovered something I didn¹t expect; the data (called tur) is appended below. p <- qplot(x = dD, y = dCu65, data = tur, color = mine) pr

Re: [R] Type-I v/s Type-III Sum-Of-Squares in ANOVA

i believe John Fox offers another solution in his book. Kevin Daniel Wollschlaeger wrote: > * On Mo, 1. Mar 2010, Ista Zahn wrote: > > > I've posted a short explanation about this at > > http://yourpsyche.org/miscellaneous that you might find helpful. I'm a > > As someone who's also stru

Re: [R] Select single column, preserve name?

On 4/19/2010 3:15 PM, Duncan Elkins wrote: Hi, list- I've got a large list of multi-column data and I'd like to filter it according to a threshold, such as, "show me the values that are above 0.4 for each line." For instance, if the line of data were: line v1 v2v3 v4

Re: [R] Select single column, preserve name?

Hi Try line[,line<=-0.7,drop=F] drop=F,keep the dimensions. On Mon, Apr 19, 2010 at 5:15 PM, Duncan Elkins wrote: > Hi, list- >I've got a large list of multi-column data and I'd like to filter > it according to a threshold, such as, "show me the values that are > above 0.4 for each line."

[R] Select single column, preserve name?

Hi, list-    I've got a large list of multi-column data and I'd like to filter it according to a threshold, such as, "show me the values that are above 0.4 for each line." For instance, if the line of data were: > line v1 v2v3 v4 v5 v6v7v8 v9 1 -0.32 0.66 -0.

Re: [R] Tinn-R

Robert Ruser gmail.com> writes: > I want to use the free distribution of R (R REvolution 3.2) [this is a little bit funny -- is the standard (CRAN) distribution not free? Or did you mean something else?] > and Tinn-R > editor as well. Unfortunately they don't cooperate. In Tinn-R > commands

Re: [R] nls for piecewise linear regression not converging to least square

On Mon, 19 Apr 2010, Karen Chang Liu wrote: > Hi R experts, > > I'm trying to use nls() for a piecewise linear regression with the first > slope constrained to 0. There are 10 data points and when it does converge > the second slope is almost always over estimated for some reason. I have > many se

[R] Overlay of barchart and xyplot

Hello R Folks, I am new to R. I have been struggling to overlay a barchart with a xyplot together on one plot but did not get this worked out. Any help and idea are greatly appreciated. I attached R scripts for barchart and xyplot below and also data I used. What I am trying to do is just to put

Re: [R] plotting RR, 95% CI as table and figure in same plot

Thanks to Thomas and Christos for helpful suggestions. The forestplot (in package rmeta) suggestion seems to work fairly well for me, though does require a bit of fiddling (no complaints, obviously using it for a different purpose than it was written). Below is an example using a slightly "h

Re: [R] How to set proxy settings for R

Also, the FAQ suggests using the alternative internet2.dll by starting R with the flag --internet2 If you start R from a desktop icon, you can add the --internet flag to the target line (right click, properties) e.g. "C:\Program Files\R\R-2.8.1\bin\Rgui.exe" --internet2 see http://cran.r-projec

Re: [R] Tinn-R

Thank you very much - it really works. Maybe it's not so useful as Tinn-R but is sufficient. 2010/4/19 Tal Galili : > Consider trying > notepad++ > with > NppToR > That's what I use (it also works with the REvolution distribution) > > > Contact > Details:-

Re: [R] plotting RR, 95% CI as table and figure in same plot

ggplot2 should work (resize to get the plot to the dimensions you need for the paper) library(Hmisc) library(pscl) library(ggplot2) ## data data("bioChemists", package = "pscl") fm_pois <- glm(art ~ ., data = bioChemists, family = poisson) summary(fm_pois) ### pull out rate-ratios and 95%

[R] Huge data sets and RAM problems

Dear all, This is the first time I am sending mail to the mailing list, so I hope I do not make a mistake... The last months I have been working on my MSc thesis project on performing data mining techniques on user logs of a software-as-a-service application. The main problem I am experiencing i

[R] nls minimum factor error

Hi, I have a small dataset that I'm fitting a segmented regression using nls on. I get a step below minimum factor error, which I presume is because residual sum of square is still "not small enough" when steps in the parameter space is already below specified/default value. However, when I l

[R] nls for piecewise linear regression not converging to least square

Hi R experts, I'm trying to use nls() for a piecewise linear regression with the first slope constrained to 0. There are 10 data points and when it does converge the second slope is almost always over estimated for some reason. I have many sets of these 10-point datasets that I need to do. Th

Re: [R] nls for piecewise linear regression not converging to least square

Try reparameterizing: nlmod2 <- nls(y2 ~ pmax(1/p, (x2 - xint)), data = dat, start = list(xint = 40.49782, p = 1), trace = TRUE, alg = "plinear") On Mon, Apr 19, 2010 at 11:32 AM, Karen Chang Liu wrote: > Hi R experts, > > I'm trying to use nls() for a piecewise linear regression with the f

[R] densCols: what are the computed densities and how to create a legend

Hi, I'm using the densCols function for a scatterplot and cannot figure out 1) how to extract the computed densities, and 2) how to create a legend based that represents the upper and lower ranges of the densities. For example: movers.den <- densCols(move$x, move$y) table(movers.den) #08306B

Re: [R] Problems with labels and scaling in star diagrams

For number 2, do the scaling yourself so that all values are between 0 and 1, then use scale=FALSE in the call to stars. For number 3 try stardata[1,,drop=FALSE] Don't have a good suggestion for 1 (though you could look at the code to see where the legend is plotted and move that code to the re

Re: [R] comparing attitudes of 2 groups / likert scales?

thanks a lot for your help!! I better get on with reading / working now! -- View this message in context: http://n4.nabble.com/comparing-attitudes-of-2-groups-likert-scales-tp2015738p2016398.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] glmer with non integer weights

Le lundi 19 avril 2010 à 03:00 -0800, Kay Cichini a écrit : > hi emmanuel, > > thanks a lot for your extensive answer. > do you think using the asin(sqrt()) transf. can be justified for publishing > prurpose or do i have to expect criticism. Hmmm ... depends of your reviewers. But if an half-asl

Re: [R] How to set proxy settings for R

Try this: setInternet2() chooseCRANmirror() On Mon, Apr 19, 2010 at 1:00 PM, danda wrote: > > Dear All, > > I would like to run R on my computer (with win xp on it) at work bu the > proxy restrictions of the university don't let me download the packages or > to connect to a cran mirror, I us

[R] nls minimum factor error

Hi, I have a small dataset that I'm fitting a segmented regression using nls on. I get a step below minimum factor error, which I presume is because residual sum of square is still "not small enough" when steps in the parameter space is already below specified/default value. However, when I look a

[R] nls for piecewise linear regression not converging to least square

Hi R experts, I'm trying to use nls() for a piecewise linear regression with the first slope constrained to 0. There are 10 data points and when it does converge the second slope is almost always over estimated for some reason. I have many sets of these 10-point datasets that I need to do. The fol

[R] How to set proxy settings for R

Dear All, I would like to run R on my computer (with win xp on it) at work bu the proxy restrictions of the university don't let me download the packages or to connect to a cran mirror, I usually get this message: > chooseCRANmirror() Warning message: In open.connection(con, "r") : unable to c

[R] help in output file

HI, Dear R-community, I AM using the following codes to grow tree and plot tree: # Classification Tree with rpart library(rpart) pdf(file="/home/cdu/changbin/dimer_tree.pdf") # grow tree fit.dimer <- rpart(outcome ~ ., method="class", data=p.dimer[,2:402]) plotcp(fit.dimer) # visualize cross-v

Re: [R] Type-I v/s Type-III Sum-Of-Squares in ANOVA

* On Mo, 1. Mar 2010, Ista Zahn wrote: > I've posted a short explanation about this at > http://yourpsyche.org/miscellaneous that you might find helpful. I'm a As someone who's also struggled with the "type X sum of squares" topic, I like the idea to completely walk through a numerical example an

Re: [R] Grouping rows of data by day

Try this: aggregate(DF[c('data1','data2')], list(gsub('\\..*', '', DF$time)), FUN = sum) On Mon, Apr 19, 2010 at 12:00 PM, jennyed wrote: > > Hi all, > > I have a set of data in hourly time steps with each row identified as > time data column1 data column2 > 1

Re: [R] How to pass a list of parameters into a function

On Mon, Apr 19, 2010 at 5:58 PM, Gene Leynes wrote: > Does anyone know how to pass a list of parameters into a function? > > > for example: > > somefun=function(x1,x2,x3,x4,x5,x6,x7,x8,x9){ >    ans=x1+x2+x3+x4+x5+x6+x7+x8+x9 >    return(ans) > } > > somefun(1,2,3,4,5,6,7,8,9) > > # I would like t

Re: [R] How to pass a list of parameters into a function

Try this: do.call(somefun, c(x1 = 1, x2 = 2, as.list(temp))) On Mon, Apr 19, 2010 at 1:58 PM, Gene Leynes > wrote: > Does anyone know how to pass a list of parameters into a function? > > > for example: > > somefun=function(x1,x2,x3,x4,x5,x6,x7,x8,x9){ >ans=x1+x2+x3+x4+x5+x6+x7+x8+x9 >

Re: [R] S4-based package failure: setGeneric example#1

Nevermind, I figured out my problem after looking at the packS4 package. I didn't realize that package.skeleton would build invalid R code. The warning hints at that, but I interpreted the warning to mean that the package.skeleton-generated code should be *edited*, and none of my edits gave me a w

Re: [R] BRugs

Perhaps a better idea is to ask on a BUGS mailing list. BRugs is just an interface to OpenBUGS and is not involved in handling the BUGS language I'd also suggest to strat trying your problem witht BRugs but in OpenBUGS directly in order to avoid confusion caused by the interface. Best wis

[R] How to pass a list of parameters into a function

Does anyone know how to pass a list of parameters into a function? for example: somefun=function(x1,x2,x3,x4,x5,x6,x7,x8,x9){ ans=x1+x2+x3+x4+x5+x6+x7+x8+x9 return(ans) } somefun(1,2,3,4,5,6,7,8,9) # I would like this to work: temp=c(x3=3,x4=4,x5=5,x6=6,x7=7,x8=8,x9=9) somefun(x1=1,x2=

Re: [R] dataframe

http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-do-I-convert-factors-to-numeric_003f -- View this message in context: http://n4.nabble.com/dataframe-tp2015650p2016230.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.o

Re: [R] plotting RR, 95% CI as table and figure in same plot

You could try the forestplot() function in rmeta, or the original grid code on which it is based, http://www.stat.auckland.ac.nz/~paul/RGraphics/chapter1.html -thomas On Mon, 19 Apr 2010, David Atkins wrote: Hi all-- I am in the process of helping colleagues write up a ms in which w

Re: [R] How to make a boxplot with exclusion of certain groups

On 4/19/2010 12:21 PM, josef.kar...@phila.gov wrote: > This seems like a simple thing, but I have been stuck for some time. My > data has 2 columns. Column 1 is the value, and column 2 is the Site where > data was collected. Column 2 contains 7 different Sites (i.e. groups). I > am only inte

[R] plotting RR, 95% CI as table and figure in same plot

Hi all-- I am in the process of helping colleagues write up a ms in which we fit zero-inflated Poisson models. I would prefer plotting the rate ratios and 95% CI (as I've found Gelman and others convincing about plotting tables...), but our journals usually like the numbers themselves. Thu

[R] Identifying names of matrix columns shared by many matrices

Greetings R-Geniuses, What is the most efficient way to handle the problem described below? Thanks Marsh Feldman Problem description: Each U.S. state has its own matrix. The rows are dates, the columns are industries, and each cell contains total statewide employment at the given time a

[R] How to make a boxplot with exclusion of certain groups

This seems like a simple thing, but I have been stuck for some time. My data has 2 columns. Column 1 is the value, and column 2 is the Site where data was collected. Column 2 contains 7 different Sites (i.e. groups). I am only interested in showing 3 groups on a single boxplot. I have trie

Re: [R] Follow up on installing formatR...

On Mon, Apr 19, 2010 at 5:28 AM, Brian Lunergan wrote: > Good morning folks: > > Made a second go at installing this and succeeded, but with some strange > behaviours along the way. First the system back story. > > My only guess is that installing RGtk2 from source is only possible if you have th

Re: [R] Tinn-R

Consider trying notepad++ with NppToR That's what I use (it also works with the REvolution distribution) Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatis

Re: [R] getting column's main

Try col=max.col(data) This should give you the index of the column with the max value. To get to the final result, combine this with names(dataframe)[col] to get the name of the column with the maximum value. HTH Jannis AuriDUL schrieb: Hello. I have data of potatoes production in EU

Re: [R] BRugs

Thanks for the reply Bob, but it still does not work, you see. I ran this model, just with the main effects and it ran fine. n=length(bi.bmi) Lgen=2 Lrace=5 Lagegp=13 Lstra=15 Lpsu=2 bi.bmi.model=function(){ # likelihood for(i in 1:n){ bi.bmi[i]~ dbern(p[i]) logit(p[i])<- a0 + a

Re: [R] stupid regexp question

Use fixed = TRUE to turn off interpretation of special characters: gsub("W.m^{-2}", "W/m2", "abc W.m^{-2} xyz", fixed = TRUE) 2010/4/19 servet ahmet çizmeli : > Hello, > > I have a stupid regexp question. I have a large data frame of strings. I would > like to convert all occurences of : > > "W.

Re: [R] stupid regexp question

On Apr 19, 2010, at 11:39 AM, servet ahmet çizmeli wrote: Hello, I have a stupid regexp question. I have a large data frame of strings. I would like to convert all occurences of : "W.m^{-2}" to "W/m2" I make the following test : gsub(glob2rx("W.m^{-2}"), "W/m2", "W.m^{-2}") Two probl

[R] Tinn-R

Hello, I want to use the free distribution of R (R REvolution 3.2) and Tinn-R editor as well. Unfortunately they don't cooperate. In Tinn-R commands: send selection, send line etc. don't work. Do you have any idea how to resolve this problem? Best, Robert _

Re: [R] var.test

On Apr 18, 2010, at 4:55 PM, anon anon wrote: Hello, I'm using var.test to do a simple F-test for equality of variances. I think I'm missing something small here: m<-rnorm(10,sd=1) n<-rnorm(5,sd=1) var.test(m,n) F test to compare two variances data: m and n F = 13.7438, num df = 9,

[R] Drawing a line with misc3d

Hi the list, I would like to draw some lines with misc3d. I find a lot of tools to draw surfaces, but nothing for simple line... Is it possible? Note that I know that it is possible to draw lines with rgl (using lines3d), but I need to do it with misc3d to export the drawing in .asy format. A

[R] stupid regexp question

Hello, I have a stupid regexp question. I have a large data frame of strings. I would like to convert all occurences of : "W.m^{-2}" to "W/m2" I make the following test : gsub(glob2rx("W.m^{-2}"), "W/m2", "W.m^{-2}") but it does not seem to work. I don't know how to do it otherwise as I coul

Re: [R] Using split and then unsplit

Hi David, do.call worked perfectly but do you have an idea why unsplit wouldn't work in that case, is that because the number of rows changed? bc when the number didn't change unsplit worked - Anna Lippel -- View this message in context: http://n4.nabble.com/Using-split-and-then-unsplit-tp2

Re: [R] Using split and then unsplit

On Apr 19, 2010, at 11:06 AM, anna wrote: Hello everyone, I use the split function splitting with the f function on a 3 columns and more than 100 000 rows data frame. Once it's split I have a list of data frames still with 3 columns and n rows. I manipulate those list elements and get a

Re: [R] Formatting data, adding column names, use reshape, a newbie question

On Mon, Apr 19, 2010 at 5:13 AM, Paul Rigor (ucla) wrote: > Hi all, > I'm an R novice. > > I have data that's already formatted as "molten" that reshape should be able > to work with. For example, the following was read in with > read.csv(filename,sep=" ", header=FALSE) > >      V1               V

Re: [R] comparing attitudes of 2 groups / likert scales?

David Winsemius wrote: > > > If you are thinking of using that quote, you might want to check the > spelling of his name. My memory is van Belle. > > Sorry, I thought I had corrected that before mailing. @BOOK{vanBelle2002, title = {Statistical rules of thumb}, publisher = {Wiley seri

Re: [R] Using split and then unsplit

here is an alternative that I just found to join my data frames with rbind: result <- do.call("rbind", myList) It worked perfectly but I still don't understand why unsplit wouldn't work... - Anna Lippel -- View this message in context: http://n4.nabble.com/Using-split-and-then-unsplit-tp20

Re: [R] comparing attitudes of 2 groups / likert scales?

On Apr 19, 2010, at 10:08 AM, Tal Galili wrote: Good luck in your work, The simple solution would be to run many non-paired wilcox on all the 20 questions (the way Dieter suggested). In which case, make sure to adjust for multiple comparisons. Read about it, and see: ?p.adjust If you ha

Re: [R] Inline Package: void vs return type functions

Many Thanks for your help Best, Sergio -- View this message in context: http://n4.nabble.com/Inline-Package-void-vs-return-type-functions-tp1838423p2015898.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

[R] Grouping rows of data by day

Hi all, I have a set of data in hourly time steps with each row identified as time data column1 data column2 1 1.042 1.083

[R] Using split and then unsplit

Hello everyone, I use the split function splitting with the f function on a 3 columns and more than 100 000 rows data frame. Once it's split I have a list of data frames still with 3 columns and n rows. I manipulate those list elements and get a list of data frames still with 3 columns but less ro

Re: [R] ecdf

The OP wrote me privately to say that the errant documantation was at: http://lib.stat.cmu.edu/S/Harrell/help/Hmisc/html/ecdf.html That is a rather old bit of information. It dates back to a time when Frank's address was at the the University of Virginia. In 2003 he moved to Vanderbilt

Re: [R] ecdf

On Apr 19, 2010, at 9:12 AM, Downey, Patrick wrote: Hi Thierry, That worked perfectly. Thanks for the suggestion. For reference, in the documentation, it What is "it"? never lists {Hmisc}'s function as starting with E instead of e. Every instance of the documentation in the r help syste

Re: [R] multiple variables pointing to single dataframe?

If you only need to retrieve x by referring to x2 and you don`t have to modify x via x2 then this works: > x <- Orange > makeActiveBinding("x2", function() x, .GlobalEnv) > x$age <- 50 > head(x2) Tree age circumference 11 5030 21 5058 31 5087 4

Re: [R] merge

Hi If the columns has the same name but different values in them then you shall either decide which one to keep yourself or you shall keep both. If they have same name and same values you could select only those which names do not match. names(data1) %in% names(data2) can select which names m

Re: [R] multiple variables pointing to single dataframe?

Hi, On Mon, Apr 19, 2010 at 10:15 AM, Alex Bryant wrote: > Hi, for example: > >> x <- Orange >> x2 <- x >> x[1,]$age <- 50 >> x2[1,] >  Tree age circumference > 1    1 118            30 > > I would like a way for x2 to also reference the modified x data frame without > having to reassign x2http:

Re: [R] merge

Maybe this what you are looking for lines1 <- "CODPROD N1 N3 N4 23 3 55 4 24 5 67 36 25 3 73 24" df1 <- read.table(textConnection(lines1),header=TRUE) lines2 <-"CODPROD N1 N2 30 34 45 45 0 78 65 0 56" df2 <- read.table(textConnection(lines2),header=TRUE) merge(df1, df2, by = inter

Re: [R] Unwanted boxes in legend

Try border=c(0,0,1,0). -tgs On Mon, Apr 19, 2010 at 4:21 AM, Steve Murray wrote: > > Dear all, > > Thanks for the response, however I'm getting the following error message > when I execute the legend command using the 'border' argument: > > Error in legend(10, par("usr")[4], c("A", "B", : > un

Re: [R] comparing attitudes of 2 groups / likert scales?

Good luck in your work, The simple solution would be to run many non-paired wilcox on all the 20 questions (the way Dieter suggested). In which case, make sure to adjust for multiple comparisons. Read about it, and see: ?p.adjust If you have some questions you can merge (by a simple mean of them)

Re: [R] multiple variables pointing to single dataframe?

Hi, for example: > x <- Orange > x2 <- x > x[1,]$age <- 50 > x2[1,] Tree age circumference 11 11830 I would like a way for x2 to also reference the modified x data frame without having to reassign x2mailto:petr.pi...@precheza.cz] Sent: Monday, April 19, 2010 3:18 AM To: Alex B

Re: [R] smart way to turn a vector into a matrix

Hi Erik, what if I do some manipulations on the new list created with split and I want to come back to its initial form? I saw you have the unsplit function but I get errors. What I am doing is that I use as the "f" function as in the help the same initial f I used to split but I get the followin

Re: [R] merge

What do you want to get? And what exactly did you do? Your question isn't very clear. Sarah On Mon, Apr 19, 2010 at 7:59 AM, n.via...@libero.it wrote: > I have a problem with the merge function. > I have to merge two big dataframes which  look like the following example.The > problems is that

[R] Extracting the coefficients of each local polynomial from loess()

Hello dear R users and Prof. Brian Ripley, I am searching for a way to extract the estimated coeffiicents of each local polynomial at given x from loess(). After searching and asking

Re: [R] ecdf

Hi Thierry, That worked perfectly. Thanks for the suggestion. For reference, in the documentation, it never lists {Hmisc}'s function as starting with E instead of e. I don't know who's in charge of documentation, but that should probably be corrected. Thanks again. -Mitch -Original Message

Re: [R] fit a deterministic function to observed data

Plotting y vs. x: plot(y ~ x) the graph seems to be flattening out at x = 0 at a level of around y = 500 so lets look at: plot(500-y ~ x) This curve is moving up rapidly so lets take the log to flatten it out: plot(log(500-y) ~ x) That looks quite linear so log(500-y)

Re: [R] ecdf

R is case sensitive. ecdf() is in the stats package, Ecdf() is in Hmisc. So you want Ecdf(x,what='1-F') Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek team Biometrie & Kwaliteitszorg Gaverstraat 4 95

[R] ecdf

Hello, I'd like to plot an empirical cumulative distribution function, except instead of the fraction of values < x, I'd like the fraction of values > x. I think this can be done using the ecdf function in {Hmisc}. I installed the package and loaded it. However, when following the example given

Re: [R] help

Hi, > > Thank you for your help. > I try your function 'combvec2' .but , it gives me an error :Erreur dans rep(1, c2) : argument 'times' incorrect > this function 'combvec ' takes any number of inputs on Matlab. So, you can > take more than two matrix. > > > The help of this function 'comb

Re: [R] help

Hi, Thank you for your help. this function 'combvec ' takes any number of inputs on Matlab. So, you can take more than two matrix. The help of this function 'combvec' is like this on Matlab: >> help combvec COMBVEC Create all combinations of vectors. Syntax combvec(a1,a2,...)

Re: [R] Comparing data frames

Thank you all for your help and suggestions. L On Sun, Apr 18, 2010 at 8:51 PM, Tal Galili wrote: > Would: > ?merge > Work for you ? > > > Contact > Details:--- > Contact me: tal.gal...@gmail.com |  972-52-7275845 > Read me: www

  1   2   >