[R] Missing data

2013-04-24 Thread Roslina Zakaria
Dear r-users, I would like to investigate about how to fill in missing data.  I started with a complete data and try to introduce missing data into the data series.  Then I would use some method to fill in the missing data and then compare with the original data how good it is.  My question is,

Re: [R] Regarding Modeling - Please! QUICK HELP

2013-04-24 Thread Rolf Turner
Sorry, this list has a "No homework" policy. Please ask your lecturer or tutor about this. cheers, Rolf Turner On 25/04/13 14:18, Andrew Cochrane wrote: I'm a student currently working with the *sleepstudy* dataset in matrix.pkg. It deals with the reaction times of sleep deprived

Re: [R] Regarding Modeling - Please! QUICK HELP

2013-04-24 Thread Patrick Coulombe
Hi Andrew, I don't know the dataset at all (and you seem to assume that your readers will), but anyway: it looks like you're trying to do an intercept-only model. If that's the case, try: > logmod11 <- lmer(log(Reaction) ~ 1 + (1|Subject),REML=FALSE) 1 is the intercept, and anything in parenthes

Re: [R] Assigning a variable value based on multiple columns

2013-04-24 Thread Patrick Coulombe
Hi Jason, I think that the easiest for you would be to keep your current elseif statements as is, but change your NA into something else (e.g., -999, or anything else). To do this in one line, you can use the package "gdata". In this code, I assume that your data are stored in the variable "datas

[R] Regarding Modeling - Please! QUICK HELP

2013-04-24 Thread Andrew Cochrane
I'm a student currently working with the *sleepstudy* dataset in matrix.pkg. It deals with the reaction times of sleep deprived students over a period of days. I am trying to model reaction times in order to describe the variation between students by days they havent slept. This is what I'm runni

[R] Bootstrapping in R

2013-04-24 Thread Preetam Pal
Hi all, 1>i have 3 vectors a,b and c, each of length 25... i want to define a new data frame z such that z[1] = (a[1] b[1] c[1]), z[2] = (a[2] b[2] c[2]) and so on...how do i do it in R 2> Then i want to draw bootstrap samples from z. Kindly suggest how i can do this in R. Thanks, Preetam

[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",

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"])),

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

[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")), main=expression(paste(i

[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

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 pat

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

[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 it

[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 http://www.R-project.org/posting-guide

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 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 give to anyone who

[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 s

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] 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 t

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" 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 back to data.f

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).

Re: [R] question

2013-04-24 Thread R. Michael Weylandt
> On Tue, Apr 23, 2013 at 2:54 PM, nafiseh hagiaghamohammadi > 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 there is 2 models,f

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 matr

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 0.

[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( digits.sec

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 i

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 wrote: > Thanks so much.. I will try it out. > Reg

[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 t

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, 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

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()` ar

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 alre

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 W

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 u

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

2013-04-24 Thread Paul Johnson
On Wed, Apr 24, 2013 at 3:11 AM, 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 following procedure without

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] [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 pat

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 res

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 an

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 wrote: > Simon, > > that was very instructive—very special thanks to you. > I already noticed that the model was bad, but it was not

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
On Apr 24, 2013, at 7: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-variabl

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 sourc

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] 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] 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] 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 Ano

[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 Chase

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 wo

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 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. t<-as.matrix((MAT1)^2+(MA

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 wrote: > In case you haven't noticed, this is making the rounds in the media, > including a handful of ref

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

[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 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 pack

[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

[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 worked

[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] 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 a

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: res1<-within(mutate(dat1,AB001A_1=1*(AB001A==AB001A[2]),AB0002A_1=1*(AB0002A==AB0002A[2]),AB362_1=1*(AB362==AB362[2]),SUM=rowSums(cbind(AB001A_1,AB0002A_1,AB362_1)),MATCH=(SUM/3)*100),{MATCH[1:2]<-NA;RANK=rank(M

[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

[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

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      

[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" wrote: >On Wed, 24 Apr 2013, meng wrote: > >> Hi all: >> For stratified count data,how to perform regression analysis? >> >> My data:

[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", lib.loc="/home/vahe/R/i686-pc-linux-gnu-libra

[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 FPE0

[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 an

[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] 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("PT

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 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

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 A&M University College Station, TX 77840-4352 From: J J [mailt

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.

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 thi

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 pro

[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 R. Michael Weylandt
On Wed, Apr 24, 2013 at 4:08 PM, Liviu Andronic wrote: > Dear Duncan, > > > On Wed, Apr 24, 2013 at 3:03 PM, Duncan Murdoch > 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 work

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 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 accidentally end up > w

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 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 loaded. (This > goes

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 ( http://www.cybaea.ne

[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 in

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

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 <- glm(

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 ( http://www.cybaea.net/Blogs/Data/A-warning-on-the-R-save-format.

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

Re: [R] News package

2013-04-24 Thread R. Michael Weylandt
On Wed, Apr 24, 2013 at 9:02 AM, Aswathy Nair 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 likely roll your

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 libra

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 di

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,

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

[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] 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:2

[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

[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 packag

[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 ~ C

[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

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 wrote: > Your example is not reproducible, so it is op

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 s