similar to: Spatstat rmh function error message

Hello, Somehow I cannot control my envelope plots on Spatstat package. I would like to plot an envelope with no legend and also I would like to label the plot as such: >plot(envelope(data, nsim=39, main="My K Plot", legend=FALSE, xlab="Distance") But somehow the "main", "legend" and "xlab" calls do not have any effect on the plot. Can someone

Hello R Users, I am trying to do a quadrat count defined by covariate properties in spatstat. I have read my elevation raster into R (from ascii) and converted to class "im" for use in spatstat. Now I have point data of class "ppp" which window is the same extent as the elevation image. I can conveniently plot the image and the points on the image as follows:

Hello R Users, I am trying to do a quadrat count defined by covariate properties in spatstat. I have read my elevation raster into R (from ascii) and converted to class "im" for use in spatstat. Now I have point data of class "ppp" which window is the same extent as the elevation image. I can conveniently plot the image and the points on the image as follows:

Hello UseRs, I need help on how to plot several raster images (such as those obtained from a kernel-smoothed intensity function) in a layout such as that obtained from the lattice package. I would like to obtain something such as obtained from using the "levelplot" or "xyplot" in lattice. I currently use: >par(mfrow=c(3,3) to set the workspace, but the resulting plots

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

I don't understand what is happening! I have a (large) object sim1, an matrix list with dim c(101,101) where each element is an 3*3 matrix. I am subsetting that with a matrix coo, of dim c(100,2), of unique indices, but the resulting object has length 99, not 100 as expected. Code reproducing the problem follows: library(RandomFields) set.seed(123) sim0 <- GaussRF(x=seq(0, 100, by=1),

Dear all, I'm working with two point patterns (ppp) in spatstat. I turned one of them into a spatial covariate (im) object. After that, I used this im object to fit a Poisson model for the second point pattern, using the covariate layer from the first one. In R, the whole thing looks somewhat like this: my_first.im <- as.im(my_first.ppp) test.ppm <- ppm (my_second.ppp, ~my_first,

Hi, I have a Excel file with three spreadsheets: PlanA, PlanB and PlanC. I'm trying to read the three spreadsheets and then adding them together. But, when I try read the PlanA there is an error message: rm(list=ls()) setwd('C:/Test/Dados/Teste') require(RODBC) Arquivo = odbcConnectExcel('T070206_1347.xls') (Geral = sqlFetch(Arquivo, 'PlanA')) (Lactacao =

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]

Hi. No biggie, but when I remove a umass device it appears that two kernel messages are being written at the same time. (da4:umass-sim1:1:0:0): lost device - 0 outstanding, 1 refs (da4:(pass6:umass-sim1:1:umass-sim1:1:0:0:0): removing device entry 0): passdevgonecb: devfs entry is gone (FreeBSD 9.1-PRERELEASE #0 r241067) Cheers Michiel

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 <-

Dear all, I use Trixbox as my PBX. Until a couple of days I've installed a GSM Gateway to communicate with our three cellular phones: 15 64227777 15 64228888 15 64229999 The GSM Gateway has just one SIM. I use the Free PBX web interface in order to set up the route and trunk parameters: Trunk: ******* Name: SIM1 Peer details: host=10.10.1.2 (IP from GSM Gateway) port=5060 type=peer

I have a dataframe with many rows like this: > df X1 X2 X3 X4 X5 X6 X7 week d sim1 FALSE TRUE TRUE TRUE TRUE TRUE TRUE 1 0.3064985 sim1 is the rowname, X1..X7,week,d are the column names. X1..X7 are factors, booleans in this case. I need to add another row, represented by the following list: list(rep(T,7),5,0.0) -- i.e, TRUE in all boolean columns,

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)

Dear R-users, I would like to show you a simple example that gives an overview of one of my current issue. Although my working setting implies a different parametric model (which cannot be framed in the glm), I guess that what I'll get from the following example it would help for the next steps. Anyway here it is. Firstly I simulated from a series of exposures, a series of deaths (given a

Hi all: I'm working on a project that I have packaged for ease of distribution. The different simulations in the package share code, so obviously I have those parts organized as functions. Now, I want to show people my code, but the structure with the internal functions might be a little confusing to follow. One thing I tried was to have the code of the functions as their own R files in the

On Friday 27 Feb 2015, ricky gutierrez wrote: > the problem is that my pbx all incoming calls using only the channel > gsm 1 , the idea is that an incoming call to channel 1 is passed to > channel 2 Ah. *Incoming* calls are not something that is within your control; they have already been routed onto a line by your telco. So you will need to speak to someone at your telco about doing

Hi Folks, I've facing an intermittent hang with a USB port which seems cam related: Event's that happen are: o USB modem (HUAWEI E220) plugged into PC ugen3.2: <HUA WEI> at usbus3 u3g0: <3G Modem> on usbus3 u3g0: Found 3 ports. umass0: <USB MASS STORAGE> on usbus3 umass0: SCSI over Bulk-Only; quirks = 0x0000 umass0:6:0:-1: Attached to scbus6 umass1: <USB

Hi, I'm trying to install Fedora Core 3 on an AMD 64-based machine using PXE (via DHCP and TFTP). I've managed to get the machine to load the kernel over the network and start booting, but instead of loading the root fs from initrd.img (that I appended as a config via pxelinux.cfg), it loads the pre-existing rootfs on my harddisk. At the top of dmesg is the message: "so old

On 26/12/2019 11:14 p.m., Spencer Graves wrote: > Hello, All: > > > ????? The default "simulate" method for lm and glm seems to ignore the > sampling variance of the parameter estimates;? see the trivial lm and > glm examples below.? Both these examples estimate a mean with formula = > x~1.? In both cases, the variance of the estimated mean is 1. That's how