similar to: Questions about deviance

My attempt to test a model using ordglm code is running into problems, and I thought if you have a moment you might illucidate the situation. Here is the data: http://www.people.fas.harvard.edu/~corr/6.4.05.RData Here is the code: # I coerce tcn8 matrix data to a vector, because ordglm will not accept matrix data. y<-as.vector(tcn8[,62]) x<-as.vector(tcn8[,60])

Dear friends. Do Polr in Mass change the sign of the coefficients ? Example (McCullagh 1980) options(contrasts=c("contr.treatment","contr.poly")) library(Mass) freq <- c(19,29,24,497,560,269) yy <- ordered(gl(3,1,6)) z4 <- polr(yy~x,weights=freq) > z4 Call: polr(formula = yy ~ x, weights = freq) Coefficients: x2 -0.6026492 Intercepts: 1|2

Dear Peter. I guess you know that Jim Lindseys code include nordr and ordglm in library gnlm - I attach the htmls which do various linear and nonlinear ordinal regressions - exemplified with just the data mentioned, McCullagh (1980) JRSS B42, 109-142. I had it work very fine. -----Oprindelig meddelelse----- Fra: Peter Malewski <p.malewski at tu-bs.de> Til: Troels Ring <tring at

Dear all, I've replicated the cheese tasting example on p175 of GLM's by McCullagh and Nelder. This is a 4 treatment (rows) by 9 ordinal response (cols) table. Here's my simple code: #### cheese library(MASS) options(contrasts = c("contr.treatment", "contr.poly")) y = c(0,0, 1, 7, 8,8,19, 8,1, 6,9,12,11, 7,6, 1, 0,0, 1,1, 6, 8,23,7,

Hi All, could someone shed some light on what the difference between the estimated dispersion parameter that is supplied with the GLM function and the one that the 'gamma.dispersion( )' function in the MASS library gives? And is there consensus for which estimated value to use? It seems that the dispersion parameter that comes with the summary command for a GLM with a Gamma dist. is

What is the current state of the R-art for ordered probit models, and more esoterically is there any available R strategy for ordered cauchit models, i.e. ordered multinomial alternatives with a cauchy link function. MCMC is an option, obviously, but for a univariate latent variable model this seems to be overkill... standard mle methods should be preferable. (??) Googling reveals that spss

I am new to the R. I managed to download and install and got some basic summary work done. I need to do ORDERED LOGIT model. I searched web-site and mail archives and I think that I need to use "ordglm" function. However, I am not sure what package has it. I tried several (like RMUTIL, ORDINAL, REPEATED) etc. but no avail. Any help is useful. I still haven't subscribed to the

I am new to generalized linear models and studying McCullagh & Nelder (1989). Especially, I have a problem resembling the \"cheese taste\" example (5.3.1. p. 109) of the book. I tried to analyse the cheese example with R but failed to do so because R allowed me to use logit link function only with binary family that supposes 0 <= y <= 1. Do I need to scale the y\'s or

i want to tag categories to its menuname. i have a csv containing menu item name and in other csv i have a column containing some strings, i want to pick that strings from categories and look into menu items if any menu item containing that string i want to create a new column next to menu item name flagged as 1 otherwise 0 and the only condition is once a menu item flagged as 1 i don't need

Hi, I am wanting to look at frequent item sets using the arules package. I need to transform my data into a "transactions" object. The data I read in from a file has 2 columns, an ID and an item. How do I convert data like this into a transactions object? I've tried class? transactions but it only confuses me. My data is like this.... basketID item 1 bread 1 cheese 1 milk 2

I'm a Gentoo amd64 user. This is what I get when I try to run the Autopano (it installed just fine): Code: szczerb at nomad ~ $ wine .wine/drive_c/Program\ Files/Kolor/Autopano\ Pro/Autopano_Win32.exe fixme:actctx:parse_depend_manifests Could not find dependent assembly L"Microsoft.VC80.CRT" err:module:import_dll Library MSVCP80.dll (which is needed by

Hi Pedro, Not too hard. Just have to watch the order of the variables: ppdf<-read.table(text="N M W I 10 106 II 124 484 III 321 874 IV 777 1140 V 896 996 VI 1706 1250 VII 635 433 VIII 1437 654 IX 693 333 X 1343 624 XI 1221 611 XII 25 15 XIII 3 NA XIV 7 8", header=TRUE) plot(rev(ppdf[,2]),1:14,col="blue",lty=1,type="l",

R-help, I have a data frame: > mydata strata nh Nh Wh fh 1 I 10 26 0.048 0.385 2 II 32 84 0.154 0.381 3 III 16 42 0.077 0.381 4 IV 4 11 0.020 0.364 5 V 10 26 0.048 0.385 7 VII 64 168 0.309 0.381 8 VIII 49 129 0.237 0.380 9 IX 22 58 0.107 0.379 91 VI 0 0 0.000 0.000 and I wish to rearrange the data are sorted according to the roman

greetings, i have run into same problem as a friend who bought a logitech c310 webcam. he told me about them being on sale, so i decided to buy one also. the 2 webcam programs that i have installed from repos are 'camorama' and 'cheese'. both work well for single photos, but i can not get 'cheese' to make smooth, continuous video recording. what is recorded is a jerky

Hello, appended is a patch that makes it possible to use PAM both for password authentication and TIS (i.e. s/key or opie or any other interactive challenge/response scheme). I have developed this starting from the patch at http://www.debian.org/Bugs/db/61/61906.html on Debian with openssh-2.1.1p4-3. After configuring ssh with --with-pam-tis, there are two PAM services, "sshd" and

Ok, so I finally got the id3/mp3 end of things worked out for the next release of QTagger and I'm back to implementing tag reading for Ogg Vorbis. I've got some code that will read tags, but I haven't been able to figure out how to write them. Here's the code that I'm using. It's a fine C/C++ mix. The code will eventually be in a C++ app, so I don't mind

Hi all, I'm trying to combine exhaustively several character arrays in R like: x=c("one","two","three") y=c("yellow","blue","green") z=c("apple","cheese") in order to get concatenation of x[1] y[1] z[1] ("one yellow apple") x[1] y[1] z[2] ("one yellow cheese") x[1] y[2] z[1]("one blue

Hello all, I've been invited to do a presentation on Asterisk for the Ohio Linuxfest in Columbus this weekend (http://www.ohiolinux.org). Rough estimates are that nearly 500 people will be attending. I've been working on an outline for a couple of weeks and I would like to have some peer review of the information presented. I am going to have to cut down the content to make it fit in

Just opened cheese in Ubuntu 10.04 (Lucid) and the camera worked fantastically. I noticed the cam listed, UVC Camera(046d:0990) is exactly the same one in icu2 except that in cheese I saw /dev/video0 after the name of the camera. What does that mean? At all related to getting the same camera to work in icu2? New to this so just asking. Thanks. Bob

Hi Pedro, melt() is probably working. The problem is I did not finish the copy and paste.? It would have been better if I had included the ggplot() command. Try ============================================================== library(reshape2) library(ggplot2) dat1? <- structure(list(N = c("I", "II", "III", "IV", "V", "VI",