similar to: plotting time series

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

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

Sir, I found your description of the dataset about nodal involvement in prostate cancer. It comes from the book biostatistics casebook. I like to use the dataset for doing logistics regression. Can you tell me where I can find the dataset. Thanks and greetings Wim van Baarle wvbaarle@wxs.nl [[alternative HTML version deleted]]

Dear, Thanks for the previous tips about 'its' for importing the following data. 5/10/1998,7 5/11/1998,5 5/12/1998,2 5/14/1998,1 5/15/1998,1 5/19/1998,1 5/20/1998,1 1. When using the following command; test<-readcsvIts('Fires98.csv',informat=its.format("%m/%d/%Y"),header=FA LSE) the function reads in the data from the csv file as; V2 05/10/1998 7

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 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, 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

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, 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

Good afternoon! I'm trying to model a time series on the following data, which represent a monthly consumption of juices: >x<-scan() 1: 2859 3613 3930 5193 4523 3226 4280 3436 3235 3379 3517 6022 13: 4465 4604 5441 6575 6092 6607 6390 6150 6488 5912 6228 10196 25: 7612 7270 8617 9535 8449 8520 9148 8077 7824 7991 7660 12130 37: 9135 9512 9631 12642

Good afternoon! I'm trying to model a ts but unfortunately i'.m very new to this kind of modeling so i 'll be very grateful if you have an advice. This was my syntax:

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,

Good afternoon! I'm trying to use the Survey package for a stratified sample which has 4 criteria on which the stratification is based. I would like to get the corrected weights and for every element i get a weight of 1 E.g: tipping design <- svydesign (id=~1, strata= ~regiune + size_loc + age_rec_hhh + size_hh, data= tabel) and then weights(design) gives

Hello Everyone I want to plot Moasic Plots, I have tried them using iplots package (using imosaic). The problem is the names dont get alligned properly, is there a way to a align the names and provide legend in Mosaic plots using R? Also I would like to know any other packages using which I can plot Mosaic Plots Thank you in advance Sunita -- View this message in context:

Hello R-users, My toy example: aa<-c(1:5) bb<-c(NA,2,NA,4,5) cc<-c(1,2,NA,4,NA) dd<-c("A","B","B","A","C") df<-data.frame(aa,bb,cc,dd=as.factor(dd)) table(unlist(df[,1:3])) Can anyone point me to what function let's me do a crosstabulation between table(unlist(df[,1:3])) and df$dd? I want to find out when dd==A (or B,

Dear all, I've been struggling to perform common operations on dates that I need to be able to correct draw date-time scales - in particular I need to be able to round/truncate/ceiling dates to arbitrary precision - e.g. to weeks, months or years (or multiples thereof). I haven't been able to find anything to do this in base R (trunc.Date only truncates to sub-day units), or in the date

Hi, I'd like to make a time series at an annual frequency. > a<-xts(x=c(2,4,5), order.by=c("1991","1992","1993")) Error in xts(x = c(2, 4, 5), order.by = c("1991", "1992", "1993")) : order.by requires an appropriate time-based object > a<-xts(x=c(2,4,5), order.by=1991:1993) Error in xts(x = c(2, 4, 5), order.by =

  Hello R-users,   I have the following question, for which my search did not really return any usable result. If I have a matrix a1, and a vector a2 like below   a1<-matrix(c(1:4),2,2) a2<-c(8,9)   is there any function like the expand.grid (or some clever calling of the function) such that it outputs a matrix or dataframe where the entire a1 matrix is repeated for each value of a2 (the