similar to: survey weights in sample with replacement

Dear helpers, I am trying to use the mda package downloaded from the R website, but the data set has missing values so I got an error message. Should I manually handle these missing values? I was trying to read the documents to specify any option related to missing values, but I did not find it. Please forgive me if I ignore something obvious. Thanks, Zhu Wang Statistical Science Department

Thanks to Brian D. Ripley and Uwe Ligges. It works now. Zhu Wang Statistical Science Department Southern Methodist University 3225 Daniel Avenue, Heroy 123 Dallas, TX 75275 Tel: (214)768-2453 Fax: (214)768-4035 > > > > -----Original Message----- > > From: Uwe Ligges [mailto:ligges at statistik.uni-dortmund.de] > > Sent: Sunday, February 02, 2003 5:22 AM >

Dear helpers, I am trying to create a library which uses some Fortran source files and Lapack and Blas subroutines. The Fortran source files from the original author contain subroutines isamax.f, sgefa.f and sgesl.f, which are part of BLAS subroutines on my Linux computer, but maybe different (old) versions. So in addition to these subroutines, there are other Lapack and Blas subroutines

Dear all, Is there any generic algorithm code for optimization implemented in R? I searched without success. Thanks, -- Zhu Wang Statistical Science Department Southern Methodist University Phone: (214)768-2453 Fax: (214)768-4035 Email: zhuw at mail.smu.edu

Dear R developers: I am trying to add a column to a data.frame. The following does the trick by expanding the name of the first data frame with the prefix foo: > data.frame(foo = as.list(data.frame(items = 1:3, bar=1:3)), items = 1:3) foo.items foo.bar items 1 1 1 1 2 2 2 2 3 3 3 3 However, the following special case produces an

Dear R Developers, For package seriation I use Fortran code. I recently got a request to add #if defined(__ICC) || defined(__INTEL_COMPILER) USE IFPORT #endif to the code since the Intel Fortran compiler otherwise has problems with rand(). However, to enable the FPP preprocessor I have to either add a compiler flag (-cpp for gFortran) which is possibly not portable or change the

Dear friends I have an array consist of r-rows and c-columns e.g. x=c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,0,0,0,0,0,0,0,0); x1=array(x, dim=c(4,6)) output is > x1 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1 2 3 4 0 0 [2,] 1 2 3 4 0 0 [3,] 1 2 3 4 0 0 [4,] 1 2 3 4 0 0 How can i ignore columns having zero sums? Help in this regard

I applied bootpred for multinomial logistic reg. (with nnet package). I used same as theta.fit and theta.predict of R for my data. but give me error. Can I do this with response vriable;7 levels predictor variables:5 (1 classifier, 4 continuous)?   Thanks alot Azam   [[alternative HTML version deleted]]

Dear friends There is a list of arrays comprising different no of rows and columns even sometimes NULL, such as [[2]] given below. How can we ignore [[2]] or others like this in the complete list. Any help in this regard is needed. Thanks [[1]] [,1] [,2] [1,] 3 1 [2,] 3 1 [3,] 3 1 [[2]] NULL [[3]] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 3 1

Dear R community I have a problem regarding which of the column in a matrix contains all of zero elements. e.g. x=c(3,3,3,3,0,0,0,0,5,5,5,5,8,8,8,8); x=matrix(x, nrow=4) the output is > x [,1] [,2] [,3] [,4] [1,] 3 0 5 8 [2,] 3 0 5 8 [3,] 3 0 5 8 [4,] 3 0 5 8 In this case the required column is second so the result should be "2".

Hi   I want to do leave-one-out cross-validation for multinomial logistic regression in R. I did multinomial logistic reg. by package nnet in R. How I do validation? by which function? response variable has 7 levels   please help me   Thanks alot Azam [[alternative HTML version deleted]]

Dear R users I have a problem regarding an addition of an extra "row" to a matrix. e.g. i have a matrix a <- matrix(1:6,2,3) > a [,1] [,2] [,3] [1,] 1 3 5 [2,] 2 4 6 I want to add a matrix having just one row. e.g. b <- matrix(7:9,1,3) > b [,1] [,2] [,3] [1,] 7 8 9 Now i want to get result like this [,1] [,2] [,3] [1,] 1 3 5

Dear R users I have a very basic question. I tried but could not find the required result. using dat <- pima f <- table(dat[,9]) > f 0 1 500 268 i want to find that class say "0" having maximum frequency i.e 500. I used >which.max(f) which provide 0 1 How can i get only the "0". Thanks and best regards Muhammad Azam Ph.D. Student Department of

Hi R-users!! I have some trouble with the survey pakage and i would be very glad if you can give me an advice. I have a sample from a survey where household were interviewed. The sample has 4 criteria on which the stratification was based: REGION, SIZE OF HOUSEHOLD, SIZE OF LOCALITY, AGE OF HEAD OF HOUSEHOLD. Since i don't have the whole information in each cell of the cross

Respected R helpers/ users I am one of the new R user. I have a problem regarding to know which of the integer in each column of the following matrix is in majority. I want to know that integer e.g. in the first column 1 is in majority. Similarly in the third column 4 is in majority. So what is the suitable way to get the desired integer for each column. I am looking for some kind reply. Thanks

Hi Martin, package arules heavily relies on ngCMatrix and uses multiplication and addition for logical operations. I think it makes sense that in a mixed operation with one dgCMatrix and one ngCMatrix the ngCMatrix should be "promoted" to a dgCMatrix. The current behavior of %*% and friends is in deed confusing: > m <- matrix(sample(c(0,1), 5*5, replace=TRUE), nrow=5) >

Dear helpers, When I use ar to fit the data with length 180, I have the following error: ar(x,method="burg") Error in acf(x, type = "covariance",lag.max=order,plot=FALSE): lag.max must be at least 1 If I use ar(x), then I have Call (x=x) order selected 0 sigma^2 estimated as 5374 Obviously I missed some points for using ar. This is R 1.7.0 under Redhat Linux

Dear R users, I am trying to convert a dataset from one format to several rectangular datasets. A consultant helped design the data entry program for our survey using Delphi/Pascal and for each household the information is stored in a file called "EA-HM-HH.TXT" where EA is the enumeration area number, HM is the homestead number and HH is the household number. Within this file the

Dear helpers, I am trying to write a function to return prediction values using package ts. I have written three different versions since I am not sure what's wrong with my func2. func and func1 return the same results.But func1 and func2 don't. In particular, the only difference between "func1" and "func2" is the function variable name being y and data, respectively.

Dear R community Hope every one be in best of his/her health. I have a situation in which there are s-sectors. Each sector is further divided into r-rows and c-columns. All it makes an array having dimension (r,c,s). e.g. x=c(1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,2,2,2,3,3,0,0,0,0,0,0,0,0,0,0,1,1,1,2,2,3,3,3,4,4,4,0,0,0,0,0,0,0,1,1,1,2,2,2,3,3,3,4,4,4,