Re: [R] fast cumulative matrix multiplication

If the matrices are not all the same size, then the order of computation will make a difference. simple example: A is 1xn, B is nx1, C is 1xn. A(BC) takes n^3 multiplies, while (AB)C requires 2n. albyn Quoting Todd Schneider todd.w.schnei...@gmail.com: Hi all, I am looking for a

Re: [R] FW: making a plot in xyplot

On Mon, Oct 26, 2009 at 12:19 PM, Ordaz, Sarah ord...@upmc.edu wrote: Hi, I'd like to clarify my previous posting.  See below - my updates are noted with *s Thanks, Sarah Ordaz ord...@upmc.edu Hello, I am a newbie to the lattice package in R, and I'm trying to make a plot using the

Re: [R] PDF too large, PNG bad quality

On Fri, Oct 23, 2009 at 1:54 AM, Jim Lemon j...@bitwrit.com.au wrote: On 10/23/2009 06:07 AM, Lasse Kliemann wrote: I wish to save a scatter plot comprising approx. 2 million points in order to include it in a LaTeX document. Using 'pdf(...)' produces a file of size about 20 MB, which is

Re: [R] Print several xyplots to the same page in a pdf file

2009/10/27 Joel Fürstenberg-Hägg joel_furstenberg_h...@hotmail.com: Hello everybody, I'm using the lattice package and the xyplot to make several graphs like below. However, I can just print the three grouped plots onto one page as I'm putting them into a pdf-file, which gives me a huge

Re: [R] R crashes

Windows XP 32 bit machine. allocated 2gb as memory for R 2.9.2 data will be 20mb and was running calculations on the data . but the r crashes with out any warning or so. -- View this message in context: http://old.nabble.com/R-crashes-tp26110355p26124740.html Sent from the R help mailing list

Re: [R] how to control the label position on panels of trellise plot

On Thu, 29-Oct-2009 at 05:22PM -0400, zrl wrote: | Hi, | | I am wondering how to control the label position on panels of trellise plot. | I have 6 panels, I want to all x label positions at the bottom of each | panel. How should I do that? | Similar, how should I control y label positions of

Re: [R] R crashes

On 10/30/2009 05:17 AM, premmad wrote: Windows XP 32 bit machine. allocated 2gb as memory for R 2.9.2 can you try again whatever you do with 2.10.0 data will be 20mb and was running calculations on the data . but the r crashes with out any warning or so. That is a bit better, but nowhere

[R] Package zelig

hello all I am using the R package Zelig for some tobit regression with robust standard errors. I have got R version 2.9.2 (2009-08-24) and Zelig Version: 3.4-5 when i do demo(robust) It ends like this way data(coalition) # Fit the model with robust standard error user.prompt() Press return

[R] Names of list members in a plot using sapply

Hi R users: I got this code to generate a graphic for each member of a lists. list1-list(A=data.frame(x=c(1,2),y=c(5,6)),B=data.frame(x=c(8,9),y=c(12,6))) names1-names(list1) sapply(1:length(list1),function(i) with(list1[[i]],plot(x,y,type=l,main=paste(Graphic of,names1[i] Is there a more

Re: [R] similarity measure for binary data

The 'Gower' metric is one that is commonly used. R On Thu, Oct 29, 2009 at 10:22 AM, karuna m m_karuna2...@yahoo.com wrote: I am doing hierarchical clustering with cluster package.  I couldnot find similarity measures like matching coefficient, Jaccard coefficient and sokal and sneath.

Re: [R] similarity measure for binary data

The 'Gower' metric is one that is commonly used. Rod On Thu, Oct 29, 2009 at 10:22 AM, karuna m m_karuna2...@yahoo.com wrote: I am doing hierarchical clustering with cluster package.  I couldnot find similarity measures like matching coefficient, Jaccard coefficient and sokal and sneath.

Re: [R] ggplot2: Histogram with negative values on x-axis doesn't work

2009/10/30 hadley wickham h.wick...@gmail.com: I read anything that mentions ggplot2 no matter where it is. ... one should hope this statement only applies to the Internet though, does it? Please do share your regexp if it's not the case. :) baptiste

Re: [R] how to test for stationarity in time series?

On 10/30/09, Luna Laurent luna.laurent@gmail.com wrote: Could anybody tell me how to test for stationarity in time series? http://www.rseek.org/?cx=010923144343702598753%3Aboaz1reyxd4q=stationarity+time-seriessa=Search+functions%2C+lists%2C+and+morecof=FORID%3A11 Liviu

Re: [R] [R-SIG-Finance] Fast optimizer

Le 30 oct. 2009 à 02:21, R_help Help a écrit : Ok. I have the following likelihood function. L - p*dpois(x,a)*dpois(y,b+c)+(1-p)*dpois(x,a+c)*dpois(y,b) where I have 100 points of (x,y) and parameters c(a,b,c,p) to estimate. Constraints are: 0 p 1 a,b,c 0 c a c b I construct a

[R] Agglomerative Information Bottleneck (AIB)

Hi; Does anybody aware of R implementation for AIB (Agglomerative Information Bottleneck) by Slomni and Tishby? in which package should i find it? Thanks _

Re: [R] Removing generating data by category

Hmm, so if read correctly you want to remove exactly duplicated rows. So maybe try the following to begin with. duplicated(newdf[ , c(id, loc, clm)]) [1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE TRUE TRUE Then you can remove the duplicated rows before proceeding

[R] insert a text in panels, always in the same position (lattice, ltext, ?prepanel?)

Dear R-users, my present problem is related to lattice. I would like to put some text in each panel, namely a p-value. I therefore wrote a simple panel function as reported here below. I'm able to write the value in each panel at the maximum value of y for each panel, but this obviously

[R] whitespace with subscript text

Readers, I have the following command: expression(A[1]B~%) but that causes an error so I changed to: expression(A[1]~B~%) but the result is too much whitespace between subscript and B. Is there a way to reduce this whitespace? Yours, rhelp at conference.jabber.org r 251 mandriva 2008

[R] Applying a function on n nearest neighbours

I'm having a problem where I have to apply a function to a subset of a variable, where the subset is defined by the n nearest neighbours of a second variable. Here's an example applied to the 'iris' dataset: $ head(iris) Sepal.Length Sepal.Width Petal.Length Petal.Width Species 1

Re: [R] similarity measure for binary data

Hi Karuna, You can try vegdist() function in the vegan package. The function computes dissimilarity indices that are useful for or popular with community ecologists. All indices use quantitative data, although they would be named by the corresponding binary index, but you can calculate the

Re: [R] Temperature Prediction Model

Hi, also interested... If you are checking for non-normal behaviour, first let us define normal behaviour: only small temperature changes and no steep ramps? If so, maybe you can make an rolling average of the last x points, and check if the following point deviates more than ... ? Or is it

Re: [R] x-y plot as an histogram

Thanks for all the suggestions. Cut + plotmeans (gplots) solved my problem. groups - cut(t[,1],breaks=5) plotmeans(t[,2]~groups) Thanks! Federico --- El jue, 29/10/09, Jim Lemon j...@bitwrit.com.au escribió: De: Jim Lemon j...@bitwrit.com.au Asunto: Re: [R] x-y plot as an histogram Para:

Re: [R] how to test for stationarity in time series?

Hi Luna, tseries package provides various stationarity and unit root tests including Augmented Dickey-Fuller, Phillips-Perron, and KPSS. Alternative implementations of the ADF and KPSS tests are in the urca package, which also includes further methods such as Elliott-Rothenberg-Stock,

Re: [R] whitespace with subscript text

e-letter wrote: Readers, I have the following command: expression(A[1]B~%) but that causes an error so I changed to: expression(A[1]~B~%) but the result is too much whitespace between subscript and B. Is there a way to reduce this whitespace? Use * rather than ~ Uwe Ligges Yours,

Re: [R] Applying a function on n nearest neighbours

On Fri, 30 Oct 2009 10:28:49 +0100 Karl Ove Hufthammer k...@huftis.org wrote: $ (pos=which(order(abs(iris$Sepal.Length-x)) %in% 2:6)) This should of course be: (pos=order(abs(iris$Sepal.Length-x))[2:6]) -- Karl Ove Hufthammer __

[R] data.frame extracting data row-wise

Dear All, I am struggling with extracting data from a data frame: x=data.frame(a=1:11,b=100:110) What I want is a list/vector in this sence: 1 100 2 101 3 102... For single rows, this works fine: as.matrix(x)[1,] For, say 2 rows, this works fine: z=c(as.matrix(x)[1,],as.matrix(x)[2,]) But

Re: [R] data.frame extracting data row-wise

On Fri, 30 Oct 2009 10:42:12 +0100 Christian Lerch t@gmx.net wrote: I am struggling with extracting data from a data frame: x=data.frame(a=1:11,b=100:110) What I want is a list/vector in this sence: 1 100 2 101 3 102... y=t(as.matrix(x)) as.vector(y) Happy to help. -- Karl Ove

Re: [R] insert a text in panels, always in the same position (lattice, ltext, ?prepanel?)

On Fri, Oct 30, 2009 at 2:10 AM, Ottorino-Luca Pantani ottorino-luca.pant...@unifi.it wrote: Dear R-users, my present problem is related to lattice. I would like to put some text in each panel, namely a p-value. I therefore wrote a simple panel function as reported here below. I'm able to

Re: [R] data.frame extracting data row-wise

Nice workaround :-). Thank you. Best, Christian Original-Nachricht Datum: Fri, 30 Oct 2009 10:54:49 +0100 Von: Karl Ove Hufthammer k...@huftis.org An: r-h...@stat.math.ethz.ch Betreff: Re: [R] data.frame extracting data row-wise On Fri, 30 Oct 2009 10:42:12 +0100 Christian

[R] Time series temporal disaggregation

  Hi, This is a newbie question. I would to be able to convert annual time series of flow data into quarterly data. I wonder if there is any existing R-function which permits to do it? In what package ?   I the archive, i found that some poeple speak about tempDis package for performing time

Re: [R] Names of list members in a plot using sapply

sapply(names(list1), function(.data.){ with(list1[[.data.]], plot(x, y, type='l', main=paste(Graphic of, .data))) }) On Fri, Oct 30, 2009 at 3:41 AM, Kenneth Roy Cabrera Torres krcab...@une.net.co wrote: Hi R users: I got this code to generate a graphic for each member of a lists.

[R] [R-pkgs] ascii package updated

Hello, ascii package has been updated to version 0.3.2. ascii is a R package for writing asciidoc, txt2tags, sphinx or org documents with embeded R code. See http://eusebe.github.com/ascii/ for some examples. News since version 0.2 are: - new sphinx driver and output (see

Re: [R] Names of list members in a plot using sapply

You could lapply over the names: junk - lapply(names(list1), function(nm) with(list1[[nm]], plot(x, y, type = l, main = paste(Graphic of, nm which is a slight improvement though that still involves dealing with a list and its names separately. On Fri, Oct 30, 2009 at 3:41 AM, Kenneth

[R] Multicore package: sharing/modifying variable accross processes

Hi, I want to parallelize some computations when it's possible on multicore machines. Each computation produces a big objects that I don't want to store if not necessary: in the end only the object that best fits my data have to be returned. In non-parallel mode, a single gloabl object is

Re: [R] Time series temporal disaggregation

See: https://stat.ethz.ch/pipermail/r-sig-finance/2009q3/004677.html On Fri, Oct 30, 2009 at 6:09 AM, Axel Leroix axel.ler...@yahoo.fr wrote: Hi, This is a newbie question. I would to be able to convert annual time series of flow data into quarterly data. I wonder if there is any existing

[R] Fractal Package - Hurst exponent

Hi, everybody. My name is Ricardo and I'm doing doctorate at USP (Brazil). I have used a lot the package 'fractal', developed by William Constantine and Donald Percival, in my research. However, I have a doubt about the execution of RoverS function. I was reading an article

[R] Back-to-back graph with lattice package

Hi, I am an [R] and lattice beginner. I an doing an epidemiological study and want to present my data in a condensed and clear manner as possible. The example below is fake, but representative for the the data I want to show. I have a data frame with different categorical data and one column of

[R] possible memory leak in predict.gbm(), package gbm ?

Dear gbm users, When running predict.gbm() on a large dataset (150,000 rows, 300 columns, 500 trees), I notice that the memory used by R grows beyond reasonable limits. My 14GB of RAM are often not sufficient. I am interpreting this as a memory leak since there should be no reason to expand memory

Re: [R] different L2 regularization behavior between lrm, glmnet, and penalized? (original question)

Dear Robert, The differences have to do with diffent scaling defaults. lrm by default standardizes the covariates to unit sd before applying penalization. penalized by default does not do any standardization, but if asked standardizes on unit second central moment. In your example: x =

[R] Multidimensional scaling

Hi there, when conducting a multidimensional scaling analysis, it is possible to do this on a correlation matrix? I calculated a dissimilarity matrix by substracting the corr.matrix from 1. Then I calculated the distance matrix with d - dist(dissimilarity matrix) and ran the isoMDS function.

[R] opacity under dispersion command under plotrix

Is there anyway to make the lines in the dispersion command come forward in a plot and allow the fill in the dispersion parameter be transparent so all of the lines I am using to note confidence intervals are shown? --- Joe King, M.A. Ph.D. Student

[R] Quarterly data in PLM package

Dear all, Does anyone know if the PLM package (to run Panel Data Analysis) accepts quarterly data? The package vignette and documentation only use annual data -and the only time index available seems to work for years. José Mr José Luis Iparraguirre Senior Research Economist

[R] Function to manage seconds?

Hi, I could not find a function that returns a time in hours and minutes given seconds as argument. I looked also in the chron package. I would need something like this: mytime(61) #61 are seconds [1] 00:01:01 mytime(7385) #7385 are seconds [1] 02:03:05 Am I missing something? or should I

[R] R strucchange question: recursive-based CUSUM

Hello R users: I'm trying now to apply the package strucchange to see whether there is a structural change in linear regression. I have noted the following problem that arises in my case with recursive-based CUSUM: generic function recresid() in efp() generates an error, since (probably) it

Re: [R] R-help Digest, Vol 80, Issue 30

Dear friends, I will be very happy if anyone tell me the way to change work directory permanently? I mean not use the function setwd() which can only change temporary, when you close the console, it will the old directory. Sys.setenv(R_USER = '') also doesn't work. [[alternative HTML

Re: [R] Function to manage seconds?

Try this: library(chron) times(100/(24*60*60)) # 100 seconds [1] 00:01:40 On Fri, Oct 30, 2009 at 5:37 AM, Stefano Leonardi stefano.leona...@unipr.it wrote: Hi, I could not find a function that returns a time in hours and minutes given seconds as argument. I looked also in the chron

[R] .Rprofile replacement function setwd() causing errors

In my .Rprofile I have the following functions which display the current directory in the main R window title bar, and modify base::setwd() to keep this up to date. I like this because I can always tell where I am in the file system. cd - function(dir) { base::setwd(dir)

[R] Help with creating some loops

Hi All: I have a data frame called all_corn. This has 31 columns. The first column is a character key. The next 15 columns (stat1,stat2,...,stat15) are the statistical forecast. The last 15 columns (sls1,sls2,...,sls5) are actual sales. I want to calculate textbook tracking signal and cuulative

Re: [R] Summing identical IDs

On Oct 29, 2009, at 5:23 PM, PDXRugger wrote: Terrific help thank you. dupbuild-aggregate(DF$Acres, list(Bldgid), sum) This line worked best. Now im going to challenge everyone (i think?) Consider the following: Acres-c(100,101,100,130,156,.5,293,300,.09,100,12.5)

Re: [R] multiple pages with ggplot2 facet_wrap?

On Wed, Oct 28, 2009 at 8:19 PM, Bill Gillespie bi...@metrumrg.com wrote: I currently use lattice functions to produce multiple pages of plots using the layout argument to specify the number of rows and columns of panels, e.g., xyplot(price ~ carat | clarity, diamonds, layout = c(2, 2))

Re: [R] fast optimizer

Date: Fri, 30 Oct 2009 09:29:06 +0100 From: Christophe Dutang duta...@gmail.com Subject: Re: [R] [R-SIG-Finance] Fast optimizer To: R_help Help rhelp...@gmail.com Cc: r-help@r-project.org Ok. I have the following likelihood function. L -

[R] Odp: Help with creating some loops

Hi r-help-boun...@r-project.org napsal dne 30.10.2009 14:40:06: Hi All: I have a data frame called all_corn. This has 31 columns. The first column is a character key. The next 15 columns (stat1,stat2,...,stat15) are the statistical forecast. The last 15 columns (sls1,sls2,...,sls5) are

Re: [R] environment or other technique to isolate some private working information

In order to instrument functions for optimization or DE solving etc., I want to set up a private ensemble of working data somewhere in my workspace that won't collide with the problem-data. I can do this already with named objects within .GlobalEnv, but I'd like to have functions like

Re: [R] Back-to-back graph with lattice package

On Oct 30, 2009, at 9:11 AM, Erik Svensson wrote: Hi, I am an [R] and lattice beginner. I an doing an epidemiological study and want to present my data in a condensed and clear manner as possible. The example below is fake, but representative for the the data I want to show. I have a data

Re: [R] opacity under dispersion command under plotrix

On Oct 30, 2009, at 7:02 AM, Joe King wrote: Is there anyway to make the lines in the dispersion command come forward in a plot and allow the fill in the dispersion parameter be transparent so all of the lines I am using to note confidence intervals are shown? Got code?

[R] R_PHP

Hi All, Once again , would like to know about your suggestions regarding the integration of PHP and R.I have read ** *4. R Web Interfaces* here: http://cran.r-project.org/doc/FAQ/R-FAQ.txt It says about many tools.Here my requirement is :-- I have a tool developed in PHP.I want to

Re: [R] Help with creating some loops

Hi: In general, how to I cast a character to the class that I am trying to change. For example, if I have a data frame df1. df1 has a column x. suppose I want to a substring of x (the first 3 chars). Then I want to do something like df1$new = substring(of x) Example Data frame df1 x abcd efgh

Re: [R] Easy method to set user-mode virtual memory space in Windows Vista and 7

On Wed, 28 Oct 2009, Kingsford Jones wrote: I thought I'd share this with the list since it appears to provide a quick fix to some memory problems, and I haven't see it discussed in relation to R. Seems you didn't look in the places mentioned in the posting guide, fo this is an actual FAQ!

Re: [R] Help with creating some loops

On Oct 30, 2009, at 11:24 AM, Vadlamani, Satish {FLNA} wrote: Hi: In general, how to I cast a character to the class that I am trying to change. For example, if I have a data frame df1. df1 has a column x. suppose I want to a substring of x (the first 3 chars). Then I want to do

[R] How to properly shade the background panels of an xyplot?

Dear R users, this is a follow up of this message http://tolstoy.newcastle.edu.au/R/e6/help/09/05/13897.html I'm reproducing the core of it for convenience. // / data(Oats, package = MEMSS) / / tp1.oats - xyplot(yield ~ nitro | Variety + Block, / / data = Oats, / /

Re: [R] .Rprofile replacement function setwd() causing errors

On 10/30/2009 9:30 AM, Michael Friendly wrote: In my .Rprofile I have the following functions which display the current directory in the main R window title bar, and modify base::setwd() to keep this up to date. I like this because I can always tell where I am in the file system. cd -

[R] how to loop thru a matrix or data frame , and append calculations to a new data frame?

How do you do a double loop through a matrix or data frame , and within each iteration , do a calculation and append it to a new, second data frame? (So, if your original matrix or data frame is 4 x 5 , then 20 calculations are done, and the new data frame, which had 0 rows to start with, now has

Re: [R] Bayesian regression stepwise function?

On Fri, Oct 23, 2009 at 6:05 PM, Ravi Varadhan rvarad...@jhmi.edu wrote: Frank, I have heard this (i.e. only head-to-head comparisons are valid) and various other folklores about AIC and BIC based model selection before, including one that these information criteria are only applicable for

Re: [R] Can I write Nothing to a variable?

OK, here's the first-cut code, which does execute properly. What I would like to do, simply for 'cleanliness', is to assign nothing to the variable l , rather than assigning seq(1,dim(data)[timedim]). Then the result of all the get(unlist(dims)) operations below will be func[j,k,] Where

Re: [R] : source R-script interactively outside R

At least in the Unix version of the ``Rscript`` front-end to R (but oddly enough not ``R CMD BATCH``) you could put something like this in test.R: cat(Save workspace image? [y/n]: ) answer - tolower( scan(stdin, , n=1) ) if (substr(answer, 1, 1) == y) { q(yes) } else { q(no) } Which

Re: [R] how to loop thru a matrix or data frame , and append calculations to a new data frame?

On Oct 30, 2009, at 11:46 AM, Robert Wilkins wrote: How do you do a double loop through a matrix or data frame , and within each iteration , do a calculation and append it to a new, second data frame? (So, if your original matrix or data frame is 4 x 5 , then 20 calculations are done, and the

Re: [R] R crashes

premmad wrote: Windows XP 32 bit machine. allocated 2gb as memory for R 2.9.2 data will be 20mb and was running calculations on the data . but the r crashes with out any warning or so. Ahh, well 20 MB is not huge data for R, so there should be no problem with processing it. The

Re: [R] deriv() to take vector of expressions as 1st arg?

Thanks a lot for the hints. The sapply method may work for me. But how can I extract just the .grad expression from the return value of the deriv function? (and secondly store in a matrix of expressions?) Thanks again. __ R-help@r-project.org mailing

[R] MatLab SimBiology

Is there any R package that implements the same capability of MatLab toolbox called SimBiology ? We are expecially interested in protein-protein interactions and network analysis. As far as I know SimBiology implements a system of ODEs reflecting the kinetic chemical reactions. We would be more

[R] polar.plot

Hi, Two questions: 1 - Say I have average speed and directions for tide and I would like to plot them on a polar plot, but with different colors so I can indicate the two directions. I'm using polar.plot from the plotrix library. How can I add a second b and dir.b series to a polar.plot?

[R] Interpreting gnls() output in comparison to nls()

Hi, I've been trying to work with the gnls() function in the nlme package. My decision to use gnls() was so that I could fit varPower and such to some of the data. However, in working with a small dataset, I've found that the results given by gnls() don't seem to make any sense and they

Re: [R] Easy method to set user-mode virtual memory space in Windows Vista and 7

Thank you for the response, Prof Ripley. Some comments below On Fri, Oct 30, 2009 at 9:27 AM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote: On Wed, 28 Oct 2009, Kingsford Jones wrote: I thought I'd share this with the list since it appears to provide a quick fix to some memory problems, and

Re: [R] R-help Digest, Vol 80, Issue 30

Which OS and how do you start R? -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On Behalf Of Huidong TIAN Sent: Friday,

[R] Multidimensional scaling on a correlation matrix

Hi there, when conducting a multidimensional scaling analysis, it is possible to do this on a correlation matrix? I calculated a dissimilarity matrix by substracting the corr.matrix from 1. Then I calculated the distance matrix with d - dist(dissimilarity matrix) and ran the isoMDS function.

Re: [R] R strucchange question: recursive-based CUSUM

Julia: I'm trying now to apply the package strucchange to see whether there is a structural change in linear regression. I have noted the following problem that arises in my case with recursive-based CUSUM: generic function recresid() in efp() generates an error, since (probably) it cannot

[R] Help with RGDAL

Hello, We (two mac users) have been attempting to install rgdal from http://www.kyngchaos.com/software:frameworks;, given that it is not available as a binary on the CRAN (binaries) of the Package Installer. I have also tried to solve this problem by looking on the net for an old question,

[R] problems with re-loading exported data

Hello, I have a data file originally as .mat file from matlab that I imported in R using the R.matlab package (readMat). This loaded a list containing amongst other things a 4x4 array of rainfall data that is quite heavy (around 40 Mb, it took 10 minutes to load on my computer). I tried

[R] identify which commands have generated graphics?

Dear all Is it possible to programmatically detect which commands have generated a graphic to the x11() graphics device? When in front of the computer, it is easy to see that after a command---say, plot(1:10)---a graphics window opens/activates and displays a graphic. But is there a way to detect

Re: [R] Help with RGDAL

On Oct 30, 2009, at 1:51 PM, Pablo Alvarez wrote: Hello, We (two mac users) have been attempting to install rgdal from http://www.kyngchaos.com/software:frameworks , given that it is not available as a binary on the CRAN (binaries) of the Package Installer. I have also tried to solve this

Re: [R] problems with re-loading exported data

Can you at least show us what you were doing. It sounds like you were trying to 'load' a text file, not something that was saved with 'save'. Here is what I get: load('trackV3.r') Error: bad restore file magic number (file may be corrupted) -- no data loaded In addition: Warning message: file

Re: [R] Efficient way to code using optim()

Hi all, I am trying to estimate a simple logit model. By using MLE, I am maximizing the log likelihood, with optim(). The thing is, each observation has different set of choice options, so I need a loop inside the objective function, which I think slows down the optimization process. The

Re: [R] identify which commands have generated graphics?

On 10/30/2009 3:13 PM, Liviu Andronic wrote: Dear all Is it possible to programmatically detect which commands have generated a graphic to the x11() graphics device? When in front of the computer, it is easy to see that after a command---say, plot(1:10)---a graphics window opens/activates and

Re: [R] opacity under dispersion command under plotrix

Heres my code: era1 - seq(1.5,5,.25) plot(era1,yhyp1$pe, xlab=ERA, ylab=Probability of Winning Cy Young,type=l, col=black, lwd=2) dispersion(era1,yhyp1$pe,y1upper,y1lower,type=l, fill=blue,arrow.cap=0.01,intervals=FALSE)#Requres plotrix dispersion(era1,yhyp2$pe,y2upper,y2lower,type=l,

[R] What should my.data.frame[c(3, 1)] - list(NULL, NULL) evaluate to?

(R 2.9.2 on Windows XP) Turns out that my.data.frame[c(1,3)] - list(NULL, NULL) and my.data.frame[c(3,1)] - list(NULL, NULL) are not the same. Is this expected? Example: d0 - structure(list(Year = c(2009L, 2009L, 2009L), Season = c(1L, 1L, 1L), Series = c(1L, 1L, 1L), Mst =

Re: [R] opacity under dispersion command under plotrix

On Oct 30, 2009, at 3:49 PM, Joe King wrote: Heres my code: era1 - seq(1.5,5,.25) plot(era1,yhyp1$pe, xlab=ERA, ylab=Probability of Winning Cy Young,type=l, col=black, lwd=2) This plots. library(plotrix) # should go here. dispersion(era1,yhyp1$pe,y1upper,y1lower,type=l,

[R] exact string match?

Dear R users: I need to compare character strings stored in 2 separate data frames. I need an exact match, so finding a in animal is no good. I've tried regexpr, match, and grepl, but to no avail. Anybody know how to accomplish this? Ben -- View this message in context:

[R] possible to increase precision

I would like to increase the precision/accuracy of R. That is, I'm dealing with numbers exceedingly close to 1, and I would like to increase the number of significant digits used in computation. In matlab, you can use vpa() to accomplish this. I'm wondering if there's an equivalent or a

[R] NA values in Standard Error for zeroinfl()

I am fitting a model using zeroinfl() and it runs without errors, returning results that are generally consistent with my hypotheses. One of my variables is percent black (pblack). This variable was highly significant in some of the other count models I ran on the way to my current

[R] possible to increase precision

I would like to increase the precision/accuracy of R. That is, I'm dealing with numbers exceedingly close to 1, and I would like to increase the number of significant digits used in computation. In matlab, you can use vpa() to accomplish this. I'm wondering if there's an equivalent or a

Re: [R] Summing identical IDs

David, You are correct. I think the frist two assumptions can be thrown out and only the latter two (c,d) can be considered. So how would i combine Acres for matching Bldgids based on assumptions c,d? David Winsemius wrote: On Oct 29, 2009, at 5:23 PM, PDXRugger wrote: Terrific

Re: [R] Efficient way to code using optim()

Unless this is a homework problem, you would be much better off using glm(). Giovanni Date: Fri, 30 Oct 2009 12:23:45 -0700 From: parkbomee bbom...@hotmail.com Sender: r-help-boun...@r-project.org Importance: Normal Precedence: list --Boundary_(ID_/D+lL9iK1qLhrkPBeoxH+Q)

[R] Obtaining illogical results from posterior LDA-classification because of too good data?

Dear list, my problem seems to be primarily a statistical one, but maybe there is a misspecification within R (and hopefully a solution). I have two groups with two measured variables as training data. According to the variables, the groups differ totally. I know that this is a very easy

Re: [R] possible to increase precision

Hi Greg, Others on this list are much more knowledgeable than I am, and I'm sure they will correct me if I'm wrong. My understanding is that R uses double precision (see http://en.wikipedia.org/wiki/Double_precision_floating-point_format), which has the following implication: x -

Re: [R] Efficient way to code using optim()

Thank you. But I'd prefer using a written function which allows me more flexible model specification. Later on, I could have random parameters. So I want to know if there is any more efficient way so that I can speed it up. Date: Fri, 30 Oct 2009 16:10:29 -0600 To: bbom...@hotmail.com CC:

Re: [R] possible to increase precision

On 30/10/2009 5:04 PM, Greg Michaelson wrote: I would like to increase the precision/accuracy of R. That is, I'm dealing with numbers exceedingly close to 1, and I would like to increase the number of significant digits used in computation. In matlab, you can use vpa() to accomplish this.

Re: [R] possible to increase precision

I think he means something like packages Rmpfr and gmp provide. However, careful numerical analysis (and working with x-1 or log(x) for numbers very near one) means that such packages are rarely necessary (although they can be very convenient). It all depends on the calculations: vague

Re: [R] Help with RGDAL

Pablo Alvarez-2 wrote: Hello, We (two mac users) have been attempting to install rgdal from http://www.kyngchaos.com/software:frameworks;, given that it is not available as a binary on the CRAN (binaries) of the Package Installer. The GDAL library contains an impressive collection

[R] Multiple linear regression with constraint (imposition)

Hi R users I want to do a multiple linear regression with R With a normal model (Y=Cste+A1*X1+A2*X2) the program would be lm(Y ~X2+X2) My model is Y=Cste+A1*X1+A2*X2 with the constraint A1=-A2 What is the program for such a model? Best regards -- View this message in context:

Re: [R] Help with RGDAL (SOLVED)

I had earlier downloaded what I thought to be the most recent version of rgdal (rgdal_0.6-18.tgz). After checking to see that all of the other dependencies in the configure files were met, I finally got around to checking the version of the rgdal and found my problem. It was different (and

Re: [R] Multiple linear regression with constraint (imposition)

Sorry there was a mistake: Hi R users I want to do a multiple linear regression with R With a normal model (Y=Cste+A1*X1+A2*X2) the program would be lm(Y~X1+X2) My model is Y=Cste+A1*X1+A2*X2 with the constraint A1=-A2 What is the program for such a model? Best regards -- View this

  1   2   >