similar to: Lack of independence in anova()

Dear R users, Hopefully someone can help me, Maybe I just misunderstand the function in the package? I am working with a logistic regression model. Until now I always worked with the basic glm function, where for the model was: ¡§ glm( disease ~ test.value + cnct , family=binomial(link=¡¦logit¡¦) ¡¨. This works fine when test .value and concentration (cnct) are continuous vairables. However,

I want to load many files in the R. The names of the files are "Sim1.txt", " Sim2.txt", "Sim3.txt", "Sim4.txt", "Sim5.txt" and so on. Can I read them at one time? What should I do? I can give the same names in R. Thanks. For example: > tst=paste("Sim",1:20,".txt",sep="") # the file names > tst [1]

Dear all, My double for loop as follows, but it is little efficient, I hope all friends can give me a "vectorized" program to replace my code. thanks x: is a matrix 202*263, that is 202 samples, and 263 independent variables num.compd<-nrow(x); # number of compounds diss.all<-0 for( i in 1:num.compd) for (j in 1:num.compd) if (i!=j) { S1<-sum(x[i,]*x[j,])

Hi, Can anyone help to write a more elegant version of my code? I am sure this can be put into a loop but I am having trouble creating the objects b1,b2,b3,...,etc. b1 <- rigamma(50,1,1) theta1 <- rgamma(50,0.5,(1/b1)) sim1 <- rpois(50,theta1) b2 <- rigamma(50,1,1) theta2 <- rgamma(50,0.5,(1/b2)) sim2 <- rpois(50,theta2) b3 <- rigamma(50,1,1) theta3 <-

Hi all, I met with a problem about the weighted least square regression. 1. I simulated a Normal vector (sim1) with mean 425906 and standard deviation 40000. 2. I simulated a second Normal vector with conditional mean b1*sim1, where b1 is just a number I specified, and variance proportional to sim1. Precisely, the standard deviation is sqrt(sim1)*50. 3. Then I run a WLS regression without the

2015-02-27 10:25 GMT-06:00 A J Stiles <asterisk_list at earthshod.co.uk>: > O.K. So what does your existing Dial() statement in extensions.conf look > like? > apology, put the gateway was sangoma but is a openvox , all my outgoing calls out for this context: [my-mobile-out] exten => _NXXXXXXX,n,Dial(SIP/1003/${EXTEN},55,rT) exten =>

I could use some advice. I've got 2 time series. Both cover approximately the same period of time (ie, 1940 to 2009). But one series has annual data and the other has monthly data. One refers to university enrollment; the other to unemployment rates. Both are currently in the same data frame. I'd like to use the monthly times series as a light grayscale background for a plot of the

To whoever this may concern I am trying to write a repeat loop and can't make out from the documentation on the website how exactly to construct the repeat, break structure of the loop. Below is the function sim2.dat that I am trying to create, in which firstly I create n random uniform(0,1) numbers. From there I assign values to k, P and F. From then on for each random number u[j] of u, I

Hello, I am a graduate student of University of Florida. I am trying to run a process variation model and I am facing the following error. ------------------------------------------------------------- Analysis of geostatistical data For an Introduction to geoR go to http://www.leg.ufpr.br/geoR geoR version 1.6-27 (built on 2009-10-15) is now loaded

On 13/09/17 13:24, K?tia Emidio wrote: > Dear Rolf, > > Thanks for your help! > > What I need is a spatial window with shape equal to the figure attached. > This figure I made using ArcGis, but it is important to me make it in R. > After having this figure I will make some analysis using spatstat among > others. The points within figure are trees... (1) I believe that

If you explain what radar track data is (for those like me who don't know) and what your needs are more specifically a number of people would probably offer pointers to R packages etc. Reid Huntsinger -----Original Message----- From: Scott_Summerill_Summerill/ACT/CNTR/FAA.FAA at faa.gov [mailto:Scott_Summerill_Summerill/ACT/CNTR/FAA.FAA at faa.gov] Sent: Thursday, March 21, 2002 3:22 PM To:

Hi, I am using the FracSim library to simulate a time series. However, the simulate function ignores my attempt to set the RNG seed I need for reproducible research. The published docs and google have not yielded an answer, so any help greatly received. Thanks, Selwyn ## Example code snippet library(FracSim) ## simulate some 1d fractal data set.seed(1234) sim1 = fracsim.1d(h=0.5,k=1000,n=5000)

Hello, I apologise if someone has already answered this but I searched and googled but didn't find anything. I have a matrix which gives me the similarity of each item to each other. I would like to turn this matrix into something like what they have in the graph package with the nodes and edges. http://cran.r-project.org/doc/packages/graph.pdf . However I cannot find a method to

Hello all, I hope something is not already posted regarding this exact problem I am trying to solve. I've read through the forums and previous postings and am still confused as to how to approach this. Basically, what I am trying to do is construct variables that utilizes an average of a variable from a grouping, or higher order, variable. For instance, in my dataset I have variables, with

I could use some advice regarding xyplot. I've got 2 time series. Both cover approximately the same period of time (ie, 1940 to 2009). But one series has annual data and the other has monthly data. One refers to university enrollment; the other to unemployment rates. Both are currently in the same data frame. I'd like to use the monthly times series as a light grayscale background for a

Dear useRs, I want to simulate a time series (stationary; the distribution of values is skewed to the right; quite a few ARMA absolute standardized residuals above 2 - about 8% of them). Is this the right way to do it? #-------------------------------- load("rdtb") #the time series > summary(rdtb) Min. 1st Qu. Median Mean 3rd Qu. Max. -1.11800 -0.65010 -0.09091

This is a generalize question: basically, say you are typing a long line of command in R console, and then you realize you forgot to add something in the beginning, is there a way to get to the beginning of the line without pressing the "left" key on the keyboard and waiting for the cursor to get to the beginning, or using the mouse? I'm using windows version of R. Thanks! -- View

Hello, This series exposes global performance counters (PCOUNTER) to the userspace through the nvif interface by reworking most of the code related to the PM engine. This interface will allow the userspace to control and monitor complex hardware events like the proprietary driver already does, for example with CUPTI and PerfKit. For now, this series only exposes performance counters on NV50,

Hi All I dont know if anyone else has noticed the same thing, but with 2 subsequent runs of the same syntax, I am getting exactly the same results. I was expecting that results differ slighlty, say in the 4th or 5th decimal place. Is this a specialty with R2WinBUGS? Does it have something to do with the seed value? Isnt the seed value reset everytime I restart winbugs? Thanks Toby

I am currently testing species co-occurrence patterns using null models and the oecosimu() function within the vegan() package. My issue is that none of the methods appear to be the ones that I want. The methods listed are r0, r1, r2, r2dtable, swap, tswap. However, I want to know how to go about implementing fixed row algorithms, as suggested in Gotelli 2000 in Ecology. Also, the null models