[R] R DCOMServer Manager 2.0

2009-01-29 Thread Smita.KUMAR
Hi All, I am from the packaging team in RBS. We use R DCOMServer Manager 2.0, this application uses RGUI 2.3.0 as one of its components. We are planning to upgrade to the latest version of RGUI 2.8.0 but found that R DCOMServer manager is not compatible with the latest version. Can you direct me

[R] plotting lines with missing data for x values

2009-01-29 Thread Gareth Campbell
I have some data (REE plots - geochemistry) where I have values 1:14 for the x axis, but have no data for some x values. Here for example, let's say that I don't have data for x=2,5,8. So x<-1:14 y<-c(4, NA, 5, 9, NA, 3.4, 8, NA, 19, 22, 12, 14, 15.3, 15) if I plot the data plot(x,y) and then

[R] How to generate a xyplot with multiple panels using an empty data frame ?

2009-01-29 Thread Daniel Kornhauser
Hi: I am coding some interactive interface using an xyplot from lattice. When a user clicks on a panel of the xyplot, a simulation is executed and the resulting data is plotted on the corresponding panel. The problem is that I start with an empty data frame and only fill it as the user request dat

Re: [R] How do I get my IT department to "bless" R?

2009-01-29 Thread Warren Young
Michael Olschimke wrote: Could you please share a link to the NY Times article? http://www.nytimes.com/2009/01/07/technology/business-computing/07program.html __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] How do I get my IT department to "bless" R?

2009-01-29 Thread Warren Young
Erik Iverson wrote: First, you should try to figure out why they would not want you to run R, so you can address those reasons specifically. Reasons imply reasoning. It's usually the case that decisions like this are made on an emotional basis, not a rational one. "All of my bus

Re: [R] How do I get my IT department to "bless" R?

2009-01-29 Thread Michael Olschimke
Could you please share a link to the NY Times article? Is it about OSS in general or specific to R? Thanks Michael Michael Bibo wrote: > Daniel Viar gmail.com> writes: > > >> I currently use R at work "under the radar", but there's a chance I >> could loose that access. I'd like to get our

Re: [R] Q about how to use Anova.mlm

2009-01-29 Thread pgseye
Thanks a lot for that John - really helpful. I generated some random numbers and seem to be able to get it to work, so that's great. One thing - it's come up with a 'Type III' test and given me a few warnings. What's the difference between Type II and Type III tests (if there's some basic guide y

Re: [R] How do I get my IT department to "bless" R?

2009-01-29 Thread Michael Bibo
Daniel Viar gmail.com> writes: > > I currently use R at work "under the radar", but there's a chance I > could loose that access. I'd like to get our company to feel > comfortable with open source and R in particular. Does anyone have > any experience with their company's IT department and man

Re: [R] Matlab inv() and R solve() differences

2009-01-29 Thread Berwin A Turlach
G'day all, On Thu, 29 Jan 2009 19:24:40 -0700 Kingsford Jones wrote: > I suppose the solution is unstable because x is ill-conditioned: While, as you show, x is ill-conditioned, I do not believe that this is serious enough to explain the differences that Pat sees between MATLAB and R. In fac

Re: [R] How do I get my IT department to "bless" R?

2009-01-29 Thread David M Smith
On Thu, Jan 29, 2009 at 2:29 PM, Daniel Viar wrote: > I currently use R at work "under the radar", but there's a chance I > could loose that access. I'd like to get our company to feel > comfortable with open source and R in particular. Does anyone have > any experience with their company's IT d

Re: [R] How do I get my IT department to "bless" R?

2009-01-29 Thread Warren Young
Daniel Viar wrote: I'd like to get our company to feel comfortable with open source Anyone still denying, here in 2009, that open source offers serious business value is a dinosaur, doomed to extinction. Their cerebella have calcified. The balance tipped a decade ago. Just like the real d

Re: [R] Q about how to use Anova.mlm

2009-01-29 Thread John Fox
Dear Paul, I noticed a typo in my response and some poor formatting in the email message; please see below: > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of John Fox > Sent: January-29-09 9:23 PM > To: 'pgseye' > Cc: r-help@r-p

Re: [R] Matlab inv() and R solve() differences

2009-01-29 Thread Kingsford Jones
I suppose the solution is unstable because x is ill-conditioned: > x [,1] [,2] [,3] [,4] [1,] 0.133 0.254 -0.214 0.116 [2,] 0.254 0.623 -0.674 0.139 [3,] -0.214 -0.674 0.910 0.011 [4,] 0.116 0.139 0.011 0.180 > cor(x) [,1] [,2] [,3] [,4] [1,] 1.000

Re: [R] Q about how to use Anova.mlm

2009-01-29 Thread John Fox
Dear Paul, First, to fit a multivariate linear model to your data, you'll have to rearrange the data from "long" format (with one observation per replicate) to "wide" format (with one observation per subject). If your data are in the data frame Data, then you'd do something like: Wide <- reshape(

Re: [R] Multidimensional scalling

2009-01-29 Thread stephen sefick
It depends on what you mean? If you would like a goodness of fit of your ordination to your distance matrix then this is doable and I would suggest that you look at the labdsv tutorial - http://ecology.msu.montana.edu/labdsv/R/ On Thu, Jan 29, 2009 at 6:50 PM, Tomek Wlodarski wrote: > Dear R d

Re: [R] How do I get my IT department to "bless" R?

2009-01-29 Thread Jim Porzak
Yes, Erik, there are all MS shops around! Ours happens to be one. However, I have absolutely no push back from IT on my use of R to do marketing analytics. The trick, Dan, is to deliver relevant and actionable results to the business. Your champions will stick up for you when, and if, you get any

Re: [R] Regression

2009-01-29 Thread stephen sefick
Provide reproducible code. And you will get an answer, probably. On Thu, Jan 29, 2009 at 8:29 PM, Edwin Wibisono wrote: > Hello, I have problem > I have a, and b in regression > then I can't plot x,y with (and) a, b lines > > can you help me ? > > thx > >[[alternative HTML version delete

[R] Matlab inv() and R solve() differences

2009-01-29 Thread Joseph P Gray
I submit the following matrix to both MATLAB and R x= 0.133 0.254 -0.214 0.116 0.254 0.623 -0.674 0.139 -0.214 -0.674 0.910 0.011 0.116 0.139 0.011 0.180 MATLAB's inv(x) provides the following 137.21 -50.68 -4.70 -46.42 -120.71 27.28 -8.94 62.19 -58.15 6.93 -7.89 36.94 8.35

[R] Regression

2009-01-29 Thread Edwin Wibisono
Hello, I have problem I have a, and b in regression then I can't plot x,y with (and) a, b lines can you help me ? thx [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Q about how to use Anova.mlm

2009-01-29 Thread pgseye
Hi, Am newish to stats and R, so I certainly appreciate any help. Basically I have 50 inidividuals whom I have 6 photos each of their optic nerve head. I want to check that the orientation of the nerve head is consistent, ie the 6 replicates show minimal or preferably no rotation differences. I'l

Re: [R] gls prediction using the correlation structure in nlme

2009-01-29 Thread Dr Carbon
On Wed, Jan 28, 2009 at 9:51 AM, Dr Carbon wrote: > How does one coerce predict.gls to incorporate the fitted correlation > structure from the gls object into predictions? In the example below > the AR(1) process with phi=0.545 is not used with predict.gls. Is > there another function that does th

Re: [R] How do I get my IT department to "bless" R?

2009-01-29 Thread Erik Iverson
This is a very broad question, and the answer is going to depend on your particular situation, which we are not privy to. I'll say two things. First, you should try to figure out why they would not want you to run R, so you can address those reasons specifically. Second, you might take a par

[R] Multidimensional scalling

2009-01-29 Thread Tomek Wlodarski
Dear R developers and users! I have calculated metric MDS by cmdscale from matrix of distances (dissimilarities). I would like to ask you how can I estimate how well this new mapping represents characteristic features of my data set? Thank you for any suggestions. Best, tomek ___

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread Gabor Grothendieck
On Thu, Jan 29, 2009 at 12:05 PM, Henrik Bengtsson wrote: > PS. About the double-letter index (e.g. ii vs. i); A few years ago > someone suggested me to use this, because it is much easier to search > for 'ii' in an editor compared with a single-letter 'i'. So true. I > made the move immediatel

Re: [R] bootstrapping in regression

2009-01-29 Thread Thomas Mang
Greg Snow wrote: What you are describing is actually a permutation test rather than a bootstrap (related concepts but with a subtle but important difference). The way to do a permutation test with multiple x's is to fit the reduced model (use all x's other than x1 if you want to test x1) on th

[R] How do I get my IT department to "bless" R?

2009-01-29 Thread Daniel Viar
I currently use R at work "under the radar", but there's a chance I could loose that access. I'd like to get our company to feel comfortable with open source and R in particular. Does anyone have any experience with their company's IT department and management that they would be willing to share?

Re: [R] convergence problem gamm / lme

2009-01-29 Thread geert aarts
problem solved: As is probably mostly the case, a convergence problem lies in the specification of the model or the data itself. Some information: I was trying to model the spatial distribution of fish of a particular age. The raw observations consisted of the number of individuals of a part

Re: [R] Optim error: initial value in 'vmmin' is not finite

2009-01-29 Thread Ben Bolker
SnowManPaddington gmail.com> writes: > > > Error in optim(method = "BFGS", c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, : > initial value in 'vmmin' is not finite > > I am running a logit model with latent class segments. I successfully got > estimates for 2 segments. However when I tried to incr

Re: [R] Power analysis for MANOVA?

2009-01-29 Thread Adam D. I. Kramer
Thanks for the response, Stephan. Really, I am trying to say, "My result is insignificant, my effect sizes are tiny, you may want to consider the possibility that there really are no meaningful differences." Computing post-hoc power makes a bit stronger of a claim in this setting. My real goal i

Re: [R] bootstrapping in regression

2009-01-29 Thread Greg Snow
What you are describing is actually a permutation test rather than a bootstrap (related concepts but with a subtle but important difference). The way to do a permutation test with multiple x's is to fit the reduced model (use all x's other than x1 if you want to test x1) on the original data and

Re: [R] Goodness of fit for gamma distributions

2009-01-29 Thread Albyn Jones
it is easy to make a qqplot for the gamma; suppose that the sample parameters are 1.101 and 2.49, the data in x: plot(qgamma(ppoints(x),1.101,2.49),sort(x)) see also lattice:qqmath albyn Quoting Dan31415 : Ah yes, that does produce a nice plot. Can i just ask what exactly it is sho

Re: [R] Importing data from clipboard on Mac OSX

2009-01-29 Thread Charles C. Berry
On Thu, 29 Jan 2009, Christian Anderson wrote: Hello R-Help, I noticed that there is a thread about importing data from the clipboard that is very poorly answered in the forum. One user suggests giving up, the other gives a solution that echoes the clipboard, but that's exactly the same as just

[R] Optim error: initial value in 'vmmin' is not finite

2009-01-29 Thread SnowManPaddington
Error in optim(method = "BFGS", c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, : initial value in 'vmmin' is not finite I am running a logit model with latent class segments. I successfully got estimates for 2 segments. However when I tried to increase the no. of segments, I got this error message at

Re: [R] bootstrapping in regression

2009-01-29 Thread Tom Backer Johnsen
Tom Backer Johnsen wrote: Thomas Mang wrote: Hi, Please apologize if my questions sounds somewhat 'stupid' to the trained and experienced statisticians of you. Also I am not sure if I used all terms correctly, if not then corrections are welcome. I have asked myself the following question r

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread SnowManPaddington
Hi ya, thanks a lot everyone!! I changed rr:ii-1 to rr:(ii-1) and the code works!!! I finally get some estimates from the optimization function (i am doing a logit model with 2 segments). Thanks thanks!!! I didn't realize rr:(ii-1) and rr:ii-1 would make such a big difference, especially because t

Re: [R] Graphic device & graphics primitives

2009-01-29 Thread Greg Snow
Sigbert, The plot2script function in the TeachingDemos package does essentially what Duncan talks about for you. Create your plot then run the function giving it a filename to save the info into (or run without arguments and then past into a script window or text editor (only tested on windows

Re: [R] tab characters

2009-01-29 Thread jim holtman
?cat > x <- '\t' > print(x) [1] "\t" > cat(x) > On Thu, Jan 29, 2009 at 3:13 PM, Nick Matzke wrote: > Hi all, > > Working at the R command line, how do I get strings to display e.g. tab or > newline characters as they should be displayed, rather than as e.g. \n or > \t? > > e.g.: >> x="

Re: [R] Welcome to the "R-help" mailing list

2009-01-29 Thread Johannes Huesing
Bishwa S Koirala [Thu, Jan 29, 2009 at 08:56:17PM CET]: >> You must know your password to change your options (including changing >> the password, itself) or to unsubscribe. It is: >> It is definitely a good idea to change your password ASAP. -- Johannes Hüsing There is somethin

[R] assign point values as labels in spplot

2009-01-29 Thread emorway
In the code to follow, I'm trying to label points with their corresponding values but have been unsuccessful after many variations to the following code. The code below creates the plot I want, I simply cannot get the black points ("+") to display the actual value. I'm guessing the problem is so

[R] RDieHarder with file input invocation

2009-01-29 Thread Tomek S.
I got hold of the 0.0.7 RDieHarder binary package as zip for Windows from CRANextras but don't know how to invoke its main method dieharder() in R with "file_input" as RNG argument - how to specify path to that input file ? the 0.0.7 as opposed to 0.1.0 doesn't seem to have ''inputfile'' arg fo

Re: [R] Stacking data

2009-01-29 Thread hadley wickham
            27193   30949    V2 > > But i have 41 columns (age column + 40 individuals) > I have the following script but an error is thrown up > can anyone help, where am i going wrong > > zz <- read.csv("Filename.csv",strip.white = TRUE) Try this: install.packages("r

[R] tab characters

2009-01-29 Thread Nick Matzke
Hi all, Working at the R command line, how do I get strings to display e.g. tab or newline characters as they should be displayed, rather than as e.g. \n or \t? e.g.: > x="\t" > x="\t" > x [1] "\t" > print(x) [1] "\t" -- Nicholas J. Matzke

Re: [R] Welcome to the "R-help" mailing list

2009-01-29 Thread Bishwa S Koirala
On Thu, 29 Jan 2009 20:55:19 +0100 r-help-requ...@r-project.org wrote: Welcome to the R-help@r-project.org mailing list! To post to this list, send your email to: r-help@r-project.org General information about the mailing list is at: https://stat.ethz.ch/mailman/listinfo/r-help If you eve

[R] Importing data from clipboard on Mac OSX

2009-01-29 Thread Christian Anderson
Hello R-Help, I noticed that there is a thread about importing data from the clipboard that is very poorly answered in the forum. One user suggests giving up, the other gives a solution that echoes the clipboard, but that's exactly the same as just ctrl-p. As I am asked this question at least once

Re: [R] ggplot2 - how to change location / position of wind rose axis labels?

2009-01-29 Thread Darren Norris
Many thanks Hadley, The solution you gave works well moving label locations within the area of the axis but does not quite provide what I was after. Apologies for the confusion - I've only been using ggplot2 it for a day. Compare the bar chart (from previous example): awind<-doh + geom_bar

Re: [R] Adding vertical line to histogram and qplot "stacked" plot

2009-01-29 Thread Duncan Murdoch
On 1/29/2009 2:03 PM, Jason Rupert wrote: R-users it appears I am leaning on your knowledge once again. Is there any way to add a vertical line to a histogram and qplot "stacked" plot? Here is my current attempt: "qplot" approach attempt: qplot(Run, data = data_dataframe, breaks = breaks, f

Re: [R] Unicode

2009-01-29 Thread Duncan Murdoch
On 1/29/2009 1:57 PM, R Heberto Ghezzo, Dr wrote: Hello, can somebody explain me why the following program does not work? It works for me, though \u2A8A and \u222C display as boxes, something like the way they do in Firefox if I look at this page: http://unicode.coeurlumiere.com/?n=8192 (bu

Re: [R] Adding vertical line to histogram and qplot "stacked" plot

2009-01-29 Thread Felipe Carrillo
Jason: Check Hadley's page, there's a few examples there. Good luck http://had.co.nz/ggplot2/geom_vline.html Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA --- On Thu, 1/29/09, Jason Rupert wrote: > From: Jason Ru

[R] Adding vertical line to histogram and qplot "stacked" plot

2009-01-29 Thread Jason Rupert
R-users it appears I am leaning on your knowledge once again.  Is there any way to add a vertical line to a histogram and qplot "stacked" plot?  Here is my current attempt:   "qplot" approach attempt: qplot(Run, data = data_dataframe, breaks = breaks, fill = Temperature, main = short_title) + sc

[R] Unicode

2009-01-29 Thread R Heberto Ghezzo, Dr
Hello, can somebody explain me why the following program does not work? Which pages of Unicode are implemented? the u22xx and 2Axx are math symbols and extensions Thanks Heberto Ghezzo Montreal plot(1) text(1.0,1.2,"a \u2A8A b \u222C c \u5222 d", cex=2) text(1.0,1.1," \u222C ", cex=2) text(1.0

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread davidr
The lines below made me understand clearly. Maybe they are already in some documentation, but if not, it might help others to avoid my misunderstanding. Thanks to all for the clarifications. -- David -Original Message- From: Duncan Murdoch [mailto:murd...@stats.uwo.ca] Sent: Thursday, Ja

[R] Problem VGAM and Predict

2009-01-29 Thread EinSchwede
Hello, since I installed the package VGAM I have problems useing the predict for othere methods. for example I have a model from glm and polr the command predict(model) I get the error: unable to find an inherited method for function "predict", for signature "polr". Has perhaps anybody a solut

[R] Accessing R as a Web service (UNCLASSIFIED)

2009-01-29 Thread Neiderer, Andrew (Civ, ARL/CISD)
Classification: UNCLASSIFIED Caveats: NONE I am new to R and somewhat to Web server programming. I am a Java programmer, however, and have done quite a bit with X3D (extensible 3D) computer graphics. A statistician that I work with is doing multidimensional scaling (MDS) and provides some x,

Re: [R] Stacking data

2009-01-29 Thread David Winsemius
What does str(zz) tell you? -- David Winsemius On Jan 29, 2009, at 10:31 AM, Amit Patel wrote: Hi I have data in the format below AgeV1 V2 V3 V4 23646 45190 50333 55166 56271 26174 35535 38227 37911 41184 27723 25691 25712 26

[R] Help

2009-01-29 Thread Alexandra Almeida
Hi everybody! I´m with a problem that probably is easy for you but I really don´t know how to solve. On the following script: for(j in 1:length(limiares)) { excessos<-limiares[j]-estacao[estacao__ R-help@r-project.org mailing list https://stat.ethz.c

Re: [R] Taking the min of each row in a matrix

2009-01-29 Thread jim holtman
?apply e.g. for the 'min' apply(yourMatrix, 1, min, na.rm=TRUE) On Thu, Jan 29, 2009 at 11:58 AM, Will Glass-Husain wrote: > Hi, > > I'm a new user. I've been reading through the manual and looking at > various examples but am still trying to make sense of the most > efficient ways to handle ma

Re: [R] bootstrapping in regression

2009-01-29 Thread Chuck Cleland
On 1/29/2009 11:43 AM, Thomas Mang wrote: > Hi, > > Please apologize if my questions sounds somewhat 'stupid' to the trained > and experienced statisticians of you. Also I am not sure if I used all > terms correctly, if not then corrections are welcome. > > I have asked myself the following quest

Re: [R] Taking the min of each row in a matrix

2009-01-29 Thread anna freni sterrantino
Hello Will, z=matrix(c(1,2,3,4,5,6), 2,byrow=T) min=apply(z,1,min) what do you need isapply ?apply and its family are very helpful. apply( matrix, columns, function) apply( matrix, rows, function) Check ?apply Cheers Anna Anna Freni Sterrantino Ph.D Student Department of Statistics U

Re: [R] Dynamic random effects model

2009-01-29 Thread Ben Bolker
I am forwarding back to the R list in hopes of getting you a better answer from someone else. This sounds like it might be logistic autoregression? RSiteSearch("logistic autoregress*") came up with http://finzi.psych.upenn.edu/R/library/repeated/html/gar.html Ben Bolker Joseph Magagnoli w

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread Wacek Kusnierczyk
Duncan Murdoch wrote: > On 1/29/2009 10:39 AM, dav...@rhotrading.com wrote: > >> And if the loop variable does not exist before the 'for', why is it >> created in the parent(?) environment at all? > > It's created in the current evaluation frame, because that's where > everything gets created un

Re: [R] Dynamic random effects model

2009-01-29 Thread Joseph Magagnoli
I am trying to estimate a model of third party intervention into civil war. The data is panel data structure. The unit of analysis is civil war year. For each civil war year my dependent variable is coded 0=no intervention and 1=intervention. I want to use a lagged dependent variable as an indep

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread Henrik Bengtsson
Thanks for straighten this out. Sorry for my misleading suggestion that "special" scoping rules comes into play; it is just about reassignments at the beginning of each loop. Here is a slightly better illustration: ii <- "start"; cat("ii:",ii,"\n"); for (ii in 1:2) { cat("Outer ii:",ii,"\n");

Re: [R] Help

2009-01-29 Thread rkevinburton
You could look into ' try' and set it up to catch errors and do the appropriate thing in your error handler. I don't have the exact syntax at hand right now but looking at ?try or ?tryCatch I think will do what you want. Kevin Alexandra Almeida wrote: > Hi everybody! > > I´m with a prob

[R] Taking the min of each row in a matrix

2009-01-29 Thread Will Glass-Husain
Hi, I'm a new user. I've been reading through the manual and looking at various examples but am still trying to make sense of the most efficient ways to handle matrices of data. If I have a 2D matrix of data, how do I get the mean, min, max value of each row? I see the "mean" function on a matri

[R] bootstrapping in regression

2009-01-29 Thread Thomas Mang
Hi, Please apologize if my questions sounds somewhat 'stupid' to the trained and experienced statisticians of you. Also I am not sure if I used all terms correctly, if not then corrections are welcome. I have asked myself the following question regarding bootstrapping in regression: Say for

[R] Standard Errors of Coefficients

2009-01-29 Thread Andreas Klein
Hello. I estimated a VAR(1) TSmodel (var_1) with estMaxLik from the dse1 package given a TSdata object (mydata). est.model <- estMaxLik(var_1,mydata) How can I obtain the standard errors for the four coefficients of the estimated model to check for significance? - Is it yet calculated and I ca

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread Duncan Murdoch
On 1/29/2009 11:06 AM, Patrick Burns wrote: Certainly not a complete description, but 'The R Inferno' talks about this on page 62. Hmmm, I don't think I agree with that description. There's only one object named i in that example (which was for (i in 1:6) { cat('\n i is', i, '\n\n') for

Re: [R] Help

2009-01-29 Thread patricia garcía gonzález
Hi, Test <- function(x, y, paired) { res <- try( ks.test(...), TRUE ) if( inherits( res, "try-error" ) ) res <- NA res } Regards patricia > Date: Thu, 29 Jan 2009 13:22:11 -0300 > From: alexandra...@gmail.com > To: r-help@r-project.org > Subject: [R] He

[R] Help

2009-01-29 Thread Alexandra Almeida
Hi everybody! I´m with a problem that probably is easy for you, but I really don´t know how to solve. On the following script: for(j in 1:length(limiares)) { excessos<-limiares[j]-estacao[estacao__ R-help@r-project.org mailing list https://stat.ethz.

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread Patrick Burns
Certainly not a complete description, but 'The R Inferno' talks about this on page 62. Patrick Burns patr...@burns-stat.com +44 (0)20 8525 0696 http://www.burns-stat.com (home of "The R Inferno" and "A Guide for the Unwilling S User") dav...@rhotrading.com wrote: I apologize for posting a wron

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread Duncan Murdoch
On 1/29/2009 10:39 AM, dav...@rhotrading.com wrote: I apologize for posting a wrong opinion; I should of course have checked before posting. Henrik's examples illustrate something I had never realized before, and it really surprised me! Where can I read the technical details of this scoping aspe

Re: [R] Inconsistency in F values from dropterm and anova

2009-01-29 Thread Prof Brian Ripley
Perhaps if you followed the posting guide and did not send HTML mail your tables would be readable. But your anova call seems wrong, as you have the models in decreasing not increasing order. The correct result is (computing F test by hand) 23.1814/93 [1] 0.2492624 1.5475/3 [1] 0.5158333

Re: [R] scoping rules for 'for' [was - Re: for/if loop ]

2009-01-29 Thread davidr
I apologize for posting a wrong opinion; I should of course have checked before posting. Henrik's examples illustrate something I had never realized before, and it really surprised me! Where can I read the technical details of this scoping aspect of 'for' as it still presents some subtle puzzles f

Re: [R] convergence problem gamm / lme

2009-01-29 Thread Simon Wood
Geert, Sorry for slow reply... I don't see any obvious problems with what you've done, so I guess it's the usual problem that PQL just doesn't *have* to converge, and the bit of extra flexibility of using a smooth is too much for it in this case. If you send me the data offline I can dig a litt

Re: [R] ggplot2 - how to change location / position of wind rose axis labels?

2009-01-29 Thread hadley wickham
Hi Darren, Sorry for taking so long to get back to you - my hard drive died and it's taking me a while to get up and running again of backups etc. One solution to your problem is to draw the axis labels yourself: labels <- data.frame(rrating = 1:10, count = 19000) awind + opts(axis.text.x =

[R] Stacking data

2009-01-29 Thread Amit Patel
Hi I have data in the format below      Age            V1       V2      V3       V4    23646         45190 50333 55166 56271    26174         35535 38227 37911 41184    27723          25691 25712 26144 26398 and would like to sort it as fo

Re: [R] Repeated measures design for GAM? - corrected question...

2009-01-29 Thread Simon Wood
Diederik, The `gamm' call looks fine provided `park' is nested in `count'. You are basically adding a random effect for each survey, and a random effect for each `park' within survey to the the smooth effects. best, Simon ps. I assume there are more than 18 observations in the full dataset --

[R] lmer for a binary dependent variable

2009-01-29 Thread Liliana Martinez
Hi,   I am trying to use the lmer function from the lme4 package in R 2.8.0. to fit a generalized mixed-effects model for a dependent variable with a binomial distribution (for more info on my experiment, look below). However, I encounter a major problem: How is it possible to find the general

Re: [R] Question On CrossTable function in gmodels package

2009-01-29 Thread Marc Schwartz
on 01/29/2009 08:01 AM eugen pircalabelu wrote: > Hi R-users, > I have the following problem with CrossTable function within “gmodels” > package: the output of the function (format “spss” and asresid=T) can not be > stored within another object. > For example: > >> library(gmodels) >> data(infer

Re: [R] t-test

2009-01-29 Thread patricia garcía gonzález
Hi, If you mean if the t.test is done as if the samples where paired, the answer is yes if you write argument paired = TRUE; and the pairs are done in order, that is 1º with 1º, 2º with 2º, etc. As you wrote, (paired = FALSE) the t.test is unpaired, and the order of elements in the vectors a

[R] Question On CrossTable function in gmodels package

2009-01-29 Thread eugen pircalabelu
Hi R-users, I have the following problem with CrossTable function within “gmodels” package: the output of the function (format “spss” and asresid=T) can not be stored within another object. For example: >library(gmodels) >data(infert, package = "datasets") > CrossTable(infert$education, infert$i

Re: [R] Plot dagger symbol in R

2009-01-29 Thread Duncan Murdoch
On 1/29/2009 8:56 AM, Rau, Roland wrote: Dear all, thank you very much. Yes, indeed, I should have included the sessionInfo(). Both examples (unicode and standard encoding) work fine; however, I prefer the suggestion by Prof. Ripley since it allows me to do things like: plot(1, main=expression(

Re: [R] Dynamic random effects model

2009-01-29 Thread Ben Bolker
Joseph Magagnoli gmail.com> writes: > > All R experts, > How do I fit a dynamic Random effects model with a binary dependent variable > in R > Thanks > JCM > You haven't given us nearly enough information to go on. If you're talking about something like a state-space model with a binary resp

[R] Inconsistency in F values from dropterm and anova

2009-01-29 Thread IMS
Hi, I'm working on fitting a glm model to my data using Gamma error structure and reciprocal link. I've been using dropterm (MASS) in the model simplification process, but the F values from analysis of deviance tables reported by dropterm and anova functions are different - sometimes significant

Re: [R] Goodness of fit for gamma distributions

2009-01-29 Thread Dan31415
Ah yes, that does produce a nice plot. Can i just ask what exactly it is showing. It seems to me to be a sort of Q-Q plot but with a different set of axes. Is this correct, if so do the same interpretation rules apply for this plot, i.e. departures from either end of the curve show poor fitting of

[R] t-test

2009-01-29 Thread Amit Patel
When doing the t-test in the below manner will r compare each element of the array with the relevant one. I.e. if i was comparing x and y would (1 and 0) and (1 and 9) be treated as separate variables. Or does it just assume one variable. # test data x <- c(1,1.1,1.15,1.2,1.21,1.23) y <- c(

Re: [R] Plot dagger symbol in R

2009-01-29 Thread Rau, Roland
Dear all, thank you very much. Yes, indeed, I should have included the sessionInfo(). Both examples (unicode and standard encoding) work fine; however, I prefer the suggestion by Prof. Ripley since it allows me to do things like: plot(1, main=expression(e^"\206")) This fails, however, with: plo

Re: [R] Multiple tables

2009-01-29 Thread jim holtman
you can use combn to create the combinations and the following will create a list of all the results: x1 <- x2 <- x3 <- x4 <- 1:10 comb <- combn(c('x1','x2', 'x3', 'x4'), 2) myTab <- lapply(seq(ncol(comb)), function(x){ table(get(comb[1, x]), get(comb[2, x])) }) # put names of the combinations

[R] ANOVA in R

2009-01-29 Thread Amit Patel
Hi I Have a very large dataset that I would like to conduct ANOVA tests on. Im not a very strong programmer so any help would be appreciated. the format is Identifier             A1       A2        B1      B2       C1   C2      Norm1         N

Re: [R] Plot dagger symbol in R

2009-01-29 Thread Prof Brian Ripley
On Thu, 29 Jan 2009, Mark Difford wrote: Hi Roland, But this is obviously not a dagger and it seems the Adobe Symbol font does not have a dagger. True, but ... Yoda was here. plot(0:1,0:1, type="n") text(x=0.5, y=0.5, labels=expression("\u2020")) text(x=0.4, y=0.6, labels=expression("\u202

Re: [R] Plot dagger symbol in R

2009-01-29 Thread baptiste auguie
Hi, If all else fails, you could consider using LaTeX itself with psfrag, or perhaps a similar idea involving eps2pgf. http://biostat.mc.vanderbilt.edu/twiki/bin/view/Main/PsFrag Hope this helps, baptiste On 29 Jan 2009, at 11:24, Rau, Roland wrote: Dear all, I would like to plot the da

Re: [R] Plot dagger symbol in R

2009-01-29 Thread Mark Difford
Hi Roland, >> But this is obviously not a dagger and it seems the Adobe Symbol font >> does not have a dagger. True, but ... Yoda was here. plot(0:1,0:1, type="n") text(x=0.5, y=0.5, labels=expression("\u2020")) text(x=0.4, y=0.6, labels=expression("\u2021")) ?plotmath, sub "Other symbols" ...

Re: [R] Plot dagger symbol in R

2009-01-29 Thread Duncan Murdoch
Rau, Roland wrote: Dear all, I would like to plot the dagger symbol in R (like LaTeX's \dagger). However, I was unable to do so. First, I thought maybe dagger actually exists just like the degree symbol: plot(0:1,0:1, type="n") text(x=0.5, y=0.5, labels=expression(degree)) plot(0:1,0:1, type=

[R] Use SOCKS proxy

2009-01-29 Thread Daniel Brewer
Hi, Is there anyway to set up R so that it uses a SOCKS proxy in Linux? I am getting some strange issues with the Institute's new web filtering system and want to be able to test whether a problem I have with the biomaRt package is caused by it. Many thanks Dan -- *

Re: [R] Graphic device & graphics primitives

2009-01-29 Thread Duncan Murdoch
Sigbert Klinke wrote: Hi, I know that some graphics devices in R store graphics primitives such that a redraw is possible (e.g. when resizing the window). Is it possible to get the current number of stored graphic primitives? Thanks in advance Sigbert Klinke ___

[R] Plot dagger symbol in R

2009-01-29 Thread Rau, Roland
Dear all, I would like to plot the dagger symbol in R (like LaTeX's \dagger). However, I was unable to do so. First, I thought maybe dagger actually exists just like the degree symbol: plot(0:1,0:1, type="n") text(x=0.5, y=0.5, labels=expression(degree)) plot(0:1,0:1, type="n") text(x=0.5, y=0.

Re: [R] Graphic device & graphics primitives

2009-01-29 Thread Prof Brian Ripley
On Thu, 29 Jan 2009, Sigbert Klinke wrote: Hi, I know that some graphics devices in R store graphics primitives such that a redraw is possible (e.g. when resizing the window). Is it possible to get the current number of stored graphic primitives? That's not what happens: the graphics engine

[R] Graphic device & graphics primitives

2009-01-29 Thread Sigbert Klinke
Hi, I know that some graphics devices in R store graphics primitives such that a redraw is possible (e.g. when resizing the window). Is it possible to get the current number of stored graphic primitives? Thanks in advance Sigbert Klinke __ R-help@

[R] Odp: Text in a character vector to indicate "ifelse" argument

2009-01-29 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 29.01.2009 10:25:13: > > Hello > > I have a data set that looks like this; > > > b2 > dato chr status PRRSvac > PRRSsanVac PRRSsanDk PRRSdk > 33 2007-12-03 090432Rød SPF > 34 200

Re: [R] convergence problem gamm / lme

2009-01-29 Thread geert aarts
Simon, thanks for your reply and your suggestions. I fitted the following glmm's gamm3<-try(glmmPQL(count~offset(offsetter)+poly(lon,3)*poly(lat,3),random=list(code_tripnr=~1),family="poisson")) Which worked OK (see summary below) I also fitted a model using quasipoisson, but that didn't

  1   2   >