similar to: How to find out the trend during a time period?

The data is in the attachment. What I wanna get is: , , Sex = Male Eye Hair Brown Blue Hazel Green Black 32 11 10 3 Brown 53 50 25 15 Red 10 10 7 7 Blond 3 30 5 8 , , Sex = Female Eye Hair Brown Blue Hazel Green Black 36 9 5 2 Brown 66 34 29 14 Red 16 7 7 7 Blond 4

Hi all: Here's a question about result of loglinear analysis. There're 2 factors:area and nation.The raw data is in the attachment. I fit the saturated model of loglinear with the command: glm_sat<-glm(fre~area*nation, family=poisson, data=data_Analysis) After that,I extract the coefficients: result_sat<-summary(glm_sat) result_coe<-result_sat$coefficients I find that all the

Hi all: My data:data(sleep) If I wanna calculate each group's extra,what I can do is: #method1 attach(sleep) mean(extra[group==1]) mean(extra[group==1]) #method2 result<-matrix(,0,2) g<-split(sleep,sleep$group) for(i in 1:length(g)) { result<-rbind(result,data.frame(unique(g[[i]]$group),mean(g[[i]]$extra))) } colnames(result)<-c("name","mean") But the above

Hi all: I have a question about lm on t-test. data(sleep) I wanna perform t-test to test the difference between the 2 groups: I can use: t.test(extra~group) The t.test result shows that:t = -1.8608; mean1=0.75,mean2=2.33 But I still wanna use: summary(lm(extra~group)) Intercept=0.75,which is mean1,just the same as t.test. group2=1.58 means the difference of the 2 groups,so

Hi all: There's a question about glht function. My data:data_ori,which inclue CD4, GROUP,time. f_GROUP<-factor(data_ori$GROUP) f_GROUP is a factor of 3 levels(0,1,2,3) result <- lme(sqrt(CD4) ~ f_GROUP*time ,random = ~time|ID,data=data_ori) glht(result, linfct = mcp(f_GROUP="Tukey") ) Error in `[.data.frame`(mf, nhypo[checknm]) : undefined columns selected I can't

Hi all: If I have a dataframe of N columns.If I wanna get the min(or max,or mean...etc)of the whole dataframe,how to do it quickly? What I can do is only: min(data[,1:ncol(data)]) Any other suggestion? Thanks! best [[alternative HTML version deleted]]

Hi all: I have a question on sample size calculation of 2 groups of data. If 2 groups of data are all normal distribution, then I can use the function "n.indep.t.test.eq" from samplesize package.But if 2 groups of data are all skewed distribution, but not normal distribution,how can I calculate the sample size then? I've tried many transformation (e.g. log arcsin…) in order to

Hi all: I have a question about linear mixed model. my linear mixed model with randomized slope and intercept with interaction of time and group(g1,g2,g3): model<- glmmPQL(log10(CD4) ~ time + factor(group)+ time:factor(group), random = ~time|id) What I get is only the main and interaction of time and group.My question is: 1. How can I get the g1,g2,g3's slope respectively?In other

time1 = 2008-03-09 time2 = 2010-9-10 How to compute how many years between time1 and time2? Thanks! best [[alternative HTML version deleted]]

Hi all,I have some questions about the covariants of regression. My target: To explore the trend of CD4 level through a period of time. Response variable: CD4 count Explanatory variable:time Also, the demology information is available,such as gender,occupation,income level... Q1,Are these variables of demology information called covariant? Q2,How can I correct the impact of

Hi,all: I met a problem of nls. My data: x y 60 0.8 80 6.5 100 20.5 120 45.9 I want to fit exp curve of data. My code: > nls(y ~ exp(a + b*x)+d,start=list(a=0,b=0,d=1)) Error in nlsModel(formula, mf, start, wts) : singular gradient matrix at initial parameter estimates I can't find out the reason for the error. Any suggesions are welcome. Many thanks. [[alternative HTML

Hi all: Which function or package can perform "competitive risk model"? Thanks a lot! My best [[alternative HTML version deleted]]

If I have 3 groups,and for each group,I get the ratio(e.g. incidence rate). Now I wanna compare 3 ratio pairwise,and get the corresponding p values,i.e: group1 vs group2 ,p value=? group1 vs group3 ,p value=? group2 vs group3 ,p value=? Which statistical test should be used? Thanks a lot for your help. [[alternative HTML version deleted]]

Hi all: I have a question about the sample size calculation. It's a pilot study,which includes 2 groups(low,high),3 time point(3,6,9 monthes).Each person has 3 results which according to the 3 time points.So it's a longitudinal study. I want to calculate the minimum sample size from the pilot study, but can't find the solution since the data is highly skewed and it's a

Dear all, I'm using R 2.6.0 Windows XP Access 2003 Lao Script for Windows 7.02 I'm trying to use Lao with R for 1) Chart Legend Labels 2) Writing out text for HTML pages With Lao Script for Windows it is possible to type Unicode (16-bit). The strings I want to use are stored in an Access database, queried trough RODBC. The result I'm getting back are only question marks. I tried

If there are many character variables,and I want to get the mosaic plot of every pair of each variable,how to do then? If the variables are numeric, I can use pairs to get paired scatter plot. But as to the character variables, how to get the "paired mosaic plot"? Many thanks. -- QQ: 1733768559 At 2015-02-07 17:04:26,"Jim Lemon" <drjimlemon at gmail.com>

Hi List, I have weekly sales observations for several products drawn via ODBC. Source data is available at https://www.dropbox.com/s/78vxae5ic8tnutf/asr.csv. This is retail sales data, so will contain seasonality and trend information. I expect to see 52 or 53 observations per year, each observation occuring on the same day of the week (Saturday). Ultimately I'm looking to feed these series

more info just 1 server webserver + appserver + dbserver On 9/14/07, Yan Meng <dreamwords at gmail.com> wrote: > > Hi, Guys > > Our team is running a web2.0 finance site in China > > http://www.caibangzi.com/ > > We use APACHE + mod_proxy_balancer + mongrel to deploy our application, it > workd super in the past few monthes. > > However, just from last

The X86 backend has this code: X86TargetMachine::X86TargetMachine(const Module &M, IntrinsicLowering *IL) : .... FrameInfo(TargetFrameInfo::StackGrowsDown, 8/*16 for SSE*/, 4), That is, it uses "4" as local area offset. Based on prior discussion this should mean that the local area starts and address ESP+4. Is this really true? On X86 stack grows down, so

Hi, all I would like to use R to perform k-means clustering on my data which included 33 samples measured with ~1000 variables. I have already used kmeans package for this analysis, and showed that there are 4 clusters in my data. However, it's really difficult to plot this cluster in 2-D format since the "huge" number of variables. One possible way is to project the