search for: pstart

Hi , I am trying to get some estimator based on lognormal distribution when we have left,interval, and right censored data. Since, there is now avalible pakage in R can help me in this, I had to write my own code using Newton Raphson method which requires first and second derivative of log likelihood but my problem after runing the code is the estimators were too high. with this email ,I provide

...[,1] [,2] [1,] -55.99805 51.68817 [2,] -56.00222 51.68911 [3,] -56.01694 51.68911 [4,] -56.03781 51.68606 [5,] -56.04639 51.68759 [6,] -56.04637 51.69445 [7,] -56.03777 51.70207 [8,] -56.02301 51.70892 [9,] -56.01317 51.71578 [10,] -56.00330 51.73481 [11,] -55.99805 51.73840 attr(,"pstart") attr(,"pstart")$from [1] 1 attr(,"pstart")$to [1] 11 attr(,"nParts") [1] 1 attr(,"shpID") [1] NA [[2]] [,1] [,2] [1,] -57.76294 50.88770 [2,] -57.76292 50.88693 [3,] -57.76033 50.88163 [4,] -57.75668 50.88091 [5,] -57.75551 50.88169 [...

Dear Colleagues, I am attempting to develop an optimization routine for a river suspended sediment mixing model. I have Calcium and Magnesium concentrations (%) for sediments from 4 potential source areas (Topsoil, Channel Banks, Roads, Drains) and I want to work out, based on the suspended sediment calcium and magnesium concentrations, what are the optimal contributions from each source area to

Hi, Is there a way, preferably with R, to read shapefiles and transform them in a format that I could then use with package PBSmapping? I have been able to read such files into R with maptools' read.shape and plot it with plot.Map, but I'd like to bring the data to PBSmapping and plot from there. I also looked at the package shapefile, but it does not seem to do what I want

...1:10 #Generate 10 Gaussian random number with mean=3 sd=1 simy <- rnorm(length(x), mean=3, sd=1) #theoretical value for each of the 10 numbers=the adjust value y1 <- function(pp,xx) {rep(pp$a, length(xx))} #Residual resid <- function(pp,observed,xx) {observed-y1(pp,xx)} #Starting parameter pStart <- list(a=0.1) #non-linear fit library(minpack.lm) nls.lm.test <- nls.lm(par=pStart, fn=resid, observed=simy, xx=x, control=nls.lm.control(nprint=1)) It works fine: It. 0, RSS = 86.2811, Par. = 0.1 It. 1, RSS = 5.69735, Par. = 2.93873 It. 2, RSS = 5.69735, Par. =...

...* 3. Set start to *restart. If *restart == end, you are done */ static INLINE unsigned -nouveau_vbuf_split(unsigned remaining, unsigned overhead, unsigned vpp, - unsigned mode, unsigned start, unsigned count, - unsigned *restart) +util_split_primitive(int max, unsigned* pmode, unsigned* pstart, unsigned end, unsigned* flags) { - int max, adj = 0; + unsigned mode = *pmode; + unsigned start = *pstart; + unsigned count = end - start; + int adj = 0; - max = remaining - overhead; if (max < 0) return 0; - max *= vpp; if (max >= count) + { + if(*flags & SPLIT_CLOSE_LOOP)...

Hi there, I have following equations to be solved for a and b: a/(a+b) = x1 ab/((a+b)^2 (a+b+1)) = x2 Is there any direct function available to solve them without disentangling them manually? Thanks for your help.

What's up fellows... I am a begginer in R and i am trying to find the parameters of one likelihood function, but when i otimize it, always appers a error or advertisement and the solve does not occur. The problem seems like that: "lMix<-function(pars,y){ beta1<-pars[1] beta2<-pars[2] beta3<-pars[3] beta4<-pars[4] beta5<-pars[5] alfa1<-pars[6]

Dear All, In the moment i'm using the map and maptools package to read shapefiles and display the maps. I'm looking for the possibility to draw points (randomly positioned or positioned according to a grid) into the polygons instead of filling the polygons with colors. For example: a map (shapefile) with 10 countries, 15 points in the polygon of country A, 20 points in the

...space to relocate guest image */ - mmu_switch_ttb(gpt & ~(0x4000 - 1)); + set_ttbr((unsigned long)(gpt) & ~(0x4000 - 1)); elf.dest = (void *)ventry; elf_load_binary(&elf); @@ -192,7 +192,7 @@ int domain_construct(struct domain *d, si->mfn_list = 0; si->first_p2m_pfn = pstart >> PAGE_SHIFT; si->flags = 0; - si->min_mfn = pstart >> PAGE_SHIFT; + //si->min_mfn = pstart >> PAGE_SHIFT; if (d->domain_id == 0) { si->flags = SIF_PRIVILEGED | SIF_INITDOMAIN; @@ -202,7 +202,7 @@ int domain_construct(struct domain *d, VCPU_REG(...

in wine/dlls/qcap/v4l.c: ReadThread() --> OutputPin_GetDeliveryBuffer((OutputPin *)capBox->pOut, &pSample, NULL, NULL, 0); --> IMediaSample_SetTime(*ppSample, tStart, tStop); in wine/dlls/quartz/memallocator.c: StdMediaSample2_GetTime(IMediaSample2 * iface, REFERENCE_TIME * pStart, REFERENCE_TIME * pEnd) --> if tstart/stop is NULL, leave pStart/pEnd not setted. in wine/dlls/quartz/videorenderer.c: VideoRenderer_Sample(LPVOID iface, IMediaSample * pSample) --> IMediaSample_GetTime(pSample, &tStart, &tStop) --> use tStart/tStop and current time to judge when...

In converting some older code: ## builds an OWIN to generate a set of points on a grid within the OWIN object: bdry <- read.shape( "../../data/boundary.shp" ) outline.verts <- bdry$Shapes[[2]]$verts[1:bdry$Shapes[[2]]$Pstart[2],] outline.order <- nrow( outline.verts ) : 2 bdry.poly <- vector( 2, mode="list" ) bdry.poly[[1]] <- list( x=outline.verts[outline.order,1], y=outline.verts[outline.order,2] ) bdry.poly[[2]] <- list( x=bdry$Shapes[[1]]$verts[1:(nrow(bdry$Shapes[[1]]$verts...

Hello, This patch series adds an initial implementation of restriper (it''s a clever name for relocation framework that allows to do selective profile changing and selective balancing with some goodies like pausing/resuming and reporting progress to the user. Profile changing is global (per-FS) so far, per-subvolume profiles require some discussion and can be implemented in future.

...= (1 << blktosec) - 1; + block_t block; + size_t nblocks = 0; + + ufs_debug("ufs_next_extent:\n"); + block = ufs_bmap(inode, lstart >> blktosec, &nblocks); + ufs_debug("blk: %u\n", block); + + if (!block) // Sparse block + inode->next_extent.pstart = EXTENT_ZERO; + else + /* + * Convert blk into sect addr and add the remaining + * sectors into pstart (sector start address). + */ + inode->next_extent.pstart = + ((sector_t) (block << frag_shift) >> SECTOR_SHIFT(fs)) | + (lstart & blkmask); + + /* + * Subt...

...= (1 << blktosec) - 1; + block_t block; + size_t nblocks = 0; + + ufs_debug("ufs_next_extent:\n"); + block = ufs_bmap(inode, lstart >> blktosec, &nblocks); + ufs_debug("blk: %u\n", block); + + if (!block) // Sparse block + inode->next_extent.pstart = EXTENT_ZERO; + else + /* + * Convert blk into sect addr and add the remaining + * sectors into pstart (sector start address). + */ + inode->next_extent.pstart = + ((sector_t) (block << frag_shift) >> SECTOR_SHIFT(fs)) | + (lstart & blkmask); + + /* + * Subt...

Hi all, Here''s a work in progress series whcih does a partial revert of the previous swiotlb changes, and does a partial replacement with Becky Bruce''s series. The most important difference is Becky''s use of phys_addr_t rather than page+offset to represent arbitrary pages. This turns out to be simpler. I didn''t replicate the map_single_page changes, since

Hi all, this patch series enables xen-swiotlb on arm and arm64. It has been heavily reworked compared to the previous versions in order to achieve better performances and to address review comments. We are not using dma_mark_clean to ensure coherency anymore. We call the platform implementation of map_page and unmap_page. We assume that dom0 has been mapped 1:1 (physical address == machine

Hi Guys Thanks for this! As I am a beginner, I bet I'm running into some really basic problems. Using the example from the Map2poly function in maptools: try <- read.shape("euadmll.shp") mappolys <- Map2poly(try) # this produces 14 warnings() of the type: # 1: level 2 circularity at 62 in: .mtlbuild(.mtafters(r1), rD) # etc # 14: From next release, default hole handling

From: Raphael S. Carvalho <raphael.scarv at gmail.com> Wrote the documentation below. I think it would be good to push the doc to the wiki as soon as the UFS support gets merged. Unix Fast File System (UFS/FFS) 1/2 on Syslinux - (usage/install) ----- There is a confusion about the name of this file system, then I decided to contact the author who replied: "The name has always been

Hi all, this little series removes the last swiotlb API exposed to modules. Diffstat: arch/x86/include/asm/xen/swiotlb-xen.h | 6 ------ arch/x86/kernel/pci-dma.c | 28 ++++------------------------ drivers/gpu/drm/nouveau/nouveau_ttm.c | 10 +++------- drivers/pci/xen-pcifront.c | 6 ------ kernel/dma/swiotlb.c | 1 - 5 files changed, 7