similar to: String replacement in an expression

Dear R-Users, I have a question concerning extraction of parameter estimates of variance function from lme fit. To fit my simulated data, we use varConstPower ( constant plus power variance function). fm<-lme(UPDRS~time,data=data.simula,random=~time,method="ML",weights=varConstPower(fixed=list(power=1))) I extract the results of this function by using the following codes:

Dear R-user, I'm having some difficulty with working PFIM 3.2, a package for implementing population PK/PD in R. I wish to evaluate the determinant of Fisher information matrix each time with successive dose from a pre defined sequence of doses and want to store those values in a vector. It's important to note that in my 'stdin.r' file, dose<-c(u) and each time u is to be

Hi The message occurred from R, when I was selected of "optimization > block diagonal Fhiser matrix" and used the attached file on PFIM. Could you please advise me about the following message? ***************************** Loading required pakage: combinat Attaching package:'combinat' The following object(s) are masked from 'package:utils':combn

Dear R users, Basically, from the following arbitrary data set: a <- data.frame(id=c(c("A1","A2","A3","A4","A5"),c("A3","A2","A3","A4","A5")),loc=c("B1","B2","B3","B4","B5"),clm=c(rep(("General"),6),rep("Life",4))) > a

Hello I would like to delete some values at random in a data frame. Does anyone know how I could do? With best regards Caroline

Dear list, I'd like to create a visual plot of a clmm() I've fitted using the 'ordinal' package in R. It's possible to do this with a glm() by using the 'effects' package. For example: library(effects) data(BEPS) mod <- lm(political.knowledge ~ age + gender + vote, data=BEPS) eff <- effect("age", mod, default.levels=100) plot(eff,

Hello, I hope this mailing list is correct. I've a cluster pacemaker with a clone resource ocfs2 with ocfs2-tools-1.4.1-25.6.x86_64 ocfs2-tools-o2cb-1.4.1-25.6.x86_64 on Opensuse 11.2 After some network problem on my switch I receive on one of 4 nodes of my cluster the following messages Aug 18 13:12:28 nodo1 openais[8462]: [TOTEM] The token was lost in the OPERATIONAL state. Aug 18 13:12:28

Dear all, I've come across a problem using strptime, can anyone explain what's going on? I'm using version 2.7.0 on Windows XP. Thank you Caroline First read in a data file using read.table alldata = read.table(file, header=F, skip=4, colClasses = c("character","numeric")) dim(alldata) [1] 223960 2 # inefficient, safe way of sorting out missing or dodgy

Solr 7 returns JSON by default but fts_solr requires XML. Would be good to have wt=xml added to the query to force it to xml all the time, this would prevent errors if solr has not had xml set as default for index. Dovecot 2.2.36 Centos 6 *Doveadm Request:* doveadm search -u caroline at covermg.com TEXT smith doveadm(caroline at covermg.com): Error: fts_solr: Invalid XML input at 1:0: not

Hello there, I’m computing the total value of an order from the price of the order items using a “for” loop and the “ifelse” function. I do this on a large dataframe (close to 1m lines). The computation of this function is painfully slow: in 1min only about 90 rows are calculated. The computation time taken for a given number of rows increases with the size of the dataset, see the example with

Dear List, I am using both the clm() and clmm() functions from the R package 'ordinal'. I am fitting an ordinal dependent variable with 5 categories to 9 continuous predictors, all of which have been normalised (mean subtracted then divided by standard deviation), using a probit link function. From this global model I am generating a confidence set of 200 models using clm() and the

I have a big problem with our cluster, if i read my mail with POP3 client (and check the option, keep a copy on server), every time I read my mail, my client retrive ALL messages an the mailbox !!! I test this with: - Apple Mail 1.0 + 2.0 - M$ Outlook 2000 + 2003 - Thunderbird We use: - Dovecot RC1 - Postfix on Debian Sarge (stable) Can you know this bug ? Or have a workaround ? This function

Dear all, I have a data set as follows: ID cycle.number cycle.result 1 2525 1 38 2 2525 2 38 3 2525 3 25 4 2525 4 25 5 2525 5 25 6 2525

Hi, Sorry if this is a repost. I searched but found no results. I am wondering if it is an easy way to construct the following matrix: r 1 0 0 0 r^2 r 1 0 0 r^3 r^2 r 1 0 r^4 r^3 r^2 r 1 where r could be any number. Thanks. Wen [[alternative HTML version deleted]]

Full_Name: Caroline Keef Version: 2.6.1 OS: Windows XP Submission from: (NULL) (195.171.203.131) I have tried to install the package rggobi which if I'm right requires the package RGtk2 If I install RGtk2 using the install.packages (I used the UK (Bristol) mirror, I haven't tried any other mirror) within R and then use library(RGtk2) I get the following message box "This

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

I have tried glm.nb in the MASS package, but many models (I have 250 models with different combinations of predictors for fish counts data) either fail to converge or even diverge. I'm attempting to use the negbin function in the AOD package, but am unsure what to use for the "random" term, which is supposed to provide a right hand formula for the overdispersion parameter.

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

I have a data set which contains two columns. The first column is a list of countries, and the second column contains their political risk ratings. I would like to create one large plot that contains 5 different sections, each with a scatter plot. To clarify, I have divided the countries into 5 groups. For each group (continent), I would like to have the name of the continent on the x-axis, and

On 6/22/2017 7:05 PM, Caroline wrote: > ##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') Since this dataset is only on your desktop it cannot help us