similar to: plot history

_______________________________________________________________________________________ Hi, I was just wondering if there was any available R code that could handle general constrained optimisation problems. At the moment I'm using nlminb and optim, both of which allow box constraints on the parameters, but ideally I'd like to be able to specify more general constraints on the solution

_______________________________________________________________________________________ Hi, I'm writing up the online help for a package I'm developing (in-house only, sorry), and I've come across an odd glitch when trying to nest a list inside the "arguments" section of the .Rd file. I was just wondering if anyone could provide some insights. I'm using R 2.2.0 on

# Hello R Users, # # I would like to fit a glm model with quasi family and # logistical link function, but this does not seam to work # with binary data. # # Please don't suggest to use the quasibinomial family. This # works out, but when applied to the true data, the # variance function does not seams to be # appropriate. # # I couldn't see in the # theory why this does not work. # Is

_______________________________________________________________________________________ Note: This e-mail is subject to the disclaimer contained at the bottom of this message. _______________________________________________________________________________________ Hi, I came across some rather unexpected behaviour the other day with assigning names and lists. Here's an example. > z

Hi list, could anyone please educate me on the following: lst<-seq(47, 239, by=12) for(n in lst) { lower=n; upper=lower+10 for(i in lower+2 : upper) { print(paste(n, " i: ", i, " lower: ",lower, " upper :", upper)) } } does not stop when i = upper A while loop fixes this but, I still don't understand why the for loop doesn't stop

Hello, I'm new to R and work with psychometrics (Item response theory). I was using a package called "cirt" to model a data that contains only 0s, 1s and 9s. The package is basically used for Monte Carlo Markov Chain estimation of parameters and I specified 5000 iterations. The model seems to run fine but when I try to recover odds ratio statistics using the function

Thomas, As F Harrel in a preceding message told me that you are the maintainer of Hmisc package, I write directly to you, with copy to the Rhelp list in case someone encountered the same problem and find some benefit in the response. I used quite often the function sas.get in Hmisc library which I is very, very useful. But, as trying to reuse it today it seems not to function anymore. I tried

Hello, I'm a grad student in the Intelligent Transportation Systems lab at Portland State Univ. in Portland, OR, USA. I'm trying to learn the basics of R to run a negative binomial in the near future, and so I ran a test regression on roadway crash data obtained from "Statistical and Econometric Methods for Transportation Data Analysis" by Washington et al (p. 250). I ran the

I am just starting with R and have the following problem: given a matrix of ones and zeroes, say mdim=4 m<-matrix(round(runif(mdim^mdim)),mdim,mdim) how to construct the matrix of those rows of m, whose elements sum to 2. Contrary to the random matrix above, the actual matrix always has at least one such row. Serguei Kaniovski

Dear List: We are generating data such that students are clustered in schools for some item response data for a simulation study. One component of our simulation is to generate measurement error from a logistic distribution with a mean of 0 and standard deviation of 1.7 to match the logistic curve of the Rasch model. We are generating an error term for each of the 40 hypothetical test items a

dear all, how can i get a vector that shows the number of the column of matrix that contains the maximum of the row ?? can?t believe in need a loop for this... i have a 100 x 3 matrix and want to get a 100 x 1 vector with values 1,2,3 . there must be a simple solution. i just cannot find it. i think am searching on the wrong end. thx for help in advance. m.

Dear Friends. I found a puzzling phenomenon in R when you use 'if' within a function: # defining a function aaa aaa=function(a) {if (a==1) {aaa=1}; if (a!=1) {aaa=2} } # using the function: > b=20 > bbb=aaa(b) > bbb [1] 2 > typeof(bbb) [1] "double" > > > c=1 > ccc=aaa(c) > ccc NULL > typeof(ccc) [1] "NULL" It seems that only the last

_______________________________________________________________________________________ Hello, Not sure if this counts as a bug or not, but I just noticed in R 2.4.0 that : and interaction are not quite equivalent. For example: > x <- factor(letters[1:4]) > y <- factor(letters[11:14]) > x:y [1] a:k b:l c:m d:n Levels: a:k a:l a:m a:n b:k b:l b:m b:n c:k c:l c:m c:n d:k d:l d:m

i've been trying to get samba 3 to join my AD domain, and have gotten stuck. when I wbinfo -t it returns "Could not check secret" log.winbindd says: nsswitch/winbindd.c: nsswitch/winbindd.c:process_loop(722) process_loop: Invalid request size from pid 5368: 1304 bytes sent, should be 1568 here's my smb.conf incase it matters. # Global parameters [global] workgroup =

Dear all, I tried to use "step" function to do model selection, but I got an error massage. What I don't understand is that data as data.frame worked well for my other programs, how come I cannot make it run this time. Could you please tell me how I can fix it? ***************************************************************************************************

What is a simple way to convert a factor into a matrix of dummy variables? fm<-lm(y~f) where f is a factor takes care of this in the estimation. I'd like to save the result of expanding f into a matrix for later use. Thanks. Charles -- Charles H. Franklin Professor, Political Science University of Wisconsin, Madison franklin at polisci.wisc.edu chfrankl at wisc.edu

I am running a two-way anova with Type III sums of squares and would like to be able to understand what the different SS mean when I use different contrasts, e.g. treatment contrasts vs helmert contrasts. I have read John Fox's "An R and S-Plus Companion to Applied Regression" approach -p. 140- suggesting that treatment contrasts do not usually result in meaningful results with Type

Hi all, I have been trying to calculate Type III SS in R for an unbalanced two-way anova. However, the Type III SS are lower for the first factor compared to type I but higher for the second factor (see below). I have the impression that Type III are always lower than Type I - is that right? And a clarification about how to fit Type III SS. Fitting model<-aov(y~a*b) in the base package and

I am running a binomial glm with response variable the no of mites of two species y->cbind(mitea,miteb) against two continuous variables (temperature and predatory mites) - see below. My model shows overdispersion as the residual deviance is 48.81 on 5 degrees of freedom. If I use quasibinomial to account for overdispersion the dispersion parameter estimate is 2501139, which seems

I am using Project fedena <http://projectfedena.org/> for School Management software a ROR application. i want to attach a biometric device for students & employee attendance , any one can suggest me how i can do it. -- You received this message because you are subscribed to the Google Groups "Ruby on Rails: Talk" group. To view this discussion on the web visit