similar to: How do I test against a simple null that two regressions coefficients are equal?

Hi there, I am looking for a function that takes a vector as input and generates all pair wise combination of the elements of the input vector. For example, the input vector is c(1,2,3,4). The output vector is c(1 and 2, 1 and 3, 1 and 4, 2 and 3, 2 and 4, 3 and 4). The representation here is generic. I can use a for loop to do it, but it gets slow when the number of elements of the input

Hi there, Maybe people who know both R and econometrics will be able to answer my questions. I want to run panel regressions in R with fixed-effect. I know two ways to do it. First, I can include factor(grouping_variable) in my regression equation. Second, I plan to subtract group mean from my variables and run OLS panel regression with function lm(). I plan to do it with the second way because

Hi there, When I extract data from SQLite to R, the data types (or modes) of the extracted data seems to be determined by the value of the first row. Please see the following example. When I put the missing values first, the column extracted is of the mode character. > str(dbGetQuery(sql.industry, + "select pya_var from annual_data3 + order by

I'd like to write the summary of a dataframe to disk, so that it looks essentially the same as what you'd see on screen; but I can't seem to do it. Can someone tell me how? Thanks! ------------------------ Philip A. Viton City Planning, Ohio State University 275 West Woodruff Avenue, Columbus OH 43210 viton.1 at osu.edu

i am sure that there is a simple solution that i am just overlooking, but there is no logical test for = there is one for != and >,<, >=, <=... what do i do if i need to check if something IS EQUAL or is there an easy way, in working with a matrix of size n x 3, to sort the matrix by one of the columns? specifically what i am trying to do is i have a matrix with individual and

Dear Sir/Madam, I am running a R program for my SNP data. There are some errors when I run glm model in Hapassoc software, sometimes it is over the memory and sometimes the matrix is singular. I want to ingore these errors and excute the next statement. Otherwise, I have a big big trouble. Do you have some idea about this problem of ingore errors. Wish to get your help assp. thanks. -- Wei

On: R 2.8.1 / Ms Windows / R-Gui-R-console I have a long-ish function, and to re-assure myself that it's actually making progress I've arranged a set of messages to the console (using cat). But what happens is that the messages don't actually display until the function has finished, which sort of negates the purpose. However, if I hit Return after starting the function, then all

I'm new to bugs, so please bear with me. Can someone tell me if the following two models are doing the same thing? The reason I ask is that with the same data, the first (based on 4 separate coeffs a1--a4) takes about 50 secs, while the second (based on a vectorized form, a[]) takes about 300. The means are about the same, though R-hat's in the second version are quite a bit better.

Can someone tell me how to make up (eg) a library's html help files by hand? I think I ought to be able to use RCMD Rdconv for this but (R-2.10.0, MS-win) when I type (in a dos session) "rdcmd rdconv --help" I get a message to the effect that a perl script rdconv can't be opened. Can I do this from within R itself? And if so, how (in particular, what is the target file to

I downloaded a US states+territories shapefile from http://www.weather.gov/geodata/catalog/national/html/us-state.htm and unzipped it into my working directory. Then in R (version 12.2.0 for ms-win32): library(rgdal) # Version: 0.6-28 library(spdep) # Version: 0.5-24 (loads sp: Version: 0.9-72) # read in the data states<-readOGR(".","s_01de10") summary(states) # now

Hello everybody, I'm trying to subset a data frame "foo", which looks like this: sent cond raw V1 10 2 A 614 7 11 2 A 724 6 12 2 A 641 5 13 2 A 665 9 33 5 B 510 7 34 5 B 659 7 35 5 B 607 6 37 5 B 467 8 38 5 B 586 9 39 5 B 523 10 57 8 C 858 8 58 8 C 513 9 59 8 C

Is there a method for creating "animations" using R. Or is there another piece of software (unix perferrably) that can generate a "movie" from many ps/eps/gif/jpg files? Thanks, Jeff. Jeff D. Hamann Hamann, Donald and Associates PO Box 1421 Corvallis, Oregon USA 97339-1421 Bus. 541-753-7333 Cell. 541-740-5988 jeff_hamann at hamanndonald.com www.hamanndonald.com

Hi all, Could you point me to a script editor in Linux that will allow me to a somewhat easy pasting of section(s) of an R script or function into R for execution? Something like WinEdt+RWinEdt macro in R for Windows? (My time is a bit pressed. I haven't succeeded in my attempts to install ESS for XEmacs. No doubt due to my confusion and lack of competence. I am abandoning it for the

Hi all, I'm building a plot of the values in tmeant (below) against positions 1 to 5, using matplot. tmeant looks like this: case1 case2 pos1 861.8466 818.5909 pos2 961.2841 976.3466 pos3 878.6080 1262.8523 pos4 950.8011 1129.6080 pos5 968.1080 1063.3920 I also have lower (object tl) and upper (object tu) bounds on the confidence intervals as follows: tl: pos1

Group, I am new to R. In my work as a program evaluator, I am regularly asked to estimate effect sizes of prevention/intervention and educational programs on various student outcomes (e.g. academic achievement). In many cases, I have access to data over three or more time periods (e.g. growth in proficiency test scores). I usually have multiple independent and dependent variables in each

Dear all: I have this: A1 B1 C1 D1 E1 A2 B2 C2 D2 E2 A3 B3 C3 D3 E3 And I want this A1 E1 B1 E1 C1 E1 D1 E1 A2 E2 B2 E2 C2 E2 D2 E2 A3 E3 B3 E3 C3 E3 D3 E3 Example: m<- matrix(1:15,nrow=3,byrow=T) m v<- unlist(list(t(m[,1:4]))) u<- rep(c(5,10,15),c(4,4,4)) data.frame(v,u) This is the result I want but I would like to learn a simpler way to do it. Any clue?

Hi R-Users I need to estimate Lyapunov exponent of my time series. After reading description of all functions available I still don't know how to determine time delay. My time series length is 4200. Is it possible to determine time delay with other function's output or I can choose any random value? Thanks for your help! -- View this message in context:

Hi All, I have a data frame with a group of variables named b1, b2, b3, ..., b18. These variables take the value 1, 2 or NA. For each observation, I want to do some computation by looping over the values for the group of variables: b1 to b18. In STATA I would do: forval i=1/18 { --- use b`i' for computation ---- } How can this be done in R? Deepankar

We have a wireless network that uses 802.1x authentication, in which domain joined computers use their machine credentials to connect. Windows machines do this automatically, and until recently Linux computers could join using wicd, wpa-supplicant, and a simple script that would retrieve the machine password with tdbdump. ( specifically tdbdump -k SECRETS/MACHINE_PASSWORD/DOMAIN

#I have a data set that looks like this. A bit more complicated actually with # three factor levels but these calculations need to be done on one factor at a #I then have a set of different rates that are applied #to it. #dataset cata <- c( 1,1,6,1,1,2) catb <- c( 1,2,3,4,5,6) doga <- c(3,5,3,6,4, 0) data1 <- data.frame(cata, catb, doga) rm(cata,catb,doga) data1 # start rates #