[R] How two compare two matrixes

Dear all I have two 10*10 matrixes and I would like to compare theirs contents. By the word content I mean to check visually (not with any mathematical formulation) how similar are the contents. I also know edit that prints my matrix in the scree but still one edit blocks the prompt to launch a

Re: [R] advice on classes/methods/extending classes

Hi Erin, One good option would be the official manual: http://cran.r-project.org/doc/manuals/R-exts.html It depends to an extent, I think, on what types of methods you would like to work with and use. FWIW, I have and really enjoy both S Programming by Venables & Ripley (mostly S3 methods and ge

[R] advice on classes/methods/extending classes

Dear R People: What is the best way to learn about classes, methods, extending classes, and namespaces, please? I know a bit about classes, but would like to learn much more. Thanks in advance for any advice! Sincerely, Erin -- Erin Hodgess Associate Professor Department of Computer and Math

[R] questions about using loop, while and next

Hello R helpers, I have a quick question about loop and next In my loop, I have some random generation of data, but if the data doesn't meet some condition, then I want it to go next, and generate data again for next round. # just an example.. # i want to generate the data again, if the sum is s

Re: [R] PCA - scores

Shari, Josh partly answered your question, but his example did not include rotation because he took out just one factor. Try: require(psych) mt.pc <- principal(mtcars,3,scores=TRUE) #this gives you the varimax rotated first 3 principal components #pc.scores <- mt.pc$scores #here are

Re: [R] how to store lme/lmer fit result

Hi All, I'm experiencing difficulties in saving a model fit by gls(). Basically if I just save() a gls object "gls.fit" in my workspace into an .RData file and later reload this .RData file, I get error when running script such as summary(gls.fit)$coefficients Dr. Bates's reply (quoted below)

Re: [R] Floating points and floor() ?

Hi Michael, In floating point calculation, 1.0-.9 is not exactly 0.1. This is easily seen by subtracting. > (1.0-.9)-0.1 [1] -2.775558e-17 > (1.0-.9)==0.1 [1] FALSE David is right, you can't "correct" this. You can only compensate by taking care that you never, ever test whether 2 FP numbers

Re: [R] Floating points and floor() ?

You are assuming that 0.1 is exactly represented. Since it is not, either 0.1 or 1-0.9 will be less than 0.1. As to what you should do about it... that depends on what you are trying to accomplish. Whatever it is, you need to use some other approach. Often the alternative involves using integer

Re: [R] R: Help center

On Mar 3, 2011, at 5:10 PM, > wrote: I also would like to stop the mailing, without unsubscribing myself from the help center. How can I proceed?Thanks All of the options for your subscription are changed on the same webpage. -- David. Lorenza __

Re: [R] Floating points and floor() ?

On Mar 3, 2011, at 8:23 PM, Folkes, Michael wrote: Perhaps somebody could clarify for me if the following is a floating point matter or otherwise, and how am I to correct for it? floor(100*.1) [1] 10 100*(1.0-.9) [1] 10 floor(100*(1-0.9)) [1] 9 Yes. It's a "floating point matter".

Re: [R] Problems with a function warning

Dear Felipe, I did not have any difficulty with it using: R version 2.12.1 (2010-12-16) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5] LC_TIME=Englis

[R] Problems with a function warning

Hello. I have the following funtion: fechasEntrega = function(FechasEntrega,fecha){ if(length(which(FechasEntrega0){ tkmessageBox(title = "Error en Fecha de Valoracion",message="Hay una fecha de entrega anterior a la fecha de valoracion. Todas las fechas de entrega deben ser posteriores

Re: [R] What am I doing wrong with this loop ?

Never mind Billgot it. Always seems to happen this way. Can't figure something out. Post to the site and wham, 5 min later (after posting), it's all clear. Oh well, thanks for the tips -- View this message in context: http://r.789695.n4.nabble.com/What-am-I-doing-wrong-with-this-loop-tp3332

Re: [R] Plotting Mean in plotting degree distribution

library(igraph) G <- erdos.renyi.game(1000, 1/1000) # a random graph dd1 = degree.distribution(G) plot(dd1, xlab = "degree", ylab="frequency") abline(h = mean(dd1)) # the mean would be a horizontal line On Thursday, March 3, 2011 at 4:43 PM, kparamas wrote: > Hi, > > I am plotting degree dis

Re: [R] What am I doing wrong with this loop ?

Bill, I addressed the first issue with the data frames and length(x). But my loops still isn't working. More importantly, you commented that I should be using if(...) ... else ... rather than ifelse(.,.,). Please help me understand the difference. I thought ifelse was just a faster way of doing if

[R] Floating points and floor() ?

Perhaps somebody could clarify for me if the following is a floating point matter or otherwise, and how am I to correct for it? > floor(100*.1) [1] 10 > 100*(1.0-.9) [1] 10 > floor(100*(1-0.9)) [1] 9 Thanks! Michael ___ Michael Folkes Salmon

Re: [R] R usage survey

No. That's not answering the question. ALL surveys are for collecting information. The substantive issue is what purpose do you have in seeking this information in the first place and what are you going to do with it when you get it? Do you have some commercial purpose in mind? If so, what i

Re: [R] lattice custom axis function -- right side margins

I went to try your suggestion, and the label appeared without the vjust argument. I usually run R within emacs using ESS; I happened to restart my emacs earlier today. That's the only thing I can think of that I changed. Next time I'll try running R outside of emacs before asking for help. Many

Re: [R] PCA - scores

Hi Shari, Yes, please look at the documentation for principal. You can access this (assuming you have loaded psych) by typing at the console: ?principal note the logical argument "scores". Here is a small example: ## require(psych) require(GPArotation) dat <- prin

Re: [R] lattice custom axis function -- right side margins

Timothy W. Hilton wrote: To clarify the trouble I'm having with ylab.right, I am not getting an error message; the right-side label just does not appear on the plot. Maybe this is mac-specific. On Windows, the label shows up just fine. You might be able to make it appear by adjusting the 'vjust

Re: [R] lattice: How to increase space between ticks and labels of z-axis?

Marius Hofert wrote: Dear expeRts, How can I increase the space between the ticks and the labels in the wireframe plot below? I tried some variations with par.settings=list(..) but it just didn't work. Marius, I tried setting the 'distance' parameter, but that was less than satisfactory. On

Re: [R] error in saved .csv

why didn't work? What was the error, or why did you assume there was an error? Can you clarify. Sent from my iPad On Mar 1, 2011, at 17:50, Tamas Barjak wrote: > Yes, the format is incorrect. I have already tried the write.table, but it > didn't work. > > > 2011/3/1 jim holtman > I am no

Re: [R] lattice custom axis function -- right side margins

To clarify the trouble I'm having with ylab.right, I am not getting an error message; the right-side label just does not appear on the plot. -Tim > On Thu, Mar 3, 2011 at 1:50 PM, Timothy W. Hilton wrote: > > > Many thanks, Richard -- the position argument does exactly what I > > needed. I'm no

[R] Plotting Mean in plotting degree distribution

Hi, I am plotting degree distribution of a graph using the function, library(igraph) dd1 = degree.distribution(G) plot(dd1, xlab = "degree", ylab="frequency") I would like to plot the mean of the distribution as a vertical line in the attached plot. Please let me know how to do this. Thanks, K

[R] R: Help center

I also would like to stop the mailing, without unsubscribing myself from the help center. How can I proceed?Thanks Lorenza Da: r-help-boun...@r-project.org [r-help-boun...@r-project.org] per conto di David Winsemius [dwinsem...@comcast.net] Inviato: giov

Re: [R] creating a count variable in R

You can probably simplify this if you can assume that the dates are in sorted order. Here is a way of doing it even if the days are in arbitrary order. The count refers to the number of times that this date has appeared so far in the sequence. con <- textConnection(" 01/01/2011 01/01/2011 02/

Re: [R] creating a count variable in R

On Mar 3, 2011, at 3:58 PM, JonC wrote: Hi R helpers, I'm trying to create a count in R , but as there is no retain function like in SAS I'm running into difficulties. Your data is not cut-pastable as presented but this should work: > dfrm$count_var <- ave(as.numeric(dfrm$Date_var), dfrm

Re: [R] Scatter plot with multiple data sets

Hi Joe, The easiest option will be to combine all 6 datasets (at least the variables you want to use in your scatter plot), and then create another variable that indicates to which group the observations belong. Here is a small example of what you might do once your data are all together (obvious

Re: [R] Regression with many independent variables

Thanks for getting back to me so quickly greg. Im not quite sure how to do what you just said, is there an example that you can show? I understand how to create the string with a formula in it but im not sure how to loop through the pairs of variables? How do I first get these 2way interaction var

Re: [R] creating a count variable in R

Use cumsum() to count the change points: > Date_var <- as.Date(rep(c("2011-02-04","2011-02-07","2011-01-29"), c(2,3,1))) > data.frame(Date_var, count=cumsum(c(TRUE, Date_var[-1]!=Date_var[-length(Date_var)]))) Date_var count 1 2011-02-04 1 2 2011-02-04 1 3 2011-02-07 2 4 2011-02-07

Re: [R] Regression with many independent variables

Thanks greg, that formula was exactly what I was looking for. Except now when I run it on my data I get the following error: "Error in model.matrix.default(mt, mf, contrasts) : cannot allocate vector of length 2043479998" I know there are probably many 2-way interactions that are zero so I thou

Re: [R] Help center

You need to unsubscribe using the same method you use to subscribe go to the webpage for the list and log in with the password you set. If all you wnat to stop is the mailing, you can do so without unsubscribing. On Mar 3, 2011, at 3:09 PM, Bobby Lee wrote: Could you please take my em

[R] PCA - scores

I am running a PCA, but would like to rotate my data and limit the number of factors that are analyzed. I can do this using the "principal" command from the psych package [principal(my.data, nfactors=3,rotate="varimax")], but the issue is that this does not report scores for the Principal

[R] Normalising proportional binomial data that has a set top value

I've got some data in proportional format which has a set maximum value of 0.5 and I want to know how to normalise it. I've calculated a lateralisation index to determine whether an organism deviates from an equal number of left and right turns within a trial (as opposed to measuring the proportio

[R] creating a count variable in R

Hi R helpers, I'm trying to create a count in R , but as there is no retain function like in SAS I'm running into difficulties. I have the following : Date_var and wish to obtain Date_var Count_var 01/01/2011

[R] Help center

Could you please take my email off the help center? Thank you very much. [[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

[R] Interpreting the coefficient of an interaction between continuous variables in a regression model

Hello, my question is triggered by an actual model I am running, but I will pose it as a very general question with a hypothetical example. Take the following regression model: I have a binomial dependent variable "Happiness", whose two values are 0 (=unhappy) and 1 (=happy). My two independent

Re: [R] Calling a function to store values

Hi, I tried this approach. cdata_1_1_3 <- readData(cData1,1,3) cdata_2_4_7 <- readData(cData2,4,7) Its still not working. The result I get is NULL. Thanks, Kumar On Thu, Mar 3, 2011 at 9:50 AM, David Winsemius wrote: > > On Mar 3, 2011, at 10:12 AM, kparamas wrote: > > Hi, >> >> I am calling

Re: [R] Ploting Histogram with Y axis is percentage of sample for each bin

On Mar 3, 2011, at 12:44 PM, leinwand wrote: I'm trying to do something very simple... I wan to plot a histogram where the y axis represent the percentage of the total sample that each bin represents. I know how to plot a histogram with the counts and density... but can't find anything

Re: [R] Probabilities greather than 1 in HIST

Hi Martin, No, I get frequencies greather than 1. See the code: #TLC uniform n=30 mi=1; mx=6 nrep=1000 xbar=rep(0,nrep) for (i in 1:nrep) {xbar[i]=mean(runif(n,min=mi,max=mx))} hist(xbar, freq=FALSE) #hist(xbar,prob=TRUE,breaks="Sturges",xlim=c(1,6),main=paste("n =",n), xlab="Média",

[R] embed latex beamer sans serif default font into R plot

Hello, I have seen instructions on how to embed Latex Computer Modern fonts into R, but these are the default serif fonts. I am trying to embed the default font used for Latex beamer (theme Warsaw), which is a sans serif font and may be the default LateX Computer Modern sans serif font.Does

Re: [R] Developing a web crawler / R "webkit" or something similar? [off topic]

On 03/03/2011 08:07 AM, Mike Marchywka wrote: Date: Thu, 3 Mar 2011 01:22:44 -0800 From: antuj...@gmail.com To: r-help@r-project.org Subject: [R] Developing a web crawler Hi, I wish to develop a web crawler in R. I have been using the functionalities available under the RCurl package. I

[R] Ploting Histogram with Y axis is percentage of sample for each bin

I'm trying to do something very simple... I wan to plot a histogram where the y axis represent the percentage of the total sample that each bin represents. I know how to plot a histogram with the counts and density... but can't find anything that gives me perenct of sample on the y axis. Any h

[R] Scatter plot with multiple data sets

Hi, I have multiple (6) data sets which I would like to plot together on one scatter graph. The reason they are all separate is that I require a different symbol to be plotted for each set. Could somebody advise on how to do this? Many thanks, Joe -- View this message in context: http://r.78969

Re: [R] Ordering several histograms

Em 3/3/2011 12:00, djbirdnerd escreveu: Hallo everyone, I want to evaluate the change of the distribution for several size classes. How can i order these separate histograms with the same y-axis along a common x-axis according to their size classes. It would like it to look a bit like this (http

Re: [R] Calling a function to store values

On Mar 3, 2011, at 1:49 PM, Kumaraguru Paramasivam wrote: > Hi, > > I tried this approach. > cdata_1_1_3 <- readData(cData1,1,3) > cdata_2_4_7 <- readData(cData2,4,7) > > Its still not working. The result I get is NULL. Then the function (which you have not provided) is faulty, but there is no

[R] 'merge' function creating duplicate columns names in the output

The "merge" command is creating duplicate column names in a dataframe that is the result of the merge. The following is the 'merge' command: x <- merge(invType , allocSlots , by.x = 'index' , by.y = 'indx' , all.x = TRUE ) The 'invType' dataframe was the result of a previous

Re: [R] plot, y-axis, uneven scale???

On Thu, Mar 3, 2011 at 4:09 PM, yan liu wrote: > Hello, > > I have a question about the y-axis of plots.  Actually I had about 60 > values.  About 80 percent of these values are less than 0.2, then the other > 20 percent values are more than 4,max is 10.  So when I plot these values > together, th

[R] lattice: How to increase space between ticks and labels of z-axis?

Dear expeRts, How can I increase the space between the ticks and the labels in the wireframe plot below? I tried some variations with par.settings=list(..) but it just didn't work. Many thanks, Marius library(lattice) u <- seq(0, 1, length.out=20) grid <- expand.grid(x=u, y=u) z <- apply(g

Re: [R] plot, y-axis, uneven scale???

You probably want to use the gap.plot function in the plotrix package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On

Re: [R] Regression with many independent variables

What you might need to do is create a character string with your formula in it (looping through pairs of variables and using paste or sprint) then convert that to a formula using the as.formula function. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@

[R] plot, y-axis, uneven scale???

Hello, I have a question about the y-axis of plots. Actually I had about 60 values. About 80 percent of these values are less than 0.2, then the other 20 percent values are more than 4,max is 10. So when I plot these values together, the y-axis's range will go 0 to 10, and my major values (80%

Re: [R] Problem on flexmix when trying to apply signature developed in one model to a new sample

Jon, if I did understand you correctly the problem is that you did not specify the newdata argument in posterior() correctly. You need to specify it in way such that evaluating the formula uses the correct object. If you have a matrix as dependent variable, you have to use a list which contai

Re: [R] Creating a weighted sample - Help

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of LouiseS > Sent: Thursday, March 03, 2011 5:21 AM > To: r-help@r-project.org > Subject: Re: [R] Creating a weighted sample - Help > > Hi > > Thanks for responses. The sample I h

[R] mailing list submission

-- "Gul dusununen gul bahcesi, diken dusunen dikenlik olur"Mevlana [[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/

Re: [R] lattice custom axis function -- right side margins

Many thanks, Richard -- the position argument does exactly what I needed. I'm not having any luck with the ylab.right argument. My R and lattice are up to date (below); is there something else I should check? Thanks for the help, Tim > sessionInfo() R version 2.12.2 (2011-02-25) Platform: i38

Re: [R] Probabilities greather than 1 in HIST

Hi: Here's an example: x <- rnorm(300) hist(x, breaks = 15, yaxt = 'n', ylab = 'Relative frequency') axis(2, at = seq(0, 50, by = 10), labels = round(seq(0, 50, by = 10)/length(x), 3)) HTH, Dennis On Thu, Mar 3, 2011 at 8:29 AM, jpmaroco wrote: > Dear David, > > Thanks for your prompt reply.

Re: [R] lattice custom axis function -- right side margins

print(my_plot(example_data, ylab.right=expression(e==mc^2)), position=c(0,0,.95,1)) You will need a recent R version for the ylab.right argument. On Thu, Mar 3, 2011 at 12:52 PM, Timothy W. Hilton wrote: > Dear R help list, > [[alternative HTML version deleted]] ___

[R] lattice custom axis function -- right side margins

Dear R help list, I have a plot with two different vertical scales that I want to display on either side of the plot. It's quite similar to the Fahrenheit-Centigrade example in the examples section of the documentation for axis.default. The right-side axis is clipped off, though, and I haven't

Re: [R] Error in model.frame.default

Hi: You need the second variable in D1 to be named fGRAZE - the variable names in the newdata data frame (D1) have to be the same as the variable names on the RHS of the model formula, in this case L.AREA and fGRAZE. HTH, Dennis On Thu, Mar 3, 2011 at 7:43 AM, Heike Schmitz wrote: > Dear R- Com

Re: [R] Probabilities greather than 1 in HIST

Hi, Does xx = rnorm(100) hist(xx,freq=FALSE) curve(dnorm,add=TRUE) give you what you want? Martyn -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of David Winsemius Sent: 03 March 2011 17:00 To: jpmaroco Cc: r-help@r-project.org Sub

Re: [R] Probabilities greather than 1 in HIST

On Mar 3, 2011, at 10:29 AM, jpmaroco wrote: Dear David, Thanks for your prompt reply. I see your point. But how can I get a histogram with relative frequencies? If I use plot(xhist,yhist) I get absolute frequencies in the Y axis. In my earlier reply I meant to type yaxt="n". -- Dav

Re: [R] Probabilities greather than 1 in HIST

On Mar 3, 2011, at 10:29 AM, jpmaroco wrote: Dear David, Thanks for your prompt reply. I see your point. But how can I get a histogram with relative frequencies? If I use plot(xhist,yhist) I get absolute frequencies in the Y axis. I do not know of any simple setting to do what you wa

Re: [R] Calling a function to store values

On Mar 3, 2011, at 10:12 AM, kparamas wrote: Hi, I am calling a function with different arguments to read different files and want the results to be stored in different matrices. Ex: cData1 = NULL cData2 = NULL readData = function(cData, start, end) { cData = //reads from the file

[R] Calling a function to store values

Hi, I am calling a function with different arguments to read different files and want the results to be stored in different matrices. Ex: cData1 = NULL cData2 = NULL readData = function(cData, start, end) { cData = //reads from the file } I am calling the functions using readData(cData

Re: [R] Probabilities greather than 1 in HIST

Dear David, Thanks for your prompt reply. I see your point. But how can I get a histogram with relative frequencies? If I use >plot(xhist,yhist) I get absolute frequencies in the Y axis. Best, Joao From: David Winsemius [via R] [mailto:ml-node+619-1272092771-215...@n4.nabble.com] S

Re: [R] Greek character and R

or even without using 'paste' plot(1,1,main=bquote(delta^13~'C Station'~.(i))) Am 03.03.2011 16:45, schrieb William Dunlap: > > > Bill Dunlap > Spotfire, TIBCO Software > wdunlap tibco.com > >> -Original Message- >> From: r-help-boun...@r-project.org >> [mailto:r-help-boun...@r-pro

Re: [R] Zero truncated Poisson distribution & R2WinBUGS

On Thu, 3 Mar 2011, Timothée wrote: Hi, I have a very similar problem... In some sites, counts data>0 In the other sites, counts = 0 I already applied zero inflated models with the zero-trick (Martin et al 2005 in Ecology letters), but I would like to use truncated distributions (Poisson and neg

Re: [R] Developing a web crawler

Hi The book whose companion website is here deals with many of the things you need for a web crawler, and assignment "other 5" on that site () is

Re: [R] Multivariate Granger Causality Tests

Beware that causality can only be inferred using information that extends far beyond the data at hand. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Multivariate-Granger-Causality-Tests-tp3332968p65

[R] Error in model.frame.default

Dear R- Community, to learn i reanalysed some data provided and analysed by Zuur et. al. in their book "Mixed effect models and Extensions in Ecology with R". When i run the last command i get a warning message i dont understand. Loyn<- read.table(file = "loyn.txt",header = TRUE) Loyn$L.AREA

Re: [R] Ordering several histograms

You could get close with the ggplot2 package using the function facet_grid or facet_wrap, but each histogram would be on a separate x-axis Scott On Thursday, March 3, 2011 at 8:00 AM, djbirdnerd wrote: > Hallo everyone, > > I want to evaluate the change of the distribution for several size clas

[R] Recodifying a factor due to results in lm

Dear community, I'm doing a lm. In the independent variables I've got a categorical one. Here is its histogram: http://r.789695.n4.nabble.com/file/n638/altitude.png I did this regression: lmeo2.52f <- lm(dat82$IncAltuDom ~ dat82$hdom2+log(dat82$CV)+ dat82$CA+ dat82$FCC+ factor(dat82$IdAlt

Re: [R] Pairwise T-Tests and Dunnett's Test (possibly using multcomp)

Hi Josh,   Thanks for your reply. You're right about letting R's method dispatch system choose the method for summary and inserting more spaces in the code. I was just messing around with the code in the former case and forgot to change it back. As far as not having enough whitespace goes, I've

Re: [R] Probabilities greather than 1 in HIST

On Mar 3, 2011, at 8:03 AM, jpmaroco wrote: Dear all, I am a newbie in R and could not find help on this problem. I am trying to plot an histogram with probabilities in the y axis. This is the code I am using: #TLC uniform n=30 mi=1; mx=6 nrep=1000 xbar=rep(0,nrep) for (i in 1:nrep) {xbar

[R] R Course***Advanced Statistical Modeling in R by XLSolutions Corp

Our New York City course: Advanced Statistical Modelling in R/S-PLUS is coming up on March 14-15 http://www.xlsolutions-corp.com/coursedetail.asp?id=13 by XLSolutions Corp Email s...@xlsolutions-corp.com Regards - Sue Turner Senior Account Manager XLSolutions Corporation North Americ

[R] Analytical Hierarchical Process (AHP) in R?

Hi Does anybody know anything about an implementation of an AHP in R or any other open source tool? I googled but could not find anything Thanks, Rainer -- NEW GERMAN FAX NUMBER!!! Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany) Centr

Re: [R] Greek character and R

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 Duncan Murdoch > Sent: Thursday, March 03, 2011 7:30 AM > To: Filoche > Cc: r-help@r-project.org > Subject: Re: [R] Gre

Re: [R] Probabilities greather than 1 in HIST

If you read ?hist, you will answer your own question. The issue in your code is the parameter prob = T, which does nothing. By default, hist reports density. -- Jonathan P. Daily Technician - USGS Leetown Science Center 11649 Leetown Road Kearneysville WV, 2543

Re: [R] R usage survey

Just out of curiosity : It is possible to get name, employer name, location, usage information and academic background details when searching for R users on LinkedIn and the many R related groups there. Does this also provide potential opportunities for misuse and "outrageous" analyses, since

Re: [R] Greek character and R

Eval it. This works at my house: plot(0) title(eval(parse(text=paste("expression(paste(delta^13,'C Station ',",i,"))" -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Filoche Sent: Thursday, March 03, 2011 9:39 AM To: r-help@r-

Re: [R] Greek character and R

On 11-03-03 8:15 AM, Filoche wrote: Dear R users. In a loop, I set the title of my graph with : mytitle = expression(paste(delta^13,'C Station ', i) title(mytitle) However, instead of using value of i, it will literally use "i" character. Any one know the way to concatenate th

Re: [R] Applying function to multiple data

Hi, It might not be the best approach, but here is what I would do. ## 1) If you have your data in 3 different data.frames: #create a named list where each element is one of your data.frame list_df <- vector(mode="list", length=3) names(list_df) <- c("Bank", "Corporate", "Sovereign")

Re: [R] Greek character and R

Hi and ty for the answer. However, it's not working. It will print "expression(d13C Station 1)". Thank for any help, Phil -- View this message in context: http://r.789695.n4.nabble.com/Greek-character-and-R-tp304p467.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] Zero truncated Poisson distribution & R2WinBUGS

Hi, I have a very similar problem... In some sites, counts data>0 In the other sites, counts = 0 I already applied zero inflated models with the zero-trick (Martin et al 2005 in Ecology letters), but I would like to use truncated distributions (Poisson and negative binomial) to model my counts in o

Re: [R] R usage survey

Hi Rex and useRs, The purpose of the survey has been mentioned on the survey link goo.gl/jw1ig but I will also reproduce it here. - Geographical distribution of R users - Application areas where R is being used - Supporting technology being used along with R - Academic background distribution of R

Re: [R] Developing a web crawler / R "webkit" or something similar?

> Date: Thu, 3 Mar 2011 01:22:44 -0800 > From: antuj...@gmail.com > To: r-help@r-project.org > Subject: [R] Developing a web crawler > > Hi, > > I wish to develop a web crawler in R. I have been using the functionalities > available under the RCurl package. > I am able to extract the html co

[R] Probabilities greather than 1 in HIST

Dear all, I am a newbie in R and could not find help on this problem. I am trying to plot an histogram with probabilities in the y axis. This is the code I am using: #TLC uniform n=30 mi=1; mx=6 nrep=1000 xbar=rep(0,nrep) for (i in 1:nrep) {xbar[i]=mean(runif(n,min=mi,max=mx))} hist(xbar,prob=TRUE

[R] Ordering several histograms

Hallo everyone, I want to evaluate the change of the distribution for several size classes. How can i order these separate histograms with the same y-axis along a common x-axis according to their size classes. It would like it to look a bit like this (http://addictedtor.free.fr/graphiques/RGraphGa

Re: [R] sqlFetch (RODBC) question

Hi Kenn, This is discussed in the package vignette, section 7. Best, Ista On Thu, Mar 3, 2011 at 2:55 PM, Kenn Konstabel wrote: > Dear all, > > I've used RODBC a lot to read in files created in MS excel and access but > found a strange problem today: a variable in my data file contained both > nu

[R] sqlFetch (RODBC) question

Dear all, I've used RODBC a lot to read in files created in MS excel and access but found a strange problem today: a variable in my data file contained both numbers and text; sqlFetch would set text within a row of numbers to NA; but if first 5 or 6 rows would be text then all numbers would be rea

Re: [R] Multivariate Granger Causality Tests

Dear Hazzard I. Petzev, you might find causality() in the package vars useful. Best, Bernhard > -Ursprüngliche Nachricht- > Von: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] Im Auftrag von hazzard > Gesendet: Donnerstag, 3. März 2011 10:07 > An: r-help@r-project

Re: [R] R usage survey

Harsh, "Suitably analyzed" for whose purposes? One man's "suitable" is another's "outrageous". That's why people want to see the gowns at the Oscars. Under what auspices are you conducting this survey? What do you intend to do with it? You don't give any assurance that the results you post w

Re: [R] Greek character and R

mytitle = parse(text=paste("expression(paste(delta^13,'C Station ',",i,"))")) title(mytitle) -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Filoche Sent: Thursday, March 03, 2011 8:16 AM To: r-help@r-project.org Subject: [R] Greek

Re: [R] Developing a web crawler

On Mar 3, 2011, at 4:22 AM, antujsrv wrote: > > I wish to develop a web crawler in R. As Rex said, there are faster languages, but R string processing got better due to the stringr package (R Journal 2010-2). When Hadley is done with it, it will be like having it all in R! -- Alexy _

[R] m out of n bootstrap

Can anyone confirm the formula for the m out of n bootstrap variance estimator? rq.boot applies a deflation factor directly to the bootstrap estimates. Presumably, the SE of the estimate of interest is then taken to be the SD of the deflated estimates. I have read Bickel's and others' papers on

Re: [R] Developing a web crawler

Perl seems like a 10x better choice for the task, but try looking at the examples in ?strsplit to get started. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of antujsrv Sent: Thursday, March 03, 2011 4:23 AM To: r-help@r-project.org

[R] vector("integer", length) : vector size specified is too large

Good day to the R community, I am interested to run the plot.count() function in the "untb" package. My script is as follows:- > library(untb) > Community1 <- > structure(c(371,167,119,78,74,53,50,31,28,25,20,19,19,17,13,12,12,10, > 9,9,8,8,7,7,7,7,6,6,6,6,5,5,5,5,4,4,4,3,3,3,2,2,2,2,2

Re: [R] Creating a weighted sample - Help

Hi Thanks for responses. The sample I have taken is a random sample from H, I, J and K. The further analysis I want to do is all around bad debt rates so it could be (H/H+I)*100 = Bad rate percentage also population stability calculations that are all related to credit scoring. I want to be abl

[R] Greek character and R

Dear R users. In a loop, I set the title of my graph with : mytitle = expression(paste(delta^13,'C Station ', i) title(mytitle) However, instead of using value of i, it will literally use "i" character. Any one know the way to concatenate the value of i to the mathematical expre

  1   2   >