Re: [R] matching observations and ranking

2013-04-24 Thread arun
Hi, It is not that clear. If VAR1 is a match between columns AB001A, AB0002A, VAR2  between AB001A, AB362 and VAR3 between AB0002A and AB362: Also, I assume row8 match would be taken as 1. dat1- read.table(text=   S.No AB001A AB0002A AB362    1   -/-    C/C   A/A       

Re: [R] Regression on stratified count data

2013-04-24 Thread peter dalgaard
On Apr 24, 2013, at 06:15 , meng wrote: Hi all: For stratified count data,how to perform regression analysis? My data: age case oc count 1 1 121 1 1 226 1 2 117 1 2 259 2 1 118 2 1 288 2 2

[R] Simulate user input in R

2013-04-24 Thread Vahe nr
Dear all, Can we simulate user input in R ? for example if we have a function which needs an input from the user to continue its work, can we automate this step (simulate the input...) Here is the sample: choose between one of the grouping factor available : c(Village, Country) we need to enter

[R] Understanding how nls() works

2013-04-24 Thread Patrick Connolly
I'm trying to understand what goes on in the process that nls() uses. This converges without much drama: rate.nls - nls(Dev ~ exp(rho * (T)) - exp(rho*Tmax - (Tmax - T)/del)+ + lam , data = xx, trace = TRUE, + start = c(rho = 0.15, del = 7, lam = .02, Tmax = 30))

Re: [R] Simulate user input in R

2013-04-24 Thread Jeff Newmiller
Your example is not reproducible, so it is open to misinterpretation. Normally the approach taken in R is not to wait for user input, but have the user call the function with the desired values. One interpretation of your question is that the function you are using is broken by design, and you

Re: [R] Simulate user input in R

2013-04-24 Thread Vahe nr
Thanks Jeff for your reply... The code is not written by me it is the NDVITS package which is used in R. TimeSeriesAnalysis {ndvits} And at that point it requires my input. Regards, Vahe On Wed, Apr 24, 2013 at 11:04 AM, Jeff Newmiller jdnew...@dcn.davis.ca.uswrote: Your example is not

[R] News package

2013-04-24 Thread Aswathy Nair
Hi, Is there any package available in R to download news content? Thanks in advance. Ashy [[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

[R] pglm package: fitted values and residuals

2013-04-24 Thread alfonso . carfora
I'm using the package pglm and I'have estimated a random probit model. I need to save in a vector the fitted values and the residuals of the model but I can not do it. I tried with the command fitted.values using the following procedure without results: library(pglm) m1_S-pglm(Feed ~

[R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Liviu Andronic
Dear all, I've bumped into the: Error in loadNamespace(name) : there is no package called ‘R.utils’ error. I've already read a bit on this ( http://www.cybaea.net/Blogs/Data/A-warning-on-the-R-save-format.html ) but I have a follow-up question. Given a workspace that automatically loads a package

[R] Looking for a better code for my problem.

2013-04-24 Thread Christofer Bogaso
Hello again, Let say I have following data: Dat - structure(list(AA = structure(c(3L, 1L, 2L, 1L, 2L, 3L, 3L, 2L, 3L, 1L, 1L, 3L, 3L, 2L, 2L, 3L, 2L, 1L, 1L, 1L), .Label = c(A, B, C), class = factor), BB = structure(c(2L, 3L, 2L, 2L, 2L, 3L, 2L, 2L, 2L, 1L, 1L, 2L, 3L, 1L, 3L, 2L, 1L, 2L, 2L, 3L

Re: [R] Looking for a better code for my problem.

2013-04-24 Thread Rui Barradas
Hello, You can easily make of the following a one-liner. Note that the order of rows is not the same as in your code, so identical() will return FALSE. idx - Dat[, 'AA'] == A | (Dat[, 'AA'] == B Dat[, 'BB'] == b) res2 - Dat[idx, ] Hope this helps, Rui Barradas Em 24-04-2013 11:21,

[R] extract function extracting only NA values

2013-04-24 Thread Fabio Berzaghi
Hello, I have five raster files in ASCII format. With four of them I have no problem extracting values based on a set of X and Y coordinates. Unfortunately with one of the files all I managed to extract is NA values. To verify the problem I have opened the raster with ArcMap and there are no

Re: [R] Regression on stratified count data

2013-04-24 Thread Achim Zeileis
On Wed, 24 Apr 2013, meng wrote: Hi all: For stratified count data,how to perform regression analysis? My data: age case oc count 1 1 121 1 1 226 1 2 117 1 2 259 2 1 118 2 1 288 2 2 1 7 2 2

Re: [R] Looking for a better code for my problem.

2013-04-24 Thread Jorge I Velez
Try subset(Dat, AA == A | (AA == B BB == b)) HTH, Jorge.- On Wed, Apr 24, 2013 at 8:21 PM, Christofer Bogaso bogaso.christo...@gmail.com wrote: Hello again, Let say I have following data: Dat - structure(list(AA = structure(c(3L, 1L, 2L, 1L, 2L, 3L, 3L, 2L, 3L, 1L, 1L, 3L, 3L, 2L, 2L,

Re: [R] extract function extracting only NA values

2013-04-24 Thread Fabio Berzaghi
So I think I might have found what is causing this problem. The values I have in this raster summary(ps0011yme) layer Min.-9458.911 1st Qu. 1955256.000 Median 10618870.000 3rd Qu. 79577730.000 Max.167780500.000 NA's0.000 From ArcMap though I get

Re: [R] Bootstrapping a 11x10 matrix

2013-04-24 Thread Adams, Jean
Herbejie Rose, You could use the boot() function in the R package boot. For example: # example data matrix m - matrix(sample(11*10), ncol=10) # function to calculate column means for indexed rows of matrix myfun - function(data, i) { apply(data[i, ], 2, mean) } # 1000 bootstrap samples

Re: [R] News package

2013-04-24 Thread R. Michael Weylandt
On Wed, Apr 24, 2013 at 9:02 AM, Aswathy Nair ashy4...@gmail.com wrote: Hi, Is there any package available in R to download news content? What news source are you looking for? You could, e.g., use the twitteR package, but to my knowledge for things like Google News or the BBC you'll need to

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Duncan Murdoch
On 13-04-24 5:46 AM, Liviu Andronic wrote: Dear all, I've bumped into the: Error in loadNamespace(name) : there is no package called ‘R.utils’ error. I've already read a bit on this ( http://www.cybaea.net/Blogs/Data/A-warning-on-the-R-save-format.html ) but I have a follow-up question. Given a

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Martin Morgan
On 04/24/2013 06:03 AM, Duncan Murdoch wrote: On 13-04-24 5:46 AM, Liviu Andronic wrote: Dear all, I've bumped into the: Error in loadNamespace(name) : there is no package called ‘R.utils’ error. I've already read a bit on this (

Re: [R] about the loglm and glm---Re:Re: Regression on stratified count data

2013-04-24 Thread Achim Zeileis
On Wed, 24 Apr 2013, meng wrote: Hi,Achim: Can all the analysis you mentioned via loglm be performed via glm(...,family=poisson)? Yes. ## transform table back to data.frame df - as.data.frame(tab) ## fit models: conditional independence, no-three way interaction, ## and saturated g1 -

Re: [R] extract function extracting only NA values

2013-04-24 Thread Fabio Berzaghi
ok so i imported the raster as a dataframe and everything is fine Min. :-69826220 Max. :167780500 so it's something with values lower than - that are interpreted as NA by the raster function On 4/24/2013 13:52, Fabio Berzaghi wrote: So I think I might have found what is causing

[R] How to make a raster image in R from my own data set

2013-04-24 Thread Kristi Glover
Hi R-user, I was trying to make a raster map with WGS84 projection in R, but I could not make it. I found one data set in Google that data is almost the same format as of mine. I wanted to make a raster map of temperature with 1 degree spatial resolution for the global scale. I could make it

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Duncan Murdoch
On 13-04-24 10:12 AM, Martin Morgan wrote: On 04/24/2013 06:03 AM, Duncan Murdoch wrote: On 13-04-24 5:46 AM, Liviu Andronic wrote: Dear all, I've bumped into the: Error in loadNamespace(name) : there is no package called ‘R.utils’ error. I've already read a bit on this (

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Liviu Andronic
Dear Duncan, On Wed, Apr 24, 2013 at 4:57 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote: What I've done sometimes in debugging is to change that error to a warning in the getNamespace() function, and add some tracing code to the serialization code to print the names of objects as they are

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Liviu Andronic
Dear Duncan, On Wed, Apr 24, 2013 at 3:03 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote: A better approach is to *never* save and load workspaces unless you know exactly what is in them. Always reply no to the question about saving your workspace (or set that as the default). If you

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread R. Michael Weylandt
On Wed, Apr 24, 2013 at 4:08 PM, Liviu Andronic landronim...@gmail.com wrote: Dear Duncan, On Wed, Apr 24, 2013 at 3:03 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote: A better approach is to *never* save and load workspaces unless you know exactly what is in them. Always reply no to

[R] Problem with raster function with values lower than -9999

2013-04-24 Thread Fabio Berzaghi
As I have reported in a previous thread, I have been having problems importing a ASCII raster that has values that go from Min. :-69826220 to Max. :167780500. The problem I am encountering is that when I use the raster function to import the ASCII file then every value smaller than - is

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Hadley Wickham
I must admit that I'm a bit surprised by this. I was always under the impression that saving/restoring workspaces was the proper workflow in R. If you use R interactively (e.g., not by running scripts), how else would you store your data, intermediary results, etc., while working on a

Re: [R] the joy of spreadsheets (off-topic)

2013-04-24 Thread peter dalgaard
In case you haven't noticed, this is making the rounds in the media, including a handful of references to R. See e.g. http://news.slashdot.org/story/13/04/17/0215211/excel-error-contributes-to-problems-with-austerity-study I suppose we can't fortune()'ify anonymous quotes, but I kind of like

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread John Kane
-Original Message- From: h.wick...@gmail.com Sent: Wed, 24 Apr 2013 10:48:01 -0500 To: landronim...@gmail.com Subject: Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’ I must admit that I'm a bit surprised by this. I was

Re: [R] Need to replicate Boltzman Signmodial Curve fit from Graph Pad

2013-04-24 Thread David Carlson
Sorry. I left out a line: pHvals - seq(min(dta$pH), max(dta$pH), length.out=100) Should be inserted just before the lines() function. - David L Carlson Associate Professor of Anthropology Texas AM University College Station, TX 77840-4352 From: J J

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Duncan Murdoch
On 13-04-24 11:08 AM, Liviu Andronic wrote: Dear Duncan, On Wed, Apr 24, 2013 at 3:03 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote: A better approach is to *never* save and load workspaces unless you know exactly what is in them. Always reply no to the question about saving your

[R] Trouble Computing Type III SS in a Cox Regression using drop1 and Anova

2013-04-24 Thread Paul Miller
Hello All, Am having some trouble computing Type III SS in a Cox Regression using either drop1 or Anova from the car package. Am hoping that people will take a look to see if they can tell what's going on. Here is my R code: cox3grp - subset(survData, Treatment %in% c(DC, DA, DO), c(PTNO,

[R] [r] Problems with Fortran calls when loaded a dll compiled with gfortran-4 Cygwin 4.5.3

2013-04-24 Thread Jens Olofsson
Dear users of R I have a subroutine in Fortran95, compiled to a DLL with gfortran in Cygwin 4.5.3. The subroutine is: subroutine MyPBP( S, p, N ) ! Expose subroutine rtest to users of this DLL !DEC$ ATTRIBUTES DLLEXPORT, C, REFERENCE, ALIAS: mypbp_ ::mypbp ! This function computes the

[R] RDA permutest envfit

2013-04-24 Thread Rémi Lesmerises
Dear all, I did a RDA and when I looked to the signification of the test with permutest, the output was non-significant. But when I used the envfit function, some of the vectors are significant. All the test's conditions are respected. What it means? Is it an error in the script? Commands

[R] Sum up column values according to row id

2013-04-24 Thread Matteo Mura
Dear All, here a problem I think many of you can solve in few minutes. I have a dataframe which contains values of plot id, diameters, heigths and basal area of trees, thus columns names are: id | dbh | h | g head(ipso, n=10)id dbh h g 1 FPE0164 36 13.62 0.10178760 2

[R] R Interactive Mode

2013-04-24 Thread Hrachya Astsatryan
Dear all, We are doing some research about the time series analysis of NDVI, and we found the NDVITS package which is a very great tool. Unfortunately when we run it, after TimeSeriesAnalysis it asks to enter Village or Country. library(ndvits,

[R] about the loglm and glm---Re:Re: Regression on stratified count data

2013-04-24 Thread meng
Hi,Achim: Can all the analysis you mentioned via loglm be performed via glm(...,family=poisson)? Many thanks. At 2013-04-24 19:37:10,Achim Zeileis achim.zeil...@uibk.ac.at wrote: On Wed, 24 Apr 2013, meng wrote: Hi all: For stratified count data,how to perform regression

Re: [R] matching observations and ranking

2013-04-24 Thread arun
Hi, May be this helps: As you wanted to match only from row3 onwards to row2, the corresponding values on row1 and row2 were set to NA. dat1- read.table(text=   S.No AB001A AB0002A AB362    P1   -/-    C/C   A/A       P2   C/C    C/C   A/A       3 

[R] R cannot find the path on my mac

2013-04-24 Thread Gitte Brinch Andersen
Hi I am really sorry for this probably quite simple question. I am new to R, and I am running a pipeline that has already been made. All I have to do is give the paths for different folders, where the pipeline can find the files with my data. But every time I try to run the pipeline it returns

[R] Regression and FMMs with flexmix

2013-04-24 Thread Robin Tviet
I am repeating this because it seems that some people think it is important to reveal your identity I don;t understand why this is so important. Hopefuly now this list will be helpful. Could someone please assist with this I am trying to understand how to use the flexmix package, I have

Re: [R] matching observations and ranking

2013-04-24 Thread arun
Just to add: If your original dataset have only few columns, then you can try this too:

[R] puzzles of Finance data programming with R

2013-04-24 Thread chenzi14
Hi This is Candice from Newcastle University. I am now coming across some problems with the programming of R. It is like this. I am asked to run three models based on a sample of CPI: Random walk Recursive AR(4) Rolling AR(4) 1.Based on the papers I find, I think the random walk may be similar to

[R] fonts not rendering during plot

2013-04-24 Thread George Dietz
Hi, I am having a problem where sometimes fonts for text in plots don't get rendered-the text just shows up as boxes. If I call R from a certain directory the problem goes away, otherwise the fonts don't render (but plots get made). In both cases, my PANGO_RC and FONTCONFIG_FILE do not change…

[R] Residuals for fracdiff

2013-04-24 Thread fmesserGmail
Hi, I am using the fracdiff package to estimate the parameters of an ARFIMA(1,d,1) model. I would also like to get the residuals of the series. I have seen another post about this (below). However, being still quite at the beginner level in terms of R, I did not quite understand how this

[R] string size limits in RCurl

2013-04-24 Thread Elmore, Ryan
Hi All, I am running into what appears to be character size limit in a JSON string when trying retrieve data from either `curlPerform()` or `getURL()`. Here is non-reproducible code [1], but it should shed some light on the problem. # Note that .base.url is the basic url for the API, q is

Re: [R] [r] Problems with Fortran calls when loaded a dll compiled with gfortran-4 Cygwin 4.5.3

2013-04-24 Thread Duncan Murdoch
On 13-04-24 1:36 PM, Jens Olofsson wrote: Dear users of R I have a subroutine in Fortran95, compiled to a DLL with gfortran in Cygwin 4.5.3. We don't support Cygwin. You should use the gfortran in Rtools, and get R to set the command line options for you, either by putting the code in a

[R] substitute rename.vars {gdata}

2013-04-24 Thread Knut Krueger
is there a equivalent to rename.vars {gdata}? As I do not use read.xls i am not using library gdata for a package except this function Knut __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] [r] Problems with Fortran calls when loaded a dll compiled with gfortran-4 Cygwin 4.5.3

2013-04-24 Thread Jens Olofsson
Dear Duncan, I know this isn't a forum for Cygwin, but for R. Pls treat me as a noob and also remember I am on Windows. How should I use R CMD SHLIB as if I write that at the prompt I get the error: unexpected symbol in R CMD. I have mango.f95 in the working directory. //Jens On 24 April 2013

Re: [R] the joy of spreadsheets (off-topic)

2013-04-24 Thread Thomas Adams
One might wonder if the Excel error was indeed THAT or perhaps a way to get the desired results, give the other issues in their analysis? On Wed, Apr 24, 2013 at 11:58 AM, peter dalgaard pda...@gmail.com wrote: In case you haven't noticed, this is making the rounds in the media, including a

[R] Distance matrices Combinations

2013-04-24 Thread eliza botto
Dear UseRs, MY PROBLEM IS A SMALL PIECE OF A REAL BIG AND A COMPLICATED PROBLEM. IF I DELIBERATE IN A VERY SIMPLE WAY THEN ALL I WANT IS TO PUT ALL THE POSSIBLE COMBINATIONS OF 75 DISTANCE MATRICES (BY TAKING 4 MATRICES, MORE COMMONLY 75C4), in the following equation.

Re: [R] Sum up column values according to row id

2013-04-24 Thread David Carlson
Something like this? mean6 - function(x) { if (length(x) 6) { mn - mean(x) } else { mn - mean(x[1:6]) } return(mn) } aggregate(g~id, ipso, mean6) - David L Carlson Associate Professor of Anthropology Texas AM

Re: [R] [r] Problems with Fortran calls when loaded a dll compiled with gfortran-4 Cygwin 4.5.3

2013-04-24 Thread Duncan Murdoch
On 13-04-24 1:51 PM, Jens Olofsson wrote: Dear Duncan, I know this isn't a forum for Cygwin, but for R. Pls treat me as a noob and also remember I am on Windows. How should I use R CMD SHLIB as if I write that at the prompt I get the error: unexpected symbol in R CMD. I have mango.f95 in the

[R] help with logkda.pari (in package untb)

2013-04-24 Thread Chase, Jennifer
Hello- I'm looking for help using the function logkda.pari. logkda.pari needs to access the program pari/gp (which of course I've installed), but I don't understand what commands I should use so that it recognizes that it is installed. Any help with this would be much appreciated. Cheers Jen

Re: [R] Trouble Computing Type III SS in a Cox Regression using drop1 and Anova

2013-04-24 Thread John Fox
Dear Paul, -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of Paul Miller Sent: Wednesday, April 24, 2013 1:18 PM To: r-help@r-project.org Subject: [R] Trouble Computing Type III SS in a Cox Regression using drop1 and Anova

Re: [R] Sum up column values according to row id

2013-04-24 Thread Rui Barradas
Hello, Note also that in the with(9 instruction there's some redundancy, th OP could simply do with(ipso, order(id, -g)) Hope this helps, Rui Barradas Em 24-04-2013 19:10, David Carlson escreveu: Something like this? mean6 - function(x) { if (length(x) 6) { mn - mean(x)

Re: [R] [r] Problems with Fortran calls when loaded a dll compiled with gfortran-4 Cygwin 4.5.3

2013-04-24 Thread Jens Olofsson
Ok. I apologise for not understanding. So, I have installed R-tools. It changed my PATH-variable. I didn't installed Cygwin dlls as stated by http://cran.r-project.org/doc/manuals/R-admin.html#The-Windows-toolset. Instead my PATH-variable contains the path to the Cygwin dlls AFTER the path to R...

Re: [R] the joy of spreadsheets (off-topic)

2013-04-24 Thread peter dalgaard
On Apr 24, 2013, at 20:01 , Thomas Adams wrote: One might wonder if the Excel error was indeed THAT or perhaps a way to get the desired results, give the other issues in their analysis? I think I'd reserve that suspicion for what they did with the NZ data: Growth for 1946-49: 7.7, 11.9,

Re: [R] R cannot find the path on my mac

2013-04-24 Thread MacQueen, Don
The why is that you haven't got the path correct yet. Assuming you are running R.app or R64.app (running the R Console), one way to find the correct path is to type file.choose() Navigate to your file and choose it. It will then tell you the correct path. You can also do things like

Re: [R] [r] Problems with Fortran calls when loaded a dll compiled with gfortran-4 Cygwin 4.5.3

2013-04-24 Thread R. Michael Weylandt michael.weyla...@gmail.com
On Apr 24, 2013, at 7:46 PM, Jens Olofsson jens.olofs...@gmail.com wrote: Ok. I apologise for not understanding. So, I have installed R-tools. It changed my PATH-variable. I didn't installed Cygwin dlls as stated by http://cran.r-project.org/doc/manuals/R-admin.html#The-Windows-toolset.

Re: [R] mgcv: how select significant predictor vars when using gam(...select=TRUE) using automatic optimization

2013-04-24 Thread Juliet Hannah
Hi Jan and Simon, If possible, could you attach the diagnostic plots. I would be curious to see them. Thanks, Juliet On Fri, Apr 19, 2013 at 4:39 AM, jholstei jan.holst...@awi.de wrote: Simon, that was very instructive—very special thanks to you. I already noticed that the model was bad,

Re: [R] [r] Problems with Fortran calls when loaded a dll compiled with gfortran-4 Cygwin 4.5.3

2013-04-24 Thread Jens Olofsson
Hello again. I may be a bit slow, but I'm learning. :-) Installed Rtools, added R to the path. Used CMD, moved to the folder of the source and wrote R CMD SHLIB Mango.f95 and voila I got a dll I loaded in R and I could call my current subroutines as intended. I am very grateful for ur help

Re: [R] RDA permutest envfit

2013-04-24 Thread Gavin Simpson
On Wed, 2013-04-24 at 10:38 -0700, Rémi Lesmerises wrote: Dear all, I did a RDA and when I looked to the signification of the test with permutest, the output was non-significant. But when I used the envfit function, some of the vectors are significant. All the test's conditions are

Re: [R] [r] Problems with Fortran calls when loaded a dll compiled with gfortran-4 Cygwin 4.5.3

2013-04-24 Thread Duncan Murdoch
On 13-04-24 2:46 PM, Jens Olofsson wrote: Ok. I apologise for not understanding. So, I have installed R-tools. It changed my PATH-variable. I didn't installed Cygwin dlls as stated by http://cran.r-project.org/doc/manuals/R-admin.html#The-Windows-toolset. Instead my PATH-variable contains the

Re: [R] Tables package - remove NAs and NaN

2013-04-24 Thread Santosh
Dear Rxperts, Sorry if I am posting a really really dumb request.. I am new to subversion and am trying to use subversion to download the tables package as suggested by Duncan. I installed subversion client(from collabnet) and tried to access tables package using the command below. svn checkout

Re: [R] pglm package: fitted values and residuals

2013-04-24 Thread Paul Johnson
On Wed, Apr 24, 2013 at 3:11 AM, alfonso.carf...@uniparthenope.it wrote: I'm using the package pglm and I'have estimated a random probit model. I need to save in a vector the fitted values and the residuals of the model but I can not do it. I tried with the command fitted.values using the

Re: [R] Tables package - remove NAs and NaN

2013-04-24 Thread Duncan Murdoch
On 13-04-24 3:23 PM, Santosh wrote: Dear Rxperts, Sorry if I am posting a really really dumb request.. I am new to subversion and am trying to use subversion to download the tables package as suggested by Duncan. I installed subversion client(from collabnet) and tried to access tables package

Re: [R] RDA permutest envfit

2013-04-24 Thread Rémi Lesmerises
Thank you! A last question: Is it still the same explanation if I remove the condition first=TRUE (and then testing for all axis) and permutest gives the same result?   Rémi Lesmerises, biol. M.Sc., Candidat Ph.D. en Biologie Université du Québec à Rimouski remilesmeri...@yahoo.ca On

Re: [R] identify object that causes Error in loadNamespace(name) : there is no package called ‘R.utils’

2013-04-24 Thread Duncan Murdoch
On 13-04-24 10:57 AM, Duncan Murdoch wrote: On 13-04-24 10:12 AM, Martin Morgan wrote: On 04/24/2013 06:03 AM, Duncan Murdoch wrote: On 13-04-24 5:46 AM, Liviu Andronic wrote: Dear all, I've bumped into the: Error in loadNamespace(name) : there is no package called ‘R.utils’ error. I've

Re: [R] RDA permutest envfit

2013-04-24 Thread Gavin Simpson
On Wed, 2013-04-24 at 13:13 -0700, Rémi Lesmerises wrote: Thank you! A last question: Is it still the same explanation if I remove the condition first=TRUE (and then testing for all axis) and permutest gives the same result? No; again `permutest(, first = FALSE)` and `envfit()` are

Re: [R] pglm package: fitted values and residuals

2013-04-24 Thread Achim Zeileis
On Wed, 24 Apr 2013, Paul Johnson wrote: On Wed, Apr 24, 2013 at 3:11 AM, alfonso.carf...@uniparthenope.it wrote: I'm using the package pglm and I'have estimated a random probit model. I need to save in a vector the fitted values and the residuals of the model but I can not do it. I tried

[R] getting started in parallel computing on a windows OS

2013-04-24 Thread Benjamin Caldwell
Dear R help, I've what I think is a fairly simple parallel problem, and am getting bogged down in documentation and packages for much more complex situations. I have a big matrix (30^5,5]. I have a function that will act on each row of that matrix sequentially and output the 'best' result from

Re: [R] Tables package - remove NAs and NaN

2013-04-24 Thread Santosh
yes, With Duncan's and Liviu's help, I was able to remove those NAs and NaNs from the tabular summary. svn .. thing has not worked for me yet.. would try this later.. Thanks so much! Regards, Santosh On Wed, Apr 24, 2013 at 1:46 PM, Santosh santosh2...@gmail.com wrote: Thanks so much.. I

Re: [R] Trouble Computing Type III SS in a Cox Regression

2013-04-24 Thread Terry Therneau
I should hope that there is trouble, since type III is an undefined concept for a Cox model. Since SAS Inc fostered the cult of type III they have recently added it as an option for phreg, but I am not able to find any hints in the phreg documentation of what exactly they are doing when you

[R] Floating point precision causing undesireable behaviour when printing as.POSIXlt times with microseconds?

2013-04-24 Thread O'Hanlon, Simon J
Dear list, When using as.POSIXlt with times measured down to microseconds the default format.POSIXlt seems to cause some possibly undesirable behaviour: According to the code in format.POSIXlt the maximum accuracy of printing fractional seconds is 1 microsecond, but if I do; options(

Re: [R] Sum up column values according to row id

2013-04-24 Thread arun
HI Not sure whether this helps: ipso- read.table(text=    id dbh h  g 1  FPE0164  36 13.62 0.10178760 2  FPE0164  31 12.70 0.07547676 21 FPE1127  57 18.85 0.25517586 13 FPE1127  39 15.54 0.11945906 12 FPE1127  34 14.78 0.09079203 6  FPE1127  32 15.12 0.08042477 5  FPE1127  28 14.13

Re: [R] getting started in parallel computing on a windows OS

2013-04-24 Thread Martin Morgan
On 04/24/2013 02:50 PM, Benjamin Caldwell wrote: Dear R help, I've what I think is a fairly simple parallel problem, and am getting bogged down in documentation and packages for much more complex situations. I have a big matrix (30^5,5]. I have a function that will act on each row of that

Re: [R] question

2013-04-24 Thread R. Michael Weylandt
On Tue, Apr 23, 2013 at 2:54 PM, nafiseh hagiaghamohammadi n_hajiaghamohammadi2...@yahoo.com wrote: Hi I fit one linear quantile regression with package quantreg and I want to khow this model is good or not.Is there method for checking it? Thanks your advice I ask this question because

Re: [R] fonts not rendering during plot

2013-04-24 Thread David Winsemius
On Apr 24, 2013, at 7:29 AM, George Dietz wrote: Hi, I am having a problem where sometimes fonts for text in plots don't get rendered-the text just shows up as boxes. If I call R from a certain directory the problem goes away, otherwise the fonts don't render (but plots get made). In

Re: [R] about the loglm and glm---Re:Re: Regression on stratified count data

2013-04-24 Thread meng
Understand. Many thanks for your detailed explaination. At 2013-04-24 22:22:55,Achim Zeileis achim.zeil...@uibk.ac.at wrote: On Wed, 24 Apr 2013, meng wrote: Hi,Achim: Can all the analysis you mentioned via loglm be performed via glm(...,family=poisson)? Yes. ## transform table

[R] Help me make faster R code for Kennard-Stone algorithm [My code is so slow from Matlab]

2013-04-24 Thread Kevin Hao
Hi all, Can you help me change my Kennard-Stone algorithm to faster one? [The original code can run fast in matlab, but when I change matlab code to R code, it is so slow.] Since my code so crude and too many loops (changed from matlab code), it is too slow. I hope that you can help to improve

Re: [R] question

2013-04-24 Thread Roger Koenker
good usually means good relative to something else, in this case the comparison seems, as Michael has already said, f0 - rq(y ~ 1, tau = ?) and then one can compute the R1 version that I originally suggested. But since there is still no explicit way to evaluate this, it is all a bit pointless.

[R] Stochastic Frontier: Finding the optimal scale/scale efficiency by frontier package

2013-04-24 Thread jpm miao
Hi, I am trying to find out the scale efficiency and optimal scale of banks by stochastic frontier analysis given the panel data of bank. I am free to choose any model of stochastic frontier analysis. The only approach I know to work with R is to estimate a translog production function by

Re: [R] the joy of spreadsheets (off-topic)

2013-04-24 Thread Ross Boylan
On 4/17/2013 5:18 AM, Kevin Wright wrote: On Tue, Apr 16, 2013 at 4:33 PM, Jim Lemon j...@bitwrit.com.au wrote: On 04/17/2013 03:25 AM, Sarah Goslee wrote: The final point does relate to Excel and any application that hides what is going on to the casual observer. I will treasure this URL to

[R] How to make a raster image in R

2013-04-24 Thread Kristi Glover
Hi R-User [[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

[R] creating raster image in R

2013-04-24 Thread Kristi Glover
Hi R-User I was trying to make a raster map with WGS84 projection in R, but I could not make it. I found one data set in Google that data is almost the same format as of mine. I wanted to make a raster map of temperature with 1 degree spatial resolution for the global scale. I could make

Re: [R] Distance matrices Combinations

2013-04-24 Thread arun
HI, Your code is not very clear. mata-m[,c(u)]       # assume that m is mat1 or w (as m was not defined) #also assume that 124 as nrow  of each matrix # For 10 distance matrices (#10C4) set.seed(25) lst1- lapply(1:10,function(i) dist(matrix(sample(1:50,5*124,replace=TRUE),nrow=124)))

Re: [R] fonts not rendering during plot

2013-04-24 Thread George Dietz
I apologize for not going in to too much detail-I didn't want to confuse people with the gory details of what I'm trying to do (trying to make a portable R installation bundled in to a mac app, among other things). I am on a mac, Mountain Lion. Anyhow, I figured out the problem was with the

[R] Assigning a variable value based on multiple columns

2013-04-24 Thread Jason Stout, M.D.
Hi All, I'm hoping someone can help me with a relatively simple problem. Take the following dataset: IDDiabetesESRDHIVContact 100NA0 210NA0 3NA 100 40

[R] Loop for main title in a plot

2013-04-24 Thread Eva Günther
Hi all, I have a problem in including my plot in a loop. Here is a simple example for one plot: # Plot simple graph with super- and subscript a-c(1,2,3,4) b-c(1,2,3,4) plot(x=a,y=b, ylab=expression(paste(Apple[P])), xlab=expression(paste(Banana^th)),

Re: [R] fonts not rendering during plot

2013-04-24 Thread David Winsemius
On Apr 24, 2013, at 8:42 PM, George Dietz wrote: I apologize for not going in to too much detail-I didn't want to confuse people with the gory details of what I'm trying to do You've slighted us. Kept us from our due. We LOVE gory OS and code details. Chop that chicken up. Throw dem bones.

Re: [R] Loop for main title in a plot

2013-04-24 Thread David Winsemius
On Apr 24, 2013, at 9:22 PM, Eva Günther wrote: Hi all, I have a problem in including my plot in a loop. Here is a simple example for one plot: # Plot simple graph with super- and subscript a-c(1,2,3,4) b-c(1,2,3,4) plot(x=a,y=b, ylab=expression(paste(Apple[P])),

[R] Linear Interpolation : Missing rates

2013-04-24 Thread Katherine Gobin
Dear R forum I have data.frame as df = data.frame(rate_name = c(USD_1w, USD_1w, USD_1w, USD_1w, USD_1m, USD_1m, USD_1m, USD_1m, USD_2m, USD_2m, USD_2m, USD_2m,  GBP_1w, GBP_1w, GBP_1w, GBP_1w, GBP_1m, GBP_1m, GBP_1m, GBP_1m, GBP_2m, GBP_2m, GBP_2m, GBP_2m, EURO_1w, EURO_1w, EURO_1w, EURO_1w,