similar to: Can you create a MySQL database with RMySQL?

Hi all: My data has 3 variables: age(3levels : <30y=1 30-50y=2, >50y=3) gender(Male=0, Female=1) CD4 cell count(raw lab measurement) y(1:death 0:alive) I perform logistic regression to find out the factors that influence y. result<-glm(y ~ factor(age) + factor(gender) + CD4,family = binomial) >From the result,I can get OR(Odds Ratio) of gender via exp(Estimate of Female,

Dear R experts, I would really appreciate if you had an idea on how to use more efficiently the aggregate method: More specifically, I would like to calculate the mean of certain values on a data frame,? grouped by various attributes, and then create a new column in the data frame that will have the corresponding mean for every row. I attach part of my code: matchMean <-

Consider: > x <- matrix(1:6, 2,3) > rownames(x) <- c("ID1", "ID2") > colnames(x) <- c("Attr1", "Attr2", "Attr3") > x Attr1 Attr2 Attr3 ID1 1 3 5 ID2 2 4 6 > write.csv(x,file="x.csv") "","Attr1","Attr2","Attr3" "ID1",1,3,5

I would like to selectively remove rows from a table. I had hoped that I could create a table and selectively add rows with something like > NewTable<-table(nrow=100, ncol=4) > NewTable[1,]<-OldTable[10,] but that doesn't work. The former call gives > NewTable ncol nrow 4 100 1 while the latter call gives a table the length of OldTable. Making a matrix, m, with the

Hi, Suppose I have the following data name     score Abel        88 Baker      54 Charlie    77 stored a  table called myData. I want to write a function that will create a table which is a subset of myData containing those have a score > 75. I know I can do this with the following command: subset(myData, score > 75) But I would like to do this via a function, something like:

Hi, I am running R2.8.1 under Linux, and I am having trouble using the variance functions in nlme My basic model was something like: model0 <- lme( log(growth) ~ light * species.group , data=data, random=~light|species ) # with 20 odd species divided in 2 groups Following the methods in Pinheiro&Bates I tried to put a variance function in the model: model1 <- update(model0,

In R-devel (svn 49628) and back to at least R 2.7.0 we get inconsistent results when attempting to attach attributes to a variable with the value NULL. If we use attributes<- it finishes but changes the value to list(). > a<-NULL > attributes(a)<-list(attr1="First attribute", attr2=2+2i) > a list() attr(,"attr1") [1] "First

Hello, I have a problem concerning logistic regressions. When I add a quadratic term to my linear model, I cannot draw the line through my scatterplot anymore, which is no problem without the quadratic term. In this example my binary response variable is "incidence", the explanatory variable is "sun": > model0<-glm(incidence~1,binomial) >

Hey there I am estimating a multilevel model using lmer. I have 5 imputed datasets so I am using mitools to pool the estimates from the 5 > > datasets. Everything seems to work until I try to use > MIcombine to produced pooled estimates. Does anyone have any suggestions? The betas and the standard errors were extracted with no problem so everything seems to work smoothly up until

Hi all, This should be easy, but I can't seem to figure it out. I have a table like this named newtable a1 a2 a3 a4 Cnts Score 1 1 0 0 4 3.28 1 0 1 1 2 2.63 I want the following: a1 a2 a3 a4 Cnts Score 1 1 0 0 4 3.28 1 1 0 0 4 3.28 1 1 0 0 4 3.28 1 1 0 0

Hi, i am unable to install the rugarch package. More than that i do not even find this package in my list of possible packages. Its possible than the name has changed, or the package is not longer availiable? Is there a similar package avaliable for garch modelling except the fGarch what i am using now? many Thanks Roland -- View this message in context:

Hi, I have 100000 units with 10 attributes (attr1, attr2, attr3, etc...) For instance: unit attr1 attr2 attr3 ... 1 a ww 12 2 a re 11 3 b ww 09 4 c yt 02 5 a qw 02 ... I'd like to answer to the question: a) what are the most frequent combinations of attributes? b) How could I describe the relations

Hi, I''ve started looking at ticket 90<http://merb.devjavu.com/projects/merb/ticket/90>and there are a couple of situations that the current implementation doesn''t cater for. 1. Tag Attributes are wiped out. - <tag1 attr1=''1''>Content</tag1> #=> { "tag1" => "Content" } - I think it should return -

Ok, So, I''ve got two data tables associated through a foreign key. class Child < ActiveRecord::Base belongs_to :parent, :class => Parent, :foreign_key => "parent_id" # also has attr1, attr2 end class Parent < ActiveRecord::Base has_many :children, :class => Child, :foreign_key => "parent_id" end So, I''m trying to build the new/edit

Ok, so I'm slowly figuring out what a factor is, and was able to follow the related thread about finding a mode by using constructs like my_mode = as.numeric(names(table(x))[which.max(table(x))]) Now, suppose I want to keep looking for other modes? For example, Rgames> sample(seq(1,10),50,replace=TRUE)->bag Rgames> bag [1] 2 8 8 10 7 3 2 9 8 3 8 9 6 6 10 10 7 1

Here's my contribution to R. When R interacts with external programs (MySQL, cURL, etc.), it often requires a pasted string that is sent to these programs. For readability reasons, it is often preferable to have complex commands (SQL for example) spread on several lines. However, the normal paste function requires to add additional ' ", ' at the end of each line and another '

Hello all, I''ve got a site that will be used by people (hopefully) all over the world. One of the requirements of registration is country and province of residence. I have a list of all the countries in the world and all their subnational entities. Currently, the two lists are divided between a pair of comma-delimited text files. If this were something less dynamic and cool, I''d

Thanks for your feedback David.... Yes, the intention is to support duck-typing. In my code example, a.attrB1 = 5 .... is probably more accurately stated as ..... a. attr1.attrB1 = 5 Certainly this would be an unsafe call because we would not know that attrB1 is an attribute of attr1 until runtime. I believe this type of call is valid in a language like Python. Are there any

When I use rbind() or rbind.data.frame() to add a row to an existing dataframe, it appears that attributes for the column of type "factor" are dropped. See the sample example below to reproduce the problem. Please suggest How I can fix this. Thanks, Sammy a=c("Male", "Male", "Female", "Male") b=c(1,2,3,4) c=c("great", "bad",

Dear all, I am still fairly new to R and try to analyze large tables of data generated from genomic experiment. Currently, I am trying to plot pair of experiments coming from different file, trying to look at the behavior of individual feature in pair of experiment. My problem is that I have independent list from different source and I would like to plot the pair of value using a common key. As