Novichok+: R Нужна функция которая берет каталог файлов данных и порог для полных случаев и вычисляет корреляцию
Novichok+
unregistered
Автор темы
Ответить
N
Novichok+ unregisteredАвтор темы
9 лет назад / 22 октября 2015 00:27
Создать функцию которая, которая принимает каталог файлов данных и порог для полных случаев и вычисляет корреляцию между сульфатом и нитратом для мониторов местах, где число полностью наблюдаемых случаев (на всех переменных) больше, чем порог. Функция должна возвращать вектор корреляций для мониторов, которые отвечают требованиям порогового. Если мониторы не отвечают требованиям порога, то функция должна возвращать числовой вектор длины 0.

Содержимое файла - кактинка 1-0
Результат работы функции - картинка 3-1

Прототип функции:
01corr <- function(directory, threshold = 0) {
02  ## 'directory' is a character vector of length 1 indicating
03  ## the location of the CSV files
04  ## 'threshold' is a numeric vector of length 1 indicating the
05  ## number of completely observed observations (on all
06  ## variables) required to compute the correlation between
07  ## nitrate and sulfate; the default is 0
08  ## Return a numeric vector of correlations
09  ## NOTE: Do not round the result!
10}
К сообщению приложены файлы:
loyolla
Member
542/569 ответов
19 лет на iXBT, с декабря 2005
Чаще пишет РІ "Программирование" (95%)
Инфо Ответить
l
loyolla Member
9 лет назад / 22 октября 2015 01:05
Все, что необходимо для написания этой функции, содержится в моих ответах : тут и тут

Мне кажется, не должно составить труда сделать имплементацию действуя по образу и подобию.

Функция для вычисления корелляции называется cor.
Novichok+
unregistered
Автор темы
Ответить
N
Novichok+ unregisteredАвтор темы
9 лет назад / 22 октября 2015 02:07
loyolla:
Все, что необходимо для написания этой функции, содержится в моих ответах : тут и тут
Мне кажется, не должно составить труда сделать имплементацию действуя по образу и подобию.
Функция для вычисления корелляции называется cor.
Спасибо,
первый ТУТ ведет к предыдущей теме а не к новой функции.
Novichok+
unregistered
Автор темы
Ответить
N
Novichok+ unregisteredАвтор темы
9 лет назад / 22 октября 2015 22:27
Запускаю код (ниже), получаю ошибку:

> cr <- corr("R/specdata", 150)
Error in df["nobs"] : object of type 'closure' is not subsettable
> head(cr)
Error in head(cr) : object 'cr' not found

Как исправить?


01corr <- function(directory, threshold = 0) {
02   
03  ids = df[df["nobs"] > threshold, ]$id
04  corrr = numeric()
05  for (i in ids) {
06    newRead = read.csv(paste(directory, "/", formatC(i, width = 3, flag = "0"),
07                             ".csv", sep = ""))
08    dff = newRead[complete.cases(newRead), ]
09    corrr = c(corrr, cor(dff$sulfate, dff$nitrate))
10  }
11  return(corrr)
12}
13  
14cr <- corr("R/specdata", 150)
15head(cr)

Добавление от 22.10.2015 22:37:

Нашел решение, не нужно.
Ваш ответ:

Нет значка Нет значка Р’РѕС‚ тут! Лампочка Восклицание Р’РѕРїСЂРѕСЃ Класс! Улыбка Злость Огорчение РџРѕРіРѕРІРѕСЂРёРј? Краснею Подмигивание Ругаю РћРґРѕР±СЂСЏСЋBIUdelSxsupxsuboffsp spoilerqurlimgvideo• list1. list1 codeprecenter-hr-rusQWE→ЙЦУ
файлыочистить
Ваше имя: Авторизуйтесь Предпросмотр В полную форму
вставить выделенную цитату в окно ответа
Если Вы считаете это сообщение ценным для дискуссии (не обязательно с ним соглашаться), Вы можете поблагодарить его автора, а также перечислить ему на счет некоторую сумму со своего баланса (при отзыве благодарности перечисленная сумма не будет вам возвращена).
Также вы можете оценить сообщение как неудачное.
В течение суток можно 20 раз оценить сообщения разных участников (купите Premium-аккаунт, либо оплачивайте оценки сверх лимита).
Если Вы считаете это сообщение ценным для дискуссии (не обязательно с ним соглашаться), Вы можете поблагодарить его автора, а также перечислить ему на счет некоторую сумму со своего баланса (при отзыве благодарности перечисленная сумма не будет вам возвращена).
Также вы можете оценить сообщение как неудачное.
В течение суток можно 20 раз оценить сообщения разных участников (купите Premium-аккаунт, либо оплачивайте оценки сверх лимита).
Последние обсуждения в Конференции
04:29Пандемия COVID-2019 (было - Миру угрожает вирус пандемического гриппа, может погибнуть МИЛЛИАРД человек.) Общий
04:28Электромобили - грядет революция? Авто
04:25РћС‚ Windows Рє Linux Unix
04:25Будет ли война России с Украиной? Политика
04:21ИИ - осознание себя? ИИ
04:17Люблю ругалики Флуд
04:15FF БЗК Canon EOS R (RP, R5, R3, R1 и т.д.) и оптика к ним Фото
04:05Выбор посудомоечной машины Бытовая техника
03:55Сисадминская курилка Администрирование
03:50Я плакаль... (картинка прикол) Юмор
03:47Ресиверы Yamaha - пообщаемся? FAQ в первом посте ДК аудио
03:44Вопросы знатокам o S.M.A.R.T. Накопители
03:40Отзывы о PrivalSystems (аналог skype) Интернет
03:40Выбор цветного струйного МФУ формата А3 в офис Периферия
03:28Как выразить АНАЛИТИЧЕСКИ значения корней алгебраического уравнения 5-й степени через значения корней его резольвенты 6-й степени (которые найдены, например, численно) Наука
03:28Поговорим о припоях и флюсах Эл. устройства
03:18Path of Exile. Врываемся в приятно выпуклый эндгейм PoE 2 6 декабря 2024-го и ждём 3.26 в первой части когда-нибудь! Игры
02:52Выбор бюджетного "цифрового" (китайского) усилителя Стерео
02:38Какие триллеры посоветуете? Кино
02:35Выбор матраса Ремонт
Последние темы без ответов
вчераАрхаичные выражения, подзабытые слова, названия, а также уходящие в прошлое культурные традиции XX века. Культура
вчераBattlefield 2042 drweb Игры
вчераМышь лагает и плавает при запуске любого приложения Тех. поддержка
вчераLLM Qwen – возможности, опыт использования, ограничения ИИ
вчераНе могу прошить UUID на материнской плате ASRock H55M-LE ни с помощью amidedos ни ее версии под Windows - осталось только слить дамп биоса, редактировать в HEX-редакторе и шить обратно? Тех. поддержка
2 дня[Решено] Установил и успешно запустил пакет motion под OpenWRT, но не могу зайти на его админку и получить видеопоток с помощью VLC-плеера Unix
2 дняDenon dcd-3000 ремонт Стерео
3 дняИИ шутит. ИИ
4 дняThecus 4560/2520 поиск прошивки. НАС
4 дняРоссийская инди-поп-музыка Культура
23:59вчераНовый смартфон Samsung за 120 долларов. Представлен Samsung Galaxy A06 5G: 50 Мп, 5000 мАч, IP54 и 4 года обновлений
23:51вчераMicrosoft создала квантовый чип Majorana 1: миллион кубитов на ладони
22:46вчераРоссийские астрономы обнаружили редкую двойную звезду с магнитным белым карликом Gaia 19cwm
22:35вчераПока видеокарты GeForce отсутствуют либо дорожают, пользователи активно переходят на Radeon. Как минимум так говорит статистика для Южной Кореи
22:16вчераЕсли хочется GeForce RTX 50, но максимально дёшево. Существование RTX 5050 Laptop подтверждено
22:06вчераУчёные «просветили стену» из титана в поисках аксионов на ускорителе в Гамбурге
21:59вчераGeForce RTX 5090 оказалась в семь раз слабее RTX 4090. Nvidia убрала поддержку PhysX из свежего драйвера для видеокарт RTX 50
21:48вчера800 л.с., рама и полный привод. Опубликованы фото супермощного Tank 400
21:45вчераTesla строит первую общедоступную зарядную станцию для электрогрузовиков в Лос-Анджелесе
21:39вчераNASA завершило сборку первого двигателя RS-25 для миссии «Артемида V»
20+ лет назад в этот день обсуждали
2002Citrix via Internet Сети
2002Minolta RD175 - НУЖНА ПОМОЩЬ!!!!! Фото
2002W2k - проблема с флоппи-диском. OС и сист. ПО
2002Получил по почте вирус. Лежит в базе OUTLOOKa. Как лучше поступить? Интернет
2002Как можно разогнать 500 Селерон на bx Разгон
2002Дождемся KT333A или нет? Сист. платы
2004unmountable_boot_volume Накопители
2004Как сделать, чтобы у письма обратный адрес был один, а отвечали на него на другой адрес? Интернет (архив)
2005Компьютер перезапускается во время POST Тех. поддержка
2005Выгонят ли Россию из большой восьмерки(G-8) Политика