[R] Error when running coeftest in plm

2012-04-07 Thread Hien Nguyen
Dear R-helpers, I am using plm package to run panel data analysis for different models with different dependent variables. After running the PLM, I run the COEFTEST(plm1, vcovBK) to get the panel corrected standard errors. I got the PCSEs for some of the models but some tests produce error mes

Re: [R] .Rhistory in R.app

2012-04-07 Thread Berend Hasselman
On 08-04-2012, at 03:31, nkschmidt wrote: > Hi Maria, you need to name your history file somefilename.Rapp.history, > example > > savehistory(file = "history.Rapp.history") > > /nkschmidt You replied to a message almost two years old, which to my mind is pretty useless. Original message is h

Re: [R] newbie question: strategy

2012-04-07 Thread sysot1t
yes, I would have expected documented samples for something simple.. as I said, I am not asking for anything complex but something rather simple that I can use to learn from and build upon.. not a highly complex example that provides me with little explanation of what is going on... quanstrat is gr

Re: [R] .Rhistory in R.app

2012-04-07 Thread nkschmidt
Hi Maria, you need to name your history file somefilename.Rapp.history, example savehistory(file = "history.Rapp.history") /nkschmidt Maria Gouskova wrote > > Dear R users, > > I am having a minor but annoying issue with R.app. It doesn't retain > the history information from the previous se

Re: [R] Drawing a line in xyplot

2012-04-07 Thread David Winsemius
On Apr 7, 2012, at 10:29 PM, wcheckle wrote: Thank you David, the bwplot option does what I need: x11(height=8,width=11) bwplot ( mortality~ attend|type, pch=95,cex=5,col=2, par.settings=list( box.rectangle = list(col = "transparent"), box.umbrella = list(col = "transparent"), plot.symbol =

Re: [R] Drawing a line in xyplot

2012-04-07 Thread wcheckle
Thank you David, the bwplot option does what I need: x11(height=8,width=11) bwplot ( mortality~ attend|type, pch=95,cex=5,col=2, par.settings=list( box.rectangle = list(col = "transparent"), box.umbrella = list(col = "transparent"), plot.symbol = list(col = "transparent") ), panel=functi

Re: [R] newbie question: strategy

2012-04-07 Thread Joshua Ulrich
On Apr 7, 2012 7:55 PM, "sysot1t" wrote: > > I am not asking anyone to deliver anything, to the best of my understanding You asked for documented code examples of the items you listed in your original email. > all that I requested in terms of functionality is more than likely done in > quantstra

Re: [R] r ols help

2012-04-07 Thread kebrab67
I see what you mean as I go trgouh my courses (chapters related to type i and II errors) Is there any way in the therory for me to choose which variables in a set of variables explain more my Y. Or should I base my reflexion solely on economic theory and discriminate between variables? -- View thi

Re: [R] newbie question: strategy

2012-04-07 Thread sysot1t
I am not asking anyone to deliver anything, to the best of my understanding all that I requested in terms of functionality is more than likely done in quantstrat as you point out. Issues is: quantstrat has extremely poor documentation. Yes, there are "samples", but they are certainly not simple and

Re: [R] assigment operator question

2012-04-07 Thread Duncan Murdoch
On 12-04-07 4:51 PM, Henrik Bengtsson wrote: On Sat, Apr 7, 2012 at 1:30 PM, Mark Heckmann wrote: Hello, using the<<- assignment operator I do not understand why the following does not work. l<<- list() l list() l$arg1<<- "test" error in l$arg1<<- "test" : Objekt 'l' not found ?"<<-" says:

Re: [R] Imputing missing values using "LSmeans" (i.e., population marginal means) - advice in R?

2012-04-07 Thread Søren Højsgaard
Dear Jenn, Could you please provide a reproducible example. In your case, you could for example provide the data (using dput(dat)). Regards Søren -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Jenn Barrett Sent: 6. april 2012 1

Re: [R] quadratic model with plateau

2012-04-07 Thread Gabor Grothendieck
Yes, I had fixed that on my end but didn't notice I had copied the old version. I think its obvious enough that anyone will fix it themselves. Regards. On Sat, Apr 7, 2012 at 7:23 PM, Rolf Turner wrote: > > One other minor glitch I just noticed: > >    x <- seq_along(x) > > should read > >    x

Re: [R] quadratic model with plateau

2012-04-07 Thread Rolf Turner
One other minor glitch I just noticed: x <- seq_along(x) should read x <- seq_along(y) cheers, Rolf On 08/04/12 10:08, Gabor Grothendieck wrote: On Sat, Apr 7, 2012 at 3:58 PM, Gabor Grothendieck wrote: On Fri, Apr 6, 2012 at 9:42 PM, help ly wrote: Dear All, I wo

Re: [R] Bayesian 95% Credible interval

2012-04-07 Thread Greg Snow
The "emp.hpd" function in the TeachingDemos package will do this (assumes a single interval result, either unimodal or multimodes but the valleys between don't drop far enough to split the interval). I am sure there are similar functions in other packages as well. On Fri, Apr 6, 2012 at 12:39 PM,

Re: [R] r ols help

2012-04-07 Thread David Winsemius
On Apr 7, 2012, at 3:58 PM, kebrab67 wrote: The eleventh variable is log(wage) and the 10 other variables are actually all at least theorically significant to the regression. like wether the person is black, immmigrant, it´s experience, age, if he works in a city... But i need only take a

Re: [R] how to compute a vector of min values ?

2012-04-07 Thread Greg Snow
Peter showed how to get the minimums from a list or data frame using sapply, here is a way to copy your 1440 vectors into a single list (doing this and keeping your data in a list instead of separate vectors will make your life easier in general): my.list <- lapply( 1:1440, function(x) get( sprint

[R] r ols help

2012-04-07 Thread kebrab67
The eleventh variable is log(wage) and the 10 other variables are actually all at least theorically significant to the regression. like wether the person is black, immmigrant, it´s experience, age, if he works in a city... But i need only take a few variables out of this ten. Actually taking only

Re: [R] assigment operator question

2012-04-07 Thread Gabor Grothendieck
On Sat, Apr 7, 2012 at 5:01 PM, Mark Heckmann wrote: > Thanks! I'll try to stick to that advice! > Maybe there is a better way... Here is what I want: > > I want to save some default settings for a package. > The user can change these using a function similar to par(). > I do not want to use optio

Re: [R] quadratic model with plateau

2012-04-07 Thread Gabor Grothendieck
On Sat, Apr 7, 2012 at 3:58 PM, Gabor Grothendieck wrote: > On Fri, Apr 6, 2012 at 9:42 PM, help ly wrote: >> Dear All, >> >> I would like to make a quadratic with a plateau model in R. Is there a >> package in R doing this? The bentcableAR package seems won't work. >> >> The link below describes

Re: [R] assigment operator question

2012-04-07 Thread Henrik Bengtsson
On Sat, Apr 7, 2012 at 2:01 PM, Mark Heckmann wrote: > Thanks! I'll try to stick to that advice! > Maybe there is a better way... Here is what I want: > > I want to save some default settings for a package. > The user can change these using a function similar to par(). > I do not want to use optio

Re: [R] How do Sweave users collaborate with Word users?

2012-04-07 Thread mlell08
would it be possible to separate the parts generated by R from the actually written paths? As creating small doc's by sweave which contain the generated reports and embedding them via OLE into the static written text. However, this could require sending around multiple documents around each time...

Re: [R] How do Sweave users collaborate with Word users?

2012-04-07 Thread Yihui Xie
I cannot reply in detail at the moment, but we have a proposal to the Google Summer of Code this year which will address the collaboration issue. The basic idea is to write everything in plain text with markdown (including R code), and compile the file with the knitr package, then convert it to Wor

Re: [R] assigment operator question

2012-04-07 Thread Mark Heckmann
Thanks! I'll try to stick to that advice! Maybe there is a better way... Here is what I want: I want to save some default settings for a package. The user can change these using a function similar to par(). I do not want to use options() here as it will be quite a lot of parameters. I was thinking

[R] How do Sweave users collaborate with Word users?

2012-04-07 Thread Alexander Shenkin
Hello All, I'm getting my workflow switched over to Sweave, which is very cool. However, I collaborate with folks (as many of you must as well) who use Word to Track Changes amongst a group while crafting a paper. In the simplest case, there will just be two people (one Sweave user and one Word u

Re: [R] assigment operator question

2012-04-07 Thread Henrik Bengtsson
On Sat, Apr 7, 2012 at 1:30 PM, Mark Heckmann wrote: > Hello, > > using the <<- assignment operator I do not understand why the following does > not work. > > l <<- list() > l > list() > l$arg1 <<- "test" > error in l$arg1 <<- "test" : Objekt 'l' not found > > ?"<<-" says:  "The operators <<- and

[R] assigment operator question

2012-04-07 Thread Mark Heckmann
Hello, using the <<- assignment operator I do not understand why the following does not work. l <<- list() l list() l$arg1 <<- "test" error in l$arg1 <<- "test" : Objekt 'l' not found ?"<<-" says: "The operators <<- and ->> cause a search to made through the environment for an existing defini

Re: [R] quadratic model with plateau

2012-04-07 Thread Gabor Grothendieck
On Fri, Apr 6, 2012 at 9:42 PM, help ly wrote: > Dear All, > > I would like to make a quadratic with a plateau model in R. Is there a > package in R doing this? The bentcableAR package seems won't work. > > The link below describes what I am looking for in R exactly: > http://support.sas.com/docum

Re: [R] R help

2012-04-07 Thread Peter Ehlers
On 2012-04-07 08:27, kebrab67 wrote: I have a 11*1562 set of data. I want to regress using ols all combination of the ten variables on the eleventh. Store the results of the t test of each regression and the r2 to compare which combination is better. the combination can use any amount o variables

Re: [R] Drawing a line in xyplot

2012-04-07 Thread ilai
On Sat, Apr 7, 2012 at 11:16 AM, David Winsemius wrote: > xyplot(mortality ~ type, data=xdat, >               panel=function(x,y){ >                   panel.xyplot(x,y, jitter.x=TRUE) >                   panel.segments(x0=c(.9, 1.9, 2.9), >                                   x1=c(1.1,2.1,3.1), >  

Re: [R] Uniroot error

2012-04-07 Thread Berend Hasselman
On 07-04-2012, at 19:48, jeka12386 wrote: > Dear All > > I am trying to find a uniroot of a function within another function (see > example) but I am getting an error message (f()values at end points not of > opposite sign). I was wondering if you would be able to advise how redefine > my functi

[R] Uniroot error

2012-04-07 Thread jeka12386
Dear All I am trying to find a uniroot of a function within another function (see example) but I am getting an error message (f()values at end points not of opposite sign). I was wondering if you would be able to advise how redefine my function so that I can find the solution. In short my first fu

Re: [R] Order sapply

2012-04-07 Thread Rui Barradas
Hello, > >I have the following code, but it seems that something must be doing > wrong, because it is giving the results I want. > The idea is to create segments while the value of Commutation is less than > 1000. > for example, from the small set of data below > I believe that code was min

Re: [R] Drawing a line in xyplot

2012-04-07 Thread David Winsemius
On Apr 7, 2012, at 11:35 AM, wcheckle wrote: here is the data (fyi this is made-up data) x = as.data.frame ( cbind (c (5,8,7,5,8,10,11,6,4,5,20,25,27,30,35,32,28,21,20,34,11,15,18,12,15,12,10,15,19,20 ), c (1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3 ),c(1,0,1,1,1,1,0,0,0,1,0,1

Re: [R] R help

2012-04-07 Thread R. Michael Weylandt
This is, more often that not, statistically speaking a bad idea: you're likely to get a false positive (https://xkcd.com/882/) But it sounds like something along these lines should work for you to get the models : x <- data.frame(rnorm(50), rnorm(50), rnorm(50), rnorm(50)) m <- vector("list", nco

Re: [R] regression on combinations of variables; was: R help

2012-04-07 Thread Uwe Ligges
On 07.04.2012 17:27, kebrab67 wrote: I have a 11*1562 set of data. I want to regress using ols all combination of the ten variables on the eleventh. Store the results of the t test of each regression and the r2 to compare which combination is better. the combination can use any amount o variabl

Re: [R] rgeos - gBuffer, width by z-value

2012-04-07 Thread Omphalodes Verna
Thanks for help... do anybody also have an idea how to export "object l" (which is "list" class) to .*shp file. I use SpatialPolygons, but output in console is : Error in SpatialPolygons(list(l)) :   cannot get a slot ("area") from an object of type "list" Cheers, OV ___

[R] R help

2012-04-07 Thread kebrab67
I have a 11*1562 set of data. I want to regress using ols all combination of the ten variables on the eleventh. Store the results of the t test of each regression and the r2 to compare which combination is better. the combination can use any amount o variables. How can I do that in R? -- View this

[R] Rsonlp package about how to use optional " ineqfun"

2012-04-07 Thread marco1728
I was trying to put some inequality constraints in the solnp package. But whenever I put them into the package. It just showed that there is unused arguments in it. So can anyone give me a simple example with inequality constraints ??? i want to know the right form of writing the inequalit

Re: [R] Drawing a line in xyplot

2012-04-07 Thread wcheckle
here is the data (fyi this is made-up data) x = as.data.frame ( cbind(c(5,8,7,5,8,10,11,6,4,5,20,25,27,30,35,32,28,21,20,34,11,15,18,12,15,12,10,15,19,20), c(1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3),c(1,0,1,1,1,1,0,0,0,1,0,1,0,0,0,0,1,0,0,0,1,1,1,0,1,0,1,0,0,0))) names(x)=c("mo

Re: [R] pause code in R code?

2012-04-07 Thread Hurr
Yes, that is the question. So is it the readline() ? -- View this message in context: http://r.789695.n4.nabble.com/pause-code-in-R-code-tp4523209p4539615.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

Re: [R] Systemfit with structural equations and cross equation parameter interaction

2012-04-07 Thread Arne Henningsen
Dear Waseem On 7 April 2012 09:20, waseem ahmad wrote: > I want to estimate simultaneous equation model with panel data. > The model looks as follows > > Y1=a0+a1*X1+a2*X2 > Y2=b0+b1*X2+b2*X1 > X1=Z1-(Y1/a1) > X2=Z2-(Y2/b1) > I > In this model Y1, Y2, X1 and X2 are endogenous variables; > Z1, Z2

Re: [R] pause code in R code?

2012-04-07 Thread Bert Gunter
On Sat, Apr 7, 2012 at 6:18 AM, Hurr wrote: > readline("press return to continue") > worked as a pause until I used it in 4 places in a different program where > none of them wait. > What problem do I look for? Your bad coding. cf the posting guide -- especially "reproducible examples" -- to pos

Re: [R] How do I get a rough quick utility plot of a time series?

2012-04-07 Thread R. Michael Weylandt
One last try -- though if you cannot figure out how to take the name of your object and put it in dput( NameOfObjectGoesHere ) after I gave you a worked example of how to do so, it seems unlikely that this will help: # Attempt to recreate your data library(zoo) x <- read.zoo(textConnection("monoMn

Re: [R] pause code in R code?

2012-04-07 Thread R. Michael Weylandt
Are you sure you are executing that line of code? Michael On Sat, Apr 7, 2012 at 9:18 AM, Hurr wrote: > readline("press return to continue") > worked as a pause until I used it in 4 places in a different program where > none of them wait. > What problem do I look for? > > > -- > View this me

Re: [R] pause code in R code?

2012-04-07 Thread Hurr
readline("press return to continue") worked as a pause until I used it in 4 places in a different program where none of them wait. What problem do I look for? -- View this message in context: http://r.789695.n4.nabble.com/pause-code-in-R-code-tp4523209p4539399.html Sent from the R help mailing l

Re: [R] Drawing a line in xyplot

2012-04-07 Thread David Winsemius
On Apr 6, 2012, at 9:41 PM, wcheckle wrote: i am trying to replicate the following graph using xyplot : attach(x) plot ( jitter(type), mortality, pch=16, xlim = c(0.25, 3.75)) lines ( c(1-0.375,1.375) , c ( median(mortality[type==1]), median(mortality[type==1])), lwd=5,col=2) lines ( c(2-0.37

Re: [R] rgeos - gBuffer, width by z-value

2012-04-07 Thread Michael Sumner
Just use a loop and collect the results in a list. BTW, R-Sig-Geo is a special interest group mailing list specifically for these sorts of questions regarding sp and rgeos packages (which you did not mention). library(sp) library(rgeos) point <- data.frame(x=c(10,20), y=c(10, 10), z = c(2,7)) p

Re: [R] newbie question: strategy

2012-04-07 Thread R. Michael Weylandt
So you want us to deliver an entire backtesting architecture for you? That's a pretty hefty request... But being R, it's already been done. This has all basically been made available in the quantstrat project (google it) -- the documentation is online and you can see lots of worked examples in the

Re: [R] Building R on Solaris (OpenIndiana) with gcc 4.6.2 for amd64 target - relocation problems solved

2012-04-07 Thread Michael Figiel
On 04/ 7/12 08:13 AM, Prof Brian Ripley wrote: On 06/04/2012 19:38, Michael Figiel wrote: Hello, the "R Installation and Administration" handbook states in Section C 5.1: For ‘amd64’ the builds have failed to complete in several different ways, currently with relocation errors for libRblas.so

[R] Systemfit with structural equations and cross equation parameter interaction

2012-04-07 Thread waseem ahmad
Hi there, I want to estimate simultaneous equation model with panel data. The model looks as follows Y1=a0+a1*X1+a2*X2 Y2=b0+b1*X2+b2*X1 X1=Z1-(Y1/a1) X2=Z2-(Y2/b1) I In this model Y1, Y2, X1 and X2 are endogenous variables; Z1, Z2 are exogenous variables and a0, a1, a2, b0, b1 and b2 are param

Re: [R] Sum of sd between matrix cols vs spearman correlation between them

2012-04-07 Thread ali_protocol
Hi Petr, Thank you. Yes, normalized.matrix" and "normalized" are the same, my mistakes. Operators: like Cyclic loess normalization and quantile Normalization, from limma or affy packages. Would cv= sd/mean be a better replacement for sd? On 4/5/12, Petr Savicky [via R] wrote: > > > On Thu, Apr

[R] rgeos - gBuffer, width by z-value

2012-04-07 Thread Omphalodes Verna
Dear list! I have problem with buffer size (width) in package rgeos. I would like to expand given geometry (points) to specified width based on the z value from attribute table. Here is example: point <- data.frame(x=c(10,20), y=c(10, 10), z = c(2,7)) point_shp <- SpatialPointsDataFrame(point[