similar to: setCoefTemplate

Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Rejected)~Gender,data=berkeley,family="binomial") berk2 <-

Hi everyone, Does anyone know how to get "black and white theme" (grey scale,, I would say) graphs using lattice or ggplot2, as it is shown in this webpage: http://lmdvr.r-forge.r-project.org/figures/figures.html? I am using Sweave, and I cannot get that color configuration. I have added the following option: trellis.device(color=FALSE) but I got a pdf file with color graphs. Thank

Hi everyone, I have just got different samples from a dataframe (independent and exclusive, there aren't common elements among them). I want to create a variable that indicate the sampling selection of the elements in the original dataframe (for example, 0 = no selected, 1= sample 1, 2=sample 2, etc.). I have tried to do it with ifelse command, but the problem is that the second line

Hello everyone, I am new in R and I am trying to create a id number (a correlative sequence of numbers) by group, and a correlative sequence of numbers inside each group (my idea is to get statistics by group without having to aggregate the database). Here an example: group id_group correlative_group A 1 1 A 1 2 A 1 3 A 1 4 B 2 1 B 2 2 B 2 3 C 3 1 C 3 2 C 3 3 C 3 4 C 3

Hi everyone, I have a dataset of friendship with this format: ego alter 4746 1 2 9742 1 3 14738 1 NA 4747 2 NA 9743 2 3 14739 2 1 4748 3 13 9744 3 5 14740 3 14 4749 4 NA 9745 4 NA 14741 4 NA 4750 5 NA 9746 5 13 14742 5 10 4751 6 12 9747 6 7 ... NA means that individuals don't select any

Hi everyone, Does anyone know if there is a package to compute raking weights using R? What I need is to create a variable with weights base in some demographic variables (e.g. sexo, age group, area) using the raking procedure. Thank you in advance! -- Sebasti?n Daza

Hi everyone, I need to get a between-component variance (e.g. random effects Anova), but using lmer I don't get the same results (variance component) than using random effects Anova. I am using a database of students, clustered on schools (there is not the same number of students by school). According to the ICC1 command, the interclass correlation is .44 > ICC1(anova1) [1] 0.4414491

Hi everyone, I am just trying to figure out how to do a xyplot where in addition to dots and lines I can change dots' colors according to an individual variable (e.g., marital disruption across time, a dummy 0/1). When I use "groups" specification (see below), I get two different lines for each individual based on groups, and what I want is to get one line connecting dots, and

Hello, I have got a cluster sample using an election dataset where I already had the final results of a county-specific election. I am trying to figure out what would be the best sampling design for my data. The structure of the dataset is: 1) polling station (in general schools where people vote, for a county, for example, there are 15 polling stations) 2) inside each polling station, there

Hi guys, I have been setting up a samba server to share some directorys around my office network. Now I have been successful in making these shared directorys so that everyone can access them, but what I really need to do is to make the shares to have limited accessibility. What I need is to have 1 folder that is completely Private to everyone else but 1 user and everytime I have tried this

<mailto:r-help@r-project.org>Hi every one, I want to calculate the combination function in R, the value not all the possible choices. I mean cmbn(5,2)=10. Is there any function unless using factorial? Regards, Amir [[alternative HTML version deleted]]

Dear Mrs Rachel Pearce, I am looking for a function "cutpt-coxph" in R - like you did some years ago. How have you solved the problem? Have you found it or a similar function? thank you, Sincerely, Friederike "The title says it all really; I am looking for a function along the lines of cutpt.coxph as described in "Survival Analysis Using S" (Tableman and Kim), Chapter

Dear R users, I would like to perform Latent Class Logit Models for the analysis of choice experiments in environmental valuation. This kind of analysis is usually performed with NLogit Software (http://www.limdep.com). I attach the results I usually obtain using NLogit and NLogit model specifications. For Random parameter models and Logit Models I usually perform my analysis with the package

i have a time Series of IBM closing px from 1/1/2000 to today I want to graph the time serie by dividing the graph by year and month all the monthly graphs with the same year will go to one page. so from 1/1/2000 to 11/19/2010. i will have 11 pages, and each page will have 12 graphs (jan to dec) except for 2010. I am able to do it in R, but when i use sweave, I can only print the last page.

Dear R Users, R Core Team, I have a two dimensional space where I measure a numerical value in two situations at different points. I have measured the change and I would like to test if there are areas in this 2D-space where there is a different amount of change (no change, increase, decrease). I don´t know if it´s better to analyse the data just with the coordinates or if its better to group them

Hello, I have a dataframe (Lx) with 5 Lb, and 5 Lw variables. I want to create several variables according to the loop presented below (data attached). Lx <- read.csv("Lx.csv", header=T, sep=",") for (i in 1:20) { Lx$sb1[i] <- Lx$Lb1[i+1]/Lx$Lb1[i] Lx$sb2[i] <- Lx$Lb2[i+1]/Lx$Lb2[i] Lx$sb3[i] <- Lx$Lb3[i+1]/Lx$Lb3[i] Lx$sb4[i] <-

Hi everyone! Does anyone know how to install the rj package in Windows (7)? I am trying to set up eclipse, but I got [INFO] The R package 'rj' is not available, R-StatET tools cannot be initialized. Thank you in advance!

Hi there, I am having trouble subsetting a data frame by a conditional via one column (of many). I read the file into R through "read.fwf," where I specified column widths. Original data is .DAT. I then utilized "names" function to read in column headings. For one column, PRVDR_NUM, I wish to further amend the entire data set, but only have PRVDR_NUM == 050108. This is

Hello, I have setup a cluster system with drbd and heartbeat. I have samba running as a PDC. Everything works great when I am on server1 But when failover ocuurs I cannot seem to log on the domain anymore. I Have symlinks to my home share : /var/cache/samba , /etc/samba/ who are synced between 2 servers. Do I have to sync another files ?? greets John

?? > z <- array(1:24,dim=2:4) > all.equal(f(z,3),f2(z,3)) [1] "Attributes: < Component ?dim?: Mean relative difference: 0.4444444 >" [2] "Mean relative difference: 0.6109091" In fact, > dim(f(z,3)) [1] 2 3 4 > dim(f2(z,3)) [1] 3 4 2 Have I made some sort of stupid error here? Or have I misunderstood what was wanted? Cheers, Bert Bert Gunter