similar to: Removing header from a matrix

Hi all, When I run a regression, I got the following response.  your help is appreciated.  Error in constants (lhs, cnames_symb) :   The hypothesis "Cost_PKL-QK.X_(Intercept)=0" is not well formed: contains bad coefficient/variable names.  In addition: Warning message: In constants(lhs, cnames_symb) : NAs introduced by coercion Dereje [[alternative HTML version deleted]]

Hi  I want to fit two equations simultaneously EQ1<-Y1~X1+X2 EQ2<-Y2~X1+X2 eqsystem<-list(Y1HAT=EQ1,Y2HAT=EQ2) fitols<-systemfit(eqsystem, method="OLS", data=BB) How do I get coefficients for the first equation? R code How do I restrict coefficient of X2 in the first equation (say , restrict it to less than  zero). R code Your help is appreciated. Dereje [[alternative

Dear R-users, I am facing problem with integrating in R a likelihood function which is a function of four parameters. It's giving me the result at the end but taking more than half an hour to run. I'm wondering is there any other efficient way deal with. The following is my code. I am ready to provide any other description of my function if you need to move forward.

Here is my model with interaction terms and control variables (I changed variables names for easy read): reg1 <- lm(y ~ x1*x2*x3 +control1 + control2 + control3) x1 ranges from 0 to 6; x2 from 0 to 5; and x3 from 0 to 4. All three are discrete ordinal variables; but I will treat them as continuous variables. (a) How can I see the predicted values of y for each of these scenarios (210

All, I am new to using R and know some basics. I wish to use kriging in R to do the following: given data Y =f(X1,X2,X3,.....,Xn) --1000+ irregular measured data set. I would like to be able to get a single value y given sinle input set (x1,x2,x3,...xn) A google search on this takes me lierally to the same example on involving analysis with soil sampling and I cannot figure out how to

Dear All, I'm wondering if there is a R function could give me all the combinations of the grouping/cluster result, given the number of the groups. e.g. 3 objects: x1 x2 x3, number of groups is 2 so the result will be group1:x1,x2; group2: x3 group1: x1;group2: x2,x3 group1: x1,x3;group2: x2 many thanks yan

Dear All, I need some help arranging data that was imported. The imported data frame looks something like this (the actual file is huge, so this is example data) DF: IDKey X1 Y1 X2 Y2 X3 Y3 X4 Y4 Name1 21 15 25 10 Name2 15 18 35 24 27 45 Name3 17 21 30 22 15 40 32 55 I would like to create a new data frame with the following NewDF: IDKey X Y Name1 21 15 Name1

Hi users. I'm new user in R. I'm workiing with Time series and I would like to know how can I do to undo the command DIFF(X), for exemple: If I have the model: m=arima(X, order=c(0,1,1), seasonal=list(order=c(0,0,1))) (note that have d=1 one difference), to find, in the same scale, the original numbers (like one "unDiff"), after the forecast, I need to develop some function or in

On Sat, 15 Jul 2017, Michael Reed via R-help wrote: > Dear All, > > I need some help arranging data that was imported. It would be helpful if you were to use dput to give us the sample data since you say you have already imported it. > The imported data frame looks something like this (the actual file is > huge, so this is example data) > > DF: > IDKey X1 Y1 X2 Y2

Suppose there are 4 variables d is a function of a , b and c I want to know how a, b and c change will make d change It will be straightforward to see it if we can graph the d surface if d is only a function of a and b, I can use 'persp' to see the surface of d. I can easily see at what values of a and b, d will get the maxium or minium or multiple modes, etc But for 4 dimention graph,

Dear R users/developers I requested help to solve the problem of formulating Multivariate Sample selection model by using Full Information Maximum Likelihood (FIML)estimation method. I could not get any response. I formulated the following code of FIML to analyse univariate sample selection problem. Would you please advise me where is my problem library (sem) library(nrmlepln) Selection

I would like to obtain Confidence Intervals for the estimates (unstandardized beta weights) of each predictor in a WLS regression: m1 = lm(x~ x1+x2+x3, weights=W, data=D) SPSS offers that output by default, and I am not able to find a way to do this in R. I read through predict.lm, but I do not find a way to get the CIs for multiple independent variables. Thank you Torvon [[alternative HTML

I'm failing to get a for loop working. I'm sure it's something simple, and I have found some posts relating to it, but I'm just not understanding why this isn't working. I have a data frame and would like to loop through specific column names, using aggregate() within a for loop. There are NA's scattered throughout the data frame and I'm thinking it has something to

Dear R-Users, I have a problem with an algorithm that iteratively goes over a data.frame and exclude n-columns each step based on a statistical criterion. So that the 'column-space' gets smaller and smaller with each iteration (like when you do stepwise regression). The problem is that in every round I use a new subset of my data.frame. However, as soon as I "generate" this

Hi how do we get intercepts standard error. I'm using the package pls. i got the coefficient but not able to get the stabdard error -- View this message in context: http://r.789695.n4.nabble.com/How-to-get-intecerpt-standard-error-in-PLS-tp3932104p3932104.html Sent from the R help mailing list archive at Nabble.com.

I would like to run a loop in R. I have never done this before, so I would be very grateful for your help ! 1. I have a sample set: 25 objects. I would like to draw 1 object from it and use it as a test set for my future external validation. The remaining 24 objects I would like to use as a training set (to select a model). I would like to repeat this process until all 25 objects are used as a

I have a matrix like this Name 1 2 3 4 5 NM_001039514 1.033557047 0.7469879518 0.9004524887 0.8613861386 0.7952499048 NM_001039723 1.0759493671 1.2315789474 0.8666666667 1.1142857143 0.9428011471 NM_001042605 0.9897435897 0.8870431894

Please forgive me if someone has seen this duplicated email, because I sent the email to the wrong address just now. Hi, all, I'm estimating an inter-rater coefficient, Aickin (1990)'s alpha. It's been a long time for me to figure out how to make it work, but it's still of no avail. The problem seems to be my code with the optim() function. Can anyone help me take a look at the

Hi guys, I have 1000 rows of a dataset. In my analysis, I need 70% of the data, run my analysis and then use the remaining 30% to test my model. Could anybody kindly help me on this? Cheers

Hi, I deal with 3D array say: , , 1 [,1] [,2] [,3] [,4] [,5] [1,] 1 4 7 10 13 [2,] 2 5 8 11 14 [3,] 3 6 9 12 15 , , 2 [,1] [,2] [,3] [,4] [,5] [1,] 16 19 22 25 28 [2,] 17 20 23 26 29 [3,] 18 21 24 27 30 , , 3 [,1] [,2] [,3] [,4] [,5] [1,] 31 34 37 40 43 [2,] 32 35 38 41 44 [3,] 33 36