similar to: Intercept in model formulae

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =

To anyone who can help: I have two brief questions concerning sapply. Following below is the code for my example. The two problems are described at the end of the code: site <- rep(2:6, each = 12) tillage <- rep(c(1,-1), each = 6, times = 5) carbon <- c(18.23, 16.06, 17.81, 16.07, 17.26, 17.08, 14.92, 15.88, 12.11, 14.23, 16.99, 13.57, 20.34, 20.3,

Dear List, i am trying to produce a 3d plot using wireframe using the code: wireframe(Residuals_FD ~ Elevation * Temperature, data = data2, scales = list(arrows = FALSE), drape = TRUE, colorkey = TRUE) As you can see when the code (using the data below) is run the plot area is set-up correctly but the actual surface is missing? Any help would be greatly appreciated. Chris #data Elevation

We upgraded to NSD 3.2.9 (from 3.2.8) because we encountered the problem "Fix denial of existence response for empty non-terminal that looks like a NSEC3-only domain (but has data below it)." (a nasty problem with DNSSEC). But we now have IXFR issues. On one name server, NSD 3.2.9 works fine, zones are IXFRed and work. On another name server, with much more zones (and big ones), we

Hello I have a time-comsuming program which need to simplify, I have tested the annotated program as follow: > #define function which will be call > calsta <- function(c, n=100000) + { + i <- seq(from=0, length=c) + logx <- lchoose(NT-n, CT-i) + lchoose(n, i) + logmax <- max(logx) + logmax + log(sum(exp(logx - logmax))) + } > CT=6000 #assignment to CT >

Hello R users, For the last few days I am struggling with the following task: my data.frame: A1 A2 A3 B1 B2 B3 58.81 53.292 54.501 13.013 17.39 19.407 56.02 56.251 54.033 20.099 13.15 10.411 55.376 53.099 57.625 13.396 21.031 13.22 58.584 53.194 54.218 13.038 16.854 19.289 55.7 55.921 53.847 19.942 13.153 9.828 55.093 52.934

Dear R members, I have been trying to find out whether one can use multivariate regression trees (for example mvpart) to analyze repeated measures data. As a non-parametric technique, CART is insensitive to most of the assumptions of parametric regression, but repeated measures data raises the issue of the independence of several data points measured on the same subject, or from the same plot

Hello, I'm EDWIN, I create (make) GUI, with call many function but I don't know why when I call function I can't. if without function, Yes I can.. can you help me ? can you make this, become true with full code? Can you help me to create data.entry with interface LM - BETA1.HAT - BETA2.HAT SD.BETA1.HAT HAT SD.BETA2.HAT RLM - BETA1.HAT - BETA2.HAT SD.BETA1.HAT HAT SD.BETA2.HAT 2 10

Is there a way to get confidence intervals around an ordered regression like polr() in the MASS package? ------------------------------------------- Joe King, M.A. Ph.D. Student University of Washington - Seattle 206-913-2912 <mailto:jp@joepking.com> jp@joepking.com ------------------------------------------- Ad astra per aspera - "Through hardships to the stars"

I got two .RData file e.g data.2005.RData & data.2006.RData I would like to combine these two different data set and make single RData file. in both file there are some NULL files are also available and I would like to clear this NULL files also. $ : NULL $ : NULL $ : num [1:43285, 1:8] -21.1 -21.1 -24.9 -24.9 -24.9 ... ..- attr(*, "dimnames")=List of 2 .. ..$ : NULL ..

I am trying to find out if someone has implemented a (McFadden-type) multivariate binary logistic regresssion package for R? From what I can tell, this is not available for R. Thank you, Lynne Baker [[alternative HTML version deleted]]

Hello. I am trying to do a non-linear fit using the 'nls' command. The data that I'm using is as follows pH k 1 3.79 34.21 2 4.14 25.85 3 4.38 20.45 4 4.57 15.61 5 4.74 12.42 6 4.92 9.64 7 5.11 7.30 8 5.35 5.15 9 5.67 3.24 with a transformation of pH to H <- 10^-pH When using the nls command for a set of parameters - a, b and c, I receive two sets of errors: >

Hi every one, I am using the package 'drc' to model root elongation using dose response data. I don't know which model I should use. Though I don't know which model I should use, I tried the following codes given below. But it produced the error messages.Can any one tell me the code in 'drc' package to find out the EC (Effective Concentration) values and Confidence

Hello, all R-users! I am working on fitting a survival analysis model using the coxph function for Cox proportional hazards regression model. Data look like usual: ========================== group block death censor Group1 1 4 1 Group1 1 12 1 ... Group2 30 4 1 Group2 30 4 1 ... Group3 57 16

I've been using R to do logistic regresssion, and that's working well, but there are two things I haven't figured out how to do. (1) Is there some pre-existing function that will let you compute the odds ratios and confidence intervals for them for a specific fit. I know how to do this manually or even write a function that I can call with the coefficients and se, but

Hello, I am trying to filter a data set like below so that the peaks in the Phase value are more obvious and can be identified by a peak finding function following the useful advise of Carl Witthoft. I have written the following for(i in length(data$Phase)){ newphase=if(abs(data$Phase[i+1]-data$Phase[i])>6){ data$Phase[i+1] }else{data$Phase[i] } } I get the following error which I have not

On 5 Feb 2017, at 19.49, Timo Sirainen <tss at iki.fi> wrote: > > On 4 Feb 2017, at 20.03, Ruga <ruga at protonmail.com> wrote: >> >> https://github.com/libressl-portable/portable/issues/278 > > I've no idea why that would happen. The only idea I had got rejected by someone. Oh, that's with OSX. I think that's the reason. Nothing to do with

A couple of questions: How to be sure that gls allowes errors to be correlated and/or have unequal variances? (is this on auto or is there a switch?) How to calculate confidence limits for a linear regresssion? -------------- next part -------------- A non-text attachment was scrubbed... Name: dthompson.vcf Type: text/x-vcard Size: 303 bytes Desc: Card for David Thompson Url :

Dear all: Does anyone has harmonic regresssion analysis package written in R (to be used in Windows platform) ? Thanks _______________________________________ YM - Â÷½u°T®§ ´Nºâ§A¨S¦³¤Wºô¡A§AªºªB¤Í¤´¥i¥H¯d¤U°T®§µ¹§A¡A·í§A¤Wºô®É´N¯à¥ß§Y¬Ý¨ì¡A¥ô¦ó»¡¸Ü³£ÉN¨«¥¢¡C http://messenger.yahoo.com.hk [[alternative HTML version deleted]]

> On 13 Nov 2018, at 21.06, Arkadiusz Mi?kiewicz <arekm at maven.pl> wrote: > > On 13/11/2018 15:54, Arkadiusz Mi?kiewicz wrote: >> >> 2.2.36 (not migrated to 2.3 yet) reports such problem: >> >>> Nov 13 15:50:58 mbox dovecot: lmtp(xxx): session=<ACYRCtLk6ltiEQAALZVUYQ>, Error: open(/var/mail/xxx/mailboxes.lock1bf6ad16b7b8b703) failed: Disk quota