similar to: how to overlook the zero in the denominator

Dear R users: I have a problem about catch the value from function. I have following two functions (part): fbolus1 <- function() {......... par<-data.frame(Parameter=c("kel","Vd"),Initial=c(0)) check(par) .....} check<-function(par) { if (par[ ,2] <= 0){ cat("\nEnter again (y/n) ?\n\n") ans<-readline() if (ans == "n" ){

Dear R users: I encountered difficulties in michaelis-menten equation. I found that when I use right model definiens, I got wrong Km vlaue, and I got right Km value when i use wrong model definiens. The value of Vd and Vmax are correct in these two models. #-----right model definiens-------- PKindex<-data.frame(time=c(0,1,2,4,6,8,10,12,16,20,24),

Dear R users: I wonder if it is possible to use "data.frame" and "edit" to show the form like: parameter upper lower A B C when I type the command edit(abc), assume the data.frame named "abc". And then I want to extract the number form it, like abc[1,2] which mean the upper level of A, to do something. The following is my example and the warning

Dear R users: I have a problem about catch the value from function. I have following two functions (part): sbolus1 <- function() { ....... for( i in 1:Subject) { kel<-par1 Vd<-par2 PKindex<-sbolus1.out(PKtime,kel,Vd,defun,par1,par2,Dose,i) } savefile(PKindex) } sbolus1.out<-function(PKtime,kel,Vd,defun,par1,par2,Dose,i) { time<-PKtime$time

Dear R user: I bulid a package, and in the package I use the function "nls" to solve some questions. If I have two sets of data, and I want to summary these two data's nls output, I write the command in the package source code like: { .......... summary(fm1) summary(fm2) } then i compiler the package and use "Install package(s) from local zip files" to install my

Hi all I tried to reproduce an example with lme and used the Orthodont dataset. library(nlme) fm2a.1 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1 | Subject) anova(fm2a.1) > numDF denDF F-value p-value > (Intercept) 1 80 4123.156 <.0001 > age 1 80 114.838 <.0001 > Sex 1 25 9.292 0.0054 or alternatively

I am really new to R and would appreciate some help to sort out a problem with data extraction from a large file. I have entered the following command to filter data from a large data set called "upanddown" such that those records in upanddown corresponding to blockType=3'UTR are transferred to a new file called "new": new<-sqldf("select * from upanddown where

Hi I just started using RGui.exe under widnows. I have a text file containing date arranged in columns and rows, each column has the same format, each row with different formats. 3 of the columns are something like this 1/12/2006 3:59:45 PM I need to calculate the different in seconds between 2 selected periods using their row’s index My solution: Read the file in a data frame and

Hello all, I got some severe iscsi connection loss on my dom0 (Gentoo 2.6.20-xen-r6, xen 3.1.1). Happening several times a day. open-iscsi version is 2.0.865.12. Target iscsi is the open-e DSS product. Here is a snip of my messages log file: May 5 16:52:50 ying connection226:0: iscsi: detected conn error (1011) May 5 16:52:51 ying iscsid: connect failed (111) May 5 16:52:51 ying iscsid:

I have a dataset called "results" that looks like this: arrive depart intercept 1 1 1 1 2 1 1 3 1 1 2 2 1 3 2 1 3 3 2 2 2 2 3 2 3 3 3 where arrive is the period of arrival, depart is the period of departure, and intercept

Hello :) I have two data frames, one has properties taken on a piece by piece basis and the other has performance on a lot by lot basis. I wish to combine these two data frames but the problem is that each lot has multiple pieces and hence i need to take a mean of the properties of multiple pieces and match it to the row having data about the lot. I was wondering if there is a simple commmand,

Hi friends - I'm plotting a ratio of bicarbonates i ggplot2 and ylab(expression(paste(frac("additive BIC","true BIC")))) worked OK - but now I have been asked to put the chemistry instead - so I wrote ?ylab(expression(paste(frac("additive",HCO[3]^"-","true",HCO[3]^"-")))) - and frac saw that as additive = numerator and HCO3- =

Dear R-core There is a different output for the apply function due to the change of unlist as mentioned in the R news. Newly, applying as.factor() (or factor()) in str(dat <- data.frame(x = 1:10, f1 = gl(2,5,labels = c("A", "B")))) (d1 <- apply(dat,2,as.factor)) newly returns a character matrix while in R-2.3.1 the same command resulted in an integer matrix that was

Hello, I am struggling to understand how denominator degrees of freedom and subsequent significance testing based upon them works in nlme models. I have a data set of 736 measurements (weight), taken within 3 different age groups, on 497 individuals who fall into two morphological catagories (horn types). My model is: Y ~ weight + horn type / age group, random=~1|individual I am modeling

Just giving the survey package a spin... I'm accustomed to stata, and it seems very similar in many respects. One thing is throwing me, however. I've gotten my data in, and specified the design. Looks like the weighting is right (based on published population estimates from these data), but now I'd like to check my "marginal means" for proportions against those that have

hello I need to know the dfn and dfd of my Anova. But in the Anova output there is only "Df". Is this the dfn or the dfd? and how do I get both of it in R? Thanks for any answers -- View this message in context: http://r.789695.n4.nabble.com/df-of-numerator-and-denominator-tp3765526p3765526.html Sent from the R help mailing list archive at Nabble.com.

Hi, wbinfo -Y BUILTIN\group can work without idmap_rid in Samba-3.0.14a. But I'm experiencing wbinfo -Y with idmap_rid failed for SID to GID conversion of BUILTIN groups. Since idmap_rid only works in a single domain, and captures workgroup's domain sid as a real domain sid in rid_idmap_get_domains(), when running "wbinfo -Y BUILTIN/System Operators", the function

?s 10:01 de 02/02/2024, Troels Ring escreveu: > Hi friends - I'm plotting a ratio of bicarbonates i ggplot2 and > > ylab(expression(paste(frac("additive BIC","true BIC")))) worked OK - but > now I have been asked to put the chemistry instead - so I wrote > >

Hello all. I was wondering if there is any way to adjust the denominator degrees of freedom in lme(). It seems to me that there is only one method that can be used. As has been pointed out previously on the list, the denominator degrees of freedom given by lme() do not match those given by SAS Proc Mixed or HLM5. Proc Mixed, for example, offers five different options for computing the

Hello everbody, i'm currently using the anova()-test for a small data.frame of 40 rows and 2 columns. It works well, but is there any preconditions for a valid variance analysis, that i should consider? Thank you for your answer, Daniel