similar to: Power analysis for glm-function

I remember reading about post hoc statistical power on R-help. But I can't seem to find them with RSiteSearch("post hoc statistical power") and variations of it. I would like to learn more about post hoc statistical power, its meaningfulness, advantages and disadvantages. I thought the issue was settled after Tukey's 1993 paper about post hoc statistical power being

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings around each. Each tree has a TreeID, which I use as the random effect (blocking factor). There are two

Hello R-user! I am running R 2.7.0 on a Power Book (Tiger). (I am still R and statistics beginner) I try to find the most important variables to divide my dataset as given in a categorical variable using randomForest. Is randomForest() able to deal with count data? Or is there no difference because only the ranks are used in the trees? Thanks in advance Birgit Birgit Lemcke Institut f?r

Dear all, I would like to compare two distance matrices to find the correlation between those and create a XY plot.could you please kindly help me and let me know how I can do that using R.(Mantel test?). Your help is very much appreciated!! Many thanks! Regards, Sima Fakheran ............................................................... Sima Fakheran Esfahani PhD Student Institute of

Dear Sir, When requesting the summary of pglm analyses in R, I encounter repeatedly the following warning message: "Warning message: In pf(Fstat, object$k - 1, object$n - object$k, ncp = 0, lower.tail = FALSE, : full precision was not achieved in 'pnbeta'" This message appears both with the pglm estimate of Lambda and with Lambda set to 0. Could you tell me what

Hi everyone, Is it possible to use tapply(x,y,mean) if not all groups of x by y are of the same length (for example if you have one missing observation)? I tried tapply(x,y,mean,na.omit=T) but it doesn't work! Steffi -- --------------------------------- Stefanie von Felten Doktorandin ETH Z?rich Institut f?r Pflanzenwissenschaften ETH Zentrum, LFW A 2 Telefon: 044 632 85 97 Telefax: 044

Good afternoon! I would like to ask you about similarity measures and clustering in R for Binary data. Would you please kindly help me and let me know about that commands in R? Thanks in advance for your kind attentions. I look forward to hearing from you as soon as possible. Best regards, Sima

Hello R-user! I am running R 2.7.0 on a Power Book (Tiger). (I am still R and statistics beginner) I tried rfImpute (randomForest) and as far as I understood should it replace NA`s using a proximity matrix: > set.seed(100000) > Subset5Imputed<-rfImpute(Sex~., data=Subset5) ntree OOB 1 2 300: 11.78% 12.36% 11.21% ntree OOB 1 2 300: 12.07% 12.64%

Dear R-users Problem: Given the following system of ordinary differential euqations dM/dt = (-n)*M-h*M dS/dt = n*M-h*S+u*R dA/dt = h*S-q*A dI/dt = q*A-p*I dJ/dt = h*M-v*J dR/dt=p*I+v*J-u*R where M,S,A,I,J,R are state variables and n,h,u,q,p,v parameters. I'm able to calculate the likelihood value based on the solutions M,S,A,I,J,R of the ODE's given the data, but without an explicit

R version 2.6.2 PowerBook G4 Hello R User, I try to perform an ANCOVA using the glm function. I have a dataset with continuous and categorical data (explanatory variables) and my response variable is also binary categorical. Fehler: NA/NaN/Inf in externem Funktionsaufruf (arg 4) Zus?tzlich: Warning messages: 1: In Ops.factor(y, mu) : - nicht sinnvoll f?r Faktoren (makes no sense for

Hello R-user! I am running R 2.7.0 on a Power Book (Tiger). (I am still R and statistics beginner) I try to find the most important variables to divide my dataset as given in a categorical variable. code: Test.rf4<-randomForest(Sex~.,na.action=na.roughfix, data=Subset4, importance=TRUE, proximity=TRUE, ntree=10000, do.trace=1000, keep.forest=FALSE) My dataset contains also ordered

If I save a dataframe with a character-typed last column to a relational database with dbWriteTable, the values in the last column of the resulting table in the database will have a '\r' (carriage return) appended. If I read back the dataframe with dbReadTable the last column in the resulting dataframe has also '\r' appended (see protocol below). Setting or unsetting sql-mode

Hi I would like to fit a non-linear regression to the follwoing data: quantiles<-c(seq(.05,.95,0.05)) slopes<-c( 0.000000e+00, 1.622074e-04 , 3.103918e-03 , 2.169135e-03 , 9.585523e-04 ,1.412327e-03 , 4.288103e-05, -1.351171e-04 , 2.885810e-04 ,-4.574773e-04 , -2.368968e-03, -3.104634e-03, -5.833970e-03, -6.011945e-03, -7.737697e-03 , -8.203058e-03, -7.809603e-03, -6.623985e-03,

Hello R-user community! I am running R 2.7.0 on a Power Book (Tiger). (I am still R and statistics beginner) Presently I try to run the function metaMDS (vegan) using an existing dissimilarity-matrix. As I would like to start with this matrix I thought I could just give the matrix using the x= -argument Test<-metaMDS(x=Dist.Gower) Fehler in inherits(comm, "dist") :

Hi Dirk, Thanks for the explamnayion - Debian is running in a VM (and nothing really installed on it) , and I think the easiest is for me just to install Debian 10 and to use that one. Thanks a lot, Rainer > On 22 Sep 2020, at 15:14, Dirk Eddelbuettel <edd at debian.org> wrote: > > > On 22 September 2020 at 14:49, Rainer M Krug wrote: > | I know this is likely

Full_Name: Menelaos Stavrinides Version: 1.3. 1 OS: Windows 98 Submission from: (NULL) (193.129.76.90) When model simplification is used in glm (binomial errors) and anova is used two compare two competitive models one can use either an "F" or a "Chi" test. R always performs an F test (Although when test="Chi" the test is labeled as Chi, there isn't any

Hi I know this is likely documented somewhere, but I can?t find it. How can I install R > 4.0.0 on Debian 9 Stretch? Thanks a lot, Rainer -- Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany) Orcid ID: 0000-0002-7490-0066 Department of Evolutionary Biology and Environmental Studies University of Z?rich Office Y34-J-74

I can think of many brute-force ways to do this outside of R, but was wondering if there was a simple/elegant solution within R instead. I have a table that looks something like the following: Factor1 Factor2 Value A 11/11/2009 5 A 11/12/2009 4 B 11/11/2009 7 B 11/13/2009 8 >From that I need to generate all permutations of Factor1 and Factor2 and force a 0 for any combination that doesn?t

Is names(model1$coef) what you're looking for? -----Original Message----- From: Kenneth Cabrera [mailto:krcabrer at epm.net.co] Sent: 29 November 2002 10:36 Cc: R-help at stat.math.ethz.ch Subject: [R] Obtaining the variable names of a glm object Hi, R users! Suppose I make a model like this:

Hello, as far as I see, R reports type I sums of squares. I'd like to get R to print out type III sums of squares. e.g. I have the following model: vardep~factor1*factor2 to get the type III sum of squares for factor1 I've tried anova(lm(vardep~factor2+factor1:factor2),lm(vardep~factor1*factor2)) but that didn't yield the desired result. Could anyone give me a hint how to proceed?