similar to: y label after axis (4)

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 would like to know if anyone has any idea of how to run an OLS with constraints? thank you Soyoko ______________________________________ Ms. Soyoko Umeno Graduate Research Assitant for the Illinois-Missouri Biotechnology Alliance (IMBA) at http://www.imba.missouri.edu/ Ph.D. Student at the Department of Agricultural and Consumer Economics at the University of Illinois at Urbana-Champaign

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

The following warning message occurs when running the nls on some data: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Warning message: In is.na(wts) : is.na() applied to non-(list or vector) of type 'NULL' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ and I would like to know what is causing it and how I can fix it. As an example, the following, from Venables and Ripley, which

I am trying to plot a graph but the points on the graph should be different symbols and colors. It should represent what is in the legend. I tried using the points command but this does not work. Is there another command in R that would allow me to use different symbols and colors for the points? Thank you kindly. data(mtcars) plot(mtcars$wt,mtcars$mpg,xlab= "Weight(lbs/1000)",

Dear R Help Please forgive my lack of knowledge,.I would be very thankful for some help. Here is my problem: I was using optim to estimate parameters of a model and I get this error message "Error in optim(x0, fn = riskll, method = "L-BFGS-B", lower = lbs, upper = ubs, : L-BFGS-B needs finite values of 'fn'" Below is the R code I have written.

I'm running R 1.2.2 under windows 98 on a Pentium 133 laptop. I can't seem to retrieve the package MASS data sets: > library(MASS) > data(wtloss) Warning message: Data set `wtloss' not found in: data(wtloss) > data(abbey) Warning message: Data set `abbey' not found in: data(abbey) And yet all the .rda files for the MASS datasets are in D:\Program

I am resending this, as I believe it has not arrived on the mailing list when I first emailed. I have a set of labels arranged along a timeframe in a. Each label has a timestamp and marks a state until the next label. The dataframe a contains 5 such timestamps and 5 associated labels. This means, on a continious scale between 1-100, there are 5 markers. E.g. 'abc' marks the timestampls

All the usual MASS data sets are nicely sitting in D:\Program Files\rw1022\library\MASS\data on my laptop and yet: > library(MASS) > data(wtloss) Warning message: Data set `wtloss' not found in: data(wtloss) what could be wrong? I'm running 1.2.2 on Win98. Murray Jorgensen -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list --

We''ve finally finished up work on our application and are now beginning our beta. The site (SMS+mobile blogging+event planning+LBS+photos) is built using 100% Ruby on Rails. Many thanks to all the people who answered questions over the months on this list. Our experience with Rails was overwhelmingly positive! We''d love to hear any feedback you might have:

Dear useRs, How can I draw a barplot, but instead of bars, I'd get lines? Thanks, Kind regards, Rafael. ____________________________________________________________________________________ [[elided Yahoo spam]] [[alternative HTML version deleted]]

I have only just discovered the joys of the dotchart (since I am reading William Cleveland's -- Sean Carmody The Stubborn Mule http://www.stubbornmule.net http://twitter.com/seancarmody [[alternative HTML version deleted]]

Seems like this should be easy but I'm struggling a bit. How do I rearrange a data frame to go from the first one to the second shown below ? State Date lbs TX 200701 400 TX 200702 650 TX 200703 950 TX 200704 1000 FL 200701 200 FL 200702 300 FL 200703 500 FL 200704 333 NJ 200701 409 NJ 200702 308 NJ 200703 300 NJ 200704 800 Date TX FL NJ 200701 400 200 409 200702 650

I suspect this has been asked 400 times, but Google yields 2-3 rather complex gymnastics in blogs (and possibly dated), and searching list archives on ruby-forum proves kinda useless. Using standard forms, it''s fairly easy to use @object.errors[:field_name] in the form to embed an error message right next to the field that has the invalid entry. The controller uses a simple conditional

Hi I am trying to add a less than equal sign to a plot. I have previously done this using unicode but is not working in this instance. Any suggestions would be great thanks example code: library(ggplot2) df<-data.frame(vis=c(0,0,1,1) , count=c(10,15,20,10) , grp=c(0,1,0,1)) df$grp <-factor(df$grp ,levels=c(0,1) , labels =c("x \u2265 10" , "x > 10")) ggplot(df,

Dear Rxperts, Running the following code: ======================================================= twlo=10; twhi=20; wt=154; vd=0.5; cl=0.046; tau=6; t=3; F=1; wtkg <- wt/2.2 # convert lbs to kg vd.pt <- wtkg * vd # compute weight-based vd (L) cl.pt <- wtkg * cl # compute CL (L/hr) k <- cl.pt/vd.pt # compute k (hr^-1) cave <-