Технический форум
Вернуться   Технический форум > Программирование > Форум программистов > Помощь студентам


Ответ
 
Опции темы Опции просмотра
Старый 15.10.2015, 23:57   #1 (permalink)
SlowSpock
Новичок
 
Регистрация: 15.10.2015
Сообщений: 1
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию Вопрос по Wichmann-Hill RNG

Я написал эмуляцию Wichmann-Hill генератора случайных чисел в R:

Код:
wh.unif = function(n, start=0, end=1){
  x = double(n) #initialize x
  y = c(171,172,170)
  z = c(30269, 30307, 30323)
  #If seed does not exist - create it
  if (!exists(".WH.seed",envir=.GlobalEnv)){
    .WH.seed=trunc(runif(3, 1, 30269))
  }else{
    .WH.seed=get(".WH.seed",envir=.GlobalEnv)
  }
  for (i in 1:n) {
    .WH.seed = (y*.WH.seed) %% z
    x[i] = sum(.WH.seed/z) %% 1
  }
#Save the seed
  assign("wh",.WH.seed,envir=.GlobalEnv)
  start + (end-start) * x
}
Но, он выдает результат отличающийся от runif, а моя задача написать функцию которая будет выдавать такой же результат без использования .Random.seed[-1]

Код:
b=wh.unif2(5);b
set.seed(wh,"Wich")
a=runif(5);a
Как я понимаю, этот форум не совсем по программе R, но может сможете указать на ошибку в коде. Заранее спасибо!
SlowSpock вне форума   Ответить с цитированием

Старый 15.10.2015, 23:57
Helpmaster
Member
 
Аватар для Helpmaster
 
Регистрация: 08.03.2016
Сообщений: 0

Посмотрите на эти темы, наверняка там уже есть нужные ответы

Вопрос по UPS
Вопрос про БП
Hill Climb Racing - игнорируя законы физики забираемся на вершину
Вопрос по АКБ

Ads

Яндекс

Member
 
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070
Ответ

Опции темы
Опции просмотра

Ваши права в разделе
Вы не можете создавать новые темы
Вы не можете отвечать в темах
Вы не можете прикреплять вложения
Вы не можете редактировать свои сообщения

BB коды Вкл.
Смайлы Вкл.
[IMG] код Выкл.
HTML код Выкл.
Trackbacks are Вкл.
Pingbacks are Вкл.
Refbacks are Выкл.




Часовой пояс GMT +4, время: 05:48.

Powered by vBulletin® Version 6.2.5.
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.