search for: firstob

...u supply a small dataset that cause the error you are talking about? One way to do this is to use supply the results of dput(tbdata) if it is small enough. If not, maybe create a subset of the data and then use dput() Did you get the problem when you tried with the tutorial Bert suggested? > firstobs <- subset(tbdata, capture.ID == 'B1-1108') > firstobs <- firstobs[,c(1,2,2,3,4)] > colnames(firstobs) <- c('id', 'start.date','end.date','lat','long') > > ###change date format and change start date to previous 14 days > first...

##MODISTools example library(MODISTools) library(lubridate) setwd('~/Documents/Modis data') #####MODISTools with buffalo data ###Read in data rename for easier coding tbdata <- read.csv('~/Desktop/All TB data for EVI, NDVI.csv') firstobs <- subset(tbdata, capture.ID == 'B1-1108') firstobs <- firstobs[,c(1,2,2,3,4)] colnames(firstobs) <- c('id', 'start.date','end.date','lat','long') ###change date format and change start date to previous 14 days firstobs$start.date <- dmy...

1. You should always cc the list unless there is a clear reason not to. 2. You still have failed to follow the posting guide: You say you have difficulty troubleshooting your code, but you have shown us no code. You got an error message that seems explicit, but with neither code nor data, I do not know whether anyone can make sense of it. In any case, I certainly cannot. Cheers, Bert Bert

Please read up [1][2][3] on what constitutes reproducibility. A sample of data that triggers the problem is essential. [1] http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example [2] http://adv-r.had.co.nz/Reproducibility.html [3] https://cran.r-project.org/web/packages/reprex/index.html -- Sent from my phone. Please excuse my brevity. On June 23, 2017 6:54:48

I have and they have not yet replied - however that was only two or three days ago. I have included a code example. > On Jun 22, 2017, at 8:25 PM, David Winsemius <dwinsemius at comcast.net> wrote: > >> >> On Jun 22, 2017, at 4:44 PM, Caroline <gliddeca at science.oregonstate.edu> wrote: >> >> I am using the R-package MODISTools (different than

..., Linear ","\n") cat ("\n") whiten <- 1 kernel <- 3 band <- 4 out <- fm(y,x,1,1) cat ("\n") cat ("\n") if (fplot == 1){ mtit <- "Ependiture, 1986:2011" xtit <- "Figure 1" ytit <- "F stat" t <- nrow(y) firstob <- 1986 + round(.15*t)/4 inc <- .25 rf <- nrow(out$f) crit <- matrix(1,rf,1)%*%cbind(15.2,6.22,7.815) ytics <- range(cbind(out$f,crit)) xx <- as.matrix(seq(firstob,firstob+(rf-1)*inc,inc)) plot(xx,out$f,type="l",ylim=ytics,lty=1,col=1,ann=0) lines(xx,crit[,1],lty=2,col=2...

...:00|7.3||3||||||||||||| RA|I|01|073|0023|A334|5|7|017|810|20070118|00:00|3.7||3||||||||||||| RD|I|01|073|0023|68103|5|7|017|810|20070121|00:00|6.9||3||||||||||||| RC|I|01|073|0023|Quer|5|7|017|810|20070124|00:00|1.8||3||||||||||||| infile 'C:\junk\RD_501_88101_2006-0.txt' dlm='|' firstobs=3 missover; rectype $2. @; if rectype = 'RD' then do; -- Zev Ross ZevRoss Spatial Analysis 303 Fairmount Ave Ithaca, NY 14850 607-277-0004 (phone) 866-877-3690 (fax, toll-free) zev at zevross.com

...166), c(0, 0.25, 0, 0.25), c(0, 0, 0, 0)) ematrix <- rbind( c(0, 0.01, 0, 0), c(0.01, 0, 0.01,0), c(0, 0.1, 0, 0), c(0, 0, 0, 0)) 2. the model: msm_covariates_sexandage <- msm(state ~ duration, subject = SerialNo, data = Data, qmatrix = twoway4.q, ematrix = ematrix,death = 4, obstrue = firstobs, covariates = ~ sex + age, method = "BFGS", control = list(reltol = 1e-10, fnscale = 50000, trace=1,REPORT=1)) msm stops after only 100 iterations and fails to compute the Hessian (See the output below). Why is it stopping at 100 only iterations yet i have specified reltol = 1e-10? On...

To Whomever is Interested, I have spent several days searching the web, help files, the R wiki and the archives of this mailing list for a solution to this problem, but nonetheless I apologize in advance if I have missed something obvious. The problem is this; I have a 5-column data frame with about 4.2 million rows, and want to create a new (and hopefully much smaller) data frame that

...uot;, which are character. Unfortunately, since the data are confidential, I'm unable to provide an example. The original SAS v6.12 code is provided below: ------------------------------------------------------------------ /* READ DATA */ DATA emiss; INFILE 'data.tab' LRECL=8000 FIRSTOBS=2 DLM='09'X MISSOVER DSD; INPUT study $ vehicle $ thc rv t5 t9 ar ol ox su bz; /* CREATE NEW VARIABLES */ ln_thc = log (thc); new = study||vehicle; /* PERFORM ANALYSIS */ PROC MIXED DATA=emiss MAXITER=1000 CONVH=1E-8 METHOD=REML NOCLPRINT NOITPRINT; CLASS new; MODEL ln_thc = rv...

...ately, since the data are confidential, I'm unable to > provide an example. > > The original SAS v6.12 code is provided below: > > ------------------------------------------------------------------ > /* READ DATA */ > DATA emiss; > INFILE 'data.tab' LRECL=8000 FIRSTOBS=2 DLM='09'X MISSOVER DSD; > INPUT study $ vehicle $ thc rv t5 t9 ar ol ox su bz; > > /* CREATE NEW VARIABLES */ > ln_thc = log (thc); > new = study||vehicle; > > /* PERFORM ANALYSIS */ > PROC MIXED DATA=emiss MAXITER=1000 CONVH=1E-8 METHOD=REML NOCLPRINT >...

...in the R-help list. Data comes from a multiyear field experiment in which 4 levels of a treatment (2, 3, 4, 6) are compared to see the effect on yield. It is a randomized complete block design. The SAS code follows: options ls=95; data uno; infile 'data.txt' delimiter=';' firstobs=2; input year plot block treat yield; run; proc mixed data=uno; class treat year block; model yield=block year treat treat*year; lsmeans year treat /pdiff; lsmeans treat*year /slice=year pdiff; ods output diffs=dos; run; data tres; set dos; if year=_...

In the model: lm.1 <- lm(variable ~ BLOC + TIL * YEAR , data=selvanera) I found TIL*YEAR interaction significant. Then I am trying to compare means of the different levels of TIL inside every YEAR using: mc.2 <- glht(lm.1, linfct = mcp(TIL*YEAR="Tukey")) summary(mc.2, test = univariate()) but it does not work. There is any way of doing this, like the SLICE option in

as a r noob i am having another problem: i have a big dataframe where each row corresponds to one entry and each column is a field... for instance, i have the column ID and time and many more... Id like to get a dataframe where all IDs are just included once (some users with that ID might have several entries but Id like to kepp only one).. when i use unique I only get a list of the levels (or