similar to: context problem

Hi, Can anyone tell me why I am not getting the correct intervals for fixed effect terms for the following generalized linear mixed model from HPDinterval: > sessionInfo() R version 2.4.1 (2006-12-18) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

Hello: Forgive me, this is surely a simple question but I can't figure it out, having consulted the help archives and "Data Manipulation With R" (Spector). I have a list of 11 data frames with one common variable in each (prov). I'd like to use lapply to go through and recode one particular level of that common variable. I can get the recode to work, but it only returns the

I believe that I have incorporated all of the feedback given (thanks!). Changes since the 2008-11-16 version include: - ksh style & coding standards compliance in test script (Roland) - "dof_init_debug == B_FALSE" vs. "!dof_init_debug" (Adam) The updated webrev is at: http://cr.opensolaris.org/~mgerdts/6750659-2008-12-03/ -- Mike Gerdts http://mgerdts.blogspot.com/

I've posted the following to R-help on May 15. It has reproducible R code for real data -- and a real (academic, i.e unpaid) consultion background. I'd be glad for some insight here, mainly not for myself. In the mean time, we've learned that it is to be expected for anova(*, "marginal") to be contrast dependent, but still are glad for advice if you have experience. Thank

Hi All, I *think* it's ddply because the function recognizes vr1, etc, in other parts of the function. Here's some code: # create dataset PROV.PM.FBCTS <- c(0.00 ,0.00, 33205.19, 25994.56, 23351.37, 26959.56 ,27632.58, 26076.24, 0.00, 0.00 , 6741.42, 18665.09 ,18129.59 ,21468.39 ,21294.60 ,22764.82, 26076.73) FBCTS.INV.TOT <- c(0 , 0, 958612, 487990, 413344, 573347,

Hi, I am getting a warning message when I am fitting a generalized linear mixed model (m1.2 below). CHOLMOD warning: matrix not positive definite Error in objective(.par, ...) : Cholmod error `matrix not positive definite' at file:../Supernodal/t_cholmod_super_numeric.c, line 614 Any idea? Thanks for your help, Reza > sessionInfo() R version 2.4.1 (2006-12-18) i386-pc-mingw32

Hi, How are you? I have a quick question.... I have code that works perfectly with R desktop, but does not work with RWeb. Could you please tell me how to modify the code below so it will work with RWeb? #Read in txt file happyguys<-read.table("c:/test8.txt", header=TRUE, row.names=1) #Subset the txt file to only include certain values test<-subset(happyguys, GRADE == 7

Dear all, I am interested in plotting maps visualizing spatial statistics in an aggregated fashion, according to administrative boundaries. More specifically, I have fitted a cross-section model on data regarding Italian "counties" (province, for Italian readers) and I would like to visualize residual behavior on a map, in order to have a first assessment of their spatial

Hi, i am trying to read a .txt file, do a couple of select if statements on my data, and then finally use the ?table function to get frequency counts on the data. Specifically, i am looking at answering the following question: What is the frequency of Grade 7 students in the province of Alberta who are smokers? I am having some problems: 1)i cannot get the column names to show up when print

Dear list Is there a way to extract a random sample without duplicated row from a dataframe ?. a=c(1,2,3,1,1,1,2,1) b=c(1,2,3,1,2,1,2,1) c=c(1,1,1,1,1,1,1,1) d=c(1,2,3,1,1,1,2,1) prov<-data.frame(a,b,c,d) prov2<-prov[sample(1:nrow(prov),5,replace=T),] prov2 a b c d 3 3 3 1 3 6 1 1 1 1 3.1 3 3 1 3 5 1 2 1 1 8 1 1 1 1 I tryed the above code but as you can see sample function

Hi, I have the file below called happyguys. It is a subset of data. How do I apply the weight variable (WTPP) to this file? Can i just multiply each column (except the first column because it is a record id) by WTPP? If the answer is yes, how do I multiply one variable name by another? Thanks, Nat PROV REGION GRADE Y_Q10A WTPP 83 48 4 7 2 342233324020 115

Hi there I'm trying to fit a logistic regression model to data that looks very similar to the data in the sample below. I don't understand why I'm getting this error; none of the data are proportional and the weights are numeric values. Should I be concerned about the warning about non-integer successes in my binomial glm? If I should be, how do I go about addressing it? I'm

I've been noticing that my group of Polycom IP 501 phones seems to randomly reset themselves nearly every night (I guess it usually happens at night, since I've never seen it happen while I've been at work during the day).. When I say "reset", I mean, the hands free volume and ring volume are set to the default and the call logs (received calls, missed calls, placed calls)

Hi, I got this working. i had to get rid of an underscore in one of my variable names to work on RWeb. Thanks for your help though!! Natalie __________________ I think you have to be either honest or careful, since the code you submitted to Rweb is different with your former code! Why the condition was changed from "test<-subset(X, GRADE == 7 & Y_Q10A < 9)" to

Estimados compañeros tengo el siguiente problema: quiero poner etiqueta a los valores de una variable. He creado una lista que se llama geren que tiene 194 valores distintos. Creo la siguiente data.frame: dg<-data.frame(c_autónoma=names(tca)[sb7$ca],prov=names(cprov)[sb7$prov],geren=names(geren)[sb7$gerencia]) sb7$gerencia tiene 107 códigos distintos incluidos dentro de los códigos de

Hi everyone, I am putting together an application in which each part of it was tested individually. Now that it is put together, I''m having problems... The one giving me problems is my "ImportController" which parses a CSV during post and puts it into my database. I have two separate CSV upload actions (one is for a bill called "bes" and the other from the provider

Hello, We have a strange situation (asterisk 1.6.2.14), where we get a result for DIALSTATUS for BUSY and No-ANSWER, but nothing for CANCEL. This is the (relevant) test dialplan: -------------------------------- [incoming-private] exten => _X., n, Dial(SIP/1001,30) exten => _X., n, NoOp(${DIALSTATUS}) exten => _X., n, Gosub(incoming-status,s-${DIALSTATUS},1) [incoming-status] exten

Hi is there a way in R to create a new column vector with the numbers of 2 others columns simply appended? Example : I have a column with provinces codes (1 to 19) I have a second column with districts codes (1 to ##, depending on the province) I want to create a third column with a code that would be unique for each district: Province  District   District_Unique 1             1            11

Hi, I have a data file with columns separated by ";" I read this file without any problem using read.csv2( ) but I had problems trying to read it with read.table( ... sep=";"). So it is not a problem for me, but I wonder if there is a bug here. drt <- read.csv2("t.txt", header=TRUE) # ok dcs <- read.table("t.txt", header=TRUE,

Dear CentOS hive mind, I'm trying to package up a perl module into an RPM for easy deployment. I want it to be as self-contained as possible (to avoid version issues with perl modules in base or EPEL). So in my spec file, I'm doing: curl -L http://cpanmin.us | perl - App::cpanminus -L %{buildroot}/opt/zonemaster Zonemaster This way, cpanminus is installed first, and then it goes on to