search for: eugen_pircalabelu

Hello R-users, Checking the archives, I recently came across this topic: "export tables to Excel files" (http://r.789695.n4.nabble.com/export-tables-to-Excel-files-td1565679.html#a1565679), and the following interesting references have been proposed: http://rwiki.sciviews.org/doku.php?id=tips:data-io:ms_windows

Hi R users! I am trying to learn some multilevel analysis, but unfortunately i am now very confused. The reason: http://www.ats.ucla.edu/stat/hlm/seminars/hlm_mlm/mlm_hlm_seminar.htm http://www.ats.ucla.edu/stat/sas/seminars/sas_mlm/mlm_sas_seminar.htm and MlmSoftRev. pdf from mlmRev package. >From what i see, the first two links seem to declare the level one variable as a random part (i

Hi R users I want to use the lme package for a multilevel analysis on the following example: > math<-c(2, 3,2, 5, 6 ,7 , 7) > sex<-c(1, 2, 1, 2, 2, 2, 1) > school_A<-c(1,1,1,2,2,2,2) > school_B<-c(10,10,10,20,20,20,20) > mydata<-data.frame(math, sex, school_A, school_B) > mydata School_A and school_B are two different school characteristics, math is an

Hi list, Using the lmList function from ?nlme package? I get the following error message: ?Error in !unlist(lapply(sum.lst, is.null)) : invalid argument type? # this is the syntax used # fm_cream <-lmList(MULTDV~TIME|cod, data=zz, na.action=na.omit) # summary(fm_cream) When I generate a simple example the function works, but on the real data I get that message. I imagine that it is a

Hi List, Can anyone tell me how could i put the "BACK" button in the following code, just under the "AAA" menu? I want this button to go back to the previous page, and since it has nothing to do with the "1" and "2" buttons, i want it somehow separated from these two buttons, but i don't know how. I searched the web for some examples but my results

Hi List, I have the following problem: I am using the multilevel package and make.univ function for available in the package and then xyplot from lattice and I want to know how could I be able to use the ?coefficient? for the straight line that passes the data ? Example from help: library(multilevel) data(univbct) #a data set already in univariate or stacked form for job satisfaction

Good morning! I have 4 questions which trouble me: 1. I want to test the hypothesis that the 2 proportions (the mean of a binomial) which come from 2 different samples are equal. I want to use the following function z= (p1-p2)/ sqrt((p1(1-p1)/n1)+(p2(1-p2)/n2)) which is one of the standard formulas for this case. Is there such a function in R? p1=the proportion from the first sample n1=the

Hi R users! I have the following problem: how appropriate is my aov model under the violation of anova assumptions? Example: a<-c(1,1,1,1,1,1,1,1,1,1,2,2,2,3,3,3,3,3,3,3) b<-c(101,1010,200,300,400, 202, 121, 234, 55,555,66,76,88,34,239, 30, 40, 50,50,60) z<-data.frame(a, b) fligner.test(z$b, factor(z$a)) aov(z$b~factor(z$a))->ll TukeyHSD(ll) Now from the aov i found that my model

Hi List, I have a problem with using Latex and Sweave for creating a document. So I downloaded the Sweave manual from http://www.statistik.lmu.de/~leisch/Sweave/Sweave-manual.pdf and i have tried to replicate the example on pages 4-5, but i encounter the following problem: IT DOES NOT WORK. I can run Sweave, create the .tex file and then when building the .tex file in TeXnikCenter i get all

Hi R users! I've come across using kpss tests for time series analysis and i have a question that troubles me since i don't have much experience with time series and the mathematical part underlining it. x<-c(253, 252, 275, 275, 272, 254, 272, 252, 249, 300, 244, 258, 255, 285, 301, 278, 279, 304, 275, 276, 313, 292, 302, 322, 281, 298, 305, 295, 286, 327, 286, 270, 289, 293, 287,

...an help please. On use of integration function integrate9; how do I get the function to return just the value with no messages please Glenn [[alternative HTML version deleted]] ------------------------------ Message: 39 Date: Mon, 22 Dec 2008 14:57:45 -0800 (PST) From: eugen pircalabelu <eugen_pircalabelu at yahoo.com> Subject: [R] newbie question on tcltk To: R-help <r-help at stat.math.ethz.ch> Message-ID: <620247.52862.qm at web38608.mail.mud.yahoo.com> Content-Type: text/plain; charset=us-ascii Hi List, Can anyone tell me how could i put the "BACK" button in the follow...

Good afternoon! I'm trying to use the survey package to get a better point of view for my data, but i need some piece of advice: i have some data from a survey which has been stratified using 2 criteria: region(7 values), size of locality(5 values) Using the survey pakage how can i define in a correct way this design (taking into account all 4 strata not just one as in the Survey

Good afternoon! My question is more of a "is there a package for doing........?" If one has 2 independent studies which both independently measure a trait of a population, but both of them are constructed from 2 separate frames (which TOGETHER refer to the population of interest), is there a possibility to infer a mean, median.. for the population but taking into account all the

Good afternoon! I have the following data: a<-data.frame (id_hh=c(1:5), strata=c(1,1,2,2,1), Nhstrata=c(100,100,200,200,100), Nrmemb=c(2,4,2,5,4)) a$ocmemb1<-c("wk","jl","st","jl","st") a$ocmemb2<-c("wk","jl","st","wk","wk") where id_hh is a code of identification for the household (my

Good afternoon! I'm trying to learn time series but i have a bit of of a problem using R packages for this. 1. > LakeHuron > sample(500:600, 98) > sample(500:600, 98)->t > fit<-arima(LakeHuron, order=c(2,1,1), xreg=t) > fit > predict(fit, n.ahead=1, newxreg=t) Now, my problem is this: is it ok to use the same t in predict function or should my newxreg contain 99

Hi R-users! Can anyone tell me what this error message refers to exactly > a<-c('m','f','m') > b<-c(1:3) > st<-c(1:3) > data.frame(a,b,st) a b st 1 m 1 1 2 f 2 2 3 m 3 3 > data.frame(a,b,st)->l > p<-svydesign(ids=~1, data=l) > pop<-c('(intercept)'=100, st1=10,st2=20) > calibrate(p,~st, pop) Error in

Hi R-users! I have the following example: a<-data.frame(cat=c(5,10,15), dog=c(5,10, 15), mouse=c(10,10,20)) b<-data.frame(cat=c(15,10,5), dog=c(15, 10, 5), mouse=c(20,10,10)) rownames(b)<-c("scared", "happy", "sad") rownames(a)<-c("scared", "happy", "sad") Let's say that a and b are 2 contingency tables and 5,10, 15, 20

Hi R users! Is there a function that extracts the simultaneous confidence intervals for multinomial proportions as described by Sison and Glaz 1995? or anyone else for that matter? I have seen that SAS has macro for this ( http://www.jstatsoft.org/v05/i06) and i was wondering if R had something similar? Thank you and have a great day ahead! ---------------------------------

Good evening R-users! I have the following problem: i want to get a weighted crosstable along with the adjusted standardized residuals test Example: a<-c(1,1,1,1,1,2,2,2,3,3) b<-c(10,10,23,24,33,45,46,70,21,66) c<-c(3,3,2,3,4,1,1,1,3,3) d<-c("a", "b","b","c","a","a","a", "b",

Hi List, I have the following time series and i want to be able to plot it while having the x-axis running from 1998 to 2006 but in a bi-annual format So here is my example: x<-sample (10:100, 17) x<-ts(x, start=c(1998,1), frequency=2) plot(x) When I plot the ts it goes from 1998 to 2006 by 2years, when in fact i want it to go from 1998 to 2006 with the following form:1998 may, 1998