similar to: constrained OLS

Hi, I would like to know if anyone has any idea of how to run an OLS with constraints? I need to contraint a coefficient estimate in the model equal to 1, and I am not sure how to include it into the OLS estimation... I was hoping to find something like "cnsreg" in STATA.. thank you Soyoko ______________________________________ Ms. Soyoko Umeno Graduate Research Assitant for the

Hi everyone, I have been trying to compute numbers of days between two dates as follows: > x <- c("1jan1960", "2jan1960", "31mar1960", "30jul1960") > z <- format(x, "%d%b%Y") > ex <- c("1jan1961", "15jan1960", "21mar1975", "10jul1981") > ez <- format(ex, "%d%b%Y") > ez-z

< In Stata, I can just use "bs" function ...> in STATA the bs command runs a bootstrap If you want to use the bootstrap function (or you could linearize b*c/a :) look under ?boot (you may need to load the package first). Usage: boot(data, statistic, R, sim="ordinary", stype="i", strata=rep(1,n), L=NULL, m=0, weights=NULL,

Hi, I am trying to figure out how to lable the second y-axis after the following codes: plot(x, y, xlab="time", ylab="pay1" ) par(new=TRUE) plot(x,y2, ann=FALSE, xaxt="n", yaxt="n", pch=7 ) axis(4) Then, I want to label the second y axis "pay2". I tried "title(ylab="pay2")", but it put this lable on the first

Hi, I would like to know if anyone knows how to compute a variance of the non-linear combination of the coefficient estimates. Say, I get a model of y~c+ax+bz (1) where x and z are the independent variables, c is the constant estimate, and a and b are the coefficient estimates. Then, I want to know the variance of b*c/a (2). How am I going to get it? In Stata, I can just use

Hi everyone, I am having a problem reading my data file in R. It only reads part of it, and stops doing it in half way through the data. I looked at the observation where reading stops, but it is really not different from the observations read it. I expanded the matrix size, but it did not help. I am using a text format to read the data file. Could someone help me? Thank you Soyoko

Hi, I would like to know if anyone knows how to draw a plot with two Y-axises and one X-axis? When you have two sets of y values that do not have the same scale, but correspond to the same x value, I would like to plot them on one graph. Could you please help me? Thank you Soyoko

Hello R helpers, I am trying to do a linear OLS regression of y on two variables x1 and x2. I want to constrain the coefficients of x1 and x2 to sum up to 1. and therefore run a constrained OLS. Can anybody help with this? (I have seen some answers to similar questions but it was not clear to me what I need to do) - I have tried the lm function with offset but I must not have used it properly.

A collegue of I is working on an assitant for managing the pxelinux configuration. Here come a screenshot of a devel release, what do you think of such tool ? http://people.mandrakesoft.com/~aginies/drakpxe2.jpg PS: Some options are dedicated to the Mandrakelinux installer but the assistant could be used for some generic configuration. -- Erwan Velu Linux Cluster Distribution Project Manager

Hello, For my puppet workflow, i need to run "module-assistant auto-install openvswitch-datapath". Which resource type will be best to execute this command. Should i use Exec resource type. Can you please suggest what will be the best way to implement this. Appreciate your response. Thanks in advance, Mani -- You received this message because you are subscribed to the Google

I might be sorry for asking this question :-) I have two equations and I tried to estimate them individually with "lm" and "gls", and then in a system (using systemfit) with "OLS", "WLS" and "SUR". Quite surprisingly (for myself at least) the results are identical to the last digit. Could someone (please!) give a hint as to what am I

Hi all, we are facing a strange issue about trust domain user issue. In the beginning. we are using samba-winbind integrated with Window AD(Server 2012r2 and Server 2008R2), they are working pretty well. but recently， the winbind client cloud not get the correct trust domian's user group info as well, like: id A\\user only show: uid=16077216(A\user) gid=16077216(A\domain

CentOS Errata and Bugfix Advisory 2012:0364 Upstream details at : https://rhn.redhat.com/errata/RHBA-2012-0364.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 0f0e04cefb5e7fc3e70a32e8a0dc83720edfbe2627cc28d8d7e18f6ce7f9b55d selinux-policy-3.7.19-126.el6_2.10.noarch.rpm

CentOS Errata and Bugfix Advisory 2020:0364 Upstream details at : https://access.redhat.com/errata/RHBA-2020:0364 The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) x86_64: 0f7494bbc8cc7270bc6e9d0a4b27dcfcd2ed798a3e82f9d09ef22111378d2cbd net-snmp-5.7.2-43.el7_7.3.x86_64.rpm

How can I extract the adjusted R^2 value from an ols object (using rms package)? library(rms) x <- rnorm(10) y <- x + rnorm(10) ols1 <- ols(y ~ x) Typing "ols1" displays adjusted R^2 among other things, but how can I assign it to a variable? I tried str(ols1) but couldn't see where to go from there. Thanks, Mark Seeto

Hi to all, In other statistical software, such as Eviews, it is possible to regress a model with the Least Squares method, assuming that the residuals follow an AR(q) process. For example the resulting regression is something like y = 1.2154 + 0.2215 x + 0.251 AR(1) How is it possible to do the same in R? Thank you very much in advance, Constantine Tsardounis http://www.costis.name

Hi, I am trying to fit an AR model, maximum order =4, order selection criterion is aic. I wonder why these two give different results: m1<-ar.ols(x, aic=TRUE, method="ols", order.max=4) m1<-auto.arima(x,d=0, D=0, max.p=4, max.P=0, max.q=0, max.Q=0, ic="aic") Could they both use the function predict to do forecasting? Is there any function that works better?

Hello everybody! I am currently using the quantreg package on R 2.2.1. for my diploma thesis and want to create plots with plot.summary.rqs, which should work like this: Usage: plot.summary.rqs(x, nrow=3, ncol=2, alpha = 0.1, ols = TRUE, ...) Arguments: x: an object produced by rq() fitting nrow: rows in mfrow ncol: columns in mfrow alpha: alpha level of

Hi, I have 142 observations off different variables and I'm trying to do a OLS. And I get this error. Any Ideas ? > f <- ols(lnmigr3~popden78+income+modern+spareha+rain) Error in La.chol2inv(x, size) : size cannot exceed nrow(x) = 1 D. ----------------------------------------- Stay ahead of the information curve. Receive GIS news and jobs on your desktop daily. Subscribe today

Hi folks! Does anybody know how to use the OLS function in fMultivar to do dynamic regression? I've tried specifying lags in OLS using a data series created in fSeries and it doesn't seem to work. I've done dynamic regression using dyn$lm and I was wondering how to accomplish the same thing using the OLS function from fMultivar. Thanks! John [[alternative HTML version